pressure.h

Go to the documentation of this file.

#ifndef PRESSURE_H
#define PRESSURE_H
#include "state_vector_observer.h"
#include "observer.h"
#include "array_ad_with_unit.h"

namespace FullPhysics {
/****************************************************************/
class Pressure : virtual public StateVectorObserver, 
                 public Observable<Pressure> {
public:
  virtual ~Pressure() {}
  virtual void add_observer(Observer<Pressure>& Obs) 
  { add_observer_do(Obs, *this);}
  virtual void remove_observer(Observer<Pressure>& Obs) 
  { remove_observer_do(Obs, *this);}

  AutoDerivativeWithUnit<double> surface_pressure() const;

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

  double surface_pressure_value() const 
  {return surface_pressure().convert(units::Pa).value.value();}

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

  virtual ArrayAdWithUnit<double, 1> pressure_grid() const = 0;

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

  int number_layer() const {return pressure_grid().value.rows() - 1;}

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

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

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

  virtual int max_number_level() const {return number_level();}

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

  virtual boost::shared_ptr<Pressure> clone() const = 0;
};
}
#endif