stokes_coefficient_fraction.h

Go to the documentation of this file.

#ifndef STOKES_COEFFICIENT_FRACTION_H
#define STOKES_COEFFICIENT_FRACTION_H

#include "stokes_coefficient_imp_base.h"
#include <boost/lexical_cast.hpp>

namespace FullPhysics {
/****************************************************************/
class StokesCoefficientFraction : public StokesCoefficientImpBase {
public:
  StokesCoefficientFraction(const blitz::Array<double, 2>& Stokes_coeff_parallel,
                            const blitz::Array<double, 1>& Coeffs,
                            const blitz::Array<bool, 1>& Flag);
  virtual ~StokesCoefficientFraction() {}
  virtual void print(std::ostream& Os) const;
  virtual boost::shared_ptr<StokesCoefficient> clone() const;
  virtual std::string sub_state_identifier() const { return "parallel_polarization_fraction"; } 
  virtual std::string state_vector_name_i(int i) const
  { return "Parallel Polarization Band " + 
      boost::lexical_cast<std::string>(i + 1); }

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

  double f_uncertainty(int Spec_index) const 
  { 
    if(sv_cov_sub.rows() < 1)
      return 0;
    double t = sv_cov_sub(Spec_index,Spec_index);
    return (t < 0 ? 0 : sqrt(t)); 
  }
  
protected:
  virtual void calc_stokes_coeff() const;
private:
  blitz::Array<double, 2> stokes_coeff_parallel;
};
}
#endif