mwr_absorber_nh3.cpp

Go to the documentation of this file.

// canoe
#include <air_parcel.hpp>
#include <configure.hpp>
 
// climath
#include <climath/interpolation.h>
 
// application
#include <application/exceptions.hpp>
 
// opacity
#include "absorption_functions.hpp"
#include "mwr_absorbers.hpp"
 
namespace GiantPlanets {
 
void MwrAbsorberNH3::CheckFail() const {
  if (!HasPar("xHe")) {
    throw NotFoundError("MwrAbsorberNH3", "parameter 'xHe'");
  }
 
  if (!HasPar("power")) {
    throw NotFoundError("MwrAbsorberNH3", "parameter 'power'");
  }
}
 
Real MwrAbsorberNH3::GetAttenuation(Real wave1, Real wave2,
                                    AirParcel const& var) const {
  Real P = var.w[IPR] / 1.E5;      // pa -> bar
  Real P_idl = var.w[IPR] / 1.E5;  // pa -> bar
  Real T = var.w[IDN];
  Real xdry = 1.;
  for (int i = 1; i <= NVAPOR; ++i) xdry -= var.w[i];
 
  Real XHe = GetPar<Real>("xHe") * xdry;
  Real XH2 = xdry - XHe;
  Real XNH3 = var.w[GetSpeciesIndex(0)];
  Real XH2O = var.w[GetSpeciesIndex(1)];
 
  Real abs;
  Real wave = (wave1 + wave2) / 2.;
 
  if (model_name_ == "Bellotti16") {
    abs = attenuation_NH3_Bellotti(wave, P, P_idl, T, XH2, XHe, XNH3, XH2O);
  } else if (model_name_ == "BellottiSwitch16") {
    abs = attenuation_NH3_Bellotti_switch(wave, P, P_idl, T, XH2, XHe, XNH3,
                                          XH2O);
  } else if (model_name_ == "Devaraj") {
    abs = attenuation_NH3_Devaraj(wave, P, P_idl, T, XH2, XHe, XNH3, XH2O);
  } else if (model_name_ == "Radtran") {
    abs = attenuation_NH3_radtran(wave, P, T, XH2, XHe, XNH3);
  } else if (model_name_ == "Hanley09") {
    abs = attenuation_NH3_Hanley(wave, P, P_idl, T, XH2, XHe, XNH3, XH2O,
                                 GetPar<Real>("power"));
  } else {
    throw NotFoundError("MwrAbsorberNH3::GetAttenuation", model_name_);
  }
 
  return 100. * abs;  // 1/cm -> 1/m
}
 
}  // namespace GiantPlanets