def ns_fence_rng_plot_mov():
scan = rslib.ns_fence_rng_model(lat=69., lon=19., alt=150.)
plot_radar_scan_movie(scan, earth=True)
plt.show()
import scheduler_library as sch
import antenna_library as alib
sim_root = '/home/danielk/IRF/E3D_PA/SORTSpp_sim/piggyback_test'
#initialize the radar setup
e3d = rl.eiscat_3d()
e3d.set_FOV(max_on_axis=25.0,horizon_elevation=30.0)
e3d.set_SNR_limits(min_total_SNRdb=10.0,min_pair_SNRdb=0.0)
e3d.set_TX_bandwith(bw = 1.0e6)
e3d.set_beam('TX', alib.e3d_array_beam_stage1(opt='dense') )
e3d.set_beam('RX', alib.e3d_array_beam() )
#initialize the observing mode
e3d_scan = rslib.ns_fence_rng_model(min_el = 30.0, angle_step = 2.0, dwell_time = 0.2)
#3 by 3 grid at 300km
az_points = n.arange(0,360,45).tolist() + [0.0];
el_points = [90.0-n.arctan(50.0/300.0)*180.0/n.pi, 90.0-n.arctan(n.sqrt(2)*50.0/300.0)*180.0/n.pi]*4+[90.0];
e3d_ionosphere = rslib.n_const_pointing_model(az_points,el_points,len(az_points), dwell_time = 7.5)
e3d_scan.set_radar_location(e3d)
e3d.set_scan(SST=e3d_scan,secondary_list=[e3d_ionosphere])
#load the input population
pop = p.filtered_master_catalog_factor(e3d,treshhold=1e-2,seed=12345,filter_name='e3d_full_beam')
pop._objs = pop._objs[:2000,:]
sim = s.simulation( \
radar = e3d,\
SIM_TIME = 24.0 * 1.0
##########################
sim_root = '/ZFS_DATA/SORTSpp/FP_sims/E3D_scanning'
#initialize the radar setup
radar = rlib.eiscat_3d(beam='interp', stage=1)
radar.set_FOV(max_on_axis=90.0, horizon_elevation=30.0)
radar.set_SNR_limits(min_total_SNRdb=10.0, min_pair_SNRdb=1.0)
radar.set_TX_bandwith(bw=1.0e6)
scan = rslib.ns_fence_rng_model(
radar._tx[0].lat,
radar._tx[0].lon,
radar._tx[0].alt,
min_el=30,
angle_step=2.0,
dwell_time=0.2,
)
radar.set_scan(scan)
#load the input population
pop = plib.filtered_master_catalog_factor(
radar=radar,
treshhold=0.01,
min_inc=50,
prop_time=48.0,
)
sim = Simulation(radar=radar,
