# 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(mid_day) - 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
target = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]
pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.cal_source(
start_date + (target_duration_s + 2 * cal_duration_s) * ephem.second,
'HBA')
cal_fields = [momxml.simbad('3C61.1')]
ra_3c61_1 = cal_fields[0].ra_angle
dec_3c61_1 = cal_fields[0].dec_angle
ra_inc = momxml.Angle(deg=60.0)
antenna_set = 'HBA_DUAL_INNER'
band = 'HBA_LOW'
stations = momxml.station_list('nl', exclude=[])
int_s = 2.0
import ephem
import sys
# 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
target = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]
pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.cal_source(start_date+(target_duration_s+2*cal_duration_s)*ephem.second, 'HBA')
cal_fields = [momxml.simbad('3C61.1')]
ra_3c61_1 = cal_fields[0].ra_angle
dec_3c61_1 = cal_fields[0].dec_angle
ra_inc = momxml.Angle(deg = 60.0)
antenna_set = 'HBA_DUAL_INNER'
band = 'HBA_LOW'
stations = momxml.station_list('all', exclude = [])
int_s = 2.0
chan = 64
target_subbands = '77..456' # 380 subbands 115 -- 189 MHz
['B', ('+', 3, 15, 5.15), ('+', 54, 34, 43.6)],
['C', ('+', 3, 23, 46.70), ('+', 56, 49, 29.5)],
['D', ('+', 3, 24, 16.38), ('+', 52, 19, 57.7)],
['E', ('+', 3, 41, 42.22), ('+', 52, 19, 57.7)],
['F', ('+', 3, 42, 11.90), ('+', 56, 49, 29.5)],
['G', ('+', 3, 50, 53.45), ('+', 54, 34, 43.6)],
['H', ('+', 2, 40, 31.67), ('+', 61, 18, 45.6)]]
targets = [TargetSource(name,
ra_angle = momxml.Angle(shms = ra),
dec_angle = momxml.Angle(sdms = dec))
for name, ra, dec in lc0_003_fields]
cal_list = SourceCatalogue()
calibrator = momxml.simbad('3C 147')
exclude = ['CS013']
subbands = '115..131,207..223,279..295,427..443'
cal_beam = Beam(calibrator, subbands)
src_beams = [Beam(target , subbands) for target in targets[0:-1]]
src_beams += [Beam(targets[-1], '207..218')]
observations = []
current_date = start_date
while current_date < end_date:
observations.append(Observation(
beam_list = [cal_beam],
antenna_set = antenna_set,
frequency_range = 'HBA_LOW',
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
target = momxml.simbad('NCP')
target.name = 'NCP-VLBI-%4d-%02d-%02d' % start_date.tuple()[0:3]
pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.cal_source(start_date+(target_duration_s+2*cal_duration_s)*ephem.second, 'HBA')
cal_fields = [momxml.TargetSource(
name='J%02d%02d%04.1f+%02d%02d%04.1f'% (ra+dec),
ra_angle=momxml.Angle(hms=ra),
dec_angle=momxml.Angle(sdms=('+',)+dec))
for ra, dec in [
(( 1, 17, 28.7,), (89, 28, 49.1)),
((10, 57, 40.0,), (88, 58, 47.0)),
(( 1, 10, 47.3,), (87, 38, 21.1)),
((22, 36, 21.1,), (88, 14, 54.7)),
(( 6, 22, 9.8,), (87, 19, 50.0)),
(( 5, 32, 46.4,), (87, 31, 40.8)),
# 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
target = momxml.simbad('NCP')
target.name = 'NCP-VLBI-%4d-%02d-%02d' % start_date.tuple()[0:3]
pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.cal_source(
start_date + (target_duration_s + 2 * cal_duration_s) * ephem.second,
'HBA')
cal_fields = [
momxml.TargetSource(name='J%02d%02d%04.1f+%02d%02d%04.1f' % (ra + dec),
ra_angle=momxml.Angle(hms=ra),
dec_angle=momxml.Angle(sdms=('+', ) + dec))
for ra, dec in [((
1,
17,
28.7,
), (89, 28, 49.1)), ((
# 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
target = momxml.simbad('NCP')
target.name = 'NCP-%4d-%02d-%02d' % start_date.tuple()[0:3]
pre_cal = source_catalogue.cal_source(start_date, 'HBA')
post_cal = source_catalogue.cal_source(
start_date + (target_duration_s + 2 * cal_duration_s) * ephem.second,
'HBA')
cal_fields = [
momxml.simbad('3C61.1'),
momxml.TargetSource('J0121+8928', momxml.Angle(hms=(1, 17, 32.82)),
momxml.Angle(sdms=('+', 89, 28, 48.7)))
]
antenna_set = 'HBA_DUAL_INNER'
band = 'HBA_LOW'
stations = momxml.station_list('all', exclude=[])
int_s = 2.0
['D', ('+', 3, 24, 16.38), ('+', 52, 19, 57.7)],
['E', ('+', 3, 41, 42.22), ('+', 52, 19, 57.7)],
['F', ('+', 3, 42, 11.90), ('+', 56, 49, 29.5)],
['G', ('+', 3, 50, 53.45), ('+', 54, 34, 43.6)],
['H', ('+', 2, 40, 31.67), ('+', 61, 18, 45.6)]]
targets = [
TargetSource(name,
ra_angle=momxml.Angle(shms=ra),
dec_angle=momxml.Angle(sdms=dec))
for name, ra, dec in lc0_003_fields
]
cal_list = SourceCatalogue()
calibrator = momxml.simbad('3C 147')
exclude = ['CS013']
subbands = '115..131,207..223,279..295,427..443'
cal_beam = Beam(calibrator, subbands)
src_beams = [Beam(target, subbands) for target in targets[0:-1]]
src_beams += [Beam(targets[-1], '207..218')]
observations = []
current_date = start_date
while current_date < end_date:
observations.append(
Observation(beam_list=[cal_beam],
antenna_set=antenna_set,
frequency_range='HBA_LOW',
start_date=ephem.Date(current_date).tuple(),