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
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
|
/* This is a test for the NCEPLIBS-g2c project.
*
* This test is for the g2c_degrib2() function, which prints a summary
* of the GRIB2 file.
*
* Ed Hartnett 9/17/22
*/
#include "g2ctest.h"
#include "grib2_int.h"
#include <libgen.h>
#define FILE_NAME "tst_degrib2.txt"
#define WAVE_FILE "data/gdaswave.t00z.wcoast.0p16.f000.grib2"
#define REF_FILE "data/ref_gdaswave.t00z.wcoast.0p16.f000.grib2.degrib2"
#define DEGRIB2_FILE "gdaswave.t00z.wcoast.0p16.f000.degrib2"
#define REF_INDEX_FILE "data/ref_gdaswave.t00z.wcoast.0p16.f000.grb2index"
#define TEST_INDEX_FILE "tst_gdaswave.t00z.wcoast.0p16.f000.grb2index"
#define FTP_FILE "data/WW3_Regional_US_West_Coast_20220718_0000.grib2"
#define REF_FTP_INDEX_FILE "data/ref_WW3_Regional_US_West_Coast_20220718_0000.grb2index"
#define FTP_DEGRIB2_FILE "WW3_Regional_US_West_Coast_20220718_0000.degrib2"
#define REF_FTP_DEGRIB2_FILE "data/ref_WW3_Regional_US_West_Coast_20220718_0000.degrib2"
#define GDAS_FILE "data/gdas.t12z.pgrb2.1p00.anl.grib2"
#define GDAS_DEGRIB2_FILE "gdas.t12z.pgrb2.1p00.anl.degrib2"
#define REF_GDAS_DEGRIB2_FILE "data/ref_gdas.t12z.pgrb2.1p00.anl.degrib2"
#define REF_GDAS_INDEX_FILE "gdas.t12z.pgrb2.1p00.anl.grb2index"
#define GDAS_NUM_MSG 2
#define FTP_NUM_MSG 688
#define MAX_FILENAME_LEN 256
#define NUM_TESTS 2
int
main()
{
printf("Testing g2c degrib2 function.\n");
#ifdef JPEG
#ifdef FTP_TEST_FILES
printf("Testing degrib2 on file %s downloaded via FTP...\n", FTP_FILE);
{
int g2cid;
int num_msg;
unsigned char sig_ref_time, month, day, hour, minute, second;
short year;
short center, subcenter;
unsigned char master_version, local_version;
int num_fields, num_local;
unsigned char discipline;
unsigned char expected_discipline[FTP_NUM_MSG] =
{10, 10, 10, 10, 0, 10, 10, 10, 0, 10, 10, 10, 10, 0, 10, 0, 10, 10,
10, 0, 0, 10, 10, 10, 10, 10, 0, 10, 10, 10, 10, 0, 0, 10, 10, 10,
10, 10, 10, 0, 0, 10, 10, 10, 10, 10, 10, 0, 0, 10, 10, 10, 10, 0,
0, 10, 10, 10, 10, 10, 10, 10, 10, 0, 10, 10, 0, 0, 0, 0, 10, 10,
10, 10, 10, 10, 10, 10, 10, 10, 0, 10, 10, 10, 10, 0, 10, 10, 10,
10, 10, 0, 10, 0, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10, 10, 0,
10, 10, 10, 10, 0, 0, 0, 10, 0, 10, 10, 10, 10, 0, 10, 0, 10, 10,
10, 10, 10, 10, 10, 10, 0, 10, 10, 10, 10, 0, 10, 0, 10, 10, 10,
0, 10, 10, 10, 0, 10, 0, 0, 10, 0, 10, 10, 10, 10, 10, 10, 10, 10,
10, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10, 0, 10, 10, 0, 0, 10, 0,
10, 10, 10, 0, 10, 10, 10, 10, 0, 10, 0, 10, 10, 10, 10, 0, 0, 10,
10, 10, 0, 10, 10, 10, 10, 10, 10, 10, 10, 0, 10, 10, 10, 10, 10,
0, 0, 10, 10, 10, 0, 10, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10,
10, 10, 10, 10, 0, 0, 0, 0, 10, 10, 10, 0, 10, 10, 10, 10, 10, 0,
10, 10, 10, 0, 10, 10, 10, 0, 10, 10, 10, 0, 10, 10, 10, 0, 0, 0,
10, 10, 10, 10, 10, 10, 10, 10, 10, 0, 10, 0, 0, 10, 10, 10, 10,
10, 10, 10, 0, 10, 10, 10, 10, 0, 0, 10, 10, 10, 10, 10, 0, 10,
10, 10, 10, 0, 10, 0, 10, 10, 10, 10, 10, 10, 10, 0, 10, 10, 10,
10, 10, 0, 0, 0, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10, 0, 10, 0,
10, 10, 0, 10, 10, 0, 0, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10, 10,
10, 0, 10, 10, 0, 0, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10, 10, 10,
10, 0, 10, 10, 10, 0, 10, 0, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10,
10, 10, 10, 0, 10, 0, 10, 10, 10, 0, 10, 10, 0, 10, 10, 0, 10, 0,
10, 0, 10, 0, 10, 10, 10, 10, 10, 10, 10, 10, 0, 10, 10, 0, 10, 10,
0, 0, 10, 10, 10, 10, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10, 10, 0,
10, 10, 10, 0, 10, 10, 10, 10, 10, 0, 10, 10, 10, 10, 10, 10, 10,
10, 10, 0, 10, 10, 0, 10, 10, 10, 10, 10, 10, 0, 0, 0, 10, 10, 10,
10, 10, 10, 0, 10, 0, 10, 10, 10, 0, 0, 10, 10, 10, 10, 0, 0, 10,
10, 10, 10, 10, 0, 0, 10, 10, 10, 10, 10, 10, 10, 0, 0, 0, 0, 10,
10, 10, 0, 10, 10, 10, 10, 0, 0, 10, 10, 10, 10, 10, 10, 10, 10,
10, 10, 0, 10, 0, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10,
10, 10, 0, 0, 0, 10, 10, 0, 10, 0, 0, 10, 0, 10, 10, 10, 10, 0, 10,
10, 10, 10, 10, 10, 10, 10, 0, 10, 10, 10, 10, 10, 0, 10, 10, 10, 0,
10, 10, 10, 10, 10, 10, 10, 10, 0, 0, 0, 10, 10, 10, 10, 10, 10, 10,
0, 10, 10, 10, 10, 10, 0, 0, 10, 10, 10, 10, 0, 10, 10, 0, 0, 0, 10,
10, 10, 10, 0, 10, 10, 10, 0, 10, 0, 0, 0, 10, 0, 10, 10, 10, 10,
10, 10, 10, 10, 0, 10, 0, 10, 0, 10, 10, 10, 10, 10, 10, 0, 10, 0,
10, 10, 10, 10, 10, 10, 10, 10, 0, 0, 10, 10, 10, 10, 0, 10, 10,
10, 0, 10, 10, 10, 0, 10, 10, 10, 0, 10, 10, 10, 10};
int t, m;
int ret;
for (t = 0; t < NUM_TESTS; t++)
{
/* g2c_set_log_level(11); */
/* Open the data file with and without the index file. */
if (t)
{
printf("\ttesting degrib2 on file %s downloaded via FTP using index...", FTP_FILE);
if ((ret = g2c_open_index(FTP_FILE, REF_FTP_INDEX_FILE, 0, &g2cid)))
return ret;
}
else
{
printf("\ttesting degrib2 on file %s downloaded via FTP without using index...", FTP_FILE);
if ((ret = g2c_open(FTP_FILE, 0, &g2cid)))
return ret;
}
/* Check some stuff. */
if ((ret = g2c_inq(g2cid, &num_msg)))
return ret;
if (num_msg != FTP_NUM_MSG)
return G2C_ERROR;
/* Check some stuff in all messages. */
for (m = 0; m < num_msg; m++)
{
if ((ret = g2c_inq_msg(g2cid, m, &discipline, &num_fields, &num_local,
¢er, &subcenter, &master_version, &local_version)))
return ret;
if (center != 7 || subcenter != 0 || master_version != 2 || local_version != 1)
return G2C_ERROR;
if (num_fields != 1 || num_local)
return G2C_ERROR;
if (discipline != expected_discipline[m])
return G2C_ERROR;
}
/* Check some stuff about the last message. */
if ((ret = g2c_inq_msg(g2cid, 687, &discipline, &num_fields, &num_local,
¢er, &subcenter, &master_version, &local_version)))
return ret;
if (center != 7 || subcenter != 0 || master_version != 2 || local_version != 1)
return G2C_ERROR;
if (num_local || num_fields != 1 || discipline != 10)
return G2C_ERROR;
/* Inquire about the date/time. */
if ((ret = g2c_inq_msg_time(g2cid, 687, &sig_ref_time, &year, &month, &day, &hour,
&minute, &second)))
return ret;
/* Check date/time. */
if (sig_ref_time != 1 || year != 2022 || month != 7 || day != 18 ||
hour != 0 || minute != 0 || second != 0)
return G2C_ERROR;
/* Output a degrib2 file. */
if ((ret = g2c_degrib2(g2cid, FTP_DEGRIB2_FILE)))
return ret;
/* Close the file. */
if ((ret = g2c_close(g2cid)))
return ret;
/* Compare the degrib2 output to our reference file. */
if ((ret = compare_degrib2_files2(FTP_DEGRIB2_FILE, REF_FTP_DEGRIB2_FILE)))
return ret;
printf("\tok!\n");
}
}
printf("ok!\n");
printf("Testing degrib2 on file %s downloaded via FTP...\n", GDAS_FILE);
{
int g2cid;
int num_msg;
int t;
int ret;
/* for (t = 0; t < NUM_TESTS; t++) */
for (t = 0; t < 1; t++)
{
/* Open the data file with and without the index file. */
if (t)
{
printf("\ttesting degrib2 on file %s downloaded via FTP using index...", GDAS_FILE);
if ((ret = g2c_open_index(GDAS_FILE, REF_GDAS_INDEX_FILE, 0, &g2cid)))
return ret;
}
else
{
printf("\ttesting degrib2 on file %s downloaded via FTP without using index...", GDAS_FILE);
if ((ret = g2c_open(GDAS_FILE, 0, &g2cid)))
return ret;
}
/* Check some stuff. */
if ((ret = g2c_inq(g2cid, &num_msg)))
return ret;
if (num_msg != GDAS_NUM_MSG)
return G2C_ERROR;
/* Output a degrib2 file. */
if ((ret = g2c_degrib2(g2cid, GDAS_DEGRIB2_FILE)))
return ret;
/* Close the file. */
if ((ret = g2c_close(g2cid)))
return ret;
/* Compare the degrib2 output to our reference file. */
if ((ret = compare_degrib2_files2(GDAS_DEGRIB2_FILE, REF_GDAS_DEGRIB2_FILE)))
return ret;
printf("\tok!\n");
}
}
printf("ok!\n");
#endif /* FTP_TEST_FILES */
#endif /* JPEG */
printf("SUCCESS!\n");
return 0;
}
|
