LocationsQG.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 QG_MODEL_LOCATIONSQG_H_
#define QG_MODEL_LOCATIONSQG_H_

#include <ostream>
#include <string>
#include <vector>

#include "model/QgFortran.h"
#include "oops/util/ObjectCounter.h"
#include "oops/util/Printable.h"

namespace qg {

/// LocationsQG class to handle locations for QG model.

class LocationsQG : public util::Printable,
                    private util::ObjectCounter<LocationsQG> {
 public:
  static const std::string classname() {return "qg::LocationsQG";}

  explicit LocationsQG(const F90locs key) : keyLoc_(key) {}

  explicit LocationsQG(const eckit::Configuration &);

  ~LocationsQG() {qg_loc_delete_f90(keyLoc_);}

  size_t size() const {
    int nobs;
    qg_loc_nobs_f90(keyLoc_, nobs);
    return nobs;
  }

  int toFortran() const {return keyLoc_;}
 private:
  void print(std::ostream & os) const {
    int nobs;
    qg_loc_nobs_f90(keyLoc_, nobs);

    std::vector<double> xyz(3);

    for (size_t jj=0; jj < static_cast<size_t>(nobs); ++jj) {
      qg_loc_element_f90(keyLoc_, jj, &xyz[0]);
      os << "loc " << jj << std::setprecision(2) << ": x = " << xyz[0]
         << ", y = " << xyz[1] << ", z = " << xyz[2] << std::endl;
    }
  }
  F90locs keyLoc_;
};

}  // namespace qg

#endif  // QG_MODEL_LOCATIONSQG_H_