Function pecsim::gcpm_v24

Function Documentation

std::vector<double> pecsim::gcpm_v24(DateTime datetime, double r, double amlt, double alatr, double akp = -1)

c Global Core Plasma Model Original Version 1.0, January 1, 2000

modified by dlg 8/3/2007 to include the polar cap model modified by dlg 1/6/2009 to include seasonal and solar cycle variations in inner plasmasphere from C&A 1992 corrected by dlg 6/3/2009 to fix ionosphere-plasmasphere bridge code that was trying to make bridge below the F2 peak.

Parameters:

  • datatime – date and time (in UTC)

  • r – geocentric radial distance in Re

  • amlt – solar magnetic local time in hours

  • alatr – solar magnetic latitude in radians

  • akp – planetary Kp index (-1: use the Kp index from data)

Returns:

std::vector<double> [0] = total electron density in 1/cm^3 [1] = total hydrogen density in 1/cm^3 [2] = total helium density in 1/cm^3 [3] = total oxygen density in 1/cm^3] [4] = f107 index [5] = rz12 index [6] = hmf2_km (F2 peak height in km)