aerosol_property_imp_base.h

Go to the documentation of this file.

#ifndef AEROSOL_PROPERTY_IMP_BASE_H
#define AEROSOL_PROPERTY_IMP_BASE_H

#include "aerosol_property.h"
#include "sub_state_vector_array.h"
#include <boost/lexical_cast.hpp>

namespace FullPhysics {
/****************************************************************/
class AerosolPropertyImpBase: public SubStateVectorArray<AerosolProperty> {
public:
  virtual ~AerosolPropertyImpBase() {}
  virtual boost::shared_ptr<AerosolProperty> clone() const = 0;
  virtual boost::shared_ptr<AerosolProperty> 
  clone(const boost::shared_ptr<Pressure>& Press,
        const boost::shared_ptr<RelativeHumidity>& Rh) const = 0;
  virtual ArrayAd<double,1> extinction_coefficient_each_layer(double wn) 
    const = 0;
  virtual ArrayAd<double, 1> scattering_coefficient_each_layer(double wn)
    const = 0;
  virtual ArrayAd<double, 3> 
  phase_function_moment_each_layer(double wn, int nmom = -1, 
                                   int nscatt = -1) const = 0;

  virtual std::string sub_state_identifier() const { return "aerosol/property"; }

  virtual std::string state_vector_name_i(int i) const 
  {
      return "Aerosol Property Coeff " + boost::lexical_cast<std::string>(i + 1);
  }

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

  virtual void print(std::ostream& Os) const { Os << desc(); }

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

  virtual std::string desc() const { return "AerosolPropertyImpBase"; }

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

  blitz::Array<double, 1> aerosol_parameter() const
  {
    return coefficient().value();
  }

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

  blitz::Array<double, 1> aerosol_parameter_uncertainty() const
  {
    blitz::Array<double, 1> uncert(coefficient().rows());
    for(int i = 0; i < sv_cov_sub.rows(); i++)
      uncert(i) = (sv_cov_sub(i,i) > 0 ? sqrt(sv_cov_sub(i, i)) : 0.0);
    return uncert;
  }

protected:
  // Don't think we need a cache, but if we end up needing it can add
  // this like we have in AerosolExtinctionImpBase.
  // mutable bool cache_stale;

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

  void init(const blitz::Array<double, 1>& Coeff, 
            const blitz::Array<bool, 1>& Used_flag)
  { 
    SubStateVectorArray<AerosolProperty>::init(Coeff, Used_flag);
  }

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

  AerosolPropertyImpBase() { }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------
  AerosolPropertyImpBase(const blitz::Array<double, 1>& Coeff, 
                         const blitz::Array<bool, 1>& Used_flag)
    : SubStateVectorArray<AerosolProperty>(Coeff, Used_flag)
      { }
};
}
#endif