pressure_fixed_level.h

Go to the documentation of this file.

#ifndef PRESSURE_FIXED_LEVEL_H
#define PRESSURE_FIXED_LEVEL_H
#include "state_vector.h"
#include "observer.h"
#include "fp_exception.h"
#include "pressure_level_input.h"
#include "pressure_imp_base.h"

namespace FullPhysics {
/****************************************************************/
class PressureFixedLevel : public PressureImpBase {
public:
  PressureFixedLevel(bool Pressure_flag, 
                     const boost::shared_ptr<PressureLevelInput>& Press_level,
                     double Surface_pressure)
    : press_level(Press_level)
    {
      blitz::Array<double, 1> coeff(1);
      blitz::Array<bool, 1> used(1);
      coeff(0) = Surface_pressure;
      used(0) = Pressure_flag;
      init(coeff, used);
      cov.resize(1, 1);
      cov(0,0) = 1;
    }
  virtual ~PressureFixedLevel() {}

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

  double surface_pressure_uncertainty() const 
  {return (cov(0, 0) < 0 ? 0.0 : sqrt(cov(0, 0))); }

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

  void set_surface_pressure(const AutoDerivative<double>& Surface_pressure) 
  {
    coeff(0) = Surface_pressure;
    cache_stale = true;
    const blitz::Array<double, 1>& plev = press_level->pressure_level();
    range_check(Surface_pressure.value(), plev(0), 
                plev(plev.rows() - 1));
    Observable<Pressure>::notify_update_do(*this);
  }

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

  int number_active_level() const 
  { return pressure_grid().value.rows(); }

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

  int number_active_layer() const { return number_active_level() - 1; }

  static const double new_level_fractional_size;

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

  int max_number_level() const {return press_level->pressure_level().rows();}

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

  blitz::Array<double, 1> pressure_active_levels() const
  { return press_level->pressure_level()
      (blitz::Range(0, number_active_level() - 1)); }

  virtual void print(std::ostream& Os) const;

  virtual boost::shared_ptr<Pressure> clone() const;

  virtual std::string sub_state_identifier() const { return "surface_pressure"; }

  virtual std::string state_vector_name_i(int i) const
  { return "Surface Pressure (Pascals)"; }
protected:
  virtual void calc_pressure_grid() const;
private:
  boost::shared_ptr<PressureLevelInput> press_level;
};
}
#endif