lir_achem.chemistry_scheme

This is where the chemistry scheme is implemented

Functions

chemistry_eq(t, densities, coefficients, ...)

Computes the derivative of the different species densities according to the chesmitry scheme
lir_achem.chemistry_scheme.chemistry_eq(t, densities, coefficients, rad_here, n_here, compute_source, EUV_file, XR_file, today)[source]

Computes the derivative of the different species densities according to the chesmitry scheme

Parameters:

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

  • densities – Numpy array of species densities (cm-3). This is in the same (12, len(altitudes_D)). It is orders as [Ne, O2m, Xm, NOp, O4p, Yp, Np, O2p, N2p, Op4S, Op2P, Op2D]

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

  • rad_here – Radiation class instance

  • n_here – Neutrals class instance

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

  • EUV_file – Path to the EUV file of interest (from GOES EUV daily average)

  • XR_file – Path to the XR file of interest (from GOES 1s XR flux)

  • today – Datetime of the start of the computation

Returns:

Derivatives of the densities, in the same format as densities