l95/src/lorenz95/AnalyticInit.cc

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.
 */
 
#include "lorenz95/AnalyticInit.h"
 
#include <cmath>
 
#include "lorenz95/GomL95.h"
#include "lorenz95/LocsL95.h"
 
#include "oops/util/Logger.h"
 
namespace lorenz95 {
 
// -----------------------------------------------------------------------------
AnalyticInit::AnalyticInit(const eckit::Configuration & config): config_(config)
{ }
// -----------------------------------------------------------------------------
/*! GomL95 analytic initialization */
void AnalyticInit::fillGeoVaLs(const LocsL95 & locs,
                               GomL95 & geovals) const {
  oops::Log::trace() << "AnalyticInit::fillGeoVaLs " << std::endl;
  const size_t nlocs = geovals.size();
  // analytic init for testing interpolation
  // mean value; default is zero.
  const double mean = config_.getDouble("mean", 0.0);
  for (size_t jj = 0; jj < nlocs; ++jj) geovals[jj] = mean;
  // add a sinus function if specified in yaml
  if (config_.has("sinus")) {
    const double zz = config_.getDouble("sinus");
    for (size_t jj = 0; jj < nlocs; ++jj)
      geovals[jj] += zz * std::sin(2.0*M_PI*locs[jj]);
  }
  if (!config_.has("mean") && !config_.has("sinus")) {
    oops::Log::warning() << "Using default analytic init (all zeros)" << std::endl;
  }
 
  oops::Log::trace() << "GomL95::GomL95 analytic init finished" << std::endl;
}
// -----------------------------------------------------------------------------
}  // namespace lorenz95