lidort_rt.cc

Go to the documentation of this file.

#include "lidort_rt.h"
#include "ostream_pad.h"

using namespace FullPhysics;
using namespace blitz;

#ifdef HAVE_LUA
#include "register_lua.h"
REGISTER_LUA_DERIVED_CLASS(LidortRt, RadiativeTransfer)
.def(luabind::constructor<const boost::shared_ptr<RtAtmosphere>&,
                          const boost::shared_ptr<StokesCoefficient>&,
                          const blitz::Array<double, 1>&,
                          const blitz::Array<double, 1>&, 
                          const blitz::Array<double, 1>&, 
                          bool, int, int, bool>())
.enum_("constants")
[
 luabind::value("maxmoments_input", Lidort_Pars::instance().maxmoments_input)
]
REGISTER_LUA_END()
#endif

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

LidortRt::LidortRt(const boost::shared_ptr<RtAtmosphere>& Atm,
                   const boost::shared_ptr<StokesCoefficient>& Stokes_coef,
                   const blitz::Array<double, 1>& Sza,
                   const blitz::Array<double, 1>& Zen,
                   const blitz::Array<double, 1>& Azm,
                   bool Pure_nadir,
                   int Number_streams,
                   int Number_moments,
                   bool Do_multi_scatt_only,
                   bool do_solar,
                   bool do_thermal)
: SpurrRt(Atm, Stokes_coef, Sza, Zen, Azm, do_solar, do_thermal)
{   
  rt_driver_.reset(new LidortRtDriver(Number_streams, Number_moments, Do_multi_scatt_only, surface_type(), Zen, Pure_nadir, do_solar, do_thermal));
}

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

void LidortRt::print(std::ostream& Os, bool Short_form) const
{
  Os << "LidortRt\n";
  OstreamPad opad1(Os, "  ");
  SpurrRt::print(opad1, Short_form);
  opad1.strict_sync();
  OstreamPad opad(Os, "    ");
  Os << "  Number stream: " << number_stream() << "\n"
     << "  Number moments: " << number_moment() << "\n"
     << "  Use first order scattering calc: " 
     << (rt_driver()->do_multi_scatt_only() ? "True\n" : "False\n")
     << (rt_driver()->pure_nadir() ? "True\n" : "False\n");

  opad << "\n";
  opad.strict_sync();
}