rotationPhase
Return to directory of Todd's tools
Determines the rotational phase of an observation relative to a zero point
in the past.
Usage:
tt.rotationPhase(timezero, period, timelist, degrees=False)
Inputs:
- timezero: a zero point in time, such as MJD or JD
- period: the period of rotation in the same units as timezero (e.g. days)
- timelist: a list of observation times for which you want to compute the phase, in the same units as timezero
- degrees: if True, return the phases in degrees rather than 0..1
Example:
CASA <7>: mylist = [56949.39142, 56949.40422, 56949.42936, 56949.44219, 56949.46753, 56949.48041, 56949.51175, 56949.52456, 56949.54834, 56949.56119]
CASA <8>: tt.rotationPhase(50371.2416666, .3003969, mylist, degrees=True)
Out[8]:
[70.031987680815817,
85.371693246496747,
115.49983372287984,
130.87549172694091,
161.24331509679303,
176.67889381986157,
214.23720417689651,
229.58889389161067,
258.08719064880626,
273.48681694221432]
--
ToddHunter - 2015-01-13