level_1b_cache.h

Go to the documentation of this file.

#ifndef LEVEL_1B_CACHE_H
#define LEVEL_1B_CACHE_H
#include "level_1b.h"
#include <vector>
namespace FullPhysics {
/****************************************************************/
class Level1bCache: public Level1b {
public:
  Level1bCache(const Level1b& L1_in);
  virtual ~Level1bCache() {}
  virtual void print(std::ostream& Os) const {Os << "Level1bCache";}
  virtual int number_spectrometer() const { return (int) lat.size(); }
  virtual DoubleWithUnit latitude(int i) const 
  { 
    range_check(i, 0, number_spectrometer());
    return lat[i];
  }

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

  void set_latitude(int i, const DoubleWithUnit& V)
  { 
    range_check(i, 0, number_spectrometer());
    lat[i] = V;
  }
  virtual DoubleWithUnit longitude(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return lon[i];
  }

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

  void set_longitude(int i, const DoubleWithUnit& V)
  { 
    range_check(i, 0, number_spectrometer());
    lon[i] = V;
  }
  virtual DoubleWithUnit sounding_zenith(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return szen[i];
  }

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

  void set_sounding_zenith(int i, const DoubleWithUnit& V)
  { 
    range_check(i, 0, number_spectrometer());
    szen[i] = V;
  }
  virtual DoubleWithUnit sounding_azimuth(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return sazm[i];
  }

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

  void set_sounding_azimuth(int i, const DoubleWithUnit& V)
  { 
    range_check(i, 0, number_spectrometer());
    sazm[i] = V;
  }
  virtual DoubleWithUnit solar_zenith(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return solzen[i];
  }

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

  void set_solar_zenith(int i, const DoubleWithUnit& V)
  { 
    range_check(i, 0, number_spectrometer());
    solzen[i] = V;
  }
  virtual DoubleWithUnit solar_azimuth(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return solazm[i];
  }

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

  void set_solar_azimuth(int i, const DoubleWithUnit& V)
  { 
    range_check(i, 0, number_spectrometer());
    solazm[i] = V;
  }
  virtual DoubleWithUnit altitude(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return alt[i];
  }

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

  void set_altitude(int i, const DoubleWithUnit& V) 
  { 
    range_check(i, 0, number_spectrometer());
    alt[i] = V;
  }
  virtual DoubleWithUnit relative_velocity(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return rvel[i];
  }

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

  void set_relative_velocity(int i, const DoubleWithUnit& V)
  { 
    range_check(i, 0, number_spectrometer());
    rvel[i] = V;
  }
  virtual blitz::Array<double, 1> stokes_coefficient(int i) const
  {
    range_check(i, 0, number_spectrometer());
    return stk_coeff[i];
  }    

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

  void set_stokes_coefficient(int i, const blitz::Array<double, 1>& V)
  { 
    range_check(i, 0, number_spectrometer());
    stk_coeff[i].reference(V.copy());
  }

  SpectralDomain sample_grid(int i) const
  {
    range_check(i, 0, number_spectrometer());
    return samp_grid[i];
  }

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

  void set_sample_grid(int i, const SpectralDomain& V)
  {
    range_check(i, 0, number_spectrometer());
    samp_grid[i] = V;
  }

  virtual Time time(int i) const
  {
    range_check(i, 0, number_spectrometer());
    return tm[i];
  }

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

  void set_time(int i, const Time& V)
  { 
    range_check(i, 0, number_spectrometer());
    tm[i] = V;
  }
  virtual SpectralRange radiance(int i) const
  { 
    range_check(i, 0, number_spectrometer());
    return rad[i];
  }

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

  void set_radiance(int i, const SpectralRange& V)
  { 
    range_check(i, 0, number_spectrometer());
    rad[i] = V;
  }

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

 void set_radiance(int i, const SpectralRange& V,
                    const std::vector<int>& Plist)
  { 
    range_check(i, 0, number_spectrometer());
    if(V.data().rows() != (int) Plist.size())
      throw Exception("Spectral Range is not the same size as Plist");
    blitz::Array<double, 1> rnew = rad[i].data().copy();
    blitz::Array<double, 1> unew = rad[i].uncertainty().copy();
    for(int j = 0; j < (int) Plist.size(); ++j) {
      rnew(Plist[j]) = V.data()(j);
      unew(Plist[j]) = V.uncertainty()(j);
    }
    rad[i] = SpectralRange(rnew, V.units(), unew);
  }

private:
  std::vector<DoubleWithUnit> lat, lon, szen, sazm, solzen, solazm, alt,
                                                    rvel;
  std::vector<blitz::Array<double, 1> > stk_coeff;
  std::vector<SpectralDomain> samp_grid;
  std::vector<Time> tm;
  std::vector<SpectralRange> rad;
};
}
#endif