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
|
/* This is a test for the NCEPLIBS-g2c project. This test is for the
* g2c file functions on the GRIB2 file that's included in the
* tests/data directory for testion,
* gdaswave.t00z.wcoast.0p16.f000.grib2.
*
* Ed Hartnett 7/5/23
*/
#include "grib2_int.h"
#include <stdio.h>
#include <stdlib.h>
#define MRMS_FILE "data/MRMS_MultiSensor_QPE_24H_Pass2_00.00_20230621-110000.grib2"
#define LAT_LEN 3500
#define LON_LEN 7000
#define LATITUDE "Latitude"
#define LONGITUDE "Longitude"
int
main()
{
printf("Testing g2c file functions with %s.\n", MRMS_FILE);
printf("Testing g2c_inq_dim()...");
{
int g2cid;
int num_msg;
size_t len0;
char dimname0[G2C_MAX_NAME + 1];
size_t len1;
char dimname1[G2C_MAX_NAME + 1];
float lat[LAT_LEN];
float lon[LON_LEN];
/* int d; */
/* Open GRIB2 file. */
g2c_set_log_level(10);
if (g2c_open(MRMS_FILE, 0, &g2cid))
return G2C_ERROR;
/* Find the number of messages. */
if (g2c_inq(g2cid, &num_msg))
return G2C_ERROR;
printf("Number of messages: %d\n", num_msg);
if (num_msg != 1)
return G2C_ERROR;
/* Find the dimensions of the product in the first message. */
if (g2c_inq_dim(g2cid, 0, 0, 0, &len0, dimname0, lat))
return G2C_ERROR;
printf("len0 %ld dimname0 %s\n", len0, dimname0);
printf("lat[0] = %f lat[%ld] = %f\n", lat[0], len0 - 1, lat[len0 - 1]);
if (len0 != LAT_LEN || strcmp(dimname0, LATITUDE))
return G2C_ERROR;
if (lat[0] != 54995000 || lat[len0 - 1] != 20005000)
return G2C_ERROR;
if (g2c_inq_dim(g2cid, 0, 0, 1, &len1, dimname1, lon))
return G2C_ERROR;
if (len1 != LON_LEN || strcmp(dimname1, LONGITUDE))
return G2C_ERROR;
printf("len1 %ld dimname1 %s\n", len1, dimname1);
printf("lon[0] = %f lon[%ld] = %f\n", lon[0], len1 - 1, lon[len1 - 1]);
if (len1 != LON_LEN || strcmp(dimname1, LONGITUDE))
return G2C_ERROR;
if (lon[0] != 230004992 || lon[len1 - 1] != 160014992)
return G2C_ERROR;
}
printf("ok!\n");
printf("SUCCESS!\n");
return 0;
}
|
