src/lorenz95/ObservationL95.cc

Go to the documentation of this file.

/*
 * (C) Copyright 2009-2016 ECMWF.
 * 
 * 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. 
 * In applying this licence, ECMWF does not waive the privileges and immunities 
 * granted to it by virtue of its status as an intergovernmental organisation nor
 * does it submit to any jurisdiction.
 */

#include "lorenz95/ObservationL95.h"

#include <string>
#include <vector>

#include "eckit/config/Configuration.h"
#include "lorenz95/GomL95.h"
#include "lorenz95/ObsBias.h"
#include "lorenz95/ObsVec1D.h"
#include "oops/base/Variables.h"
#include "oops/util/DateTime.h"
#include "oops/util/Logger.h"

// -----------------------------------------------------------------------------
namespace lorenz95 {
// -----------------------------------------------------------------------------

ObservationL95::ObservationL95(const ObsTable & ot, const eckit::Configuration & conf)
  : obsdb_(ot), inputs_(conf)
{}

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

ObservationL95::~ObservationL95() {}

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

void ObservationL95::simulateObs(const GomL95 & gom, ObsVec1D & ovec,
                                 const ObsBias & bias) const {
  for (size_t jj = 0; jj < gom.size(); ++jj) {
    const int ii = gom.getindx(jj);
    ovec(ii) = gom[jj] + bias.value();
  }
}

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

void ObservationL95::print(std::ostream & os) const {
  os << "ObservationL95: Lorenz 95 Obs Operator";
}

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

}  // namespace lorenz95