Ground track Ex >> Ez

[JamesMCannon]

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)