GeometryQGIterator.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 "eckit/config/Configuration.h"
#include "model/GeometryQGIterator.h"
#include "model/QgFortran.h"
#include "oops/util/Logger.h"

// -----------------------------------------------------------------------------

namespace qg {

// -----------------------------------------------------------------------------

GeometryQGIterator::GeometryQGIterator(const GeometryQGIterator& iter) {
  qg_geo_iter_clone_f90(keyIter_, iter.toFortran());
}

// -----------------------------------------------------------------------------

GeometryQGIterator::GeometryQGIterator(const GeometryQG& geom, const int & index) {
  qg_geo_iter_setup_f90(keyIter_, geom.toFortran(), index);
}


// -----------------------------------------------------------------------------

GeometryQGIterator::~GeometryQGIterator() {
  qg_geo_iter_delete_f90(keyIter_);
}

// -----------------------------------------------------------------------------

bool GeometryQGIterator::operator==(const GeometryQGIterator & other) const {
  int equals = 0;
  qg_geo_iter_equals_f90(keyIter_, other.toFortran(), equals);
  return (equals == 1);
}

// -----------------------------------------------------------------------------

bool GeometryQGIterator::operator!=(const GeometryQGIterator & other) const {
  int equals = 0;
  qg_geo_iter_equals_f90(keyIter_, other.toFortran(), equals);
  return (equals == 0);
}

// -----------------------------------------------------------------------------

oops::GeoLocation GeometryQGIterator::operator*() const {
  double lat, lon;
  qg_geo_iter_current_f90(keyIter_, lat, lon);
  return oops::GeoLocation(lat, lon);
}

// -----------------------------------------------------------------------------

GeometryQGIterator& GeometryQGIterator::operator++() {
  qg_geo_iter_next_f90(keyIter_);
  return *this;
}

// -----------------------------------------------------------------------------

void GeometryQGIterator::print(std::ostream & os) const {
  os << "GeometryQGIterator, key: " <<  keyIter_ << std::endl;
}

// -----------------------------------------------------------------------------

}  // namespace qg