1 : #include <stdio.h>
2 : #include <stdlib.h>
3 : #include "grib2.h"
4 :
5 6 : g2int g2_info(unsigned char *cgrib,g2int *listsec0,g2int *listsec1,
6 : g2int *numfields,g2int *numlocal)
7 : //$$$ SUBPROGRAM DOCUMENTATION BLOCK
8 : // . . . .
9 : // SUBPROGRAM: g2_info
10 : // PRGMMR: Gilbert ORG: W/NP11 DATE: 2002-10-28
11 : //
12 : // ABSTRACT: This subroutine searches through a GRIB2 message and
13 : // returns the number of gridded fields found in the message and
14 : // the number (and maximum size) of Local Use Sections.
15 : // Also various checks are performed
16 : // to see if the message is a valid GRIB2 message.
17 : //
18 : // PROGRAM HISTORY LOG:
19 : // 2002-10-28 Gilbert
20 : //
21 : // USAGE: int g2_info(unsigned char *cgrib,g2int *listsec0,g2int *listsec1,
22 : // g2int *numfields,g2int *numlocal)
23 : // INPUT ARGUMENT:
24 : // cgrib - Character pointer to the GRIB2 message
25 : //
26 : // OUTPUT ARGUMENTS:
27 : // listsec0 - pointer to an array containing information decoded from
28 : // GRIB Indicator Section 0.
29 : // Must be allocated with >= 3 elements.
30 : // listsec0[0]=Discipline-GRIB Master Table Number
31 : // (see Code Table 0.0)
32 : // listsec0[1]=GRIB Edition Number (currently 2)
33 : // listsec0[2]=Length of GRIB message
34 : // listsec1 - pointer to an array containing information read from GRIB
35 : // Identification Section 1.
36 : // Must be allocated with >= 13 elements.
37 : // listsec1[0]=Id of orginating centre (Common Code Table C-1)
38 : // listsec1[1]=Id of orginating sub-centre (local table)
39 : // listsec1[2]=GRIB Master Tables Version Number (Code Table 1.0)
40 : // listsec1[3]=GRIB Local Tables Version Number
41 : // listsec1[4]=Significance of Reference Time (Code Table 1.1)
42 : // listsec1[5]=Reference Time - Year (4 digits)
43 : // listsec1[6]=Reference Time - Month
44 : // listsec1[7]=Reference Time - Day
45 : // listsec1[8]=Reference Time - Hour
46 : // listsec1[9]=Reference Time - Minute
47 : // listsec1[10]=Reference Time - Second
48 : // listsec1[11]=Production status of data (Code Table 1.2)
49 : // listsec1[12]=Type of processed data (Code Table 1.3)
50 : // numfields- The number of gridded fields found in the GRIB message.
51 : // That is, the number of occurences of Sections 4 - 7.
52 : // numlocal - The number of Local Use Sections ( Section 2 ) found in
53 : // the GRIB message.
54 : //
55 : // RETURN VALUES:
56 : // ierr - Error return code.
57 : // 0 = no error
58 : // 1 = Beginning characters "GRIB" not found.
59 : // 2 = GRIB message is not Edition 2.
60 : // 3 = Could not find Section 1, where expected.
61 : // 4 = End string "7777" found, but not where expected.
62 : // 5 = End string "7777" not found at end of message.
63 : // 6 = Invalid section number found.
64 : //
65 : // REMARKS: None
66 : //
67 : // ATTRIBUTES:
68 : // LANGUAGE: C
69 : // MACHINE:
70 : //
71 : //$$$
72 : {
73 :
74 6 : g2int ierr,mapsec1len=13;
75 6 : g2int mapsec1[13]={2,2,1,1,1,2,1,1,1,1,1,1,1};
76 : g2int i,j,istart,iofst,lengrib,lensec0,lensec1;
77 : g2int ipos,isecnum,nbits,lensec;
78 :
79 6 : ierr=0;
80 6 : *numlocal=0;
81 6 : *numfields=0;
82 : //
83 : // Check for beginning of GRIB message in the first 100 bytes
84 : //
85 6 : istart=-1;
86 6 : for (j=0;j<100;j++) {
87 12 : if (cgrib[j]=='G' && cgrib[j+1]=='R' &&cgrib[j+2]=='I' &&
88 6 : cgrib[j+3]=='B') {
89 6 : istart=j;
90 6 : break;
91 : }
92 : }
93 6 : if (istart == -1) {
94 0 : printf("g2_info: Beginning characters GRIB not found.");
95 0 : ierr=1;
96 0 : return(ierr);
97 : }
98 : //
99 : // Unpack Section 0 - Indicator Section
100 : //
101 6 : iofst=8*(istart+6);
102 6 : gbit(cgrib,listsec0+0,iofst,8); // Discipline
103 6 : iofst=iofst+8;
104 6 : gbit(cgrib,listsec0+1,iofst,8); // GRIB edition number
105 6 : iofst=iofst+8;
106 6 : iofst=iofst+32;
107 6 : gbit(cgrib,&lengrib,iofst,32); // Length of GRIB message
108 6 : iofst=iofst+32;
109 6 : listsec0[2]=lengrib;
110 6 : lensec0=16;
111 6 : ipos=istart+lensec0;
112 : //
113 : // Currently handles only GRIB Edition 2.
114 : //
115 6 : if (listsec0[1] != 2) {
116 0 : printf("g2_info: can only decode GRIB edition 2.");
117 0 : ierr=2;
118 0 : return(ierr);
119 : }
120 : //
121 : // Unpack Section 1 - Identification Section
122 : //
123 6 : gbit(cgrib,&lensec1,iofst,32); // Length of Section 1
124 6 : iofst=iofst+32;
125 6 : gbit(cgrib,&isecnum,iofst,8); // Section number ( 1 )
126 6 : iofst=iofst+8;
127 6 : if (isecnum != 1) {
128 0 : printf("g2_info: Could not find section 1.");
129 0 : ierr=3;
130 0 : return(ierr);
131 : }
132 : //
133 : // Unpack each input value in array listsec1 into the
134 : // the appropriate number of octets, which are specified in
135 : // corresponding entries in array mapsec1.
136 : //
137 84 : for (i=0;i<mapsec1len;i++) {
138 78 : nbits=mapsec1[i]*8;
139 78 : gbit(cgrib,listsec1+i,iofst,nbits);
140 78 : iofst=iofst+nbits;
141 : }
142 6 : ipos=ipos+lensec1;
143 : //
144 : // Loop through the remaining sections to see if they are valid.
145 : // Also count the number of times Section 2
146 : // and Section 4 appear.
147 : //
148 : for (;;) {
149 42 : if (cgrib[ipos]=='7' && cgrib[ipos+1]=='7' && cgrib[ipos+2]=='7' &&
150 6 : cgrib[ipos+3]=='7') {
151 6 : ipos=ipos+4;
152 6 : if (ipos != (istart+lengrib)) {
153 0 : printf("g2_info: '7777' found, but not where expected.\n");
154 0 : ierr=4;
155 0 : return(ierr);
156 : }
157 : break;
158 : }
159 :
160 30 : iofst=ipos*8;
161 30 : gbit(cgrib,&lensec,iofst,32); // Get Length of Section
162 30 : iofst=iofst+32;
163 30 : gbit(cgrib,&isecnum,iofst,8); // Get Section number
164 30 : iofst=iofst+8;
165 30 : ipos=ipos+lensec; // Update beginning of section pointer
166 30 : if (ipos > (istart+lengrib)) {
167 0 : printf("g2_info: '7777' not found at end of GRIB message.\n");
168 0 : ierr=5;
169 0 : return(ierr);
170 : }
171 30 : if ( isecnum>=2 && isecnum<=7 ) {
172 30 : if (isecnum == 2) // Local Section 2
173 : // increment counter for total number of local sections found
174 0 : (*numlocal)++;
175 :
176 30 : else if (isecnum == 4)
177 : // increment counter for total number of fields found
178 6 : (*numfields)++;
179 : }
180 : else {
181 0 : printf("g2_info: Invalid section number found in GRIB message: %d\n" ,isecnum);
182 0 : ierr=6;
183 0 : return(ierr);
184 : }
185 :
186 30 : }
187 :
188 6 : return(0);
189 :
190 : }
|
