source: trunk/WRF.COMMON/WRFV3/external/io_grib1/WGRIB/gds.h @ 2759

Last change on this file since 2759 was 2759, checked in by aslmd, 2 years ago

adding unmodified code from WRFV3.0.1.1, expurged from useless data +1M size

File size: 7.3 KB
Line 
1/* version 1.4.3 of grib headers  w. ebisuzaki */
2/* this version is incomplete */
3/* 5/00 - dx/dy or di/dj controlled by bit 1 of resolution byte */
4/* 8/00 - dx/dy or di/dj for polar and lambert not controlled by res. byte */
5/* Added headers for the triangular grid of the gme model of DWD
6         Helmut P. Frank, 13.09.2001 */
7/* Clean up of triangular grid properties access and added spectral information
8         Luis Kornblueh, 27.03.2002 */
9
10#ifndef INT3
11#define INT3(a,b,c) ((1-(int) ((unsigned) (a & 0x80) >> 6)) * (int) (((a & 127) << 16)+(b<<8)+c))
12#endif
13#ifndef INT2
14#define INT2(a,b)   ((1-(int) ((unsigned) (a & 0x80) >> 6)) * (int) (((a & 127) << 8) + b))
15#endif
16
17#ifndef UINT4
18#define UINT4(a,b,c,d) ((int) ((a << 24) + (b << 16) + (c << 8) + (d)))
19#endif
20
21#ifndef UINT3
22#define UINT3(a,b,c) ((int) ((a << 16) + (b << 8) + (c)))
23#endif
24
25#ifndef UINT2
26#define UINT2(a,b) ((int) ((a << 8) + (b)))
27#endif
28
29
30#define GDS_Len1(gds)           (gds[0])
31#define GDS_Len2(gds)           (gds[1])
32#define GDS_Len3(gds)           (gds[2])
33#define GDS_LEN(gds)            ((int) ((gds[0]<<16)+(gds[1]<<8)+gds[2]))
34
35#define GDS_NV(gds)             (gds[3])
36#define GDS_DataType(gds)       (gds[5])
37
38#define GDS_LatLon(gds)         (gds[5] == 0)
39#define GDS_Mercator(gds)       (gds[5] == 1)
40#define GDS_Gnomonic(gds)       (gds[5] == 2)
41#define GDS_Lambert(gds)        (gds[5] == 3)
42#define GDS_Gaussian(gds)       (gds[5] == 4)
43#define GDS_Polar(gds)          (gds[5] == 5)
44#define GDS_RotLL(gds)          (gds[5] == 10)
45#define GDS_Harmonic(gds)       (gds[5] == 50)
46#define GDS_Triangular(gds)     (gds[5] == 192)
47#define GDS_ssEgrid(gds)        (gds[5] == 201) /* semi-staggered E grid */
48#define GDS_fEgrid(gds)         (gds[5] == 202) /* filled E grid */
49#define GDS_ss2dEgrid(gds)      (gds[5] == 203) /* semi-staggered E grid 2 d*/
50
51#define GDS_has_dy(mode)        ((mode) & 128)
52#define GDS_LatLon_nx(gds)      ((int) ((gds[6] << 8) + gds[7]))
53#define GDS_LatLon_ny(gds)      ((int) ((gds[8] << 8) + gds[9]))
54#define GDS_LatLon_La1(gds)     INT3(gds[10],gds[11],gds[12])
55#define GDS_LatLon_Lo1(gds)     INT3(gds[13],gds[14],gds[15])
56#define GDS_LatLon_mode(gds)    (gds[16])
57#define GDS_LatLon_La2(gds)     INT3(gds[17],gds[18],gds[19])
58#define GDS_LatLon_Lo2(gds)     INT3(gds[20],gds[21],gds[22])
59
60#define GDS_LatLon_dx(gds)      (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
61#define GDS_LatLon_dy(gds)      (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
62#define GDS_Gaussian_nlat(gds)  ((gds[25]<<8)+gds[26])
63
64#define GDS_LatLon_scan(gds)    (gds[27])
65
66#define GDS_Polar_nx(gds)       ((gds[6] << 8) + gds[7])
67#define GDS_Polar_ny(gds)       ((gds[8] << 8) + gds[9])
68#define GDS_Polar_La1(gds)      INT3(gds[10],gds[11],gds[12])
69#define GDS_Polar_Lo1(gds)      INT3(gds[13],gds[14],gds[15])
70#define GDS_Polar_mode(gds)     (gds[16])
71#define GDS_Polar_Lov(gds)      INT3(gds[17],gds[18],gds[19])
72#define GDS_Polar_scan(gds)     (gds[27])
73#define GDS_Polar_Dx(gds)       INT3(gds[20], gds[21], gds[22])
74#define GDS_Polar_Dy(gds)       INT3(gds[23], gds[24], gds[25])
75#define GDS_Polar_pole(gds)     ((gds[26] & 128) == 128)
76
77#define GDS_Lambert_nx(gds)     ((gds[6] << 8) + gds[7])
78#define GDS_Lambert_ny(gds)     ((gds[8] << 8) + gds[9])
79#define GDS_Lambert_La1(gds)    INT3(gds[10],gds[11],gds[12])
80#define GDS_Lambert_Lo1(gds)    INT3(gds[13],gds[14],gds[15])
81#define GDS_Lambert_mode(gds)   (gds[16])
82#define GDS_Lambert_Lov(gds)    INT3(gds[17],gds[18],gds[19])
83#define GDS_Lambert_dx(gds)     INT3(gds[20],gds[21],gds[22])
84#define GDS_Lambert_dy(gds)     INT3(gds[23],gds[24],gds[25])
85#define GDS_Lambert_NP(gds)     ((gds[26] & 128) == 0)
86#define GDS_Lambert_scan(gds)   (gds[27])
87#define GDS_Lambert_Latin1(gds) INT3(gds[28],gds[29],gds[30])
88#define GDS_Lambert_Latin2(gds) INT3(gds[31],gds[32],gds[33])
89#define GDS_Lambert_LatSP(gds)  INT3(gds[34],gds[35],gds[36])
90#define GDS_Lambert_LonSP(gds)  INT3(gds[37],gds[37],gds[37])
91
92#define GDS_ssEgrid_n(gds)      UINT2(gds[6],gds[7])
93#define GDS_ssEgrid_n_dum(gds)  UINT2(gds[8],gds[9])
94#define GDS_ssEgrid_La1(gds)    INT3(gds[10],gds[11],gds[12])
95#define GDS_ssEgrid_Lo1(gds)    INT3(gds[13],gds[14],gds[15])
96#define GDS_ssEgrid_mode(gds)   (gds[16])
97#define GDS_ssEgrid_La2(gds)    UINT3(gds[17],gds[18],gds[19])
98#define GDS_ssEgrid_Lo2(gds)    UINT3(gds[20],gds[21],gds[22])
99#define GDS_ssEgrid_di(gds)     (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
100#define GDS_ssEgrid_dj(gds)     (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
101#define GDS_ssEgrid_scan(gds)   (gds[27])
102
103#define GDS_fEgrid_n(gds)       UINT2(gds[6],gds[7])
104#define GDS_fEgrid_n_dum(gds)   UINT2(gds[8],gds[9])
105#define GDS_fEgrid_La1(gds)     INT3(gds[10],gds[11],gds[12])
106#define GDS_fEgrid_Lo1(gds)     INT3(gds[13],gds[14],gds[15])
107#define GDS_fEgrid_mode(gds)    (gds[16])
108#define GDS_fEgrid_La2(gds)     UINT3(gds[17],gds[18],gds[19])
109#define GDS_fEgrid_Lo2(gds)     UINT3(gds[20],gds[21],gds[22])
110#define GDS_fEgrid_di(gds)      (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
111#define GDS_fEgrid_dj(gds)      (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
112#define GDS_fEgrid_scan(gds)    (gds[27])
113
114#define GDS_ss2dEgrid_nx(gds)     UINT2(gds[6],gds[7])
115#define GDS_ss2dEgrid_ny(gds)     UINT2(gds[8],gds[9])
116#define GDS_ss2dEgrid_La1(gds)    INT3(gds[10],gds[11],gds[12])
117#define GDS_ss2dEgrid_Lo1(gds)    INT3(gds[13],gds[14],gds[15])
118#define GDS_ss2dEgrid_mode(gds)   (gds[16])
119#define GDS_ss2dEgrid_La2(gds)    INT3(gds[17],gds[18],gds[19])
120#define GDS_ss2dEgrid_Lo2(gds)    INT3(gds[20],gds[21],gds[22])
121#define GDS_ss2dEgrid_di(gds)     (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
122#define GDS_ss2dEgrid_dj(gds)     (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
123#define GDS_ss2dEgrid_scan(gds)   (gds[27])
124
125
126#define GDS_Merc_nx(gds)        UINT2(gds[6],gds[7])
127#define GDS_Merc_ny(gds)        UINT2(gds[8],gds[9])
128#define GDS_Merc_La1(gds)       INT3(gds[10],gds[11],gds[12])
129#define GDS_Merc_Lo1(gds)       INT3(gds[13],gds[14],gds[15])
130#define GDS_Merc_mode(gds)      (gds[16])
131#define GDS_Merc_La2(gds)       INT3(gds[17],gds[18],gds[19])
132#define GDS_Merc_Lo2(gds)       INT3(gds[20],gds[21],gds[22])
133#define GDS_Merc_Latin(gds)     INT3(gds[23],gds[24],gds[25])
134#define GDS_Merc_scan(gds)      (gds[27])
135#define GDS_Merc_dx(gds)        (gds[16] & 128 ? INT3(gds[28],gds[29],gds[30]) : 0)
136#define GDS_Merc_dy(gds)        (gds[16] & 128 ? INT3(gds[31],gds[32],gds[33]) : 0)
137
138/* rotated Lat-lon grid */
139
140#define GDS_RotLL_nx(gds)       UINT2(gds[6],gds[7])
141#define GDS_RotLL_ny(gds)       UINT2(gds[8],gds[9])
142#define GDS_RotLL_La1(gds)      INT3(gds[10],gds[11],gds[12])
143#define GDS_RotLL_Lo1(gds)      INT3(gds[13],gds[14],gds[15])
144#define GDS_RotLL_mode(gds)     (gds[16])
145#define GDS_RotLL_La2(gds)      INT3(gds[17],gds[18],gds[19])
146#define GDS_RotLL_Lo2(gds)      INT3(gds[20],gds[21],gds[22])
147#define GDS_RotLL_dx(gds)       (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
148#define GDS_RotLL_dy(gds)       (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
149#define GDS_RotLL_scan(gds)     (gds[27])
150#define GDS_RotLL_LaSP(gds)     INT3(gds[32],gds[33],gds[34])
151#define GDS_RotLL_LoSP(gds)     INT3(gds[35],gds[36],gds[37])
152#define GDS_RotLL_RotAng(gds)   ibm2flt(&(gds[38]))
153
154/* Triangular grid of DWD */
155#define GDS_Triangular_ni2(gds) INT2(gds[6],gds[7])
156#define GDS_Triangular_ni3(gds) INT2(gds[8],gds[9])
157#define GDS_Triangular_ni(gds)  INT3(gds[13],gds[14],gds[15])
158#define GDS_Triangular_nd(gds)  INT3(gds[10],gds[11],gds[12])
159
160/* Harmonics data */
161#define GDS_Harmonic_nj(gds)     ((int) ((gds[6] << 8) + gds[7]))
162#define GDS_Harmonic_nk(gds)     ((int) ((gds[8] << 8) + gds[9]))
163#define GDS_Harmonic_nm(gds)     ((int) ((gds[10] << 8) + gds[11]))
164#define GDS_Harmonic_type(gds)   (gds[12])
165#define GDS_Harmonic_mode(gds)   (gds[13])
166
167/* index of NV and PV */
168#define GDS_PV(gds)             ((gds[3] == 0) ? -1 : (int) gds[4] - 1)
169#define GDS_PL(gds)             ((gds[4] == 255) ? -1 : (int) gds[3] * 4 + (int) gds[4] - 1)
170
Note: See TracBrowser for help on using the repository browser.