Hard coded solution, future = pass ec and ef array to maxt assuming > 2 populations
Set variables given from GEMINI
Convert densities and altitudes into
Call GLOW to calculate ionized and excited species, and airglow emission rates:
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed
arrows point from an interface to procedures which implement that interface.
This could include the module procedures in a generic interface or the
implementation in a submodule of an interface in a parent module.