예제 #1
0
        frequency_range=band,
        start_date=ephem.Date(current_date).tuple(),
        duration_seconds=target_duration_s,
        stations=stations,
        clock_mhz=200,
        bit_mode=8,
        backend=backend))

current_date += target_duration_s * ephem.second + 61 * ephem.second

observations.append(
    momxml.Observation(beam_list=[momxml.Beam(post_cal, target_subbands)],
                       antenna_set=antenna_set,
                       frequency_range=band,
                       start_date=ephem.Date(current_date).tuple(),
                       duration_seconds=cal_duration_s,
                       stations=stations,
                       clock_mhz=200,
                       bit_mode=8,
                       backend=backend))

date_folder = momxml.Folder(children=observations,
                            name='EoR-%s' % target.name,
                            description='EoR-%s' % target.name,
                            grouping_parent=True)
main_folder = momxml.Folder(children=[date_folder],
                            name='EoR-NCP',
                            description='EoR NCP commissioning observations',
                            mom_id=182243)
print momxml.xml([main_folder], 'LC0_019')
예제 #2
0
sys.stderr.write('MAIN: '+str(target) + '\n')
for cal in cal_fields:
    sys.stderr.write(' CAL: '+str(cal)+ '\n')


observations = []
current_date = transit_date - 0.5*target_duration_s*ephem.second

backend = momxml.BackendProcessing(channels_per_subband     = chan,
                                   integration_time_seconds = int_s)

observations.append(momxml.Observation(
        beam_list        = [momxml.Beam(target, target_subbands)] + [momxml.Beam(field, cal_subbands) for field in cal_fields],
        antenna_set      = antenna_set,
        frequency_range  = band,
        start_date       = ephem.Date(current_date).tuple(),
        duration_seconds = target_duration_s,
        stations         = stations,
        clock_mhz        = 200,
        backend = backend,
        bit_mode         = 8))

obs_folder = momxml.Folder(name        = 'EoR-3C196',
                           description = 'EoR 3C 196 observations',
                           children    = observations,
                           grouping_parent = True,
                           mom_id      = 449283)

with open('eor-3c196-%d%02d%02d.xml' % (start_search_date.tuple()[0:3]), 'w') as output:
    output.write(momxml.xml([obs_folder], 'LC3_028'))
예제 #3
0
        
        
current_date += target_duration_s*ephem.second + 61*ephem.second

        
observations.append(momxml.Observation(
    beam_list        = [momxml.Beam(post_cal, target_subbands)],
    antenna_set      = antenna_set,
    frequency_range  = band,
    start_date       = ephem.Date(current_date).tuple(),
    duration_seconds = cal_duration_s,
    stations         = stations,
    clock_mhz        = 200,
    bit_mode         = 8,
    backend          = backend))


date_folder = momxml.Folder(children=observations,
                            name='EoR-%s' % target.name,
                            description = 'EoR-%s' % target.name,
                            grouping_parent = True,
                            label = '0')

obs_folder = momxml.Folder(name        = 'EoR-ELAIS-N1',
                           description = 'ELAIS N1 observations for the EoR',
                           children    = [date_folder],
                           mom_id      = 260476)

print momxml.xml([obs_folder], 'LC0_019')
sys.stderr.write('MAIN: '+str(target) + '\n')
for field in aux_pointings:
    sys.stderr.write(' GRID: '+str(field)+ '\n')

all_fields = [target] + aux_pointings

observations = []
current_date = transit_date - 0.5*target_duration_s*ephem.second

backend = momxml.BackendProcessing(channels_per_subband     = chan,
                                   integration_time_seconds = int_s)

observations.append(momxml.Observation(
    beam_list        = [momxml.Beam(field, subbands) for field in all_fields],
    antenna_set      = antenna_set,
    frequency_range  = band,
    start_date       = ephem.Date(current_date).tuple(),
    duration_seconds = target_duration_s,
    stations         = stations,
    clock_mhz        = 200,
    backend          = backend,
    bit_mode         = 8))

obs_folder = momxml.Folder(name        = 'EoR-3C196',
                           description = 'EoR 3C 196 observations',
                           children    = observations,
                           mom_id      = 322103)

print momxml.xml([obs_folder], 'LC1_039')
예제 #5
0
    bit_mode         = 8,
    backend          = backend,
    initial_status   = 'approved'))

        
        
current_date += target_duration_s*ephem.second + 61*ephem.second

        
observations.append(momxml.Observation(
    beam_list        = [momxml.Beam(post_cal, target_subbands)],
    antenna_set      = antenna_set,
    frequency_range  = band,
    start_date       = ephem.Date(current_date).tuple(),
    duration_seconds = cal_duration_s,
    stations         = stations,
    clock_mhz        = 200,
    bit_mode         = 8,
    backend          = backend,
    initial_status   = 'approved'))


