LinearGetValues.h

Go to the documentation of this file.

/*
 * (C) Copyright 2020 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.
 */
 
#pragma once
 
#include <fstream>
#include <map>
#include <memory>
#include <ostream>
#include <string>
#include <vector>
 
#include "eckit/exception/Exceptions.h"
 
#include "oops/util/DateTime.h"
#include "oops/util/Logger.h"
#include "oops/util/ObjectCounter.h"
#include "oops/util/Printable.h"
#include "oops/util/Timer.h"
 
#include "ufo/GeoVaLs.h"
#include "ufo/Locations.h"
 
#include "mpasjedi/getvalues/LinearGetValues.interface.h"
 
// -----------------------------------------------------------------------------
 
namespace ufo {
  class GeoVaLs;
  class Locations;
}
 
namespace mpas {
  class GeometryMPAS;
  class IncrementMPAS;
  class StateMPAS;
  class VarChaModel2GeoVars;
  class LinVarChaModel2GeoVars;
 
// -----------------------------------------------------------------------------
 
class LinearGetValues : public util::Printable,
  private util::ObjectCounter<LinearGetValues> {
 public:
  static const std::string classname() {return "mpas::LinearGetValues";}
 
  LinearGetValues(const GeometryMPAS &, const ufo::Locations &,
                  const eckit::Configuration &);
  virtual ~LinearGetValues();
 
  void setTrajectory(const StateMPAS & state, const util::DateTime & t1,
                     const util::DateTime & t2, ufo::GeoVaLs & geovals);
  void fillGeoVaLsTL(const IncrementMPAS & inc, const util::DateTime & t1,
                     const util::DateTime & t2, ufo::GeoVaLs & geovals) const;
  void fillGeoVaLsAD(IncrementMPAS & inc, const util::DateTime & t1,
                     const util::DateTime & t2,
                     const ufo::GeoVaLs & geovals) const;
 
 private:
  const LinVarChaModel2GeoVars * getLinVarCha(const util::DateTime &) const;
 
  typedef std::map< util::DateTime, LinVarChaModel2GeoVars * >::iterator lvcIter;
  typedef std::map< util::DateTime, LinVarChaModel2GeoVars * >::const_iterator lvcIterCnst;
 
  void print(std::ostream &) const;
  F90lineargetvalues keyLinearGetValues_;
  ufo::Locations locs_;
  std::shared_ptr<const GeometryMPAS> geom_;
  std::map< util::DateTime, LinVarChaModel2GeoVars * > linearmodel2geovars_;
  std::unique_ptr<VarChaModel2GeoVars> model2geovars_;
};
 
// -----------------------------------------------------------------------------
 
}  // namespace mpas