ObsSpaces.h

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.
 */
 
#ifndef OOPS_BASE_OBSSPACES_H_
#define OOPS_BASE_OBSSPACES_H_
 
#include <cstddef>
#include <map>
#include <memory>
#include <ostream>
#include <string>
#include <utility>
#include <vector>
 
#include "eckit/config/LocalConfiguration.h"
#include "eckit/exception/Exceptions.h"
#include "eckit/geometry/Point2.h"
 
#include "oops/interface/ObsSpace.h"
#include "oops/mpi/mpi.h"
#include "oops/util/ConfigFunctions.h"  // for vectoriseAndFilter
#include "oops/util/DateTime.h"
#include "oops/util/Logger.h"
#include "oops/util/ObjectCounter.h"
#include "oops/util/parameters/Parameters.h"
#include "oops/util/Printable.h"
#include "oops/util/Timer.h"
 
namespace oops {
 
// -----------------------------------------------------------------------------
template <typename OBS>
class ObsSpaces : public util::Printable,
                  private util::ObjectCounter<ObsSpaces<OBS> > {
  typedef ObsSpace<OBS>                   ObsSpace_;
 
 public:
  typedef typename ObsSpace_::Parameters_ Parameters_;
 
  static const std::string classname() {return "oops::ObsSpaces";}
 
  ObsSpaces(const std::vector<Parameters_> &, const eckit::mpi::Comm &,
            const util::DateTime &, const util::DateTime &,
            const eckit::mpi::Comm & time = oops::mpi::myself());
  ObsSpaces(const eckit::Configuration &, const eckit::mpi::Comm &,
            const util::DateTime &, const util::DateTime &,
            const eckit::mpi::Comm & time = oops::mpi::myself());
  ~ObsSpaces();
 
/// Save files
  void save() const;
 
/// Access
  std::size_t size() const {return spaces_.size();}
  ObsSpace_ & operator[](const std::size_t ii) {return *spaces_.at(ii);}
  const ObsSpace_ & operator[](const std::size_t ii) const {return *spaces_.at(ii);}
 
/// Assimilation window
  const util::DateTime & windowStart() const {return wbgn_;}
  const util::DateTime & windowEnd() const {return wend_;}
 
 private:
  void print(std::ostream &) const;
 
  std::vector<std::shared_ptr<ObsSpace_> > spaces_;
  const util::DateTime wbgn_;
  const util::DateTime wend_;
};
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
ObsSpaces<OBS>::ObsSpaces(const std::vector<Parameters_> & params, const eckit::mpi::Comm & comm,
                          const util::DateTime & bgn, const util::DateTime & end,
                          const eckit::mpi::Comm & time)
  : spaces_(0), wbgn_(bgn), wend_(end)
{
  spaces_.reserve(params.size());
  for (const Parameters_ & param : params) {
    Log::debug() << "ObsSpaces::ObsSpaces : conf " << param << std::endl;
    auto tmp = std::make_shared<ObsSpace_>(param, comm, bgn, end, time);
    spaces_.push_back(std::move(tmp));
  }
  ASSERT(spaces_.size() >0);
}
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
ObsSpaces<OBS>::ObsSpaces(const eckit::Configuration & conf, const eckit::mpi::Comm & comm,
                          const util::DateTime & bgn, const util::DateTime & end,
                          const eckit::mpi::Comm & time)
  : ObsSpaces(  // Split conf into subconfigurations, extract the "obs space" section from each
                // of them, then validate and deserialize that section into a Parameters_ object
              validateAndDeserialize<Parameters_>(util::vectoriseAndFilter(conf, "obs space")),
              comm, bgn, end, time)
{}
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
ObsSpaces<OBS>::~ObsSpaces() {}
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
void ObsSpaces<OBS>::save() const {
  for (std::size_t jj = 0; jj < spaces_.size(); ++jj) {
    spaces_[jj]->save();
  }
}
 
// -----------------------------------------------------------------------------
 
template <typename OBS>
void ObsSpaces<OBS>::print(std::ostream & os) const {
  for (std::size_t jj = 0; jj < spaces_.size(); ++jj) {
    os << *spaces_[jj];
  }
}
 
// -----------------------------------------------------------------------------
 
}  // namespace oops
 
#endif  // OOPS_BASE_OBSSPACES_H_