-
Notifications
You must be signed in to change notification settings - Fork 0
/
calculate_LST.py
34 lines (25 loc) · 740 Bytes
/
calculate_LST.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
from astropy.time import Time
from astropy.coordinates import EarthLocation
from astropy import units as u
def calculate_LST(time_in,lon):
loc = EarthLocation(lat=0*u.deg,lon=lon,height=0*u.m)
time = Time(time_in,location=loc)
lst = time.sidereal_time('mean')
return lst
"""
#EXAMPLE
from calculate_LST import *
lst = calculate_LST('2017-01-02 09:06:11.5','21d25m42.3s')
print lst
#17h20m41.9048s
jd = Time(2457755.87936,format='jd')
jd
#<Time object: scale='utc' format='jd' value=2457755.87936>
lst = calculate_LST(jd,'21d25m41.0s')
lst
#<Longitude 17.34639899731436 hourangle>
print lst
#17h20m47.0364s
# Seach for Files that have GC at Transit
$ lst_select.py --ra=17.4_17.6 -C hsa7458_v001 zen.*xx.*uv
"""