AstroTime.h

Go to the documentation of this file.

// $Id: AstroTime.h,v 1.21 2012/03/02 15:50:31 mpound Exp $

#ifndef CARMA_SERVICES_ASTROTIME_H
#define CARMA_SERVICES_ASTROTIME_H

#include "carma/util/Time.h"
#include "carma/services/IERSTable.h"
#include "carma/services/Location.h"
#include <memory>

namespace carma  {

  namespace services {

    class Angle;
    class HourAngle;
    class Length;

        class AstroTime {
        public:
           AstroTime();

           AstroTime(Location location);

          virtual ~AstroTime();

          void setSite(Angle longitude, Angle latitude, Length altitude);
           //__attribute__((deprecated));


          void setSite(double longitude, double latitude, double altitude);

          void setSite(Location location);

          Location getSite() const;

          double localSiderealTime(double mjd = carma::util::Time::MJD());

          double meanSiderealTime(double mjd = carma::util::Time::MJD());

          std::string lstString(double mjd = carma::util::Time::MJD(),
                                int precision = 1);

          double eqnEqx(double mjd = carma::util::Time::MJD());

          double ut1Utc(double mjd = carma::util::Time::MJD());

          double xPolar(double mjd = carma::util::Time::MJD());

          double yPolar(double mjd = carma::util::Time::MJD());

        
          double hourAngle(double mjd, double RA);

          HourAngle hourAngle(double mjd, Angle RA);


      double iersTableAge();

          //====================================================
          // STATIC METHODS - don't require Location
          //====================================================
          static double leap(const double mjd);

        static double elapsedJulCent( double mjd = carma::util::Time::MJD() );

        static double julianDay( double mjd = carma::util::Time::MJD() );

    static double modulo24(double hours);

        static const double JULIAN_DAY_ZERO ;

        static const double JULIAN_CENTURY ;

        static const double SECONDS_PER_DAY;

        //  Yes, I realize Time.h already has this, but
        //  I want the full PER_DAY set in one place.
        //  These could in principal be replaced by 
        //  calls to Units.convert(), which knows
        //  about sidereal vs. solar, etc. Must check
        //  the accuracy of its GnuUnits though.

        static const double MINUTES_PER_DAY;
        
        static const double HOURS_PER_DAY;

        static const double HOURS_PER_SIDEREAL_DAY;

        static const double SOLAR_DAY;

        static const double SIDEREAL_DAY;

        static const double SECONDS_PER_HOUR;
  
        static const double MINUTES_PER_HOUR ;

        static const double HOURS_PER_RADIAN;

        private:
        Location siteLocation_;    // Location holding site information 
    // Can't use auto_ptr because copy
    // construction releases ownership of
    // the contents.  Perhaps auto_ptr_ref
    // would work in this case, but I don't
    // really understand its usage.
        //std::auto_ptr<IERSTable> iers_;  
    IERSTable iers_;

        void checkIERSandLoadIfNecessary(void);  // helper function to load new IERS table
        void loadIERS(void);       // Re-read IERS Table from disk
    double lastIERSLoadTime_;  // MJD when the IERS table was last read from disk.
        };
    }
}

#endif //CARMA_SERVICES_ASTROTIME_H