Usage for computeAzElFromRADecMJD
Return to directory of Todd's CASA extensions
Computes the az/el for a specified J2000 RA/Dec, MJD (or date/time string) and observatory name using the
casa measures tool. See also
computeRADecFromAzElMJD,
computeAzElFromRADecLST, and
computeAzElForMS.
The current known observatories (stored in /usr/lib64/casapy/data/geodetic/Observatories/ ) are:
- ['ALMA', 'ARECIBO', 'ATCA', 'BIMA', 'CLRO', 'DRAO', 'DWL', 'GB', 'GBT', 'GMRT', 'IRAM PDB', 'IRAM_PDB', 'JCMT', 'MOPRA', 'MOST', 'NRAO12M', 'NRAO_GBT', 'PKS', 'SAO SMA', 'SMA', 'VLA', 'VLBA', 'WSRT', 'ATF', 'ATA', 'CARMA', 'ACA', 'OSF', 'OVRO_MMA', 'EVLA', 'ASKAP', 'APEX', 'SMT', 'NRO', 'ASTE', 'LOFAR', 'MeerKAT', 'KAT-7', 'EVN', 'LWA1', 'PAPER_SA', 'PAPER_GB', 'e-MERLIN', 'MERLIN2', 'Effelsberg', 'MWA32T']
Usage:
au.computeAzElFromRADecMJD(raDec, mjd, observatory='ALMA', verbose=True, degrees=False)
- raDec must either be a tuple in radians: [ra,dec] or a string of the form "hh:mm:ss.sss -dd:mm:ss.ss"
- mjd must either be a floating point value in days, or a date string in of the following formats:
- 2011/10/15 05:00:00
- 2011/10/15-05:00:00
- 2011-10-15 05:00:00
- 2011-10-15-05:00:00
- prints the az,el in degrees (unless verbose=False)
- degrees: if False, returns the [az,el] in radians; if True, returns it in degrees
Examples:
CASA <2>: au.computeAzElFromRADecMJD('18:13:58.1 -8:54:16.7', '2013-05-31, 08:08:31', observatory='SMA')
RA Dec in radians = (4.7733373323053749, -0.15541526731224065)
MJD= 56443.33925, MJDseconds= 4876704511.0, JD= 2456443.83925
MJD = 56443.3392477
Azim = 111.066 deg Elev = 26.072 deg
Out[2]: [1.938469912199654, 0.45504278347839677]
CASA<3>: au.computeAzElFromRADecMJD('18:13:58.1 -8:54:16.7', 55640, observatory='SMA')
RA Dec in radians = (4.7733373323053749, -0.15541526731224065)
Azim = -92.544 deg Elev = -19.910 deg
Out[3]: [-1.6151962971338532, -0.3474933437311884]
--
ToddHunter - 2013-05-10