EarthOceanTide.hpp

Go to the documentation of this file.

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

#ifndef GPSTK_OCEAN_TIDE_HPP
#define GPSTK_OCEAN_TIDE_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 <string>


namespace gpstk
{

   class EarthOceanTide
   {
   public:
      EarthOceanTide()
         : minX(0.05), 
           maxN(4), 
           isLoaded(false)
      {
         fileName = "InputData\\Earth\\OT_CSRC.TID";
      }

      ~EarthOceanTide(){}

     
      struct CSR_OTIDE
      {
         double   KNMP[20];
         int      NTACT;
         int      NDOD[1200][6];
         double   CSPM[1200][4];
         int      NM[1200][2];  
      };

         

      void loadTideFile(std::string fileName, int NMAX = 4, double XMIN = 0.05);

      void getOceanTide(double mjdUtc, double dC[], double dS[] );

      void setTideFile(std::string file)
      {
         fileName = file;
         isLoaded = false;
      }

      void test();

   protected:
      //Get Doson  FUNDAMENTAL ARGUMENTS
      //void NUTARG(double mjdUtc, double BETA[6]);   

   protected:
      std::string fileName;

      int maxN;
      double minX;

      bool isLoaded;

      double FAC[41];

      int    NWAV;   // number of lines to skip
      int    NTOT;   // number of data lines
      int    NMX;      // max degree
      int    MMX;     // max order 
      double RRE;
      double RHOW;
      double XME;
      double PFCN;
      double XXX;

      CSR_OTIDE  tideData;

      static const double G;  //= 6.67259e-11;
      static const double GE; //= 9.780327;

   }; // End of class 'EarthOceanTide'

      // @}

}  // End of namespace 'gpstk'


#endif   // GPSTK_OCEAN_TIDE_HPP

---