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
|
/*
* (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 "eccodes.h"
#undef NDEBUG
#include <assert.h>
static void compare_doubles(const double d1, const double d2, const double epsilon)
{
assert(fabs(d1 - d2) < epsilon);
}
int main(int argc, char* argv[])
{
FILE* in = NULL;
codes_handle* h = NULL;
size_t values_len = 0;
int elem_index = 0;
double a_value = 0; /* value at index */
int err = 0;
long compressed = 0;
const char* infile = "../data/bufr/mhen_55.bufr";
const char* key_name = "#5#brightnessTemperature";
const double tolerance = 1e-6;
ProductKind prod_kind = PRODUCT_ANY;
in = fopen(infile, "rb");
assert(in);
while ((h = codes_handle_new_from_file(NULL, in, PRODUCT_BUFR, &err)) != NULL || err != CODES_SUCCESS) {
CODES_CHECK(codes_get_product_kind(h, &prod_kind), 0);
assert(prod_kind == PRODUCT_BUFR);
CODES_CHECK(codes_get_long(h, "compressedData", &compressed), 0);
assert(compressed == 1);
CODES_CHECK(codes_set_long(h, "unpack", 1), 0);
CODES_CHECK(codes_get_size(h, key_name, &values_len), 0);
assert(values_len == 2070);
elem_index = values_len - 1; /* last value of key */
CODES_CHECK(codes_get_double_element(h, key_name, elem_index, &a_value), 0);
printf("Key = %s, index = %d, value = %.5f\n", key_name, elem_index, a_value);
compare_doubles(a_value, 237.43, tolerance);
codes_handle_delete(h);
}
fclose(in);
grib_context_delete(grib_context_get_default());
return 0;
}
|
