GeoVaLsWriter.h

Go to the documentation of this file.

/*
 * (C) Copyright 2017-2019 UCAR
 *
 * This software is licensed under the terms of the Apache Licence Version 2.0
 * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
 */
 
#ifndef OOPS_BASE_GEOVALSWRITER_H_
#define OOPS_BASE_GEOVALSWRITER_H_
 
#include <memory>
 
#include "eckit/config/LocalConfiguration.h"
#include "oops/base/ObsFilterBase.h"
#include "oops/base/Variables.h"
#include "oops/interface/GeoVaLs.h"
#include "oops/interface/ObsDataVector.h"
#include "oops/interface/ObsDiagnostics.h"
#include "oops/interface/ObsSpace.h"
#include "oops/interface/ObsVector.h"
#include "oops/util/dot_product.h"
#include "oops/util/Logger.h"
 
namespace oops {
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
class GeoVaLsWriter : public ObsFilterBase<OBS> {
  typedef GeoVaLs<OBS>             GeoVaLs_;
  typedef ObsDiagnostics<OBS>      ObsDiags_;
  typedef ObsSpace<OBS>            ObsSpace_;
  typedef ObsVector<OBS>           ObsVector_;
  template <typename DATA> using ObsDataPtr_ = std::shared_ptr<ObsDataVector<OBS, DATA> >;
 
 public:
  GeoVaLsWriter(const ObsSpace_ &, const eckit::Configuration & conf,
                ObsDataPtr_<int>, ObsDataPtr_<float>): conf_(conf), novars_() {}
  ~GeoVaLsWriter() {}
 
  void preProcess() const override {}
 
  void priorFilter(const GeoVaLs_ & gv) const override {
    const double zz = sqrt(dot_product(gv, gv));
    Log::info() << "GeoVaLsWriter norm = " << zz << std::endl;
    gv.write(conf_);
  }
 
  void postFilter(const ObsVector_ &, const ObsDiags_ &) const override {}
 
  Variables requiredVars() const override {return novars_;};
  Variables requiredHdiagnostics() const override {return novars_;};
 
 private:
  const eckit::LocalConfiguration conf_;
  const Variables novars_;  // could be used to determine what needs saving
 
  void print(std::ostream &) const override;
};
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
void GeoVaLsWriter<OBS>::print(std::ostream & os) const {
  os << "GeoVaLsWriter: " << conf_;
}
 
// -----------------------------------------------------------------------------
 
}  // namespace oops
 
#endif  // OOPS_BASE_GEOVALSWRITER_H_