linemixing_hitran.h

Go to the documentation of this file.

 
#ifndef LINEMIXING_HITRAN_H
#define LINEMIXING_HITRAN_H
 
#include "absorptionlines.h"
#include "arts_conversions.h"
#include "linescaling.h"
#include "matpack_data.h"
#include "matpack_complex.h"
#include "mystring.h"
 
 
struct HitranRelaxationMatrixData {
  Tensor4 W0pp, B0pp;
  Tensor4 W0rp, B0rp;
  Tensor4 W0qp, B0qp;
  
  Tensor4 W0pr, B0pr;
  Tensor4 W0rr, B0rr;
  Tensor4 W0qr, B0qr;
  
  Tensor4 W0pq, B0pq;
  Tensor4 W0rq, B0rq;
  Tensor4 W0qq, B0qq;
  HitranRelaxationMatrixData() = default;
  friend std::ostream& operator<<(std::ostream& os, const HitranRelaxationMatrixData& hit) {
    return os << hit.W0pp << '\n' << hit.B0pp << '\n'
              << hit.W0rp << '\n' << hit.B0rp << '\n'
              << hit.W0qp << '\n' << hit.B0qp << '\n'
              << hit.W0pr << '\n' << hit.B0pr << '\n'
              << hit.W0rr << '\n' << hit.B0rr << '\n'
              << hit.W0qr << '\n' << hit.B0qr << '\n'
              << hit.W0pq << '\n' << hit.B0pq << '\n'
              << hit.W0rq << '\n' << hit.B0rq << '\n'
              << hit.W0qq << '\n' << hit.B0qq << '\n';
  }
};
 
 
namespace lm_hitran_2017 {
enum class calctype {
  FullVP,
  FullRosenkranz,
  FullW,
  SDVP,
  SDRosenkranz,
  SDW,
  NoneVP,
  NoneRosenkranz,
  NoneW
};
 
Vector compute(const Numeric p,
               const Numeric t,
               const Numeric xco2,
               const Numeric xh2o,
               const Vector& invcm_grid,
               const Numeric stotmax,
               const calctype type=calctype::FullW);
 
Vector compute(const HitranRelaxationMatrixData& hitran,
               const ArrayOfAbsorptionLines& bands,
               const SpeciesIsotopologueRatios& isotopologue_ratio,
               const Numeric P,
               const Numeric T,
               const Vector& vmrs,
               const Vector& f_grid);
 
ENUMCLASS(ModeOfLineMixing, unsigned char,
  VP,  // Sets LineShape::VP, will not use LineMixing code; Sets ByLTE mode
  VP_Y,  // Sets LineShape::VP, will use LineMixing code with pressure > linemixinglimit;  Sets ByRosenkranzRelmatLTE mode
  SDVP,  // Sets LineShape::SDVP, will not use LineMixing code; Sets ByLTE mode
  SDVP_Y,  // Sets LineShape::SDVP, will use LineMixing code with pressure > linemixinglimit;  Sets ByHITRANRosenkranzRelmat mode
  FullW, // Sets LineShape::Lorentz, will use LineMixing code with pressure > linemixinglimit;  Sets ByHITRANFullRelmat mode
  VP_W  // Sets LineShape::Voigt, will use LineMixing code with pressure > linemixinglimit;  Sets ByHITRANFullRelmat mode
)
 
constexpr bool typeVP(ModeOfLineMixing x)
{
  return x == ModeOfLineMixing::VP or x == ModeOfLineMixing::VP_Y or x == ModeOfLineMixing::FullW or x == ModeOfLineMixing::VP_W;
}
 
constexpr bool typeLP(ModeOfLineMixing x)
{
  return x == ModeOfLineMixing::FullW;
}
 
constexpr bool typeFull(ModeOfLineMixing x)
{
  return x == ModeOfLineMixing::FullW or x == ModeOfLineMixing::VP_W;
}
 
void read(HitranRelaxationMatrixData& hitran,
          ArrayOfAbsorptionLines& bands,
          const SpeciesIsotopologueRatios& isotopologue_ratio,
          const String& basedir,
          const Numeric linemixinglimit,
          const Numeric fmin,
          const Numeric fmax,
          const Numeric stot,
          const ModeOfLineMixing mode);
 
void hitran_lm_eigenvalue_adaptation(AbsorptionLines& band,
                                     const Vector& temperatures,
                                     const HitranRelaxationMatrixData& hitran,
                                     const Numeric P0,
                                     const Index ord);
 
Tensor5 hitran_lm_eigenvalue_adaptation_test(const AbsorptionLines& band,
                                             const Vector& temperatures,
                                             const HitranRelaxationMatrixData& hitran,
                                             const Vector& pressures);
} // namespace lm_hitran_2017
 
#endif  // LINEMIXING_HITRAN_H