absorption.cc

Go to the documentation of this file.

 
#include "absorption.h"
#include "arts.h"
#include "arts_conversions.h"
#include "auto_md.h"
#include "file.h"
#include "linescaling.h"
#include "lineshape.h"
#include "logic.h"
#include "math_funcs.h"
#include "messages.h"
#include "partfun.h"
#include <algorithm>
#include <cfloat>
#include <cmath>
#include <cstdlib>
#include <map>
 
#include "global_data.h"
 
void checkPartitionFunctions(const ArrayOfArrayOfAbsorptionLines& abs_lines_per_species) {
  for (auto& abs_lines: abs_lines_per_species) {
    for (auto& band: abs_lines) {
      ARTS_USER_ERROR_IF (not PartitionFunctions::has_partfun(band.Isotopologue()),
                          "Species: ", band.Isotopologue().FullName(), " has no partition function\n",
                          "You must recompile ARTS partition functions with data for this species to continue your calculations,\n"
                          "or exclude the species from your computation setup")
    }
  }
}
 
void checkIsotopologueRatios(const ArrayOfArrayOfAbsorptionLines& abs_lines_per_species,
                             const Species::IsotopologueRatios& isoratios) {
  
  for (auto& abs_lines: abs_lines_per_species) {
    for (auto& band: abs_lines) {
      ARTS_USER_ERROR_IF (std::isnan(isoratios[band.Isotopologue()]),
                          "Species: ", band.Isotopologue().FullName(), " has no isotopologue ratios\n",
                          "You must add its isotopologue ratios to your included data or\n"
                          "exclude the species from your computation setup")
    }
  }
}
 
Numeric wavenumber_to_joule(Numeric e) {
  return Conversion::kaycm2joule(e);
}
 
 
void set_vmr_from_first_species(Vector& vmr,
                                const String& species_name,
                                const ArrayOfArrayOfSpeciesTag& abs_species,
                                const Matrix& abs_vmrs) {
  const Index index = find_first_species(abs_species, Species::fromShortName(species_name));
 
  vmr.resize(abs_vmrs.ncols());
  if (index < 0)
    vmr = -99;
  else
    vmr = abs_vmrs(index, Range(joker));
}