SolarRadiationPressure.hpp

Go to the documentation of this file.

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

#ifndef GPSTK_SOLAR_RADIATION_PRESSURE_HPP
#define GPSTK_SOLAR_RADIATION_PRESSURE_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 "ForceModel.hpp"


namespace gpstk
{

   class SolarRadiationPressure : public ForceModel
   {
   public:
      enum ShadowModel 
      {
         SM_CYLINDRICAL,         // cylindrical
         SM_CONICAL               // conical
      };

   public:
      SolarRadiationPressure() {}

      virtual ~SolarRadiationPressure() {}


      double getShadowFunction(Vector<double> r, 
                               Vector<double> r_Sun,
                               Vector<double> r_Moon, 
                               SolarRadiationPressure::ShadowModel sm = SM_CONICAL);

         // this is the real one
      virtual void doCompute(UTCTime t, EarthBody& bRef, Spacecraft& sc);

      virtual std::string modelName() const
      { return "SolarRadiationPressure"; }

      virtual int forceIndex() const
      { return FMI_SRP; }

      Vector<double> accelSRP(Vector<double> r, Vector<double> r_Sun);


      double partial_illumination(Vector<double> r, Vector<double> r_Sun );

    
   protected:      

      // some local variables

      double crossArea;
      double srpArea;

      double dryMass;

      double reflectCoeff;      // CR

      Vector<double> dadcr;

   }; // End of class 'SolarRadiationPressure'

      // @}

}  // End of namespace 'gsptk'

#endif   // GPSTK_SOLAR_RADIATION_PRESSURE_HPP

---