temperature_imp_base.h

Go to the documentation of this file.

#ifndef TEMPERATURE_IMP_BASE_H
#define TEMPERATURE_IMP_BASE_H
#include "temperature.h"
#include "sub_state_vector_array.h"
#include <boost/function.hpp>

namespace FullPhysics {
/****************************************************************/
class TemperatureImpBase: public SubStateVectorArray<Temperature> {
public:
  virtual ~TemperatureImpBase() {}
  virtual AutoDerivativeWithUnit<double> 
  temperature(const AutoDerivativeWithUnit<double>& Press) const
  { fill_cache(); 
    return AutoDerivativeWithUnit<double>(tgrid(Press.convert(units::Pa).value),
                                          units::K); 
  }
  virtual boost::shared_ptr<Temperature> clone() const = 0;
  virtual boost::shared_ptr<Temperature> 
  clone(const boost::shared_ptr<Pressure>& Press) const = 0;
  
  virtual void update_sub_state_hook() 
  { cache_stale = true; }
  
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  virtual void print(std::ostream& Os) const { Os << desc(); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  virtual std::string desc() const { return "TemperatureImpBase"; }
protected:
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  mutable bool cache_stale;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  mutable boost::function<AutoDerivative<double>(AutoDerivative<double>)> tgrid;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  virtual void calc_temperature_grid() const = 0;

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  void init(const blitz::Array<double, 1>& Coeff, 
            const blitz::Array<bool, 1>& Used_flag,
            const boost::shared_ptr<Pressure>& Press,
            bool Mark_according_to_press = true,
            int Pdep_start = 0)
  { SubStateVectorArray<Temperature>::init(Coeff, Used_flag, Press,
                                           Mark_according_to_press,
                                           Pdep_start);
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  TemperatureImpBase() : cache_stale(true) { }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  TemperatureImpBase(const blitz::Array<double, 1>& Coeff, 
                     const blitz::Array<bool, 1>& Used_flag,
                     const boost::shared_ptr<Pressure>& Press,
                     bool Mark_according_to_press = true,
                     int Pdep_start = 0)
    : SubStateVectorArray<Temperature>(Coeff, Used_flag, Press,
                                       Mark_according_to_press, Pdep_start),
      cache_stale(true) { }
private:
  void fill_cache() const
  {
    if(cache_stale)
      calc_temperature_grid();
    cache_stale = false;
  }

};
}
#endif