import sys
from numpy import array
# 17:29 - 06:10
cal_duration_s = 5*60
source_catalogue = momxml.SourceCatalogue()
mid_day = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(momxml.next_sunset(mid_day) + 20*momxml.ephem.minute)
end_date = momxml.ephem.Date(momxml.next_sunrise(start_date) - 20*momxml.ephem.minute)
total_duration_s = (end_date - start_date)*24*3600.0
target_duration_s = total_duration_s - 2*cal_duration_s - 2*61.0
targets = [momxml.simbad('NCP')]
targets[0].name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]
# targets.append(momxml.TargetSource('NCP-32A-%4d-%02d-%02d' % start_date.tuple()[0:3],
# ra_angle = momxml.Angle(hms=(9,30,0.0)),
# dec_angle = momxml.Angle(deg=86.0)))
# targets.append(momxml.TargetSource('NCP-32B-%4d-%02d-%02d' % start_date.tuple()[0:3],
# ra_angle = momxml.Angle(hms=(14,0,0.0)),
# dec_angle = momxml.Angle(deg=86.0)))
pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.find_source('3C 196')
#source_catalogue.cal_source(start_date+(target_duration_s+2*cal_duration_s)*ephem.second,
# 'HBA')
station_set = 'all'
# 17:29 - 06:10
cal_duration_s = 5 * 60
source_catalogue = momxml.SourceCatalogue()
mid_day = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(
momxml.next_sunset(mid_day) + 20 * momxml.ephem.minute)
end_date = momxml.ephem.Date(
momxml.next_sunrise(start_date) - 20 * momxml.ephem.minute)
total_duration_s = (end_date - start_date) * 24 * 3600.0
target_duration_s = total_duration_s - 2 * cal_duration_s - 2 * 61.0
targets = [momxml.simbad('NCP')]
targets[0].name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]
s3c61_1 = momxml.simbad('3C61.1')
ra_3c61_1 = s3c61_1.ra_angle
dec_3c61_1 = s3c61_1.dec_angle
aux_dec = momxml.Angle(deg=78.75000001)
ra_inc = momxml.Angle(deg=360 / 18.0)
ra_offset = momxml.Angle(deg=360 / 18.0)
aux_fields = []
for i in range(18):
aux_fields.append(
momxml.TargetSource('NCP-' + chr(ord('A') + i),
ra_angle=ra_3c61_1 + ra_inc * (i) - ra_offset,
def find_source(self, source_name):
selection = [row for row in self.source_table
if source_name in row[0]]
if len(selection) != 1:
return simbad(source_name)
return target_source_from_row(selection[0])
def find_source(self, source_name):
selection = [row for row in self.source_table if source_name in row[0]]
if len(selection) != 1:
return simbad(source_name)
return target_source_from_row(selection[0])
import sys
from numpy import array
# 17:29 - 06:10
cal_duration_s = 5*60
source_catalogue = momxml.SourceCatalogue()
mid_day = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(momxml.next_sunset(mid_day) + 20*momxml.ephem.minute)
end_date = momxml.ephem.Date(momxml.next_sunrise(start_date) - 20*momxml.ephem.minute)
total_duration_s = (end_date - start_date)*24*3600.0
target_duration_s = total_duration_s - 2*cal_duration_s - 2*61.0
targets = [momxml.simbad('NCP')]
targets[0].name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]
targets.append(momxml.simbad('3C 220.3'))
pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.find_source('3C 196')
#source_catalogue.cal_source(start_date+(target_duration_s+2*cal_duration_s)*ephem.second,
# 'HBA')
station_set = 'all'
if len(sys.argv) >= 3:
station_set = sys.argv[2]
antenna_set = 'LBA_OUTER'
band = 'LBA_HIGH'
stations = momxml.station_list(station_set, exclude = [])
import sys
from numpy import array
# 17:29 - 06:10
cal_duration_s = 5*60
source_catalogue = momxml.SourceCatalogue()
mid_day = momxml.ephem.Date(sys.argv[1])
start_date = momxml.ephem.Date(momxml.next_sunset(mid_day) + 20*momxml.ephem.minute)
end_date = momxml.ephem.Date(momxml.next_sunrise(start_date) - 20*momxml.ephem.minute)
total_duration_s = (end_date - start_date)*24*3600.0
target_duration_s = total_duration_s - 2*cal_duration_s - 2*61.0
targets = [momxml.simbad('NCP')]
targets[0].name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]
s3c61_1 = momxml.simbad('3C61.1')
ra_3c61_1 = s3c61_1.ra_angle
dec_3c61_1 = s3c61_1.dec_angle
aux_dec = momxml.Angle(deg=82.50000001)
ra_inc = momxml.Angle(deg=360/12.0)
ra_offset = momxml.Angle(deg=360/12.0/2.0)
aux_fields = []
for i in range(12):
aux_fields.append(momxml.TargetSource('NCP-'+chr(ord('A')+i),
ra_angle = ra_3c61_1 + ra_inc*(i) - ra_offset,
dec_angle = aux_dec))