ionosphere_tables.h

Go to the documentation of this file.

/*
 * This file is part of Vlasiator.
 * Copyright 2010-2016 Finnish Meteorological Institute
 *
 * For details of usage, see the COPYING file and read the "Rules of the Road"
 * at http://www.physics.helsinki.fi/vlasiator/
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License along
 * with this program; if not, write to the Free Software Foundation, Inc.,
 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
 */

#pragma once
 
// Tabulated chapman function for atmospheric absorption of EUV values, from Laundal et al. 2024
const static Real chapman_euv_table[1201] = {
   1.00000e+00, 9.99998e-01, 9.99994e-01, 9.99986e-01, 9.99976e-01, 9.99962e-01, 9.99946e-01, 9.99926e-01, 9.99904e-01, 9.99878e-01,
   9.99850e-01, 9.99818e-01, 9.99784e-01, 9.99746e-01, 9.99706e-01, 9.99662e-01, 9.99616e-01, 9.99566e-01, 9.99514e-01, 9.99458e-01,
   9.99400e-01, 9.99338e-01, 9.99274e-01, 9.99206e-01, 9.99136e-01, 9.99062e-01, 9.98986e-01, 9.98906e-01, 9.98824e-01, 9.98738e-01,
   9.98650e-01, 9.98558e-01, 9.98464e-01, 9.98366e-01, 9.98266e-01, 9.98162e-01, 9.98056e-01, 9.97946e-01, 9.97834e-01, 9.97719e-01,
   9.97600e-01, 9.97479e-01, 9.97354e-01, 9.97227e-01, 9.97096e-01, 9.96963e-01, 9.96827e-01, 9.96687e-01, 9.96545e-01, 9.96399e-01,
   9.96251e-01, 9.96100e-01, 9.95945e-01, 9.95788e-01, 9.95628e-01, 9.95464e-01, 9.95298e-01, 9.95129e-01, 9.94957e-01, 9.94781e-01,
   9.94603e-01, 9.94422e-01, 9.94238e-01, 9.94051e-01, 9.93860e-01, 9.93667e-01, 9.93471e-01, 9.93272e-01, 9.93070e-01, 9.92865e-01,
   9.92657e-01, 9.92446e-01, 9.92232e-01, 9.92015e-01, 9.91795e-01, 9.91572e-01, 9.91346e-01, 9.91117e-01, 9.90885e-01, 9.90651e-01,
   9.90413e-01, 9.90172e-01, 9.89928e-01, 9.89682e-01, 9.89432e-01, 9.89179e-01, 9.88924e-01, 9.88665e-01, 9.88404e-01, 9.88139e-01,
   9.87872e-01, 9.87601e-01, 9.87328e-01, 9.87051e-01, 9.86772e-01, 9.86490e-01, 9.86205e-01, 9.85916e-01, 9.85625e-01, 9.85331e-01,
   9.85034e-01, 9.84734e-01, 9.84431e-01, 9.84125e-01, 9.83817e-01, 9.83505e-01, 9.83190e-01, 9.82872e-01, 9.82552e-01, 9.82228e-01,
   9.81901e-01, 9.81572e-01, 9.81240e-01, 9.80904e-01, 9.80566e-01, 9.80225e-01, 9.79880e-01, 9.79533e-01, 9.79183e-01, 9.78830e-01,
   9.78474e-01, 9.78116e-01, 9.77754e-01, 9.77389e-01, 9.77021e-01, 9.76651e-01, 9.76277e-01, 9.75901e-01, 9.75522e-01, 9.75139e-01,
   9.74754e-01, 9.74366e-01, 9.73975e-01, 9.73581e-01, 9.73184e-01, 9.72784e-01, 9.72382e-01, 9.71976e-01, 9.71567e-01, 9.71156e-01,
   9.70742e-01, 9.70324e-01, 9.69904e-01, 9.69481e-01, 9.69055e-01, 9.68626e-01, 9.68195e-01, 9.67760e-01, 9.67322e-01, 9.66882e-01,
   9.66438e-01, 9.65992e-01, 9.65543e-01, 9.65091e-01, 9.64636e-01, 9.64178e-01, 9.63718e-01, 9.63254e-01, 9.62787e-01, 9.62318e-01,
   9.61846e-01, 9.61371e-01, 9.60893e-01, 9.60412e-01, 9.59928e-01, 9.59441e-01, 9.58952e-01, 9.58460e-01, 9.57964e-01, 9.57466e-01,
   9.56965e-01, 9.56461e-01, 9.55955e-01, 9.55445e-01, 9.54933e-01, 9.54417e-01, 9.53899e-01, 9.53378e-01, 9.52855e-01, 9.52328e-01,
   9.51798e-01, 9.51266e-01, 9.50731e-01, 9.50193e-01, 9.49652e-01, 9.49108e-01, 9.48561e-01, 9.48012e-01, 9.47460e-01, 9.46904e-01,
   9.46347e-01, 9.45786e-01, 9.45222e-01, 9.44656e-01, 9.44087e-01, 9.43514e-01, 9.42940e-01, 9.42362e-01, 9.41781e-01, 9.41198e-01,
   9.40612e-01, 9.40023e-01, 9.39431e-01, 9.38837e-01, 9.38239e-01, 9.37639e-01, 9.37036e-01, 9.36430e-01, 9.35822e-01, 9.35210e-01,
   9.34596e-01, 9.33979e-01, 9.33359e-01, 9.32737e-01, 9.32111e-01, 9.31483e-01, 9.30852e-01, 9.30219e-01, 9.29582e-01, 9.28943e-01,
   9.28301e-01, 9.27656e-01, 9.27009e-01, 9.26358e-01, 9.25705e-01, 9.25049e-01, 9.24391e-01, 9.23729e-01, 9.23065e-01, 9.22398e-01,
   9.21729e-01, 9.21056e-01, 9.20381e-01, 9.19703e-01, 9.19022e-01, 9.18339e-01, 9.17653e-01, 9.16964e-01, 9.16273e-01, 9.15578e-01,
   9.14881e-01, 9.14181e-01, 9.13479e-01, 9.12774e-01, 9.12066e-01, 9.11355e-01, 9.10642e-01, 9.09925e-01, 9.09207e-01, 9.08485e-01,
   9.07761e-01, 9.07034e-01, 9.06304e-01, 9.05572e-01, 9.04837e-01, 9.04099e-01, 9.03359e-01, 9.02616e-01, 9.01870e-01, 9.01121e-01,
   9.00370e-01, 8.99616e-01, 8.98860e-01, 8.98100e-01, 8.97338e-01, 8.96574e-01, 8.95807e-01, 8.95037e-01, 8.94264e-01, 8.93489e-01,
   8.92711e-01, 8.91930e-01, 8.91147e-01, 8.90361e-01, 8.89573e-01, 8.88782e-01, 8.87988e-01, 8.87191e-01, 8.86392e-01, 8.85591e-01,
   8.84786e-01, 8.83979e-01, 8.83170e-01, 8.82357e-01, 8.81542e-01, 8.80725e-01, 8.79905e-01, 8.79082e-01, 8.78257e-01, 8.77429e-01,
   8.76598e-01, 8.75765e-01, 8.74929e-01, 8.74091e-01, 8.73250e-01, 8.72406e-01, 8.71560e-01, 8.70712e-01, 8.69860e-01, 8.69006e-01,
   8.68150e-01, 8.67291e-01, 8.66429e-01, 8.65565e-01, 8.64698e-01, 8.63829e-01, 8.62957e-01, 8.62082e-01, 8.61205e-01, 8.60326e-01,
   8.59444e-01, 8.58559e-01, 8.57672e-01, 8.56782e-01, 8.55890e-01, 8.54995e-01, 8.54097e-01, 8.53198e-01, 8.52295e-01, 8.51390e-01,
   8.50483e-01, 8.49573e-01, 8.48660e-01, 8.47745e-01, 8.46828e-01, 8.45907e-01, 8.44985e-01, 8.44060e-01, 8.43132e-01, 8.42202e-01,
   8.41270e-01, 8.40335e-01, 8.39397e-01, 8.38457e-01, 8.37515e-01, 8.36570e-01, 8.35622e-01, 8.34672e-01, 8.33720e-01, 8.32765e-01,
   8.31808e-01, 8.30848e-01, 8.29886e-01, 8.28921e-01, 8.27954e-01, 8.26984e-01, 8.26012e-01, 8.25038e-01, 8.24061e-01, 8.23082e-01,
   8.22100e-01, 8.21116e-01, 8.20129e-01, 8.19140e-01, 8.18149e-01, 8.17155e-01, 8.16159e-01, 8.15160e-01, 8.14159e-01, 8.13156e-01,
   8.12150e-01, 8.11141e-01, 8.10131e-01, 8.09118e-01, 8.08102e-01, 8.07085e-01, 8.06064e-01, 8.05042e-01, 8.04017e-01, 8.02990e-01,
   8.01960e-01, 8.00928e-01, 7.99894e-01, 7.98857e-01, 7.97818e-01, 7.96777e-01, 7.95733e-01, 7.94687e-01, 7.93638e-01, 7.92588e-01,
   7.91535e-01, 7.90479e-01, 7.89421e-01, 7.88361e-01, 7.87299e-01, 7.86234e-01, 7.85167e-01, 7.84098e-01, 7.83027e-01, 7.81953e-01,
   7.80877e-01, 7.79798e-01, 7.78717e-01, 7.77634e-01, 7.76549e-01, 7.75462e-01, 7.74372e-01, 7.73280e-01, 7.72185e-01, 7.71089e-01,
   7.69990e-01, 7.68889e-01, 7.67785e-01, 7.66680e-01, 7.65572e-01, 7.64462e-01, 7.63350e-01, 7.62235e-01, 7.61118e-01, 7.59999e-01,
   7.58878e-01, 7.57755e-01, 7.56629e-01, 7.55501e-01, 7.54371e-01, 7.53239e-01, 7.52104e-01, 7.50968e-01, 7.49829e-01, 7.48688e-01,
   7.47545e-01, 7.46399e-01, 7.45252e-01, 7.44102e-01, 7.42950e-01, 7.41796e-01, 7.40640e-01, 7.39482e-01, 7.38321e-01, 7.37159e-01,
   7.35994e-01, 7.34827e-01, 7.33658e-01, 7.32487e-01, 7.31313e-01, 7.30138e-01, 7.28961e-01, 7.27781e-01, 7.26599e-01, 7.25415e-01,
   7.24229e-01, 7.23041e-01, 7.21851e-01, 7.20659e-01, 7.19465e-01, 7.18268e-01, 7.17070e-01, 7.15869e-01, 7.14667e-01, 7.13462e-01,
   7.12256e-01, 7.11047e-01, 7.09836e-01, 7.08623e-01, 7.07408e-01, 7.06191e-01, 7.04972e-01, 7.03751e-01, 7.02528e-01, 7.01303e-01,
   7.00076e-01, 6.98847e-01, 6.97616e-01, 6.96383e-01, 6.95148e-01, 6.93911e-01, 6.92672e-01, 6.91431e-01, 6.90188e-01, 6.88943e-01,
   6.87696e-01, 6.86447e-01, 6.85196e-01, 6.83943e-01, 6.82689e-01, 6.81432e-01, 6.80173e-01, 6.78913e-01, 6.77650e-01, 6.76386e-01,
   6.75119e-01, 6.73851e-01, 6.72581e-01, 6.71308e-01, 6.70034e-01, 6.68758e-01, 6.67480e-01, 6.66201e-01, 6.64919e-01, 6.63635e-01,
   6.62350e-01, 6.61063e-01, 6.59774e-01, 6.58482e-01, 6.57190e-01, 6.55895e-01, 6.54598e-01, 6.53300e-01, 6.51999e-01, 6.50697e-01,
   6.49393e-01, 6.48088e-01, 6.46780e-01, 6.45470e-01, 6.44159e-01, 6.42846e-01, 6.41531e-01, 6.40215e-01, 6.38896e-01, 6.37576e-01,
   6.36254e-01, 6.34930e-01, 6.33604e-01, 6.32277e-01, 6.30948e-01, 6.29617e-01, 6.28284e-01, 6.26950e-01, 6.25614e-01, 6.24276e-01,
   6.22936e-01, 6.21595e-01, 6.20252e-01, 6.18907e-01, 6.17561e-01, 6.16212e-01, 6.14862e-01, 6.13511e-01, 6.12157e-01, 6.10802e-01,
   6.09446e-01, 6.08087e-01, 6.06727e-01, 6.05366e-01, 6.04002e-01, 6.02637e-01, 6.01271e-01, 5.99902e-01, 5.98532e-01, 5.97161e-01,
   5.95787e-01, 5.94413e-01, 5.93036e-01, 5.91658e-01, 5.90278e-01, 5.88897e-01, 5.87514e-01, 5.86129e-01, 5.84743e-01, 5.83356e-01,
   5.81966e-01, 5.80576e-01, 5.79183e-01, 5.77789e-01, 5.76394e-01, 5.74997e-01, 5.73598e-01, 5.72198e-01, 5.70796e-01, 5.69393e-01,
   5.67988e-01, 5.66582e-01, 5.65174e-01, 5.63765e-01, 5.62354e-01, 5.60942e-01, 5.59528e-01, 5.58113e-01, 5.56696e-01, 5.55278e-01,
   5.53859e-01, 5.52437e-01, 5.51015e-01, 5.49591e-01, 5.48165e-01, 5.46739e-01, 5.45310e-01, 5.43881e-01, 5.42449e-01, 5.41017e-01,
   5.39583e-01, 5.38148e-01, 5.36711e-01, 5.35273e-01, 5.33833e-01, 5.32392e-01, 5.30950e-01, 5.29507e-01, 5.28062e-01, 5.26615e-01,
   5.25168e-01, 5.23719e-01, 5.22268e-01, 5.20817e-01, 5.19364e-01, 5.17909e-01, 5.16454e-01, 5.14997e-01, 5.13539e-01, 5.12079e-01,
   5.10619e-01, 5.09157e-01, 5.07693e-01, 5.06229e-01, 5.04763e-01, 5.03296e-01, 5.01828e-01, 5.00358e-01, 4.98887e-01, 4.97416e-01,
   4.95942e-01, 4.94468e-01, 4.92992e-01, 4.91516e-01, 4.90038e-01, 4.88559e-01, 4.87078e-01, 4.85597e-01, 4.84114e-01, 4.82630e-01,
   4.81146e-01, 4.79660e-01, 4.78172e-01, 4.76684e-01, 4.75195e-01, 4.73704e-01, 4.72213e-01, 4.70720e-01, 4.69226e-01, 4.67731e-01,
   4.66235e-01, 4.64738e-01, 4.63240e-01, 4.61741e-01, 4.60241e-01, 4.58740e-01, 4.57237e-01, 4.55734e-01, 4.54230e-01, 4.52725e-01,
   4.51219e-01, 4.49711e-01, 4.48203e-01, 4.46694e-01, 4.45184e-01, 4.43673e-01, 4.42161e-01, 4.40648e-01, 4.39134e-01, 4.37619e-01,
   4.36103e-01, 4.34587e-01, 4.33069e-01, 4.31551e-01, 4.30031e-01, 4.28511e-01, 4.26990e-01, 4.25468e-01, 4.23945e-01, 4.22422e-01,
   4.20897e-01, 4.19372e-01, 4.17846e-01, 4.16319e-01, 4.14791e-01, 4.13263e-01, 4.11734e-01, 4.10204e-01, 4.08673e-01, 4.07142e-01,
   4.05609e-01, 4.04076e-01, 4.02543e-01, 4.01008e-01, 3.99473e-01, 3.97937e-01, 3.96401e-01, 3.94864e-01, 3.93326e-01, 3.91787e-01,
   3.90248e-01, 3.88708e-01, 3.87168e-01, 3.85627e-01, 3.84085e-01, 3.82543e-01, 3.81000e-01, 3.79457e-01, 3.77913e-01, 3.76369e-01,
   3.74824e-01, 3.73278e-01, 3.71732e-01, 3.70185e-01, 3.68638e-01, 3.67091e-01, 3.65542e-01, 3.63994e-01, 3.62445e-01, 3.60896e-01,
   3.59346e-01, 3.57795e-01, 3.56245e-01, 3.54694e-01, 3.53142e-01, 3.51590e-01, 3.50038e-01, 3.48485e-01, 3.46932e-01, 3.45379e-01,
   3.43826e-01, 3.42272e-01, 3.40718e-01, 3.39163e-01, 3.37608e-01, 3.36053e-01, 3.34498e-01, 3.32943e-01, 3.31387e-01, 3.29831e-01,
   3.28275e-01, 3.26719e-01, 3.25163e-01, 3.23606e-01, 3.22050e-01, 3.20493e-01, 3.18936e-01, 3.17379e-01, 3.15822e-01, 3.14265e-01,
   3.12708e-01, 3.11150e-01, 3.09593e-01, 3.08036e-01, 3.06479e-01, 3.04922e-01, 3.03365e-01, 3.01808e-01, 3.00251e-01, 2.98694e-01,
   2.97137e-01, 2.95580e-01, 2.94024e-01, 2.92467e-01, 2.90911e-01, 2.89355e-01, 2.87799e-01, 2.86244e-01, 2.84689e-01, 2.83134e-01,
   2.81579e-01, 2.80024e-01, 2.78470e-01, 2.76916e-01, 2.75363e-01, 2.73810e-01, 2.72257e-01, 2.70705e-01, 2.69153e-01, 2.67601e-01,
   2.66050e-01, 2.64500e-01, 2.62950e-01, 2.61400e-01, 2.59851e-01, 2.58303e-01, 2.56755e-01, 2.55208e-01, 2.53661e-01, 2.52116e-01,
   2.50570e-01, 2.49026e-01, 2.47482e-01, 2.45939e-01, 2.44397e-01, 2.42855e-01, 2.41315e-01, 2.39775e-01, 2.38236e-01, 2.36698e-01,
   2.35161e-01, 2.33624e-01, 2.32089e-01, 2.30555e-01, 2.29022e-01, 2.27489e-01, 2.25958e-01, 2.24428e-01, 2.22899e-01, 2.21372e-01,
   2.19845e-01, 2.18320e-01, 2.16795e-01, 2.15273e-01, 2.13751e-01, 2.12231e-01, 2.10712e-01, 2.09194e-01, 2.07678e-01, 2.06163e-01,
   2.04650e-01, 2.03139e-01, 2.01628e-01, 2.00120e-01, 1.98613e-01, 1.97108e-01, 1.95604e-01, 1.94102e-01, 1.92602e-01, 1.91104e-01,
   1.89607e-01, 1.88112e-01, 1.86619e-01, 1.85129e-01, 1.83640e-01, 1.82153e-01, 1.80668e-01, 1.79185e-01, 1.77704e-01, 1.76226e-01,
   1.74750e-01, 1.73276e-01, 1.71804e-01, 1.70334e-01, 1.68867e-01, 1.67403e-01, 1.65941e-01, 1.64481e-01, 1.63024e-01, 1.61569e-01,
   1.60118e-01, 1.58668e-01, 1.57222e-01, 1.55778e-01, 1.54338e-01, 1.52900e-01, 1.51465e-01, 1.50033e-01, 1.48604e-01, 1.47178e-01,
   1.45755e-01, 1.44336e-01, 1.42919e-01, 1.41506e-01, 1.40097e-01, 1.38691e-01, 1.37288e-01, 1.35889e-01, 1.34493e-01, 1.33101e-01,
   1.31713e-01, 1.30328e-01, 1.28948e-01, 1.27571e-01, 1.26198e-01, 1.24829e-01, 1.23465e-01, 1.22104e-01, 1.20748e-01, 1.19396e-01,
   1.18048e-01, 1.16704e-01, 1.15365e-01, 1.14031e-01, 1.12701e-01, 1.11376e-01, 1.10056e-01, 1.08740e-01, 1.07430e-01, 1.06124e-01,
   1.04823e-01, 1.03528e-01, 1.02237e-01, 1.00952e-01, 9.96726e-02, 9.83982e-02, 9.71294e-02, 9.58661e-02, 9.46084e-02, 9.33564e-02,
   9.21103e-02, 9.08700e-02, 8.96357e-02, 8.84075e-02, 8.71853e-02, 8.59694e-02, 8.47597e-02, 8.35565e-02, 8.23596e-02, 8.11693e-02,
   7.99857e-02, 7.88087e-02, 7.76385e-02, 7.64752e-02, 7.53189e-02, 7.41695e-02, 7.30274e-02, 7.18924e-02, 7.07647e-02, 6.96444e-02,
   6.85316e-02, 6.74263e-02, 6.63287e-02, 6.52388e-02, 6.41568e-02, 6.30826e-02, 6.20164e-02, 6.09582e-02, 5.99083e-02, 5.88665e-02,
   5.78331e-02, 5.68081e-02, 5.57915e-02, 5.47836e-02, 5.37842e-02, 5.27936e-02, 5.18118e-02, 5.08389e-02, 4.98749e-02, 4.89200e-02,
   4.79742e-02, 4.70376e-02, 4.61102e-02, 4.51922e-02, 4.42836e-02, 4.33844e-02, 4.24948e-02, 4.16148e-02, 4.07445e-02, 3.98839e-02,
   3.90330e-02, 3.81920e-02, 3.73610e-02, 3.65398e-02, 3.57287e-02, 3.49277e-02, 3.41367e-02, 3.33559e-02, 3.25853e-02, 3.18250e-02,
   3.10749e-02, 3.03351e-02, 2.96056e-02, 2.88865e-02, 2.81778e-02, 2.74795e-02, 2.67917e-02, 2.61143e-02, 2.54473e-02, 2.47909e-02,
   2.41449e-02, 2.35094e-02, 2.28844e-02, 2.22698e-02, 2.16658e-02, 2.10721e-02, 2.04890e-02, 1.99163e-02, 1.93539e-02, 1.88020e-02,
   1.82604e-02, 1.77292e-02, 1.72082e-02, 1.66975e-02, 1.61970e-02, 1.57066e-02, 1.52263e-02, 1.47561e-02, 1.42959e-02, 1.38456e-02,
   1.34051e-02, 1.29745e-02, 1.25536e-02, 1.21423e-02, 1.17405e-02, 1.13483e-02, 1.09654e-02, 1.05919e-02, 1.02275e-02, 9.87231e-03,
   9.52609e-03, 9.18877e-03, 8.86026e-03, 8.54045e-03, 8.22922e-03, 7.92647e-03, 7.63207e-03, 7.34591e-03, 7.06786e-03, 6.79782e-03,
   6.53564e-03, 6.28121e-03, 6.03439e-03, 5.79506e-03, 5.56309e-03, 5.33834e-03, 5.12068e-03, 4.90998e-03, 4.70609e-03, 4.50889e-03,
   4.31824e-03, 4.13400e-03, 3.95602e-03, 3.78418e-03, 3.61834e-03, 3.45835e-03, 3.30408e-03, 3.15538e-03, 3.01214e-03, 2.87420e-03,
   2.74142e-03, 2.61369e-03, 2.49085e-03, 2.37278e-03, 2.25934e-03, 2.15041e-03, 2.04585e-03, 1.94554e-03, 1.84934e-03, 1.75713e-03,
   1.66879e-03, 1.58420e-03, 1.50323e-03, 1.42578e-03, 1.35171e-03, 1.28092e-03, 1.21329e-03, 1.14872e-03, 1.08710e-03, 1.02832e-03,
   9.72276e-04, 9.18868e-04, 8.67998e-04, 8.19568e-04, 7.73485e-04, 7.29657e-04, 6.87995e-04, 6.48410e-04, 6.10819e-04, 5.75138e-04,
   5.41288e-04, 5.09191e-04, 4.78772e-04, 4.49958e-04, 4.22677e-04, 3.96862e-04, 3.72446e-04, 3.49365e-04, 3.27557e-04, 3.06963e-04,
   2.87525e-04, 2.69188e-04, 2.51898e-04, 2.35604e-04, 2.20256e-04, 2.05808e-04, 1.92213e-04, 1.79428e-04, 1.67411e-04, 1.56122e-04,
   1.45521e-04, 1.35573e-04, 1.26243e-04, 1.17495e-04, 1.09299e-04, 1.01624e-04, 9.44401e-05, 8.77197e-05, 8.14363e-05, 7.55646e-05,
   7.00806e-05, 6.49614e-05, 6.01853e-05, 5.57318e-05, 5.15813e-05, 4.77152e-05, 4.41161e-05, 4.07673e-05, 3.76532e-05, 3.47588e-05,
   3.20701e-05, 2.95738e-05, 2.72576e-05, 2.51095e-05, 2.31184e-05, 2.12740e-05, 1.95663e-05, 1.79861e-05, 1.65247e-05, 1.51739e-05,
   1.39260e-05, 1.27738e-05, 1.17107e-05, 1.07302e-05, 9.82645e-06, 8.99392e-06, 8.22742e-06, 7.52213e-06, 6.87351e-06, 6.27736e-06,
   5.72974e-06, 5.22700e-06, 4.76572e-06, 4.34272e-06, 3.95505e-06, 3.59996e-06, 3.27491e-06, 2.97753e-06, 2.70562e-06, 2.45715e-06,
   2.23022e-06, 2.02309e-06, 1.83415e-06, 1.66190e-06, 1.50496e-06, 1.36205e-06, 1.23199e-06, 1.11370e-06, 1.00619e-06, 9.08519e-07,
   8.19848e-07, 7.39395e-07, 6.66443e-07, 6.00333e-07, 5.40459e-07, 4.86265e-07, 4.37244e-07, 3.92929e-07, 3.52892e-07, 3.16743e-07,
   2.84124e-07, 2.54710e-07, 2.28201e-07, 2.04325e-07, 1.82836e-07, 1.63505e-07, 1.46128e-07, 1.30516e-07, 1.16500e-07, 1.03924e-07,
   9.26468e-08, 8.25417e-08, 7.34923e-08, 6.53936e-08, 5.81503e-08, 5.16763e-08, 4.58936e-08, 4.07318e-08, 3.61273e-08, 3.20225e-08,
   2.83657e-08, 2.51101e-08, 2.22136e-08, 1.96383e-08, 1.73501e-08, 1.53184e-08, 1.35157e-08, 1.19171e-08, 1.05006e-08, 9.24625e-09,
   8.13625e-09, 7.15467e-09, 6.28724e-09, 5.52121e-09, 4.84520e-09, 4.24904e-09, 3.72366e-09, 3.26099e-09, 2.85383e-09, 2.49576e-09,
   2.18110e-09, 1.90477e-09, 1.66227e-09, 1.44962e-09, 1.26328e-09, 1.10011e-09, 9.57322e-10, 8.32472e-10, 7.23383e-10, 6.28134e-10,
   5.45029e-10, 4.72575e-10, 4.09452e-10, 3.54500e-10, 3.06696e-10, 2.65142e-10, 2.29047e-10, 1.97718e-10, 1.70546e-10, 1.46997e-10,
   1.26604e-10, 1.08958e-10, 9.36992e-11, 8.05157e-11, 6.91340e-11, 5.93153e-11, 5.08516e-11, 4.35617e-11, 3.72876e-11, 3.18921e-11,
   2.72559e-11
};
 
