sc_atmJ70.h

//-----------------------------------------------------------------------------------------
//     J70 Atmosphere model header.
//-----------------------------------------------------------------------------------------
//     Copyright (c) 2002 Princeton Satellite Systems.  All rights reserved.
//-----------------------------------------------------------------------------------------
 
 
#ifndef __AtmJ70__
#define __AtmJ70__
 
#include <matrixlib/matrixlib.h>
 
class AtmJ70
{
                
        private:
                double  rho;            // Density (g/cm^3)
                double  nHe;            // Number density of helium
                double  nN2;            // Number density of nitrogen
                double  nO2;            // Number density of oxygen
                double  nAr;            // Number density of argon
                double  nO;             // Number density of monatomic oxygen
                double  tZ;             // Temperature
                double  eM;             // Mean molecular mass
                
                double  aP;             // Geomagnetic index 6.7 hours before the computation
                double  f;                      // Daily 10.7 cm solar flux (e-22 watts/m^2/cycle/sec)
                double  fHat;           // 81-day mean of f (e-22 watts/m^2/cycle/sec)
                double  fHat400;        // fHat 400 days before computation date
 
                // kSigma  = 0, 1 or 2 (nominal +2 sig -2 sig)
                // kTIming = 0, 1 or 2 (nominal, early, late)
                bool SolarFluxPredictions( double jD, int kTiming, int kSigma, const char fileName[] );
 
        public:
                                AtmJ70( void );
                                AtmJ70( double aP, double  f, double fHat, double fHat400 );
                                AtmJ70( double jD, int kTiming, int kSigma, const char fileName[] );
                bool    Initialize( double jD, int kTiming, int kSigma, const char fileName[] );
                void    Update( double jD, ml_matrix rECI );
                double  GetDensity( void ){ return this->rho; };
                double  GetNumberDensityOfHelium( void ){ return this->nHe; };
                double  GetNumberDensityOfArgon( void ){ return this->nAr; };
                double  GetNumberDensityOfNitrogen( void ){ return this->nN2; };
                double  GetNumberDensityOfOxygen( void ){ return this->nO2; };
                double  GetNumberDensityOfMonatomicOxygen( void ){ return this->nO; };
                double  GetTemperature( void ){ return this->tZ; };
                double  GetMeanMolecularMass( void ){ return this->eM*0.001; };
};
 
 
#endif