StaticBInit.h

Go to the documentation of this file.

/*
 * (C) Copyright 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_RUNS_STATICBINIT_H_
#define OOPS_RUNS_STATICBINIT_H_

#include <string>

#include "eckit/config/LocalConfiguration.h"
#include "oops/base/instantiateCovarFactory.h"
#include "oops/base/ModelSpaceCovarianceBase.h"
#include "oops/base/Variables.h"
#include "oops/interface/Geometry.h"
#include "oops/interface/State.h"
#include "oops/runs/Application.h"

namespace oops {

  template <typename MODEL> class StaticBInit : public Application {
    typedef ModelSpaceCovarianceBase<MODEL>  Covariance_;
    typedef Geometry<MODEL>                  Geometry_;
    typedef State<MODEL>                     State_;

   public:
    // -----------------------------------------------------------------------------
    StaticBInit() {
      instantiateCovarFactory<MODEL>();
    }
    // -----------------------------------------------------------------------------
    virtual ~StaticBInit() {}
    // -----------------------------------------------------------------------------
    int execute(const eckit::Configuration & fullConfig) const {
      //  Setup resolution
      const eckit::LocalConfiguration resolConfig(fullConfig, "Geometry");
      const Geometry_ resol(resolConfig);

      //  Setup variables
      const eckit::LocalConfiguration varConfig(fullConfig, "Variables");
      const Variables vars(varConfig);

      //  Setup background state
      const eckit::LocalConfiguration bkgconf(fullConfig, "State");
      State_ xx(resol, vars, bkgconf);

      //  Initialize static B matrix
      const eckit::LocalConfiguration covarconf(fullConfig, "Covariance");
      boost::scoped_ptr< Covariance_ >
       Bmat(CovarianceFactory<MODEL>::create(covarconf, resol, vars, xx, xx));

      return 0;
    }
    // -----------------------------------------------------------------------------
   private:
    std::string appname() const {
      return "oops::StaticBInit<" + MODEL::name() + ">";
    }
    // -----------------------------------------------------------------------------
  };

}  // namespace oops

#endif  // OOPS_RUNS_STATICBINIT_H_