CENamespace.h

Go to the documentation of this file.

/***************************************************************************
 *  CENamespace.h: CppEphem                                                *
 * ----------------------------------------------------------------------- *
 *  Copyright © 2016-2019 JCardenzana                                      *
 * ----------------------------------------------------------------------- *
 *                                                                         *
 *  This program is free software: you can redistribute it and/or modify   *
 *  it under the terms of the GNU General Public License as published by   *
 *  the Free Software Foundation, either version 3 of the License, or      *
 *  (at your option) any later version.                                    *
 *                                                                         *
 *  This program 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 General Public License for more details.                           *
 *                                                                         *
 *  You should have received a copy of the GNU General Public License      *
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.  *
 *                                                                         *
 ***************************************************************************/
 
#ifndef CENamespace_h
#define CENamespace_h
 
#include <cmath>
#include <map>
#include <sofa.h>
#include <memory>
#include "CECorrections.h"
 
 
// Create a namespace with some useful stuff
namespace CppEphem {
    // Define an angle type so that we can differentiate between degrees or radians
 
    // Approximate temperature at sea-level in Kelvin
    inline double SeaLevelTemp_K() {return 288.2 ;}
    inline double SeaLevelTemp_C() {return SeaLevelTemp_K()-273.15 ;}
    inline double SeaLevelTemp_F() {return (1.8*SeaLevelTemp_C())+32.0 ;}
 
    // Temperature conversion methods
    inline double Temp_C2F(const double& temp_C) {return (1.8 * temp_C) + 32.0;}
    inline double Temp_C2K(const double& temp_C) {return temp_C+273.15;}
    inline double Temp_F2C(const double& temp_F) {return (temp_F - 32.0)/1.8;}
    inline double Temp_F2K(const double& temp_F) {return Temp_C2K(Temp_F2C(temp_F));}
    inline double Temp_K2C(const double& temp_K) {return temp_K - 273.15;}
    inline double Temp_K2F(const double& temp_K) {return 1.8*temp_K-459.67;}
    
    // Some constants taken from the sofa dictionary, but with some more discriptive names
    inline double julian_date_J2000() {return DJ00 ;}           
    inline double c() {return DC * DAU / DAYSEC ;}              
    inline double c_au_per_day() {return DC ;}                  
    inline double m_per_au() {return DAU ;}                     
    inline double sec_per_day() {return DAYSEC;}                
 
    /*********************************************
     * Return the 'dut1' value which represents 
     * 'UT1-UTC' for a given MJD or it's error
     *********************************************/
    double dut1(const double& mjd) ;
    double dut1Error(const double& mjd=51544.5) ;
    double dut1Calc(const double& mjd) ;
    
    /*********************************************
     * Polar Motion methods
     * For the moment, 0.0 should be sufficient. This
     * assumption may need to be revisited
     *********************************************/
    double xp(const double& mjd);
    double yp(const double& mjd);
 
    /*********************************************
     * Earth longitude & obliquity correction
     *********************************************/
    double deps(const double& mjd);
    double dpsi(const double& mjd);
    
    /*********************************************
     * TT-UT1 correction
     *********************************************/
    double ttut1(const double& mjd);
 
    inline double EstimateAltitude_m(double pressure_hPa)
        {return -29.3 * SeaLevelTemp_K() * std::log(pressure_hPa/1013.25) ;}
    inline double EstimatePressure_hPa(double elevation_m)
        {return 1013.25 * std::exp(-elevation_m / (29.3*SeaLevelTemp_K() )) ;}
 
    // Methods for getting the corrections values
    std::string NutationFile(void);
    std::string TtUt1HistFile(void);
    std::string TtUt1PredFile(void);
    void        SetNutationFile(const std::string& filename);
    void        SetTtUt1HistFile(const std::string& filename);
    void        SetTtUt1PredFile(const std::string& filename);
    void        CorrectionsInterp(bool set_interp);
    static      CECorrections corrections;
 
    namespace StrOpt {
        // Method for splitting a string based on some delimiter into a vector of strings
        void                     split(const std::string&        s, 
                                       const char&               delim, 
                                       std::vector<std::string>* elems);
        std::vector<std::string> split(const std::string& s, 
                                       const char&        delim);
 
        // Join strings using a set delimiter
        template<typename T>
        std::string join(const std::vector<T>& values,
                         const char&           delim);
        std::string join_angle(const std::vector<double>& values,
                               const char&                delim);
    }
}
 
#endif /* CENamespace_h */