// Ionosoheric Sigma calculation function coefficients from
// Juusola et al. 2025
// Note: MLT is in hours
const static Real c1p = 0.351;
const static Real c1h = 0.720;
const static Real c2p = 0.697;
const static Real c2h = 0.617;
const static Real c3p = 0.707;
const static Real c3h = 0.846;
const static Real c4P_values[] = {
   0.272, // 00
   0.212, // 01
   0.268, // 02
   0.357, // 03
   0.241, // 04
   0.225, // 05
   0.283, // 06
   0.621, // 07
   1.185, // 08
   1.643, // 09
   1.756, // 10
   1.562, // 11
   1.421, // 12
   1.327, // 13
   1.001, // 14
   0.959, // 15
   0.672, // 16
   0.300, // 17
   0.150, // 18
   0.241, // 19
   0.253, // 20
   0.340, // 21
   0.391, // 22
   0.413  // 23
};
const static Real c4H_values[] = {
   0.184, // 00
   0.248, // 01
   0.256, // 02
   0.343, // 03
   0.344, // 04
   0.341, // 05
   0.285, // 06
   0.676, // 07
   1.617, // 08
   2.610, // 09
   2.818, // 10
   2.962, // 11
   2.579, // 12
   2.049, // 13
   1.367, // 14
   0.701, // 15
   0.277, // 16
   0.141, // 17
   0.081, // 18
   0.172, // 19
   0.178, // 20
   0.223, // 21
   0.322, // 22
   0.373  // 23
};
const static Real c5P_values[] = {
   0.564, // 00
   0.602, // 01
   0.565, // 02
   0.507, // 03
   0.578, // 04
   0.598, // 05
   0.538, // 06
   0.338, // 07
   0.175, // 08
   0.075, // 09
   0.011, // 10
   0.000, // 11
   0.000, // 12
   0.037, // 13
   0.164, // 14
   0.222, // 15
   0.314, // 16
   0.470, // 17
   0.645, // 18
   0.581, // 19
   0.587, // 20
   0.533, // 21
   0.504, // 22
   0.498  // 23
};
const static Real c5H_values[] = {
   0.757, // 00
   0.702, // 01
   0.705, // 02
   0.649, // 03
   0.652, // 04
   0.662, // 05
   0.683, // 06
   0.453, // 07
   0.236, // 08
   0.115, // 09
   0.079, // 10
   0.042, // 11
   0.035, // 12
   0.030, // 13
   0.126, // 14
   0.306, // 15
   0.529, // 16
   0.653, // 17
   0.811, // 18
   0.710, // 19
   0.753, // 20
   0.721, // 21
   0.657, // 22
   0.635  // 23
};