SphericalHarmonicGravity.hpp

Go to the documentation of this file.

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

#ifndef GPSTK_SPHERICAL_HARMONIC_GRAVITY_HPP
#define GPSTK_SPHERICAL_HARMONIC_GRAVITY_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 "ForceModel.hpp"
#include "EarthSolidTide.hpp"
#include "EarthOceanTide.hpp"
#include "EarthPoleTide.hpp"

namespace gpstk
{
   using namespace StringUtils;



   class SphericalHarmonicGravity : public ForceModel
   {
   public:

      SphericalHarmonicGravity(int n, int m);


      virtual ~SphericalHarmonicGravity() {};


      virtual void initialize() = 0;

   
      Vector<double> gravity(Vector<double> r, Matrix<double> E);


      Matrix<double> gravityGradient(Vector<double> r, Matrix<double> E);
      

      virtual void doCompute(UTCTime utc, EarthBody& rb, Spacecraft& sc);


      SphericalHarmonicGravity& setDesiredDegree(const int& n, const int& m)
      { desiredDegree = n; desiredOrder = m; return (*this); }



      SphericalHarmonicGravity& enableSolidTide(bool b = true)
      { correctSolidTide = b; return (*this);}

      SphericalHarmonicGravity& enableOceanTide(bool b = true)
      { correctOceanTide = b; return (*this);}
      
      SphericalHarmonicGravity& enablePoleTide(bool b = true)
      { correctPoleTide = b; return (*this);}

      virtual std::string modelName() const
      {return "SphericalHarmonicGravity";}

      virtual int forceIndex() const
      { return FMI_GEOEARTH; }


      virtual void test();

   protected:

      void computeVW(Vector<double> r, Matrix<double> E);

      void correctCSTides(UTCTime t,bool solidFlag = false, bool oceanFlag = false, bool poleFlag = false);

      double normFactor(int n, int m);

   protected:

      struct GravityModelData
      {
         std::string modelName;

         double GM;
         double refDistance;

         bool includesPermTide;

         double refMJD;

         double dotC20;
         double dotC21;
         double dotS21;

         int maxDegree;
         int maxOrder;

         Matrix<double> unnormalizedCS;

      } gmData;

      Matrix<double> V, W;

      int desiredDegree, desiredOrder;
         
      bool   correctSolidTide;
      bool   correctPoleTide;
      bool   correctOceanTide;

      EarthSolidTide   solidTide;
      EarthPoleTide   poleTide;
      EarthOceanTide  oceanTide;


   }; // End of namespace 'gpstk'

      // @}

}  // End of namespace 'gpstk'

#endif   // GPSTK_SPHERICAL_HARMONIC_GRAVITY_HPP

---