computeUTForElevation
Return to directory of Todd's CASA extensions
For a specified RA/Dec, date/MJD, and elevation, compute the UT time that the
object will cross that point (by repeatedly calling
computeAzElFromRADecMJD). It also converts
this time to local Chilean time using the python pytz module, which accounts for daylight savings time changes in Chile. See also
computeAzElFromRADecLST.
Usage:
au.computeUTForElevation(elevation, radec=None, date='', mjd=0, observatory='ALMA', intervalMinutes=5, vis=None, field=None, useJPL=False)
- elevation: in degrees
- raDec: if specified, it must either be a tuple in radians: [ra,dec], or a string of the form "hh:mm:ss.sss -dd:mm:ss.ss"
- mjd: modified Julian date (optional, default = today)
- date: string date (optional), Either of these formats is valid: 2011/10/15, 2011-10-15
- observatory: location to run the calculation for
- intervalMinutes: grid size to use in the search (reduce for higher accuracy & longer computation time)
- vis: the measurement set from which to read the field coordinates
- field: the field ID or name in the measurement set to use, or the name of a planet
- useJPL: if True, query JPL Horizons for planets instead of casa ephemerides
Example:
CASA <4>: au.computeUTForElevation(30, '09:47:33 +11:25:43', date='2013-06-01')
rising: 18:15 = 14:14 Chile time
setting: 01:00 = 20:59 Chile time
CASA <11>: au.computeUTForElevation(30, field='Uranus',date='2014-04-26')
Since vis has not been specified, I will assume that Uranus is a planet.
Converted MJD=56773.000000 to 2014-04-26 00:00:00 UT
Using casa's ephemerides for Uranus. These are in J2000 coordinates in this casa version (>=4.0.0)
J2000 Position: 00:51:04.25760, +004:46:11.49600 Azim, Elev: -104.483, -41.407
rising: 11:25 UT = 08:24 Chile time
setting: 18:40 UT = 15:39 Chile time
--
ToddHunter - 2013-11-06