1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
|
/*
* (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;
ECCODES_ASSERT(argc == 3);
option = argv[1];
filename = argv[2];
if (strcmp(option, "-n")==0) {
mode = NO_VALUES;
}
fin = fopen(filename, "rb");
ECCODES_ASSERT(fin);
h = grib_handle_new_from_file(0, fin, &err);
ECCODES_ASSERT(!err);
ECCODES_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);
ECCODES_ASSERT(!err);
ECCODES_ASSERT(iter);
ECCODES_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) {
ECCODES_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++;
}
ECCODES_ASSERT(n == numberOfDataPoints);
ECCODES_ASSERT(grib_iterator_has_next(iter) == 0);
{
// Test getting the previous value from a geoiterator.
// Note: Only implemented for regular grids
char gridType[128] = {0,};
size_t len = sizeof(gridType);
GRIB_CHECK(grib_get_string(h, "gridType", gridType, &len), 0);
if (strstr(gridType, "regular_")) {
double last_lat = 0, last_lon = 0, last_value = 0;
int result = grib_iterator_previous(iter, &last_lat, &last_lon, &last_value);
ECCODES_ASSERT(result == 1);
ECCODES_ASSERT(lat == last_lat);
ECCODES_ASSERT(lon == last_lon);
ECCODES_ASSERT(value == last_value);
}
}
grib_iterator_delete(iter);
grib_handle_delete(h);
return 0;
}
|
