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 "complex.h"
#include "exceptions.h"
#include "matpackI.h"

using std::ostringstream;
using std::runtime_error;

extern const Numeric PI;
extern const Numeric DEG2RAD;
extern const Numeric RAD2DEG;

#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);
}