GlobalTropModel.hpp

Go to the documentation of this file.

//==============================================================================
//
//  This file is part of GNSSTk, the ARL:UT GNSS Toolkit.
//
//  The GNSSTk 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 3.0 of the License, or
//  any later version.
//
//  The GNSSTk 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 GNSSTk; if not, write to the Free Software Foundation,
//  Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin.
//  Copyright 2004-2022, The Board of Regents of The University of Texas System
//
//==============================================================================
 
//==============================================================================
//
//  This software was developed by Applied Research Laboratories at the
//  University of Texas at Austin, under contract to an agency or agencies
//  within the U.S. Department of Defense. The U.S. Government retains all
//  rights to use, duplicate, distribute, disclose, or release this software.
//
//  Pursuant to DoD Directive 523024
//
//  DISTRIBUTION STATEMENT A: This software has been approved for public
//                            release, distribution is unlimited.
//
//==============================================================================
 
#ifndef GLOBAL_TROP_MODEL_HPP
#define GLOBAL_TROP_MODEL_HPP
 
#include "CommonTime.hpp"
#include "TropModel.hpp"
 
namespace gnsstk
{
 
   class GlobalTropModel : public TropModel
   {
   public:
      GlobalTropModel();
 
      GlobalTropModel(const double& ht, const double& lat, const double& lon,
                      const double& mjd)
      {
         validCoeff = validHeight = validLat = validLon = validDay = valid =
            false;
         setReceiverHeight(ht);
         setReceiverLatitude(lat);
         setReceiverLongitude(lon);
         setTime(mjd);
      }
 
      GlobalTropModel(const Position& RX, const CommonTime& time)
      {
         validCoeff = validHeight = validLat = validLon = validDay = valid = false;
         setReceiverHeight(RX.getAltitude());
         setReceiverLatitude(RX.getGeodeticLatitude());
         setReceiverLongitude(RX.getLongitude());
         setTime(time);
      }
 
      virtual std::string name()
      { return std::string("Global"); }
 
      virtual double correction(double elevation) const;
 
      virtual double correction(const Position& RX, const Position& SV);
 
      virtual double correction(const Position& RX,
                                const Position& SV,
                                const CommonTime& tt)
      {
         setTime(tt);
         return correction(RX,SV);
      }
 
      virtual double dry_zenith_delay() const;
 
      virtual double wet_zenith_delay() const;
 
      virtual double dry_mapping_function(double elevation) const;
 
      virtual double wet_mapping_function(double elevation) const;
 
      void getGPT(double& P, double& T, double& U);
 
      virtual void setWeather(const double& T, const double& P, const double& H)
      { setHumidity(H); }
 
      virtual void setWeather(const WxObservation& wx)
      { setHumidity(wx.humidity); }
 
      void setHumidity(const double& rh)
      {
         if(rh < 0.0 || rh > 100.)
            GNSSTK_THROW(InvalidParameter("Invalid humidity (%)"));
         humid = rh;
      }
 
      virtual void setReceiverLatitude(const double& lat);
 
      virtual void setReceiverLongitude(const double& lon);
 
      virtual void setReceiverHeight(const double& ht);
 
      virtual void setDayOfYear(const int& doy);
 
      virtual void setParameters(const CommonTime& time, const Position& rxPos);
 
      double getHeightLimit()
      { return HEIGHT_LIMIT; }
 
   private:
      void setTime(const double& mjd);
 
      void setTime(const CommonTime& time);
 
      // GMF
      static const double ADryMean[55];
      static const double BDryMean[55];
      static const double ADryAmp[55];
      static const double BDryAmp[55];
      static const double AWetMean[55];
      static const double BWetMean[55];
      static const double AWetAmp[55];
      static const double BWetAmp[55];
      // GPT
      static const double Ageoid[55];
      static const double Bgeoid[55];
      static const double APressMean[55];
      static const double BPressMean[55];
      static const double APressAmp[55];
      static const double BPressAmp[55];
      static const double ATempMean[55];
      static const double BTempMean[55];
      static const double ATempAmp[55];
      static const double BTempAmp[55];
 
      static const double Factorial[19];
 
      GNSSTK_EXPORT
      static const double HEIGHT_LIMIT;
 
      double height, latitude, longitude, dayfactor, undul;
      double P[10][10], aP[55], bP[55];
      bool validHeight, validLat, validLon, validDay, validCoeff;
 
      void updateGTMCoeff();
 
      void testValidity() const;
 
      void setValid()
      {
         try{
            valid = validHeight && validLat && validLon && validDay;
            if(valid && !validCoeff) {
               updateGTMCoeff();
               validCoeff = true;
               getGPT(press,temp,undul);
            }
         } catch(Exception& e) { GNSSTK_RETHROW(e); }
      }
 
   }; // end class GlobalTropModel
 
}  // end namespace
 
#endif   // #define GLOBAL_TROP_MODEL_HPP