inv_z_transform_tau=10**(3 / (max_range / antenna.wavelength / 2)))
pade_task = TroposphericRadioWaveSSPadePropagator(antenna=antenna,
env=environment,
max_range_m=max_range,
comp_params=pade_params)
pade_field = pade_task.calculate()
petool_task = PETOOLPropagationTask(antenna=antenna,
env=environment,
two_way=False,
max_range_m=max_range,
dx_wl=1,
n_dx_out=1,
dz_wl=0.05,
n_dz_out=20)
petool_field = petool_task.calculate()
trm = TwoRayModel(src=antenna, env=environment)
x_grid_m = np.linspace(1, max_range, 5000)
z_grid_m = np.linspace(1, environment.z_max, 500)
trm_f = trm.calculate(x_grid_m, z_grid_m)
trm_field = Field(x_grid=x_grid_m, z_grid=z_grid_m, freq_hz=antenna.freq_hz)
trm_field.field[:, :] = trm_f
pade_vis = FieldVisualiser(pade_field,
trans_func=lambda v: 10 * cm.log10(1e-16 + abs(v)),
label='Pade + nonlocal BC',
x_mult=1)
plt = pade_vis.plot2d(min=-40, max=0)
plt.xlabel('Range (m)')
grid_optimizator_abs_threshold=5e-3,
z_order=4,
exp_pade_order=(10, 11),
dx_wl=500,
dz_wl=0.8,
storage=PickleStorage()))
pade_field_4 = pade_task_4.calculate()
petool_task = PETOOLPropagationTask(antenna=ant,
env=env,
two_way=False,
max_range_m=150e3,
dx_wl=500,
n_dx_out=1,
dz_wl=3)
petool_field = petool_task.calculate()
env.z_max = 3000
petool_task_m = PETOOLPropagationTask(antenna=ant,
env=env,
two_way=False,
max_range_m=150e3,
dx_wl=500,
n_dx_out=1,
dz_wl=3,
n_dz_out=2)
petool_field_m = petool_task_m.calculate()
pade_vis_4 = FieldVisualiser(
pade_field_4,
env=env,