rt_atmosphere.cc

Go to the documentation of this file.

#include "rt_atmosphere.h"

using namespace FullPhysics;
using namespace blitz;

#ifdef HAVE_LUA
#include "register_lua.h"
REGISTER_LUA_CLASS(RtAtmosphere)
REGISTER_LUA_END()
#endif

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

AccumulatedTimer RtAtmosphere::timer("RtAtmosphere calculations");

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

std::string RtAtmosphere::timer_info() const
{
  std::ostringstream os;
  os << timer;
  return os.str();
}


//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  
ArrayAd<double, 1> 
RtAtmosphere::optical_depth_wrt_state_vector(double wn, int spec_index) const
{
  firstIndex i1; secondIndex i2; thirdIndex i3;
  Range ra(Range::all());
  ArrayAd<double, 1> od(optical_depth_wrt_iv(wn, spec_index));
  Array<double, 3> ivjac(intermediate_variable(wn, spec_index).jacobian());
  ArrayAd<double, 1> res(od.rows(), ivjac.depth());
  res.value() = od.value();
  for(int i = 0; i < res.rows(); ++i)
    res.jacobian()(i, ra) = 
      sum(od.jacobian()(i, ra)(i2) * ivjac(i, ra, ra)(i2,i1),i2);
  return res;
}

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  
ArrayAd<double, 1> 
RtAtmosphere::single_scattering_albedo_wrt_state_vector
(double wn, int spec_index) const
{
  firstIndex i1; secondIndex i2; thirdIndex i3;
  Range ra(Range::all());
  ArrayAd<double, 1> ss(single_scattering_albedo_wrt_iv(wn, spec_index));
  Array<double, 3> ivjac(intermediate_variable(wn, spec_index).jacobian());
  ArrayAd<double, 1> res(ss.rows(), ivjac.depth());
  res.value() = ss.value();
  for(int i = 0; i < res.rows(); ++i)
    res.jacobian()(i, ra) = 
      sum(ss.jacobian()(i, ra)(i2) * ivjac(i, ra, ra)(i2,i1),i2);
  return res;
}

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

ArrayAd<double, 3>
RtAtmosphere::scattering_moment_wrt_state_vector
(double wn, int spec_index, int nummom, int numscat) const
{
  firstIndex i1; secondIndex i2; thirdIndex i3; fourthIndex i4; fifthIndex i5;
  ArrayAd<double, 3> 
    pf(scattering_moment_wrt_iv(wn, spec_index, nummom, numscat));
  Array<double, 3> ivjac(intermediate_variable(wn, spec_index).jacobian());
  ArrayAd<double, 3> res(pf.rows(), pf.cols(), pf.depth(), ivjac.depth());
  res.value() = pf.value();
  res.jacobian() = sum(pf.jacobian()(i1, i2, i3, i5) * ivjac(i2,i5, i4), i5);
  return res;
}