MoonPosition.hpp

Go to the documentation of this file.

#pragma ident "$Id: MoonPosition.hpp 2741 2011-06-22 16:37:02Z nwu $"

#ifndef MOONPOSITION_HPP
#define MOONPOSITION_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
//
//  Dagoberto Salazar - gAGE ( http://www.gage.es ). 2007
//
//============================================================================


#include <cmath>
#include <string>

#include "DayTime.hpp"
#include "Triple.hpp"
#include "icd_200_constants.hpp"
#include "AstronomicalFunctions.hpp"
#include "Vector.hpp"
#include "Matrix.hpp"


namespace gpstk
{


   class MoonPosition
   {
   public:

      MoonPosition() throw() {}

      virtual ~MoonPosition() {}


      Triple getPosition(const DayTime& t) const
         throw(InvalidRequest);


      Triple getPositionCIS(const DayTime& t) const
         throw(InvalidRequest);


      DayTime getInitialTime() const throw(InvalidRequest)
      { return initialTime; }


      DayTime getFinalTime() const throw(InvalidRequest)
      { return finalTime; }


   private:

      static const DayTime initialTime;

      static const DayTime finalTime;

         // Coefficients for fundamental arguments
         // Units are degrees for position and Julian centuries for time

      static const double ELP0, ELP1, ELP2, ELP3;

      static const double EM0, EM1, EM2, EM3;

      static const double EMP0, EMP1, EMP2, EMP3;

      static const double D0, D1, D2, D3;

      static const double F0, F1, F2, F3;

      static const double OM0, OM1, OM2, OM3;

      static const double E1, E2;

      static const double PAC, PA0, PA1;
      static const double PBC;
      static const double PCC;
      static const double PDC;
      static const double PEC, PE0, PE1, PE2;
      static const double PFC;
      static const double PGC;
      static const double PHC;
      static const double cPIC;
      static const double PJC, PJ0, PJ1;
      static const double CW1;
      static const double CW2;

         // Coefficients for Moon position
         //      Tx(N): coefficient of L, B or P term (deg)
         //      ITx(N,0-4): coefficients of M, M', D, F, E**n in argument
         //
      static const size_t NL, NB, NP;
      static Vector<double> TL, TB, TP;
      static Matrix<int> ITL, ITB, ITP;

   }; // end class MoonPosition



} // namespace gpstk
#endif  // MOONPOSITION_HPP

---