ObsBiasCorrection.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/ObsBiasCorrection.h"
 
#include <iostream>
#include <string>
 
#include "eckit/config/Configuration.h"
#include "eckit/config/LocalConfiguration.h"
#include "lorenz95/ObsBias.h"
#include "oops/util/Logger.h"
 
// -----------------------------------------------------------------------------
namespace lorenz95 {
// -----------------------------------------------------------------------------
ObsBiasCorrection::ObsBiasCorrection(const ObsTable &, const Parameters_ & params)
  : bias_(0.0), active_(false)
{
  if (params.covariance.value() != boost::none &&
      params.covariance.value()->standardDeviation.value() != boost::none) {
    active_ = true;
  }
  if (active_) {oops::Log::trace() << "ObsBiasCorrection::ObsBiasCorrection created." << std::endl;}
}
// -----------------------------------------------------------------------------
ObsBiasCorrection::ObsBiasCorrection(const ObsBiasCorrection & other,
                                     const bool copy)
  : bias_(0.0), active_(other.active_)
{
  if (active_ && copy) bias_ = other.bias_;
}
// -----------------------------------------------------------------------------
void ObsBiasCorrection::diff(const ObsBias & b1, const ObsBias & b2) {
  if (active_) bias_ = b1.value() - b2.value();
}
// -----------------------------------------------------------------------------
void ObsBiasCorrection::zero() {
  bias_ = 0.0;
}
// -----------------------------------------------------------------------------
ObsBiasCorrection & ObsBiasCorrection::operator=(const ObsBiasCorrection & rhs) {
  if (active_) {
    bias_ = rhs.bias_;
  } else {
    bias_ = 0.0;
  }
  return *this;
}
// -----------------------------------------------------------------------------
ObsBiasCorrection & ObsBiasCorrection::operator+=(const ObsBiasCorrection & rhs) {
  if (active_) bias_ += rhs.bias_;
  return *this;
}
// -----------------------------------------------------------------------------
ObsBiasCorrection & ObsBiasCorrection::operator-=(const ObsBiasCorrection & rhs) {
  if (active_) bias_ -= rhs.bias_;
  return *this;
}
// -----------------------------------------------------------------------------
ObsBiasCorrection & ObsBiasCorrection::operator*=(const double fact) {
  if (active_) bias_ *= fact;
  return *this;
}
// -----------------------------------------------------------------------------
void ObsBiasCorrection::axpy(const double fact, const ObsBiasCorrection & rhs) {
  if (active_) bias_ += fact * rhs.bias_;
}
// -----------------------------------------------------------------------------
double ObsBiasCorrection::dot_product_with(const ObsBiasCorrection & rhs) const {
  double zz = 0.0;
  if (active_) zz = bias_ * rhs.bias_;
  return zz;
}
// -----------------------------------------------------------------------------
size_t ObsBiasCorrection::serialSize() const {
  return 1;
}
// -----------------------------------------------------------------------------
void ObsBiasCorrection::serialize(std::vector<double> & vect) const {
  vect.push_back(bias_);
}
// -----------------------------------------------------------------------------
void ObsBiasCorrection::deserialize(const std::vector<double> & vect, size_t & index) {
  if (!vect.empty()) bias_ = vect.at(index);
  ++index;
}
// -----------------------------------------------------------------------------
void ObsBiasCorrection::print(std::ostream & os) const {
  if (active_) {os << std::endl << "ObsBiasCorrection = " << bias_;}
}
// -----------------------------------------------------------------------------
}  // namespace lorenz95