chapman_boa_rt.h

Go to the documentation of this file.

#ifndef CHAPMAN_BOA_RT_H
#define CHAPMAN_BOA_RT_H

#include "radiative_transfer.h"
#include "atmosphere_oco.h"
#include "chapman_boa.h"
#include "spectral_bound.h"

namespace FullPhysics {

/****************************************************************/
class ChapmanBoaRT : public RadiativeTransfer,
                     public Observer<RtAtmosphere>,
                     public boost::noncopyable  {
public:
  ChapmanBoaRT(const boost::shared_ptr<AtmosphereOco>& Atm,
               const blitz::Array<double, 1>& Sza);

  ChapmanBoaRT(const boost::shared_ptr<AtmosphereOco>& Atm,
               const blitz::Array<double, 1>& Sza, 
               const SpectralBound& Spec_bound);

  virtual ~ChapmanBoaRT() {}

  //-----------------------------------------------------------------------
  //-----------------------------------------------------------------------
  void notify_update(const RtAtmosphere& updated_atm) { chapman_cache_stale = true; }

  virtual int number_stokes() const { return 1; }

  virtual int number_spectrometer() const
  { return sza.extent(blitz::firstDim);}

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

  const boost::shared_ptr<AtmosphereOco>& atmosphere_ptr() const {return atm;}

  // See description in base class
  virtual Spectrum reflectance
  (const SpectralDomain& Spec_domain, int Spec_index, 
   bool Skip_jacobian = false) const;
  virtual blitz::Array<double, 2> stokes(const SpectralDomain& Spec_domain, int Spec_index) const;
  virtual ArrayAd<double, 2> stokes_and_jacobian(const SpectralDomain& Spec_domain, int Spec_index) const;
 
  virtual void print(std::ostream& Os, bool Short_form = false) const;

private:

  void compute_chapman_factors(const int Spec_idx) const;

  mutable blitz::Array<bool, 1> chapman_cache_stale;

  mutable std::vector<boost::shared_ptr<ChapmanBOA> > chapman_boa;

  SpectralBound spec_bound;

  boost::shared_ptr<AtmosphereOco> atm;

  blitz::Array<double, 1> sza;

};
}

#endif