ils_table.h

Go to the documentation of this file.

#ifndef ILS_TABLE_H
#define ILS_TABLE_H
#include "ils_function.h"
#include "hdf_file.h"
#include "linear_interpolate.h"
#include "log_interpolate.h"

namespace FullPhysics {
/****************************************************************/
template <class Lint>
class IlsTable : public IlsFunction {
public:
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  IlsTable(const blitz::Array<double, 1>& Wavenumber, 
           const blitz::Array<double, 2>& Delta_lambda, 
           const blitz::Array<double, 2>& Response,
           const std::string& Band_name, const std::string& Hdf_band_name,
           bool Interpolate_wavenumber = false)
    : band_name_(Band_name), 
      hdf_band_name_(Hdf_band_name), 
      interpolate_wavenumber(Interpolate_wavenumber),
      from_hdf_file(false)
  { 
    create_delta_lambda_to_response(Wavenumber, Delta_lambda, Response); 
  }
  IlsTable(const HdfFile& Hdf_static_input, int Spec_index, 
           const std::string& Band_name, const std::string& Hdf_band_name,
           const std::string& Hdf_group = "Instrument/ILS");

  virtual ~IlsTable() {}
  virtual void ils
  (const AutoDerivative<double>& wn_center,
   const blitz::Array<double, 1>& wn,
   ArrayAd<double, 1>& res) const;

  virtual void print(std::ostream& Os) const;
  virtual std::string band_name() const {return band_name_;}
  virtual std::string hdf_band_name() const { return hdf_band_name_;}
  
  virtual blitz::Array<double, 1> wavenumber() const { return wavenumber_; }
  virtual blitz::Array<double, 2> delta_lambda() const { return delta_lambda_; }
  virtual blitz::Array<double, 2> response() const { return response_; }
  
  void create_delta_lambda_to_response(const blitz::Array<double, 1>& Wavenumber, 
                                       const blitz::Array<double, 2>& Delta_lambda, 
                                       const blitz::Array<double, 2>& Response);
private:
  // Save some typing by giving these typedefs
  typedef ArrayAd<double, 1> arrad;
  typedef AutoDerivative<double> ad;
  typedef typename std::map<double, Lint>::const_iterator it;

  std::map<double, Lint> delta_lambda_to_response;

  // Store the input arrays so we can re-build the ils map for testing
  // purposes, or for whatever reason you would need to do this mid program
  blitz::Array<double, 1> wavenumber_;
  blitz::Array<double, 2> delta_lambda_;
  blitz::Array<double, 2> response_;

  std::string band_name_, hdf_band_name_;
  bool interpolate_wavenumber;
  
  // Stash some information for print out with print
  bool from_hdf_file;
  std::string hdf_file_name;
  std::string hdf_group;
};

typedef IlsTable<LinearInterpolate<AutoDerivative<double>, AutoDerivative<double> > > IlsTableLinear;
typedef IlsTable<LinearLogInterpolate<AutoDerivative<double>, AutoDerivative<double> > > IlsTableLog;

}
#endif