Using the GroundSampler to set up equivalent samplers for Ex, Ey, and Ez, for a variety of ionospheres and grounds, Ex dominates several orders of magnitude above Ez. Ey behaves predictably with much lower amplitude and Ez behaves appropriately when compared to Lehtninen's Full Wave Method. Code used to create the figure is included below.
using LongwaveModePropagator
using LongwaveModePropagator: QE, ME
using Plots
using LaTeXStrings
using Colors, ColorSchemes
const LMP = LongwaveModePropagator
txpower = 1e6
# "standard" vertical dipole transmitter at 24 kHz
tx = Transmitter(VerticalDipole(),24e3,txpower)
# sample vertical electric field every 5 km out to 2000 km from tx
rx_z = GroundSampler(0:1e3:1500e3, Fields.Ez)
rx_y = GroundSampler(0:1e3:1500e3, Fields.Ey)
rx_x = GroundSampler(0:1e3:1500e3, Fields.Ex)
# vertical magnetic field
bfield = BField(50e-6, π/2, 0)
# "typical" earth ground
ground = Ground(10, 1e-4)
h=75
electrons = Species(QE, ME, z->waitprofile(z, h, 0.4), electroncollisionfrequency)
waveguide = HomogeneousWaveguide(bfield, electrons, ground)
# return the complex electric field, amplitude (dB uV/m), and phase (radians)
Ez, az, pz = propagate(waveguide, tx, rx_z)
Ey, ay, py = propagate(waveguide, tx, rx_y)
Ex, ax, px = propagate(waveguide, tx, rx_x)
# transform units to V/m
az_V = 10 .^((az.-120) ./20)
ay_V = 10 .^((ay.-120) ./20)
ax_V = 10 .^((ax.-120) ./20)
#----------------------------------------------------------------------------
#Plot
plot(rx_x.distance/1000, [ax_V ay_V az_V],label=["x" "y" "z"], xlabel="distance (km)", ylabel="amplitude (E V/m)",title="E, h' = 75, \$ \\beta \$ = 0.4, Ground",dpi=400,ylims = (1e-8, 1e4),xlims = (0,1500),yscale=:log10,minorgrid=true)
Using the GroundSampler to set up equivalent samplers for Ex, Ey, and Ez, for a variety of ionospheres and grounds, Ex dominates several orders of magnitude above Ez. Ey behaves predictably with much lower amplitude and Ez behaves appropriately when compared to Lehtninen's Full Wave Method. Code used to create the figure is included below.