date_folder = momxml.Folder(children=observations, name='EoR-%s' % target.name,
                            description = 'EoR-%s' % target.name,
                            grouping_parent = True)
main_folder = momxml.Folder(children=[date_folder], name='EoR-NCP', description='EoR NCP observations')


with open('eor-ncp-%d%02d%02d.xml' % (start_date.tuple()[0:3]), 'w') as output:
    output.write(momxml.xml([main_folder], 'LT5_009', description='The LOFAR EoR project'))
예제 #6
0
start_date = (2012, 12, 16, 6, 0, 0)

all_fields = TARGETS_8bit
field_sets = [MSSSFieldSet(obs[0], obs[1:]) for obs in all_fields]

date_tagged_field_sets = []
for field_set in field_sets:
    for ha in HOUR_ANGLES:
        date_tagged_field_sets.append(
            (field_set,
             momxml.next_date_with_lofar_lst(field_set.mean_ra_rad() + ha,
                                             start_date=start_date)))

observations = []
for field_set, date in sorted(date_tagged_field_sets,
                              key=lambda x: str(ephem.Date(x[1]))):
    stderr('%20s: %s' % (field_set.name, str(ephem.Date(date))))
    observations += field_set.obs_list(date)

subfolder = momxml.Folder(name='MSSS-HBA-2012-12-16',
                          description='8 bit MSSS HBA experiment',
                          children=observations)

folder = momxml.Folder(name='MSSS commissioning',
                       description='MSSS experimental observations',
                       mom_id=173637,
                       children=[subfolder])

print momxml.as_xml_mom_project([folder], 'Commissioning2012')
예제 #7
0
                            ra_angle=ra_angle,
                            dec_angle=dec_angle))

sys.stderr.write('MAIN: ' + str(target) + '\n')
for cal in cal_fields:
    sys.stderr.write(' CAL: ' + str(cal) + '\n')

observations = []
current_date = transit_date - 0.5 * target_duration_s * ephem.second

observations.append(
    momxml.Observation(
        beam_list=[momxml.Beam(target, target_subbands)] +
        [momxml.Beam(field, cal_subbands) for field in cal_fields],
        antenna_set=antenna_set,
        frequency_range=band,
        start_date=ephem.Date(current_date).tuple(),
        duration_seconds=target_duration_s,
        stations=stations,
        clock_mhz=200,
        integration_time_seconds=int_s,
        channels_per_subband=chan,
        bit_mode=8))

obs_folder = momxml.Folder(name='EoR-3C295',
                           description='3C 295 observations for the EoR',
                           children=observations,
                           mom_id=227443)

print momxml.xml([obs_folder], 'LC0_019')
예제 #8
0
                                             stokes_downsampling_steps=4,
                                             polarizations='I')),
        bit_mode=8))

current_date += target_duration_s * ephem.second + 61 * ephem.second

observations.append(
    momxml.Observation(beam_list=[momxml.Beam(post_cal, target_subbands)],
                       antenna_set=antenna_set,
                       frequency_range=band,
                       start_date=ephem.Date(current_date).tuple(),
                       duration_seconds=cal_duration_s,
                       stations=stations,
                       clock_mhz=200,
                       backend=momxml.BackendProcessing(
                           integration_time_seconds=int_s,
                           channels_per_subband=chan),
                       bit_mode=8))

date_folder = momxml.Folder(children=observations,
                            name='%s' % target.name,
                            description='%s' % target.name)

obs_folder = momxml.Folder(
    name='ELAIS-N1',
    description='ELAIS N1 observations for Surveys, EoR, and Transients',
    children=[date_folder],
    mom_id=266709)

print momxml.as_xml_mom_project([obs_folder], 'DDT_003')
예제 #9
0
    sys.stderr.write(' CAL: ' + str(cal) + '\n')

observations = []
current_date = transit_date - 0.5 * target_duration_s * ephem.second

backend = momxml.BackendProcessing(channels_per_subband=chan,
                                   integration_time_seconds=int_s)

observations.append(
    momxml.Observation(
        beam_list=[momxml.Beam(target, target_subbands)] +
        [momxml.Beam(field, cal_subbands) for field in cal_fields],
        antenna_set=antenna_set,
        frequency_range=band,
        start_date=ephem.Date(current_date).tuple(),
        duration_seconds=target_duration_s,
        stations=stations,
        clock_mhz=200,
        backend=backend,
        bit_mode=8))

start_tuple = ephem.Date(current_date).tuple()

folder = momxml.Folder(name='3C196:11h%04d%02d%02d-%02d%02d' % start_tuple[:5],
                       children=observations,
                       grouping_parent=True)

with open('cobalt-3c196-%04d%02d%02d-%02d%02d.xml' % start_tuple[:5],
          mode='w') as out:
    out.write(momxml.xml([folder], project='COBALT'))