pressure_level_input.h

Go to the documentation of this file.

#ifndef PRESSURE_LEVEL_INPUT_H
#define PRESSURE_LEVEL_INPUT_H
#include "printable.h"
#include "hdf_file.h"
#include <blitz/array.h>

namespace FullPhysics {
/****************************************************************/
class PressureLevelInput : public Printable<PressureLevelInput> {
public:
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  PressureLevelInput(const blitz::Array<double, 1>& Press_level)
    : press_level(Press_level.copy()) {}

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

  PressureLevelInput(const HdfFile& Hdf_file, 
                     const std::string& Hdf_group = "Pressure")
    : press_level(Hdf_file.read_field<double, 1>(Hdf_group + "/Pressure")) {}

  virtual ~PressureLevelInput() {}

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

  const blitz::Array<double, 1>& pressure_level() const {return press_level; }

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

  void print(std::ostream& Os) const;
private:
  blitz::Array<double, 1> press_level;
};
}
#endif