fastem.cc

Go to the documentation of this file.

/* Copyright (C) 2002-2012
   Sreerekha Ravi<rekha@sat.physik.uni-bremen.de>
   Stefan Buehler <sbuehler@ltu.se>                  
 
   This program is free software; you can redistribute it and/or modify it
   under the terms of the GNU General Public License as published by the
   Free Software Foundation; either version 2, or (at your option) any
   later version.
 
   This program is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the Free Software
   Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307,
   USA. 
*/
 
/*===========================================================================
  === External declarations
  ===========================================================================*/
 
#include <cmath>
#include <stdexcept>
#include "arts_constants.h"
#include "arts_conversions.h"
#include "matpack_complex.h"
#include "exceptions.h"
#include "matpackI.h"
 
using std::ostringstream;
using std::runtime_error;
 
inline constexpr Numeric PI=Constant::pi;
inline constexpr Numeric DEG2RAD=Conversion::deg2rad(1);
inline constexpr Numeric RAD2DEG=Conversion::rad2deg(1);
 
#ifdef ENABLE_FASTEM
extern "C" {
#endif
 
void rttov_fastem5_(const Index& fastem_version,
                    const Numeric& frequency,
                    const Numeric& za,
                    const Numeric& temperature,
                    const Numeric& salinity,
                    const Numeric& wind_speed,
                    Numeric* emissivity,
                    Numeric* reflectivity,
                    const Numeric& transmittance,
                    const Numeric& rel_azimuth);
 
#ifdef ENABLE_FASTEM
}
#endif
 
// Define dummy function that throws a runtime error if ARTS is
// compiled without FASTEM support.
#ifndef ENABLE_FASTEM
void rttov_fastem5_(const Index&,
                    const Numeric&,
                    const Numeric&,
                    const Numeric&,
                    const Numeric&,
                    const Numeric&,
                    Numeric*,
                    Numeric*,
                    const Numeric&,
                    const Numeric&) {
  throw std::runtime_error(
      "This version of ARTS was compiled without FASTEM support.");
}
 
#endif
 
 
void fastem(  // Output:
    Vector& emissivity,
    Vector& reflectivity,
    // Input:
    const Numeric frequency,
    const Numeric za,
    const Numeric temperature,
    const Numeric salinity,
    const Numeric wind_speed,
    const Numeric transmittance,
    const Numeric rel_azimuth,
    const Index fastem_version) {
  emissivity.resize(4);
  reflectivity.resize(4);
 
  rttov_fastem5_(fastem_version,
                 frequency / 1e9,
                 180 - za,
                 temperature,
                 salinity * 1e3,
                 wind_speed,
                 emissivity.get_c_array(),
                 reflectivity.get_c_array(),
                 transmittance,
                 rel_azimuth);
}