#pragma ident "$Id: MoonForce.hpp 3140 2012-06-18 15:03:02Z susancummins $"

/**
 * @file MoonForce.hpp
 * This class calculate the gravitational effect of the moon.
 */

#ifndef GPSTK_MOON_FORCE_HPP
#define GPSTK_MOON_FORCE_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., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//
//  Wei Yan - Chinese Academy of Sciences . 2009, 2010
//
//============================================================================

#include "ForceModel.hpp"

namespace gpstk
{
      /** @addtogroup GeoDynamics */
      //@{

      /**
       * This class calculate the gravitational effect of the moon.
       * 
       *
       *  This Model is checked on Sep 25th,2009, OK!!!
       * 
       */
   class MoonForce : public ForceModel
   {
   public:
         /// Default constructor
      MoonForce();
         
         /// Default destructor
      virtual ~MoonForce(){};

         /** Call the relevant methods to compute the acceleration.
          * @param utc Time reference class
          * @param rb  Reference Body class
          * @param sc  Spacecraft parameters and state
          * @return the acceleration [m/s^s]
          */
      virtual void doCompute(UTCTime utc, EarthBody& rb, Spacecraft& sc);
      
         /// Return force model name
      virtual std::string modelName() const
      { return "MoonForce"; }

         /// return the force model index
      virtual int forceIndex() const
      { return FMI_GEOMOON; }

   protected:

         /// Gravitational constant for the body. [m^3/s^2]
      double mu;


   }; // End of class 'MoonForce'

      // @}

}  // End of namespace 'gpstk'

#endif   // GPSTK_MOON_FORCE_HPP




