HarrisPriesterDrag.hpp

Go to the documentation of this file.

#pragma ident "$Id: HarrisPriesterDrag.hpp 2457 2010-08-18 14:20:12Z coandrei $"

#ifndef   GPSTK_HARRIS_PRIESTER_DRAG_HPP
#define   GPSTK_HARRIS_PRIESTER_DRAG_HPP


//============================================================================
//
//  This file is part of GPSTk, the GPS Toolkit.
//
//  The GPSTk is free software; you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published
//  by the Free Software Foundation; either version 2.1 of the License, or
//  any later version.
//
//  The GPSTk 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 Lesser General Public License for more details.
//
//  You should have received a copy of the GNU Lesser General Public
//  License along with GPSTk; if not, write to the Free Software Foundation,
//  Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
//  Wei Yan - Chinese Academy of Sciences . 2009, 2010
//
//============================================================================


#include "AtmosphericDrag.hpp"
#include "Matrix.hpp"
#include <map>

namespace gpstk
{

   class HarrisPriesterDrag : public AtmosphericDrag
   {
   public:
      HarrisPriesterDrag();

      virtual ~HarrisPriesterDrag() {};

      void test();

      void updateF107(double f107 = 157);


      virtual double computeDensity(UTCTime utc, EarthBody& rb, Vector<double> r, Vector<double> v);

   protected:

      void init();

      Matrix<double> getDensityCoeficentsByF107(double f107 = 157);

   protected:
      static const double upper_limit;        

      static const double lower_limit;
         
      static const double ra_lag;       
    
      double workingF107;

      Matrix<double> workingDens;

      std::map<int, Matrix<double> > hpCoeficient;
         
      static const double hpcoef[10][177];

   }; // End of class 'HarrisPriesterDrag'

      // @}

}  // End of namespace 'gpstk'

#endif   // GPSTK_HARRIS_PRIESTER_DRAG_HPP

---