CosineOfLatitudeTimesOrbitNode.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 "ufo/predictors/CosineOfLatitudeTimesOrbitNode.h"
 
#include "ioda/ObsSpace.h"
 
#include "ufo/utils/Constants.h"
 
namespace ufo {
 
static PredictorMaker<CosineOfLatitudeTimesOrbitNode>
       makerFuncCosineOfLatitudeTimesOrbitNode_("cosine_of_latitude_times_orbit_node");
 
// -----------------------------------------------------------------------------
 
CosineOfLatitudeTimesOrbitNode::CosineOfLatitudeTimesOrbitNode(
                                const eckit::Configuration & conf, const std::vector<int> & jobs)
  : PredictorBase(conf, jobs) {
  // override the preconditioner from options
  if (conf.has("predictor.options"))
    precond_ = conf.getDouble("predictor.options.preconditioner");
}
 
// -----------------------------------------------------------------------------
void CosineOfLatitudeTimesOrbitNode::compute(const ioda::ObsSpace & odb,
                                             const GeoVaLs &,
                                             const ObsDiagnostics &,
                                             ioda::ObsVector & out) const {
  const std::size_t nlocs = odb.nlocs();
 
  // assure shape of out
  ASSERT(out.nlocs() == nlocs);
 
  // retrieve the sensor view angle
  std::vector<float> cenlat(nlocs, 0.0);
  std::vector<float> node(nlocs, 0.0);
  odb.get_db("MetaData", "latitude", cenlat);
  odb.get_db("MetaData", "sensor_azimuth_angle", node);
 
  const std::size_t njobs = jobs_.size();
  for (std::size_t jl = 0; jl < nlocs; ++jl) {
    for (std::size_t jb = 0; jb < njobs; ++jb) {
      out[jl*njobs+jb] = node[jl] * cos(cenlat[jl] * Constants::deg2rad);
    }
  }
}
 
// -----------------------------------------------------------------------------
 
}  // namespace ufo