Function pecsim::ne_iri_cap

Function Documentation

double pecsim::ne_iri_cap(double r, double alatr, double amlt, const std::array<int, 2> &itime)

This subroutine provides polar cap densities based on IRI and on the GCPM polar cap model.

powern = -2.8618 refalt = reference altitude where polar cap profile and IRI are made to agree

The mathematical form of the density model comes from Persoon et al. and Chandler et al, 1991. It approximates Alouette/ISIS, DE 1 RIMS, and DE 1 PWI observations, while allowing for the density to rise and fall with the IRI model through the refn parameter.

Parameters:

  • r – radial distance (RE)

  • alatr – geomagnetic latitude in radians

  • amlt – geomagnetic local time in hours

  • itime –

Returns:

double c The density bridge has the form: density = exp(-2.8618*log(r)+refn) where r is height in km this is obtained from the IRI density at 350km altitude and at the latitude, L-shell, magnetic local time, and local time of the point provided by the empirical model user. It is units of m-3, like the IRI total electron density.