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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
|
/*
* This is test code from the NCEPLIBS-g2c project.
*
* This test is for decoding full grib2 message.
*
* Dusan Jovic, July, 2021
* Ed Hartnett
*/
#include "grib2_int.h"
int
main()
{
/* This is the GRIB2 message. */
unsigned char cgrib[] = {
0x47, 0x52, 0x49, 0x42, 0x00, 0x00, 0x02, 0x02, 0x00, 0x00, 0x00, 0x00,
0x00, 0x00, 0x00, 0xc3, 0x00, 0x00, 0x00, 0x15, 0x01, 0x00, 0x07, 0x00,
0x00, 0x02, 0x01, 0x01, 0x07, 0xe5, 0x07, 0x0e, 0x06, 0x00, 0x00, 0x00,
0x01, 0x00, 0x00, 0x00, 0x48, 0x03, 0x00, 0x00, 0x00, 0x00, 0x79, 0x00,
0x00, 0x00, 0x00, 0x06, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x0b, 0x00,
0x00, 0x00, 0x0b, 0x00, 0x00, 0x00, 0x00, 0xff, 0xff, 0xff, 0xff, 0x01,
0xc9, 0xc3, 0x80, 0x01, 0x31, 0x2d, 0x00, 0x30, 0x02, 0x62, 0x5a, 0x00,
0x01, 0xc9, 0xc3, 0x80, 0x00, 0x0f, 0x42, 0x40, 0x00, 0x0f, 0x42, 0x40,
0x40, 0x00, 0x00, 0x00, 0x22, 0x04, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
0x02, 0x00, 0x60, 0x00, 0x00, 0x00, 0x01, 0x00, 0x00, 0x00, 0x06, 0x01,
0x00, 0x00, 0x00, 0x00, 0x00, 0xff, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00,
0x00, 0x00, 0x15, 0x05, 0x00, 0x00, 0x00, 0x79, 0x00, 0x00, 0x00, 0x00,
0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x01, 0x00, 0x00, 0x00, 0x00, 0x06,
0x06, 0xff, 0x00, 0x00, 0x00, 0x15, 0x07, 0xff, 0xef, 0xf7, 0xe0, 0x00,
0x00, 0x00, 0x08, 0x00, 0x01, 0x83, 0x38, 0xee, 0x3f, 0xa7, 0x80, 0x37,
0x37, 0x37, 0x37};
/* This is the data encoded in the GRIB2 message. */
float fld_ok[121] = {/* 11x11 grid */
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1,
1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1,
0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1,
0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1,
1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1};
int i;
int ret;
printf("Testing decoding full grib2 message.\n");
printf("Testing g2_unpack7()...");
{
g2int iofst = 1360;
g2int igdsnum = 0;
g2int *igdstmpl = NULL;
g2int idrsnum = 0;
g2int idrstmpl[5] = {0, 0, 0, 1, 0};
g2int ndpts = 121;
float *fld;
/* g2c_set_log_level(10); */
/* Call g2_unpack7() on our message. */
if ((ret = g2_unpack7(cgrib, &iofst, igdsnum, igdstmpl, idrsnum, idrstmpl, ndpts, &fld)))
return ret;
/* Check the data. */
for (i = 0; i < ndpts; i++)
if (fld[i] != fld_ok[i])
return G2C_ERROR;
/* Free memory allocated by g2_unpack7(). */
free(fld);
}
printf("ok!\n");
printf("Testing g2c_unpack7()...");
{
int igdsnum = 0;
long long int *igdstmpl = NULL;
int idrsnum = 0;
long long int idrstmpl[5] = {0, 0, 0, 1, 0};
int ndpts = 121;
float *fld;
/* g2c_set_log_level(10); */
/* Allocate memory to hold data. */
if (!(fld = malloc(sizeof(float) * ndpts)))
return G2C_ERROR;
/* Call g2c_unpack7() on our message. */
if ((ret = g2c_unpack7(&cgrib[1360 / 8], igdsnum, 0, igdstmpl, idrsnum, 5, idrstmpl, ndpts, fld)))
return ret;
/* Check the data. */
for (i = 0; i < ndpts; i++)
if (fld[i] != fld_ok[i])
return G2C_ERROR;
/* Free memory allocated by g2_unpack7(). */
free(fld);
}
printf("Testing g2_info() and g2_getfld()...");
{
g2int listsec0[3];
g2int listsec0_ok[3] = {2, 2, 195};
g2int listsec1[13];
g2int listsec1_ok[13] = {7, 0, 2, 1, 1, 2021, 7, 14, 6, 0, 0, 0, 1};
g2int numfields;
g2int numlocal;
gribfield *gfld = NULL;
/* Call g2_info() on our message. */
/* g2c_set_log_level(10); */
if ((ret = g2_info(cgrib, listsec0, listsec1, &numfields, &numlocal)))
return ret;
/* Check results. */
for (i = 0; i < 3; i++)
if (listsec0[i] != listsec0_ok[i])
return G2C_ERROR;
for (i = 0; i < 13; i++)
if (listsec1[i] != listsec1_ok[i])
return G2C_ERROR;
if (numfields != 1)
return G2C_ERROR;
if (numlocal != 0)
return G2C_ERROR;
/* Now read the data field. */
if ((ret = g2_getfld(cgrib, 1, 1, 1, &gfld)) != 0)
return G2C_ERROR;
/* Check results. */
if (gfld->version != 2)
return G2C_ERROR;
if (gfld->ndpts != 121) /* 11x11 grid */
return G2C_ERROR;
for (i = 0; i < gfld->ndpts; i++)
if (gfld->fld[i] != fld_ok[i])
return G2C_ERROR;
/* Free resources. */
g2_free(gfld);
}
printf("ok!\n");
printf("SUCCESS!\n");
return 0;
}
|
