vlf4ions.compute_electron_density

Created on Wed Aug 21 10:41:33 2024

@author: pteysseyre

Functions

LMP_findminimum(station, pathtoCSVfiles)

Gives the time profile of H, B. NOTE : Make sure that the relevant csv files are stored on the computer AND that they have been detrended

compute_Ne(z, H, B)

Applies Wait's profile to H and B to obtain Ne at altitude z

get_quiet_hb(sza)

Returns an estimate of the quiet values for h' and beta (Appendice A, MacRae & Thomson, 2000)
vlf4ions.compute_electron_density.get_quiet_hb(sza)[source]

Returns an estimate of the quiet values for h’ and beta (Appendice A, MacRae & Thomson, 2000)

Parameters:

sza – Solar zenith angle (in rad)

Returns:

  • H_quiet: quiet value of h’

  • B_quiet: quiet value of beta

vlf4ions.compute_electron_density.LMP_findminimum(station, pathtoCSVfiles)[source]
Gives the time profile of H, B. NOTEMake sure that the relevant csv files are stored on the computer AND that they have been detrended

They should contain A/P values, with h and beta values in their column/rows. Lastly, because of the detrend, the modelled phase may be shifted by a multiple of 90°

Parameters:

  • station – station class instance

  • sza – Solar zenith angle at the middle of the path

  • pathtoCSVfiles – path to the files were the LMP results are stored

Returns:

H, B arrays of h’ and beta values

History : written in march 2024 by P.Teysseyre In May 2025, reworked to use DA/DP Based on golden ratio method

vlf4ions.compute_electron_density.compute_Ne(z, H, B)[source]

Applies Wait’s profile to H and B to obtain Ne at altitude z

Parameters:

  • z – Height of interest in the ionosphere (in km)

  • H – h’ values (in km)

  • B – beta value (in km-1)

Returns:

Electron density (in m^-3) at altitude z