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,\
population = pop,\
sim_root = sim_root,\
simulation_name = s.auto_sim_name('piggyback_test')
)
sim.calc_observation_params(\
duty_cycle=0.01, \
SST_fraction=0.1, \
tracking_fraction=1.0, \
SST_time_slice=0.2, \
interleaving_time_slice = 7.5, \
import space_object as so
import population as p
import plothelp
import numpy as n
import matplotlib.pyplot as plt
t=n.linspace(0,5*24*3600,num=10000,dtype=n.float)
# test propagation
#load the input population
m = p.filtered_master_catalog_factor(None,treshhold=1e-2,seed=12345,filter_name='e3d_planar_beam')
o = m.get_object(15866)
print(o)
ecefs=o.get_state(t)
fig = plt.figure(figsize=(15,15))
ax = fig.add_subplot(111, projection='3d')
ax.view_init(15, 5)
plothelp.draw_earth_grid(ax)
ax.plot(ecefs[0,:],ecefs[1,:],ecefs[2,:],".",alpha=0.5,color="black")
plt.title("Orbital propagation test ID {}".format(o.oid))
plt.show()