hitran_absorber.hpp

Go to the documentation of this file.

#ifndef SRC_OPACITY_HITRAN_ABSORBER_HPP_
#define SRC_OPACITY_HITRAN_ABSORBER_HPP_
 
// C/C++
#include <string>
#include <vector>
 
// harp
#include "absorber.hpp"
 
class HitranAbsorber : public Absorber {
  friend std::ostream& operator<<(std::ostream& os, HitranAbsorber const& ab);
  // friend HitranAbsorber const& MakeCKAbsorber<>(HitranAbsorber const& albl,
  //   int const *ck_axis, Real const *ck_wave, int nbins);
 public:
  HitranAbsorber(std::string name) : Absorber(name) {}
 
  virtual ~HitranAbsorber() {}
  void LoadCoefficient(std::string fname, size_t bid = 0) override;
  Real GetAttenuation(Real wave1, Real wave2,
                      AirParcel const& var) const override;
 
 protected:
  size_t len_[3];            
  std::vector<Real> axis_;   
  std::vector<Real> kcoeff_; 
  AthenaArray<Real> refatm_; 
  Real getRefTemp(Real pres) const;
};
 
#endif  // SRC_HARP_HITRAN_ABSORBER_HPP_