/* version 1.4.3 of grib headers w. ebisuzaki */
/* this version is incomplete */
/* 5/00 - dx/dy or di/dj controlled by bit 1 of resolution byte */
/* 8/00 - dx/dy or di/dj for polar and lambert not controlled by res. byte */
/* Added headers for the triangular grid of the gme model of DWD
Helmut P. Frank, 13.09.2001 */
/* Clean up of triangular grid properties access and added spectral information
Luis Kornblueh, 27.03.2002 */
#ifndef INT3
#define INT3(a,b,c) ((1-(int) ((unsigned) (a & 0x80) >> 6)) * (int) (((a & 127) << 16)+(b<<8)+c))
#endif
#ifndef INT2
#define INT2(a,b) ((1-(int) ((unsigned) (a & 0x80) >> 6)) * (int) (((a & 127) << 8) + b))
#endif
#ifndef UINT4
#define UINT4(a,b,c,d) ((int) ((a << 24) + (b << 16) + (c << 8) + (d)))
#endif
#ifndef UINT3
#define UINT3(a,b,c) ((int) ((a << 16) + (b << 8) + (c)))
#endif
#ifndef UINT2
#define UINT2(a,b) ((int) ((a << 8) + (b)))
#endif
#define GDS_Len1(gds) (gds[0])
#define GDS_Len2(gds) (gds[1])
#define GDS_Len3(gds) (gds[2])
#define GDS_LEN(gds) ((int) ((gds[0]<<16)+(gds[1]<<8)+gds[2]))
#define GDS_NV(gds) (gds[3])
#define GDS_DataType(gds) (gds[5])
#define GDS_LatLon(gds) (gds[5] == 0)
#define GDS_Mercator(gds) (gds[5] == 1)
#define GDS_Gnomonic(gds) (gds[5] == 2)
#define GDS_Lambert(gds) (gds[5] == 3)
#define GDS_Gaussian(gds) (gds[5] == 4)
#define GDS_Polar(gds) (gds[5] == 5)
#define GDS_RotLL(gds) (gds[5] == 10)
#define GDS_Harmonic(gds) (gds[5] == 50)
#define GDS_Triangular(gds) (gds[5] == 192)
#define GDS_ssEgrid(gds) (gds[5] == 201) /* semi-staggered E grid */
#define GDS_fEgrid(gds) (gds[5] == 202) /* filled E grid */
#define GDS_ss2dEgrid(gds) (gds[5] == 203) /* semi-staggered E grid 2 d*/
#define GDS_has_dy(mode) ((mode) & 128)
#define GDS_LatLon_nx(gds) ((int) ((gds[6] << 8) + gds[7]))
#define GDS_LatLon_ny(gds) ((int) ((gds[8] << 8) + gds[9]))
#define GDS_LatLon_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_LatLon_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_LatLon_mode(gds) (gds[16])
#define GDS_LatLon_La2(gds) INT3(gds[17],gds[18],gds[19])
#define GDS_LatLon_Lo2(gds) INT3(gds[20],gds[21],gds[22])
#define GDS_LatLon_dx(gds) (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
#define GDS_LatLon_dy(gds) (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
#define GDS_Gaussian_nlat(gds) ((gds[25]<<8)+gds[26])
#define GDS_LatLon_scan(gds) (gds[27])
#define GDS_Polar_nx(gds) ((gds[6] << 8) + gds[7])
#define GDS_Polar_ny(gds) ((gds[8] << 8) + gds[9])
#define GDS_Polar_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_Polar_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_Polar_mode(gds) (gds[16])
#define GDS_Polar_Lov(gds) INT3(gds[17],gds[18],gds[19])
#define GDS_Polar_scan(gds) (gds[27])
#define GDS_Polar_Dx(gds) INT3(gds[20], gds[21], gds[22])
#define GDS_Polar_Dy(gds) INT3(gds[23], gds[24], gds[25])
#define GDS_Polar_pole(gds) ((gds[26] & 128) == 128)
#define GDS_Lambert_nx(gds) ((gds[6] << 8) + gds[7])
#define GDS_Lambert_ny(gds) ((gds[8] << 8) + gds[9])
#define GDS_Lambert_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_Lambert_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_Lambert_mode(gds) (gds[16])
#define GDS_Lambert_Lov(gds) INT3(gds[17],gds[18],gds[19])
#define GDS_Lambert_dx(gds) INT3(gds[20],gds[21],gds[22])
#define GDS_Lambert_dy(gds) INT3(gds[23],gds[24],gds[25])
#define GDS_Lambert_NP(gds) ((gds[26] & 128) == 0)
#define GDS_Lambert_scan(gds) (gds[27])
#define GDS_Lambert_Latin1(gds) INT3(gds[28],gds[29],gds[30])
#define GDS_Lambert_Latin2(gds) INT3(gds[31],gds[32],gds[33])
#define GDS_Lambert_LatSP(gds) INT3(gds[34],gds[35],gds[36])
#define GDS_Lambert_LonSP(gds) INT3(gds[37],gds[37],gds[37])
#define GDS_ssEgrid_n(gds) UINT2(gds[6],gds[7])
#define GDS_ssEgrid_n_dum(gds) UINT2(gds[8],gds[9])
#define GDS_ssEgrid_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_ssEgrid_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_ssEgrid_mode(gds) (gds[16])
#define GDS_ssEgrid_La2(gds) UINT3(gds[17],gds[18],gds[19])
#define GDS_ssEgrid_Lo2(gds) UINT3(gds[20],gds[21],gds[22])
#define GDS_ssEgrid_di(gds) (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
#define GDS_ssEgrid_dj(gds) (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
#define GDS_ssEgrid_scan(gds) (gds[27])
#define GDS_fEgrid_n(gds) UINT2(gds[6],gds[7])
#define GDS_fEgrid_n_dum(gds) UINT2(gds[8],gds[9])
#define GDS_fEgrid_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_fEgrid_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_fEgrid_mode(gds) (gds[16])
#define GDS_fEgrid_La2(gds) UINT3(gds[17],gds[18],gds[19])
#define GDS_fEgrid_Lo2(gds) UINT3(gds[20],gds[21],gds[22])
#define GDS_fEgrid_di(gds) (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
#define GDS_fEgrid_dj(gds) (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
#define GDS_fEgrid_scan(gds) (gds[27])
#define GDS_ss2dEgrid_nx(gds) UINT2(gds[6],gds[7])
#define GDS_ss2dEgrid_ny(gds) UINT2(gds[8],gds[9])
#define GDS_ss2dEgrid_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_ss2dEgrid_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_ss2dEgrid_mode(gds) (gds[16])
#define GDS_ss2dEgrid_La2(gds) INT3(gds[17],gds[18],gds[19])
#define GDS_ss2dEgrid_Lo2(gds) INT3(gds[20],gds[21],gds[22])
#define GDS_ss2dEgrid_di(gds) (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
#define GDS_ss2dEgrid_dj(gds) (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
#define GDS_ss2dEgrid_scan(gds) (gds[27])
#define GDS_Merc_nx(gds) UINT2(gds[6],gds[7])
#define GDS_Merc_ny(gds) UINT2(gds[8],gds[9])
#define GDS_Merc_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_Merc_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_Merc_mode(gds) (gds[16])
#define GDS_Merc_La2(gds) INT3(gds[17],gds[18],gds[19])
#define GDS_Merc_Lo2(gds) INT3(gds[20],gds[21],gds[22])
#define GDS_Merc_Latin(gds) INT3(gds[23],gds[24],gds[25])
#define GDS_Merc_scan(gds) (gds[27])
#define GDS_Merc_dx(gds) (gds[16] & 128 ? INT3(gds[28],gds[29],gds[30]) : 0)
#define GDS_Merc_dy(gds) (gds[16] & 128 ? INT3(gds[31],gds[32],gds[33]) : 0)
/* rotated Lat-lon grid */
#define GDS_RotLL_nx(gds) UINT2(gds[6],gds[7])
#define GDS_RotLL_ny(gds) UINT2(gds[8],gds[9])
#define GDS_RotLL_La1(gds) INT3(gds[10],gds[11],gds[12])
#define GDS_RotLL_Lo1(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_RotLL_mode(gds) (gds[16])
#define GDS_RotLL_La2(gds) INT3(gds[17],gds[18],gds[19])
#define GDS_RotLL_Lo2(gds) INT3(gds[20],gds[21],gds[22])
#define GDS_RotLL_dx(gds) (gds[16] & 128 ? INT2(gds[23],gds[24]) : 0)
#define GDS_RotLL_dy(gds) (gds[16] & 128 ? INT2(gds[25],gds[26]) : 0)
#define GDS_RotLL_scan(gds) (gds[27])
#define GDS_RotLL_LaSP(gds) INT3(gds[32],gds[33],gds[34])
#define GDS_RotLL_LoSP(gds) INT3(gds[35],gds[36],gds[37])
#define GDS_RotLL_RotAng(gds) ibm2flt(&(gds[38]))
/* Triangular grid of DWD */
#define GDS_Triangular_ni2(gds) INT2(gds[6],gds[7])
#define GDS_Triangular_ni3(gds) INT2(gds[8],gds[9])
#define GDS_Triangular_ni(gds) INT3(gds[13],gds[14],gds[15])
#define GDS_Triangular_nd(gds) INT3(gds[10],gds[11],gds[12])
/* Harmonics data */
#define GDS_Harmonic_nj(gds) ((int) ((gds[6] << 8) + gds[7]))
#define GDS_Harmonic_nk(gds) ((int) ((gds[8] << 8) + gds[9]))
#define GDS_Harmonic_nm(gds) ((int) ((gds[10] << 8) + gds[11]))
#define GDS_Harmonic_type(gds) (gds[12])
#define GDS_Harmonic_mode(gds) (gds[13])
/* index of NV and PV */
#define GDS_PV(gds) ((gds[3] == 0) ? -1 : (int) gds[4] - 1)
#define GDS_PL(gds) ((gds[4] == 255) ? -1 : (int) gds[3] * 4 + (int) gds[4] - 1)