def test_simobserve_subcommands_syntax(self):
"""
Check all the commands to configure a simulated observation will run OK.
(Not the same as ensuring they do the right thing!)
"""
script = []
sim.open_sim(script, self.output_ms_path)
sim.setpb(script,
telescope_name='VLA',
primary_beam_hwhm=1.5 * u.degree,
frequency=2.5 * u.GHz)
sim.setconfig(script,
telescope_name='VLA',
antennalist_path=test_data.vla_antennalist)
sim.setspwindow(script,
freq_start=2.5 * u.GHz,
freq_resolution=125 * u.kHz,
freq_delta=125 * u.kHz,
n_channels=1,
)
sim.setfeed(script, )
sim.setfield(script, SkyCoord(192, -35, unit='deg'))
sim.setlimits(script)
sim.setauto(script)
ref_time = Time('2014-05-01T19:55:45', format='isot', scale='tai')
sim.settimes(script, integration_time=10 * u.s, reference_time=ref_time)
sim.observe(script, stop_delay=60 * u.s)
sim.set_simplenoise(script, noise_std_dev=0.1 * u.Jy)
sim.corrupt(script)
sim.close_sim(script)
self.casa.run_script(script)
def test_setconfig(self):
script = []
sim.setpb(script,
telescope_name='VLA',
primary_beam_hwhm=1.5 * u.degree,
frequency=2.5 * u.GHz)
self.casa.run_script(script)
def test_simobserve_subcommands_syntax(self):
"""
Check all the commands to configure a simulated observation will run OK.
(Not the same as ensuring they do the right thing!)
"""
script = []
script.append("sm.open('{}')".format(self.output_ms_path))
sim.setpb(script,
telescope_name='VLA',
primary_beam_hwhm=1.5 * u.degree,
frequency=2.5 * u.GHz)
sim.setconfig(script,
telescope_name='VLA',
antennalist_path=test_data.vla_antennalist)
sim.setspwindow(script,
freq_start=2.5 * u.GHz,
freq_resolution=125 * u.kHz,
freq_delta=125 * u.kHz,
n_channels=1,
)
sim.setfeed(script, )
sim.setfield(script, SkyCoord(192, -35, unit='deg'))
sim.setlimits(script)
sim.setauto(script)
ref_time = Time('2014-05-01T19:55:45', format='isot', scale='tai')
sim.settimes(script, integration_time=10 * u.s, reference_time=ref_time)
sim.observe(script, stop_delay=60 * u.s)
sim.set_simplenoise(script, noise_std_dev=0.1 * u.Jy)
sim.corrupt(script)
sim.close_sim(script)
self.casa.run_script(script)
def test_setconfig(self):
script = []
sim.setpb(script,
telescope_name='VLA',
primary_beam_hwhm=1.5 * u.degree,
frequency=2.5 * u.GHz)
self.casa.run_script(script)
srclist = [
(SkyCoord(centre_ra, centre_dec), 1. * u.mJy, obs_central_frequency),
(SkyCoord(centre_ra, centre_dec + primary_beam_fwhm * 0.5), 1. * u.mJy,
obs_central_frequency),
]
# Convert the sources to a CASA 'componentlist'
component_list_path = sim.make_componentlist(script, srclist,
component_list_path)
# Open the visibility file
sim.open_sim(script, output_visibility)
# Configure the virtual telescope
sim.setpb(script,
telescope_name='VLA',
primary_beam_hwhm=primary_beam_fwhm * 0.5,
frequency=obs_central_frequency)
sim.setconfig(script, telescope_name='VLA', antennalist_path='./vla.c.cfg')
sim.setspwindow(
script,
freq_start=obs_central_frequency - 0.5 * obs_frequency_bandwidth,
freq_resolution=obs_frequency_bandwidth,
freq_delta=obs_frequency_bandwidth,
n_channels=1,
)
sim.setfeed(script, )
sim.setfield(script, pointing_centre)
sim.setlimits(script)
sim.setauto(script)
ref_time = Time('2014-05-01T19:55:45', format='isot', scale='tai')
sim.settimes(script, integration_time=10 * u.s, reference_time=ref_time)
srclist = [
(SkyCoord(centre_ra, centre_dec), 1. * u.mJy, obs_central_frequency),
(SkyCoord(centre_ra, centre_dec + primary_beam_fwhm * 0.5), 1. * u.mJy,
obs_central_frequency),
]
# Convert the sources to a CASA 'componentlist'
component_list_path = sim.make_componentlist(script, srclist,
component_list_path)
# Open the visibility file
sim.open_sim(script, output_visibility)
# Configure the virtual telescope
sim.setpb(script,
telescope_name='VLA',
primary_beam_hwhm=primary_beam_fwhm * 0.5,
frequency=obs_central_frequency)
sim.setconfig(script,
telescope_name='VLA',
antennalist_path='./vla.c.cfg')
sim.setspwindow(script,
freq_start=obs_central_frequency - 0.5 * obs_frequency_bandwidth,
freq_resolution=obs_frequency_bandwidth,
freq_delta=obs_frequency_bandwidth,
n_channels=1,
)
sim.setfeed(script, )
sim.setfield(script, pointing_centre)
sim.setlimits(script)
sim.setauto(script)
ref_time = Time('2014-05-01T19:55:45', format='isot', scale='tai')
