def test6():
tel = TelescopeAnalysis()
tel.add_uniform_core(400, 500)
tel.add_log_spiral(12 * 3, 535, 1500, 0.5, 3)
tel.lon_deg = 116.631289
tel.lat_deg = -26.697024
tel.dec_deg = tel.lat_deg
# tel.plot_layout(filename='zzz_layout.png')
# tel.obs_length_h = 4
# tel.num_times = int((tel.obs_length_h * 3600) / 30)
tel.obs_length_h = 0
tel.num_times = 1
tel.grid_cell_size_m = tel.station_diameter_m / 2
tel.gen_uvw_coords()
# tel.grid_uvw_coords()
# tel.plot_grid('zzz_uv_grid_%02.1fh_%i.png' %
# (tel.obs_length_h, tel.num_times))
tel.uv_hist(100)
def test7():
tel = TelescopeAnalysis()
# Settings
tel.lon_deg = 116.631289
tel.lat_deg = -26.697024
tel.dec_deg = tel.lat_deg
tel.obs_length_h = 0
tel.num_times = 1
tel.grid_cell_size_m = tel.station_diameter_m / 2
# Telescope 1: ska v5 < 1500m
tel.add_ska1_v5(r_max=1500)
tel.plot_layout(filename='zzz_layout_ref.png', plot_decorations=False,
x_lim=[-1.5e3, 1.5e3], y_lim=[-1.5e3, 1.5e3],
plot_radii=[500, 1500])
tel.gen_uvw_coords()
tel.grid_uvw_coords()
tel.plot_grid(filename='zzz_uvw_ref_%02.1fh.png' % tel.obs_length_h,
plot_radii=[1500],
x_lim=[-3e3, 3e3], y_lim=[-3e3, 3e3])
tel.obs_length_h = 4
tel.num_times = int((tel.obs_length_h * 3600) / 60)
tel.gen_uvw_coords()
tel.grid_uvw_coords()
tel.plot_grid(filename='zzz_uvw_ref_%02.1fh.png' % tel.obs_length_h,
plot_radii=[1500],
x_lim=[-3e3, 3e3], y_lim=[-3e3, 3e3])
# Clear the telescope model
tel.clear()
# Telescope 2
tel.add_ska1_v5(r_max=500)
# tel.add_log_spiral_clusters(3 * 3, 3, 535, 1500, 0.5, 6, 70)
tel.add_log_spiral(3 * 6 * 3, 535, 1500, 0.5, 3)
tel.plot_layout(plot_decorations=False, filename='zzz_layout_2.png',
x_lim=[-1.5e3, 1.5e3], y_lim=[-1.5e3, 1.5e3],
plot_radii=[500, 1500])
tel.gen_uvw_coords()
tel.grid_uvw_coords()
tel.plot_grid(filename='zzz_uvw_2_%02.1fh.png' % tel.obs_length_h,
plot_radii=[1500],
x_lim=[-3e3, 3e3], y_lim=[-3e3, 3e3])
tel.obs_length_h = 4
tel.num_times = int((tel.obs_length_h * 3600) / 60)
tel.gen_uvw_coords()
tel.grid_uvw_coords()
tel.plot_grid(filename='zzz_uvw_2_%02.1fh.png' % tel.obs_length_h,
plot_radii=[1500],
x_lim=[-3e3, 3e3], y_lim=[-3e3, 3e3])
def test5():
tel = TelescopeAnalysis()
tel.add_ska1_v5(r_min=0, r_max=1000)
# tel.add_uniform_core(368, 5000)
fig = tel.plot_layout(filename='zzz_layout.png')
tel.lon_deg = 116.631289
tel.lat_deg = -26.697024
tel.dec_deg = tel.lat_deg
tel.obs_length_h = 4
tel.num_times = int((tel.obs_length_h * 3600) / 30)
tel.grid_cell_size_m = tel.station_diameter_m / 2
tel.gen_uvw_coords()
tel.grid_uvw_coords()
fig = tel.plot_grid()
fig.savefig('zzz_uv_grid_%02.1fh_%i.png' % (tel.obs_length_h, tel.num_times))
plt.close(fig)
