lir_achem.mitra_rowe_scheme

This is where the Mitra-Rowe scheme is solve

Functions

chemistry_mr_eq(t, densities, coefficients, ...)

Computes the derivative of the different species densities according to the Mitra-Rowe chemistry scheme
lir_achem.mitra_rowe_scheme.chemistry_mr_eq(t, densities, coefficients, rad_here, n_here, Phi_EUV, compute_source, today, compute_sza=False, HXR_bins=True)[source]

Computes the derivative of the different species densities according to the Mitra-Rowe chemistry scheme

Parameters:

  • t – Time (in s) since the start of the computation

  • densities – Numpy array of species densities (cm-3). This is in the shape (7, len(altitudes_D)). It is ordered as [Ne, O2m, Xm, NOp, Yp, O2p, O4p, Phi_SXR, Phi_HXR]

  • coefficients – Chemistry coefficients, as returned by the get_mr_coefficients function

  • rad_here – Radiation class instance

  • n_here – Neutrals class instance

  • Phi_EUV – Absorbed EUV flux at D-region altitudes

  • compute_source – Boolean. If True, the ionisation source will be recomputed at each time-step

  • today – Datetime of the start of the computation (in UT)

  • compute_sza – Boolean. If True, the solar zenith angle will be recomputed at each time step and the absorption, Ch and H recomputed. Default:False

  • HXR_bins – Bollean, Default=True. If True, the ionisation will be computed using the bins in rad_here. If False, it will only use the average GOES HXR flux.

Returns:

Derivatives of the densities, in the same format as densities