mirror of https://github.com/ecmwf/eccodes.git
82 lines
2.0 KiB
C++
82 lines
2.0 KiB
C++
/*
|
|
* (C) Copyright 2005- 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 <stdio.h>
|
|
#include <stdlib.h>
|
|
|
|
#include "grib_api_internal.h"
|
|
|
|
typedef enum IterMode
|
|
{
|
|
WITH_VALUES,
|
|
NO_VALUES
|
|
} IterMode;
|
|
|
|
int main(int argc, char** argv)
|
|
{
|
|
int err = 0;
|
|
double lat = 0, lon = 0, value = 0;
|
|
long n = 0, numberOfDataPoints = 0;
|
|
char* filename = NULL;
|
|
FILE* fin = NULL;
|
|
grib_handle* h = NULL;
|
|
grib_iterator* iter = NULL;
|
|
char* option = NULL;
|
|
IterMode mode = WITH_VALUES;
|
|
int flags = 0;
|
|
double* pValue = NULL;
|
|
bool verbose = false;
|
|
|
|
Assert(argc == 3);
|
|
option = argv[1];
|
|
filename = argv[2];
|
|
if (strcmp(option, "-n")==0) {
|
|
mode = NO_VALUES;
|
|
}
|
|
|
|
fin = fopen(filename, "rb");
|
|
Assert(fin);
|
|
h = grib_handle_new_from_file(0, fin, &err);
|
|
Assert(!err);
|
|
Assert(h);
|
|
|
|
GRIB_CHECK(grib_get_long(h,"numberOfDataPoints", &numberOfDataPoints),0);
|
|
|
|
flags = (mode == NO_VALUES) ? GRIB_GEOITERATOR_NO_VALUES : 0;
|
|
iter = grib_iterator_new(h, flags, &err);
|
|
Assert(!err);
|
|
Assert(iter);
|
|
|
|
Assert(grib_iterator_has_next(iter));
|
|
n = 0;
|
|
|
|
pValue = (mode == NO_VALUES) ? NULL : &value;
|
|
while (grib_iterator_next(iter, &lat, &lon, pValue)) {
|
|
if (n < numberOfDataPoints - 1) {
|
|
Assert(grib_iterator_has_next(iter));
|
|
}
|
|
if (verbose) {
|
|
printf("%9.3f %9.3f ",lat, lon);
|
|
if (pValue) printf(" %.10e\n", *pValue);
|
|
else printf("\n");
|
|
}
|
|
n++;
|
|
}
|
|
Assert(n == numberOfDataPoints);
|
|
|
|
Assert(grib_iterator_has_next(iter) == 0);
|
|
|
|
grib_iterator_delete(iter);
|
|
|
|
grib_handle_delete(h);
|
|
|
|
return 0;
|
|
}
|