GeoVaLsWriter.h

Go to the documentation of this file.

/*
 * (C) Copyright 2017-2018 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 "eckit/config/LocalConfiguration.h"
#include "oops/base/FilterBase.h"
#include "oops/interface/GeoVaLs.h"
#include "oops/interface/ObservationSpace.h"
#include "oops/interface/ObsVector.h"
#include "oops/util/dot_product.h"
#include "oops/util/Logger.h"

namespace oops {

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

template <typename MODEL>
class GeoVaLsWriter : public FilterBase<MODEL> {
  typedef GeoVaLs<MODEL>             GeoVaLs_;
  typedef ObservationSpace<MODEL>    ObsSpace_;
  typedef ObsVector<MODEL>           ObsVector_;

 public:
  explicit GeoVaLsWriter(const eckit::Configuration & conf) : conf_(conf) {}
  ~GeoVaLsWriter() {}

  void priorFilter(const ObsSpace_ &) const override {}
  void postFilter(const GeoVaLs_ & gv, const ObsVector_ &, const ObsSpace_ &) const override {
    const double zz = sqrt(dot_product(gv, gv));
    Log::debug() << "GeoVaLsWriter norm = " << zz << std::endl;
    gv.write(conf_);
  }

 private:
  const eckit::LocalConfiguration conf_;
  void print(std::ostream &) const override;
};

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

template <typename MODEL>
void GeoVaLsWriter<MODEL>::print(std::ostream & os) const {
  os << "GeoVaLsWriter " << conf_;
}

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

}  // namespace oops

#endif  // OOPS_BASE_GEOVALSWRITER_H_