Finding passes over sites

Out:

a    : 7.2000e+06   x : -1.2630e+06
e    : 1.0000e-01   y : 2.4221e+06
i    : 7.5000e+01   z : 6.3519e+06
omega: 0.0000e+00   vx: -2.1341e+03
Omega: 7.9000e+01   vy: -6.8303e+03
anom : 7.2000e+01   vz: 2.9544e+03
Temporal points: 71071
/home/danielk/IRF/IRF_GITLAB/pyant/pyant/coordinates.py:89: RuntimeWarning: invalid value encountered in greater
  proj[proj > 1.0] = 1.0
/home/danielk/IRF/IRF_GITLAB/pyant/pyant/coordinates.py:90: RuntimeWarning: invalid value encountered in less
  proj[proj < -1.0] = -1.0
/home/danielk/IRF/IRF_GITLAB/SORTS/sorts/passes.py:283: RuntimeWarning: invalid value encountered in less
  check = zenith_ang < 90.0-station.min_elevation
tx-0 passes: 2
rx-1 passes: 2
Pass 0 length tx-0: 148
Pass 0 length rx-1: 148
Full tx-0 to rx-1 pass: 148
Using the predefined paired passes function
/home/danielk/IRF/IRF_GITLAB/SORTS/sorts/passes.py:341: RuntimeWarning: invalid value encountered in less
  check_st = zenith_st_ang < 90.0-station.min_elevation
tx-0 and rx-1 passes: 2
Full tx-0 to rx-1 pass: 148
------------------------ Performance analysis -----------------------
 Name                     |   Executions | Mean time     | Total time
--------------------------+--------------+---------------+--------------
 equidistant_sampling     |            1 | 2.80380e-04 s | 0.05 %
 propagate                |            1 | 5.49081e-01 s | 94.52 %
 find_passes              |            2 | 7.47430e-03 s | 2.57 %
 sim_passes               |            1 | 1.21832e-04 s | 0.02 %
 find_simultaneous_passes |            1 | 1.56384e-02 s | 2.69 %
 total                    |            1 | 5.80912e-01 s | 100.00 %
---------------------------------------------------------------------

import numpy as np
import matplotlib.pyplot as plt
import pyorb

import sorts
eiscat3d = sorts.radars.eiscat3d
from sorts.profiling import Profiler

from sorts.propagator import SGP4
Prop_cls = SGP4
Prop_opts = dict(
    settings = dict(
        out_frame='ITRF',
    ),
)
prop = Prop_cls(**Prop_opts)

p = Profiler()
p.start('total')


orb = pyorb.Orbit(M0 = pyorb.M_earth, direct_update=True, auto_update=True, degrees=True, a=7200e3, e=0.1, i=75, omega=0, Omega=79, anom=72, epoch=53005.0)
print(orb)

p.start('equidistant_sampling')
t = sorts.equidistant_sampling(
    orbit = orb,
    start_t = 0,
    end_t = 3600*24*1,
    max_dpos=1e4,
)
p.stop('equidistant_sampling')

print(f'Temporal points: {len(t)}')

p.start('propagate')

states = prop.propagate(t, orb.cartesian[:,0], orb.epoch, A=1.0, C_R = 1.0, C_D = 1.0)
p.stop('propagate')

# fig = plt.figure(figsize=(15,15))
# ax = fig.add_subplot(111, projection='3d')
# ax.plot(states[0,:], states[1,:], states[2,:],"-b")
# plt.show()


p.start('find_passes')
passes_tx0 = sorts.find_passes(t, states, eiscat3d.tx[0])
print(f'tx-0 passes: {len(passes_tx0)}')
p.stop('find_passes')

p.start('find_passes')
passes_rx1 = sorts.find_passes(t, states, eiscat3d.rx[1])
print(f'rx-1 passes: {len(passes_rx1)}')
p.stop('find_passes')

p.start('sim_passes')
#finding simultaneous passes
chtx0 = np.full((len(t),), False, dtype=np.bool)
chtx0[passes_tx0[0].inds] = True
chrx1 = np.full((len(t),), False, dtype=np.bool)
chrx1[passes_rx1[0].inds] = True

inds = np.where(np.logical_and(chtx0, chrx1))[0]
p.stop('sim_passes')

print(f'Pass 0 length tx-0: {len(passes_tx0[0].inds)}')
print(f'Pass 0 length rx-1: {len(passes_rx1[0].inds)}')
print(f'Full tx-0 to rx-1 pass: {len(inds)}')

print('Using the predefined paired passes function')

p.start('find_simultaneous_passes')
tx0rx1_passes = sorts.find_simultaneous_passes(t, states, [eiscat3d.tx[0], eiscat3d.rx[1]])
p.stop('find_simultaneous_passes')

print(f'tx-0 and rx-1 passes: {len(tx0rx1_passes)}')
print(f'Full tx-0 to rx-1 pass: {len(tx0rx1_passes[0].inds)}')


p.stop('total')
print(p.fmt(normalize='total'))