twostream_rt.cc

Go to the documentation of this file.

#include "twostream_rt.h"
#include "ostream_pad.h"

using namespace FullPhysics;
using namespace blitz;

#ifdef HAVE_LUA
#include "register_lua.h"
REGISTER_LUA_DERIVED_CLASS(TwostreamRt, 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>&>())
.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>())
REGISTER_LUA_END()
#endif

// For debugging purposes, it can be useful to dump out the input
// used by this class. We'll leave this code in place, in case we
// need it again, but normally this turned off.
const bool dump_data = false;

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

TwostreamRt::TwostreamRt(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 do_fullquadrature, bool do_solar, bool do_thermal)
: SpurrRt(Atm, Stokes_coef, Sza, Zen, Azm, do_solar, do_thermal)
{   
  rt_driver_.reset(new TwostreamRtDriver(atm->number_layer(), surface_type(), do_fullquadrature, do_solar, do_thermal));
  if(dump_data)
    std::cout << "# Nlayer:\n" << atm->number_layer() << "\n"
              << "# Surface type:\n" << surface_type() << "\n"
              << "# do_fullquadrature:\n" << do_fullquadrature << "\n";
}

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

void TwostreamRt::print(std::ostream& Os, bool Short_form) const
{
  Os << "TwostreamRt\n";
  OstreamPad opad1(Os, "  ");
  SpurrRt::print(opad1, Short_form);
  opad1.strict_sync();
  Os << "do_full_quadrature = " << (rt_driver()->do_full_quadrature() ? "true" : "false") << "\n";
  opad1.strict_sync();
}