absorber_vmr_fixed_level.h

Go to the documentation of this file.

#ifndef ABSORBER_VMR_FIXED_LEVEL_H
#define ABSORBER_VMR_FIXED_LEVEL_H

#include "absorber_vmr_imp_base.h"
#include "pressure.h"
#include "pressure_level_input.h"
#include <boost/lexical_cast.hpp>

namespace FullPhysics {
/****************************************************************/
class AbsorberVmrFixedLevel : public AbsorberVmrImpBase {
public:
  AbsorberVmrFixedLevel(const boost::shared_ptr<Pressure>& Press,
        const boost::shared_ptr<PressureLevelInput>& Press_level,          
        const blitz::Array<bool, 1>& Used_flag, 
        const blitz::Array<double, 1>& Vmr,
        const std::string& Gas_name);
  virtual ~AbsorberVmrFixedLevel() {}
  virtual void print(std::ostream& Os) const;

  virtual std::string sub_state_identifier() const { return "absorber_levels/" + gas_name(); }

  virtual std::string state_vector_name_i(int i) const
  { return gas_name() + " VMR for Press Lvl " + 
      boost::lexical_cast<std::string>(i + 1); }
  virtual boost::shared_ptr<AbsorberVmr> clone() const;
  virtual boost::shared_ptr<AbsorberVmr> 
  clone(const boost::shared_ptr<Pressure>& Press) const;

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

  blitz::Array<double, 1> volume_mixing_ratio_level() const 
  {return coeff.value();}

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

  blitz::Array<double, 1> volume_mixing_ratio_active_level() const 
  {return coeff.value()(blitz::Range(0, press->number_level() - 1));}
protected:
  virtual void calc_vmr() const;
private:
  boost::shared_ptr<PressureLevelInput> press_level;
};
}
#endif