stream_gribapi.c 79.2 KB
Newer Older
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1
2
3
4
5
6
#if defined (HAVE_CONFIG_H)
#  include "config.h"
#endif

#include <stdio.h>

Uwe Schulzweida's avatar
Uwe Schulzweida committed
7
#include "dmemory.h"
Uwe Schulzweida's avatar
Uwe Schulzweida committed
8
9
10
11
12
13
#include "cdi.h"
#include "stream_int.h"
#include "file.h"
#include "varscan.h"
#include "datetime.h"
#include "vlist.h"
Uwe Schulzweida's avatar
Uwe Schulzweida committed
14
#include "stream_grb.h"
15
#include "calendar.h"
Uwe Schulzweida's avatar
Uwe Schulzweida committed
16
17
18
19


#if  defined  (HAVE_LIBGRIB_API)
#  include "cgribex.h"      /* gribGetSize, gribRead, gribGetZip */
Uwe Schulzweida's avatar
Uwe Schulzweida committed
20
#  include "gribapi.h"
Uwe Schulzweida's avatar
Uwe Schulzweida committed
21
22
23
24
25
26
27
#  include "grib_api.h"
#endif


extern int cdiInventoryMode;

typedef struct {
Uwe Schulzweida's avatar
Uwe Schulzweida committed
28
  int param;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
29
30
31
32
33
34
35
36
  int level1;
  int level2;
  int ltype;
} compvar2_t; 


#if  defined  (HAVE_LIBGRIB_API)
static
37
int gribapiGetGridType(grib_handle *gh)
Uwe Schulzweida's avatar
Uwe Schulzweida committed
38
{
39
40
  int gridtype = GRID_GENERIC;
  int gribgridtype;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
41
42
  long lpar;

43
44
45
46
47
48
    {
      GRIB_CHECK(grib_get_long(gh, "gridDefinitionTemplateNumber", &lpar), 0);
      gribgridtype = (int) lpar;

      switch (gribgridtype)
	{
49
	case  GRIB2_GTYPE_LATLON:     { GRIB_CHECK(grib_get_long(gh, "Ni", &lpar), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
50
	                                if ( lpar == (long) GRIB_MISSING_LONG ) break;
51
                                      }
52
53
	case  GRIB2_GTYPE_LATLON_ROT: { gridtype = GRID_LONLAT;    break; }
	case  GRIB2_GTYPE_LCC:        { gridtype = GRID_LCC;       break; }
54
	case  GRIB2_GTYPE_GAUSSIAN:   { GRIB_CHECK(grib_get_long(gh, "Ni", &lpar), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
55
	                                if ( lpar == (long) GRIB_MISSING_LONG )
56
57
58
59
60
					  gridtype = GRID_GAUSSIAN_REDUCED;
					else
					  gridtype = GRID_GAUSSIAN;
					break;
                                      }
61
62
63
	case  GRIB2_GTYPE_SPECTRAL:   { gridtype = GRID_SPECTRAL;  break; }
	case  GRIB2_GTYPE_GME:        { gridtype = GRID_GME;       break; }
	case  GRIB2_GTYPE_NUMBER:     { gridtype = GRID_REFERENCE; break; }
64
	}
Uwe Schulzweida's avatar
Uwe Schulzweida committed
65
66
67
68
69
70
    }

  return (gridtype);
}
#endif

Uwe Schulzweida's avatar
Uwe Schulzweida committed
71
#if  defined  (HAVE_LIBGRIB_API)
72
static
73
int gribapiGetIsRotated(grib_handle *gh)
74
75
{
  int isRotated = 0;
76
77
  int gribgridtype;
  long lpar;
78
79

    {
80
81
82
83
      GRIB_CHECK(grib_get_long(gh, "gridDefinitionTemplateNumber", &lpar), 0);
      gribgridtype = (int) lpar;

      if ( gribgridtype == GRIB2_GTYPE_LATLON_ROT ) isRotated = 1;
84
85
86
87
    }

  return (isRotated);
}
Uwe Schulzweida's avatar
Uwe Schulzweida committed
88
#endif
89

Uwe Schulzweida's avatar
Uwe Schulzweida committed
90
91
#if  defined  (HAVE_LIBGRIB_API)
static
92
int gribapiGetZaxisType(long editionNumber, int grib_ltype)
Uwe Schulzweida's avatar
Uwe Schulzweida committed
93
{
Uwe Schulzweida's avatar
Uwe Schulzweida committed
94
95
96
  int zaxistype = ZAXIS_GENERIC;

  if ( editionNumber <= 1 )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
97
    {
98
      zaxistype = grib1ltypeToZaxisType(grib_ltype);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
99
100
101
    }
  else
    {
102
      zaxistype = grib2ltypeToZaxisType(grib_ltype);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
103
104
105
106
107
108
    }

  return (zaxistype);
}
#endif

109
110
111
112
113
#if  defined  (HAVE_LIBGRIB_API)
static
int gribapiGetTimeUnits(grib_handle *gh)
{
  int timeunits = -1;
114
  long unitsOfTime;
115
  long lpar;
116
117
  size_t len = 8;
  char stepunits[8];
118
119
  static int lprint = TRUE;

120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
  GRIB_CHECK(grib_get_long(gh, "indicatorOfUnitOfTimeRange", &unitsOfTime), 0);  

  switch (unitsOfTime)
    {
    case 13:  timeunits = TUNIT_SECOND;  break;
    case  0:  timeunits = TUNIT_MINUTE;  break;
    case  1:  timeunits = TUNIT_HOUR;    break;
    case 10:  timeunits = TUNIT_3HOURS;  break;
    case 11:  timeunits = TUNIT_6HOURS;  break;
    case 12:  timeunits = TUNIT_12HOURS; break;
    case  2:  timeunits = TUNIT_DAY;     break;
    default:  timeunits = TUNIT_HOUR;    break;
    }

  /*
135
136
  GRIB_CHECK(grib_get_string(gh, "stepUnits", stepunits, &len), 0);

Uwe Schulzweida's avatar
Uwe Schulzweida committed
137
138
139
140
141
142
143
144
145
146
147
  len--;

  if      ( memcmp(stepunits, "s",   len) == 0 ) timeunits = TUNIT_SECOND;
  else if ( memcmp(stepunits, "m",   len) == 0 ) timeunits = TUNIT_MINUTE;
  else if ( memcmp(stepunits, "h",   len) == 0 ) timeunits = TUNIT_HOUR;
  else if ( memcmp(stepunits, "3h",  len) == 0 ) timeunits = TUNIT_3HOURS;
  else if ( memcmp(stepunits, "6h",  len) == 0 ) timeunits = TUNIT_6HOURS;
  else if ( memcmp(stepunits, "12h", len) == 0 ) timeunits = TUNIT_12HOURS;
  else if ( memcmp(stepunits, "D",   len) == 0 ) timeunits = TUNIT_DAY;
  else if ( memcmp(stepunits, "M",   len) == 0 ) timeunits = TUNIT_MONTH;
  else if ( memcmp(stepunits, "Y",   len) == 0 ) timeunits = TUNIT_YEAR;
148
  else if ( lprint )
149
    {
150
      Message("Step units >%s< unsupported!", stepunits);
151
      lprint = FALSE;
152
    }
153
  */
154
155
156
157
158
159
160
161
162

  return (timeunits);
}
#endif

#if  defined  (HAVE_LIBGRIB_API)
static
int gribapiTimeIsFC(grib_handle *gh)
{
163
  long editionNumber;
164
165
  int isFC = TRUE;

166
  GRIB_CHECK(grib_get_long(gh, "editionNumber", &editionNumber), 0);
167

168
169
170
171
172
173
174
175
  if ( editionNumber == 2 )
    {
      long sigofrtime;

      GRIB_CHECK(grib_get_long(gh, "significanceOfReferenceTime", &sigofrtime), 0);

      if ( sigofrtime == 3 ) isFC = FALSE;
    }
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191

  return (isFC);
}
#endif

#if  defined  (HAVE_LIBGRIB_API)
static
int gribapiGetTsteptype(grib_handle *gh)
{
  int tsteptype = 0;
  int timerange;
  long lpar;
  static int lprint = TRUE;

  if ( gribapiTimeIsFC(gh) )
    {
192
193
194
195
196
      int status;
      status = grib_get_long(gh, "stepType", &lpar);
      if ( status == 0 )
	{
	  timerange = (int) lpar;
197

198
	  // printf("timerange %d\n", timerange);
199

200
	  switch ( timerange )
201
	    {
202
203
204
205
206
207
208
209
210
211
212
	    case  0:  tsteptype = TSTEP_AVG;    break;
	    case  1:  tsteptype = TSTEP_ACCUM;  break;
	    case  2:  tsteptype = TSTEP_MIN;    break;
	    case  3:  tsteptype = TSTEP_MAX;    break;
	    case  4:  tsteptype = TSTEP_DIFF;   break;
	    default:
	      if ( lprint )
		{
		  Message("Time range %d unsupported", timerange);
		  lprint = FALSE;
		}
213
214
215
216
217
218
219
220
221
	    }
	}
    }

  return (tsteptype);
}
#endif

#if  defined  (HAVE_LIBGRIB_API)
222
static
223
224
void gribapiGetValidityDateTime(grib_handle *gh, int *vdate, int *vtime)
{
225
226
  int rdate, rtime;
  int timeUnits, timePeriod;
227
  long lpar;
228
229
230
231
232
233
234
235
236
  long sigofrtime = 3;
  long editionNumber;

  GRIB_CHECK(grib_get_long(gh, "editionNumber", &editionNumber), 0);

  if ( editionNumber == 2 )
    {
      GRIB_CHECK(grib_get_long(gh, "significanceOfReferenceTime", &sigofrtime), 0);
    }
237
238
239
240
241
242
243
244
245
246

  if ( sigofrtime == 3 )
    {
      GRIB_CHECK(grib_get_long(gh, "dataDate", &lpar), 0);
      *vdate = (int) lpar;
      GRIB_CHECK(grib_get_long(gh, "dataTime", &lpar), 0);
      *vtime = (int) lpar*100;
    }
  else
    {
247
      /*
248
249
250
251
      GRIB_CHECK(grib_get_long(gh, "validityDate", &lpar), 0);
      *vdate = (int) lpar;
      GRIB_CHECK(grib_get_long(gh, "validityTime", &lpar), 0);
      *vtime = (int) lpar*100;
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
      */
      GRIB_CHECK(grib_get_long(gh, "dataDate", &lpar), 0);
      rdate = (int) lpar;
      GRIB_CHECK(grib_get_long(gh, "dataTime", &lpar), 0);
      rtime = (int) lpar*100;    
      GRIB_CHECK(grib_get_long(gh, "forecastTime", &lpar), 0);
      timePeriod = (int) lpar;
      timeUnits = gribapiGetTimeUnits(gh);
      
      {
	static int lprint = TRUE;
	extern int grib_calendar;
	int ryear, rmonth, rday, rhour, rminute, rsecond;
	int time_period = timePeriod;
	int julday, secofday, addsec;

	cdiDecodeDate(rdate, &ryear, &rmonth, &rday);
	cdiDecodeTime(rtime, &rhour, &rminute, &rsecond);

	encode_caldaysec(grib_calendar, ryear, rmonth, rday, rhour, rminute, rsecond, &julday, &secofday);

	addsec = 0;
	switch ( timeUnits )
	  {
	  case TUNIT_SECOND:  addsec =         time_period; break;
	  case TUNIT_MINUTE:  addsec =    60 * time_period; break;
	  case TUNIT_HOUR:    addsec =  3600 * time_period; break;
	  case TUNIT_3HOURS:  addsec = 10800 * time_period; break;
	  case TUNIT_6HOURS:  addsec = 21600 * time_period; break;
	  case TUNIT_12HOURS: addsec = 43200 * time_period; break;
	  case TUNIT_DAY:     addsec = 86400 * time_period; break;
	  default:
	    if ( lprint )
	      {
	        Warning("Time unit %d unsupported", timeUnits);
		lprint = FALSE;
	      }
	  }

	julday_add_seconds(addsec, &julday, &secofday);

	decode_caldaysec(grib_calendar, julday, secofday, &ryear, &rmonth, &rday, &rhour, &rminute, &rsecond);
	/*
	  printf("new %d/%d/%d %d:%d\n", ryear, rmonth, rday, rhour, rminute);
	*/
	*vdate = cdiEncodeDate(ryear, rmonth, rday);
	*vtime = cdiEncodeTime(rhour, rminute, rsecond);
      }
300
301
302
303
    }
}
#endif

Uwe Schulzweida's avatar
Uwe Schulzweida committed
304
#if  defined  (HAVE_LIBGRIB_API)
305
static
306
void gribapiGetGrid(grib_handle *gh, grid_t *grid)
Uwe Schulzweida's avatar
Uwe Schulzweida committed
307
{
308
  long editionNumber;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
309
  int gridtype;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
310
  size_t datasize;
311
312
  long numberOfPoints;
  long lpar;
313
314
315

  GRIB_CHECK(grib_get_long(gh, "editionNumber", &editionNumber), 0);

316
  gridtype = gribapiGetGridType(gh);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
317
318
319
320
321
322
323
324
  /*
  if ( streamptr->unreduced && gridtype == GRID_GAUSSIAN_REDUCED )
    {
      gridtype = GRID_GAUSSIAN;
      ISEC2_NumLon = 2*ISEC2_NumLat;
      ISEC4_NumValues = ISEC2_NumLon*ISEC2_NumLat;
    }
  */
325
  memset(grid, 0, sizeof(grid_t));
Uwe Schulzweida's avatar
Uwe Schulzweida committed
326

Uwe Schulzweida's avatar
Uwe Schulzweida committed
327
328
329
  GRIB_CHECK(grib_get_size(gh, "values", &datasize), 0);
  GRIB_CHECK(grib_get_long(gh, "numberOfPoints", &numberOfPoints), 0);

Uwe Schulzweida's avatar
Uwe Schulzweida committed
330
331
332
333
334
335
336
337
338
339
340
341
  switch (gridtype)
    {
    case GRID_LONLAT:
    case GRID_GAUSSIAN:
      {
	int nlon, nlat;

	GRIB_CHECK(grib_get_long(gh, "Ni", &lpar), 0);
	nlon = lpar;
	GRIB_CHECK(grib_get_long(gh, "Nj", &lpar), 0);
	nlat = lpar;

342
343
344
345
346
347
	if ( gridtype == GRID_GAUSSIAN )
          {
            GRIB_CHECK(grib_get_long(gh, "numberOfParallelsBetweenAPoleAndTheEquator", &lpar), 0);
            grid->np = lpar;
          }

Uwe Schulzweida's avatar
Uwe Schulzweida committed
348
	if ( numberOfPoints != nlon*nlat )
349
	  Error("numberOfPoints (%d) and gridSize (%d) differ!",
Uwe Schulzweida's avatar
Uwe Schulzweida committed
350
		(int)numberOfPoints, nlon*nlat);
351

352
353
354
355
356
357
358
359
360
361
362
	grid->size  = numberOfPoints;
	grid->xsize = nlon;
	grid->ysize = nlat;
	grid->xinc  = 0;
	grid->yinc  = 0;
	grid->xdef  = 0;
	GRIB_CHECK(grib_get_double(gh, "longitudeOfFirstGridPointInDegrees", &grid->xfirst), 0);
	GRIB_CHECK(grib_get_double(gh, "longitudeOfLastGridPointInDegrees",  &grid->xlast), 0);
	GRIB_CHECK(grib_get_double(gh, "latitudeOfFirstGridPointInDegrees",  &grid->yfirst), 0);
	GRIB_CHECK(grib_get_double(gh, "latitudeOfLastGridPointInDegrees",   &grid->ylast), 0);
	GRIB_CHECK(grib_get_double(gh, "iDirectionIncrementInDegrees", &grid->xinc), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
363
	if ( gridtype == GRID_LONLAT )
364
	  GRIB_CHECK(grib_get_double(gh, "jDirectionIncrementInDegrees", &grid->yinc), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
365

366
	if ( IS_EQUAL(grid->xinc, GRIB_MISSING_DOUBLE) ) grid->xinc = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
367

368
	/* if ( IS_NOT_EQUAL(grid->xfirst, 0) || IS_NOT_EQUAL(grid->xlast, 0) ) */
Uwe Schulzweida's avatar
Uwe Schulzweida committed
369
	  {
370
	    if ( grid->xsize > 1 )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
371
	      {
372
		if ( (grid->xfirst >= grid->xlast) && (grid->xfirst >= 180) ) grid->xfirst -= 360;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
373

374
375
376
		if ( editionNumber <= 1 )
		  {
		    /* correct xinc if necessary */
377
		    if ( IS_EQUAL(grid->xfirst, 0) && grid->xlast > 354 )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
378
		      {
379
			double xinc = 360. / grid->xsize;
380

381
			if ( fabs(grid->xinc-xinc) > 0.0 )
382
			  {
383
384
			    grid->xinc = xinc;
			    if ( CDI_Debug ) Message("set xinc to %g", grid->xinc);
385
			  }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
386
387
388
		      }
		  }
	      }
389
	    grid->xdef = 2;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
390
	  }
391
	grid->ydef = 0;
392
        /* if ( IS_NOT_EQUAL(grid->yfirst, 0) || IS_NOT_EQUAL(grid->ylast, 0) ) */
Uwe Schulzweida's avatar
Uwe Schulzweida committed
393
	  {
394
	    if ( grid->ysize > 1 )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
395
	      {
396
397
398
		if ( editionNumber <= 1 )
		  {
		  }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
399
	      }
400
	    grid->ydef = 2;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
401
402
403
404
405
	  }
	break;
      }
    case GRID_GAUSSIAN_REDUCED:
      {
406
407
408
409
	int nlat, i;
	size_t dummy;
	long *pl;

410
411
412
        GRIB_CHECK(grib_get_long(gh, "numberOfParallelsBetweenAPoleAndTheEquator", &lpar), 0);
        grid->np = lpar;

413
414
415
	GRIB_CHECK(grib_get_long(gh, "Nj", &lpar), 0);
	nlat = lpar;

416
	grid->size   = numberOfPoints;
417

418
        grid->rowlon = (int *) malloc(nlat*sizeof(int));
419
420
421
        pl          = (long *) malloc(nlat*sizeof(long));
	dummy       = nlat;
	GRIB_CHECK(grib_get_long_array(gh, "pl", pl, &dummy), 0);
422
	for ( i = 0; i < nlat; ++i ) grid->rowlon[i] = pl[i];
423
424
	free(pl);

425
426
427
428
429
430
431
432
433
	grid->ysize  = nlat;
	grid->xinc   = 0;
	grid->yinc   = 0;
	grid->xdef   = 0;
	GRIB_CHECK(grib_get_double(gh, "longitudeOfFirstGridPointInDegrees", &grid->xfirst), 0);
	GRIB_CHECK(grib_get_double(gh, "longitudeOfLastGridPointInDegrees",  &grid->xlast), 0);
	GRIB_CHECK(grib_get_double(gh, "latitudeOfFirstGridPointInDegrees",  &grid->yfirst), 0);
	GRIB_CHECK(grib_get_double(gh, "latitudeOfLastGridPointInDegrees",   &grid->ylast), 0);
	GRIB_CHECK(grib_get_double(gh, "iDirectionIncrementInDegrees", &grid->xinc), 0);
434

435
	if ( IS_EQUAL(grid->xinc, GRIB_MISSING_DOUBLE) ) grid->xinc = 0;
436

437
	/* if ( IS_NOT_EQUAL(grid->xfirst, 0) || IS_NOT_EQUAL(grid->xlast, 0) ) */
Uwe Schulzweida's avatar
Uwe Schulzweida committed
438
	  {
439
	    if ( grid->xsize > 1 )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
440
	      {
441
		if ( (grid->xfirst > grid->xlast) && (grid->xfirst >= 180) ) grid->xfirst -= 360;
442
443
444
445

		if ( editionNumber <= 1 )
		  {
		    /* correct xinc if necessary */
446
		    if ( IS_EQUAL(grid->xfirst, 0) && grid->xlast > 354 )
447
		      {
448
			double xinc = 360. / grid->xsize;
449

450
			if ( fabs(grid->xinc-xinc) > 0.0 )
451
			  {
452
453
			    grid->xinc = xinc;
			    if ( CDI_Debug ) Message("set xinc to %g", grid->xinc);
454
455
456
			  }
		      }
		  }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
457
	      }
458
	    grid->xdef = 2;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
459
	  }
460
461
	grid->ydef  = 0;
        /* if ( IS_NOT_EQUAL(grid->yfirst, 0) || IS_NOT_EQUAL(grid->ylast, 0) ) */
Uwe Schulzweida's avatar
Uwe Schulzweida committed
462
	  {
463
	    if ( grid->ysize > 1 )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
464
	      {
465
466
467
		if ( editionNumber <= 1 )
		  {
		  }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
468
	      }
469
	    grid->ydef = 2;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
470
471
472
	  }
	break;
      }
473
      /*
Uwe Schulzweida's avatar
Uwe Schulzweida committed
474
475
476
    case GRID_LCC:
      {
	if ( ISEC4_NumValues != ISEC2_NumLon*ISEC2_NumLat )
477
	  Error("numberOfPoints (%d) and gridSize (%d) differ!",
Uwe Schulzweida's avatar
Uwe Schulzweida committed
478
479
		ISEC4_NumValues, ISEC2_NumLon*ISEC2_NumLat);

480
481
482
	grid->size  = ISEC4_NumValues;
	grid->xsize = ISEC2_NumLon;
	grid->ysize = ISEC2_NumLat;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
483

484
485
486
487
488
489
490
491
492
	grid->lcc_xinc      = ISEC2_Lambert_dx;
	grid->lcc_yinc      = ISEC2_Lambert_dy;
	grid->lcc_originLon = ISEC2_FirstLon * 0.001;
	grid->lcc_originLat = ISEC2_FirstLat * 0.001;
	grid->lcc_lonParY   = ISEC2_Lambert_Lov * 0.001;
	grid->lcc_lat1      = ISEC2_Lambert_LatS1 * 0.001;
	grid->lcc_lat2      = ISEC2_Lambert_LatS2 * 0.001;
	grid->lcc_projflag  = ISEC2_Lambert_ProjFlag;
	grid->lcc_scanflag  = ISEC2_ScanFlag;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
493

494
	grid->xdef   = 0;
495
	grid->ydef   = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
496
497
498

	break;
      }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
499
      */
Uwe Schulzweida's avatar
Uwe Schulzweida committed
500
501
    case GRID_SPECTRAL:
      {
502
503
	size_t len = 256;
	char typeOfPacking[256];
504
	GRIB_CHECK(grib_get_string(gh, "packingType", typeOfPacking, &len), 0);
505
506
	grid->lcomplex = 0;
	if ( strncmp(typeOfPacking, "spectral_complex", len) == 0 ) grid->lcomplex = 1;
507

508
	grid->size  = datasize;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
509
	GRIB_CHECK(grib_get_long(gh, "J", &lpar), 0);
510
	grid->trunc = lpar;
511

Uwe Schulzweida's avatar
Uwe Schulzweida committed
512
513
514
515
	break;
      }
    case GRID_GME:
      {
516
517
518
519
520
	grid->size  = numberOfPoints;
	if ( grib_get_long(gh, "nd", &lpar) == 0 ) grid->nd  = lpar;
	if ( grib_get_long(gh, "Ni", &lpar) == 0 ) grid->ni  = lpar;
	if ( grib_get_long(gh, "n2", &lpar) == 0 ) grid->ni2 = lpar;
	if ( grib_get_long(gh, "n3", &lpar) == 0 ) grid->ni3 = lpar;
521
522
523

	break;
      }
524
    case GRID_REFERENCE:
525
      {
526
        char uuid[17];
527
528
529
530
	char reference_link[8192];
	size_t len = sizeof(reference_link);
	reference_link[0] = 0;

531
	grid->size  = numberOfPoints;
532
533
	if ( grib_get_long(gh, "numberOfGridUsed", &lpar) == 0 )
	  {
534
535
	    grid->number   = lpar;
	    if ( grib_get_long(gh, "numberOfGridInReference", &lpar) == 0 ) grid->position = lpar;
536
537
	    if ( grib_get_string(gh, "gridDescriptionFile", reference_link, &len) == 0 )
	      {
538
		if ( strncmp(reference_link, "file://", 7) == 0 )
539
		  grid->reference = strdupx(reference_link);
540
	      }
541
542
543
544
545
            len = (size_t) 16;
            if ( grib_get_string(gh, "uuidOfHGrid", uuid, &len) == 0)
              {
                strncpy(grid->uuid, uuid, 16);
              }
546
	  }
547

Uwe Schulzweida's avatar
Uwe Schulzweida committed
548
549
550
551
	break;
      }
    case GRID_GENERIC:
      {
552
553
554
555
556
	int nlon = 0, nlat = 0;

	if ( grib_get_long(gh, "Ni", &lpar) == 0 ) nlon = lpar;
	if ( grib_get_long(gh, "Nj", &lpar) == 0 ) nlat = lpar;

557
	grid->size  = numberOfPoints;
558

559
	if ( nlon > 0 && nlat > 0 && nlon*nlat == grid->size )
560
	  {
561
562
	    grid->xsize = nlon;
	    grid->ysize = nlat;
563
564
565
	  }
	else
	  {
566
567
	    grid->xsize = 0;
	    grid->ysize = 0;
568
569
	  }

Uwe Schulzweida's avatar
Uwe Schulzweida committed
570
571
572
573
	break;
      }
    default:
      {
574
	Error("Unsupported grid type: %s", gridNamePtr(gridtype));
Uwe Schulzweida's avatar
Uwe Schulzweida committed
575
576
577
578
	break;
      }
    }

579
  grid->isRotated = FALSE;
580
  if ( gribapiGetIsRotated(gh) )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
581
    {
582
583
584
585
      grid->isRotated = TRUE;
      GRIB_CHECK(grib_get_double(gh, "latitudeOfSouthernPoleInDegrees",  &grid->ypole), 0);
      GRIB_CHECK(grib_get_double(gh, "longitudeOfSouthernPoleInDegrees", &grid->xpole), 0);
      GRIB_CHECK(grib_get_double(gh, "angleOfRotation", &grid->angle), 0);
586
      /* change from south to north pole */
587
588
589
590
591
592
593
594
595
596
597
598
      grid->ypole = -grid->ypole;
      grid->xpole =  grid->xpole - 180;
    }

  grid->xvals = NULL;
  grid->yvals = NULL;
  grid->type  = gridtype;
}
#endif

#if  defined  (HAVE_LIBGRIB_API)
static
599
int gribapiGetZaxisHasBounds(grib_handle *gh)
600
{
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
  int lbounds = 0;
  long editionNumber;

  GRIB_CHECK(grib_get_long(gh, "editionNumber", &editionNumber), 0);

  if ( editionNumber <= 1 )
    {
      int status, grb1_ltype;
      long lpar;
      status = grib_get_long(gh, "indicatorOfTypeOfLevel", &lpar);
      if ( status == 0 ) grb1_ltype = (int) lpar;

      switch (grb1_ltype)
	{
	case GRIB1_LTYPE_SIGMA_LAYER:
	case GRIB1_LTYPE_HYBRID_LAYER:
	case GRIB1_LTYPE_LANDDEPTH_LAYER:
	  {
	    lbounds = 1;
	    break;
	  }
	}
    }
  else
    {
      int status;
      long lpar, ltype1 = -1, ltype2 = -2;
628

629
630
631
632
633
634
635
      status = grib_get_long(gh, "typeOfFirstFixedSurface", &lpar);
      if ( status == 0 ) ltype1 = lpar;
      status = grib_get_long(gh, "typeOfSecondFixedSurface", &lpar);
      if ( status == 0 ) ltype2 = lpar;

      if ( ltype1 == ltype2 && ltype1 != 255 ) lbounds = 1;
    }
636

637
  return (lbounds);
638
639
640
}

static
641
void grib1GetLevel(grib_handle *gh, int leveltype, int lbounds, int *level1, int *level2)
642
643
644
{
  double dlevel;

645
  if ( lbounds == 0 )
646
    {
647
648
649
650
651
      GRIB_CHECK(grib_get_double(gh, "level", &dlevel), 0);
      if ( leveltype == 100 ) dlevel *= 100;
      if ( dlevel < -2.e9 || dlevel > 2.e9 ) dlevel = 0;
      *level1 = (int) dlevel;
      *level2 = 0;
652
653
654
    }
  else
    {
655
656
657
658
659
      long lpar;
      GRIB_CHECK(grib_get_long(gh, "topLevel", &lpar), 0);
      *level1 = lpar;
      GRIB_CHECK(grib_get_long(gh, "bottomLevel", &lpar), 0);
      *level2 = lpar;
660
    }
661
662
663
664
665
666
667
}

static
void grib2GetLevel(grib_handle *gh, int leveltype, int lbounds, int *level1, int *level2)
{
  double dlevel;
  long factor;
668

669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
  if ( lbounds == 0 )
    {
      if ( leveltype == GRIB2_LTYPE_LANDDEPTH )
	{
	  GRIB_CHECK(grib_get_long(gh, "scaleFactorOfFirstFixedSurface", &factor), 0);
	  GRIB_CHECK(grib_get_double(gh, "scaledValueOfFirstFixedSurface", &dlevel), 0);
	  if      ( factor == 0 ) dlevel *= 100;   //  m to cm
	  else if ( factor == 1 ) dlevel *=  10;   // dm to cm
	  else if ( factor == 3 ) dlevel *=   0.1; // mm to cm
	}
      else
	{
	  GRIB_CHECK(grib_get_double(gh, "level", &dlevel), 0);
	  if ( leveltype == GRIB2_LTYPE_ISOBARIC ) dlevel *= 100;
	  if ( dlevel < -2.e9 || dlevel > 2.e9 ) dlevel = 0;
	}

      *level1 = (int) dlevel;
      *level2 = 0;
    }
  else
    {
      if ( leveltype == GRIB2_LTYPE_LANDDEPTH )
	{
	  GRIB_CHECK(grib_get_long(gh, "scaleFactorOfFirstFixedSurface", &factor), 0);
	  GRIB_CHECK(grib_get_double(gh, "scaledValueOfFirstFixedSurface", &dlevel), 0);
	  if      ( factor == 0 ) dlevel *= 100;   //  m to cm
	  else if ( factor == 1 ) dlevel *=  10;   // dm to cm
	  else if ( factor == 3 ) dlevel *=   0.1; // mm to cm
	  *level1 = (int) dlevel;
	  GRIB_CHECK(grib_get_long(gh, "scaleFactorOfSecondFixedSurface", &factor), 0);
	  GRIB_CHECK(grib_get_double(gh, "scaledValueOfSecondFixedSurface", &dlevel), 0);
	  if      ( factor == 0 ) dlevel *= 100;   //  m to cm
	  else if ( factor == 1 ) dlevel *=  10;   // dm to cm
	  else if ( factor == 3 ) dlevel *=   0.1; // mm to cm
	  *level2 = (int) dlevel;
	}
      else
	{
	  long lpar;
	  GRIB_CHECK(grib_get_long(gh, "topLevel", &lpar), 0);
	  *level1 = lpar;
	  GRIB_CHECK(grib_get_long(gh, "bottomLevel", &lpar), 0);
	  *level2 = lpar;
	}
    }
715
716
717
718
}

static
void gribapiAddRecord(int streamID, int param, grib_handle *gh,
719
		      long recsize, off_t position, int datatype, int comptype)
720
721
722
723
724
725
{
  long editionNumber;
  int zaxistype;
  int gridID = CDI_UNDEFID, varID;
  int levelID = 0;
  int tsID, recID;
726
  int level1 = 0, level2 = 0;
727
728
729
730
731
732
733
734
735
736
737
  int numavg;
  int tsteptype;
  int lbounds = 0;
  record_t *record;
  grid_t grid;
  int vlistID;
  stream_t *streamptr;
  int leveltype;
  long lpar;
  int status;
  char name[256], longname[256], units[256];
738
  size_t vlen;
739
740
741
742
743
744
745
746
747
748
749
750
751
752

  streamptr = stream_to_pointer(streamID);

  vlistID = streamInqVlist(streamID);
  tsID    = streamptr->curTsID;
  recID   = recordNewEntry(streamID, tsID);
  record  = &streamptr->tsteps[tsID].records[recID];

  tsteptype = gribapiGetTsteptype(gh);
  // numavg  = ISEC1_AvgNum;
  numavg  = 0;

  GRIB_CHECK(grib_get_long(gh, "editionNumber", &editionNumber), 0);

753
754
  lbounds = gribapiGetZaxisHasBounds(gh);

755
756
757
758
759
760
  if ( editionNumber <= 1 )
    {
      status = grib_get_long(gh, "indicatorOfTypeOfLevel", &lpar);
      if ( status == 0 )
	{
	  leveltype = (int) lpar;
761
762
	  grib1GetLevel(gh, leveltype, lbounds, &level1, &level2);
	  if ( leveltype == 99 ) leveltype = 100;
763
	}
764
      else
765
766
767
768
769
770
771
772
773
774
	{
	  leveltype = 0;
	}
    }
  else
    {
      status = grib_get_long(gh, "typeOfFirstFixedSurface", &lpar);
      if ( status == 0 )
	{
	  leveltype = (int) lpar;
775
	  grib2GetLevel(gh, leveltype, lbounds, &level1, &level2);
776
777
	  if ( leveltype == 99 ) leveltype++;
	}
778
      else
779
780
781
	{
	  leveltype = 0;
	}
Uwe Schulzweida's avatar
Uwe Schulzweida committed
782
    }
783

784
785
786
787
788
789
790
791
792
793
  // fprintf(stderr, "param %d %d %d %d\n", param, level1, level2, leveltype);

  (*record).size     = recsize;
  (*record).position = position;
  (*record).param    = param;
  (*record).ilevel   = level1;
  (*record).ilevel2  = level2;
  (*record).ltype    = leveltype;

  gribapiGetGrid(gh, &grid);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
794
795
796

  gridID = varDefGrid(vlistID, grid, 0);

Uwe Schulzweida's avatar
Uwe Schulzweida committed
797
  zaxistype = gribapiGetZaxisType(editionNumber, leveltype);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
798

799
  switch (zaxistype)
Uwe Schulzweida's avatar
Uwe Schulzweida committed
800
    {
801
802
803
804
805
806
    case ZAXIS_HYBRID:
    case ZAXIS_HYBRID_HALF:
      {
        int vctsize;
        size_t dummy;
        double *vctptr;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
807

808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
        GRIB_CHECK(grib_get_long(gh, "NV", &lpar), 0);
        vctsize = lpar;
        if ( vctsize > 0 )
          {
            vctptr = (double *) malloc(vctsize*sizeof(double));
            dummy = vctsize;
            GRIB_CHECK(grib_get_double_array(gh, "pv", vctptr, &dummy), 0);
            varDefVCT(vctsize, vctptr);
            free(vctptr);
          }
        break;
      }
    case ZAXIS_REFERENCE:
      {
        size_t len;
        char uuid[17];
        long nlev, nvgrid;

        GRIB_CHECK(grib_get_long(gh, "NV", &lpar), 0);
827
        if ( lpar != 3 )
828
829
830
831
832
833
834
835
836
837
838
839
          {
            fprintf(stderr, "Warning ...\n");
          }
        GRIB_CHECK(grib_get_long(gh, "nlev", &nlev), 0);
        GRIB_CHECK(grib_get_long(gh, "numberOfVGridUsed", &nvgrid), 0);

        len = (size_t) 16;
        uuid[16] = 0;
        GRIB_CHECK(grib_get_string(gh, "uuidOfVGrid", uuid, &len), 0);
        varDefZAxisReference((int) nlev, (int) nvgrid, uuid);
        break;
      }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
840
    }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
841

842
  // if ( datatype > 32 ) datatype = DATATYPE_PACK32;
843
  if ( datatype <  0 ) datatype = DATATYPE_PACK;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
844

845
846
847
848
  name[0] = 0;
  longname[0] = 0;
  units[0] = 0;

Uwe Schulzweida's avatar
Uwe Schulzweida committed
849
850
  vlen = 256;
  GRIB_CHECK(grib_get_string(gh, "shortName", name, &vlen), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
851
852
853
854
855
856
857
858
859
860
861
862
  if      ( vlen == 8 && memcmp(name, "unknown", vlen) == 0 ) name[0] = 0;
  else if ( vlen == 2 && memcmp(name, "~", vlen)       == 0 ) name[0] = 0;

  if ( name[0] != 0 )
    {
      vlen = 256;
      GRIB_CHECK(grib_get_string(gh, "name", longname, &vlen), 0);
      if ( vlen == 8 && memcmp(longname, "unknown", vlen) == 0 ) longname[0] = 0;
      vlen = 256;
      GRIB_CHECK(grib_get_string(gh, "units", units, &vlen), 0);
      if ( vlen == 8 && memcmp(units, "unknown", vlen) == 0 ) units[0] = 0;
    }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
863
  // fprintf(stderr, "param %d name %s %s %s\n", param, name, longname, units); 
Uwe Schulzweida's avatar
Uwe Schulzweida committed
864

Uwe Schulzweida's avatar
Uwe Schulzweida committed
865
  varAddRecord(recID, param, gridID, zaxistype, lbounds, level1, level2,
866
	       datatype, &varID, &levelID, 0, numavg, leveltype,
Uwe Schulzweida's avatar
Uwe Schulzweida committed
867
	       name, longname, units);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
868
869
870
871

  (*record).varID   = varID;
  (*record).levelID = levelID;

872
  varDefCompType(varID, comptype);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
873
874
875

  if ( varInqInst(varID) == CDI_UNDEFID )
    {
Uwe Schulzweida's avatar
Uwe Schulzweida committed
876
877
878
879
880
      long center, subcenter;
      int instID;
      GRIB_CHECK(grib_get_long(gh, "centre", &center), 0);
      GRIB_CHECK(grib_get_long(gh, "subCentre", &subcenter), 0);
      instID    = institutInq((int)center, (int)subcenter, NULL, NULL);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
881
      if ( instID == CDI_UNDEFID )
Uwe Schulzweida's avatar
Uwe Schulzweida committed
882
	instID = institutDef((int)center, (int)subcenter, NULL, NULL);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
883
884
885
886
887
888
      varDefInst(varID, instID);
    }

  if ( varInqModel(varID) == CDI_UNDEFID )
    {
      int modelID;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
889
      long processID;
890
891
892
893
894
895
896
897
      status = grib_get_long(gh, "generatingProcessIdentifier", &processID);
      if ( status == 0 )
	{
	  modelID = modelInq(varInqInst(varID), processID, NULL);
	  if ( modelID == CDI_UNDEFID )
	    modelID = modelDef(varInqInst(varID), processID, NULL);
	  varDefModel(varID, modelID);
	}
Uwe Schulzweida's avatar
Uwe Schulzweida committed
898
899
900
901
    }

  if ( varInqTable(varID) == CDI_UNDEFID )
    {
Uwe Schulzweida's avatar
Uwe Schulzweida committed
902
      int pdis, pcat, pnum;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
903

Uwe Schulzweida's avatar
Uwe Schulzweida committed
904
      cdiDecodeParam(param, &pnum, &pcat, &pdis);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
905

Uwe Schulzweida's avatar
Uwe Schulzweida committed
906
907
908
909
      if ( pdis == 255 )
	{
	  int tableID;
	  int tabnum = pcat;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
910

Uwe Schulzweida's avatar
Uwe Schulzweida committed
911
	  tableID = tableInq(varInqModel(varID), tabnum, NULL);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
912

Uwe Schulzweida's avatar
Uwe Schulzweida committed
913
914
915
	  if ( tableID == CDI_UNDEFID )
	    tableID = tableDef(varInqModel(varID), tabnum, NULL);
	  varDefTable(varID, tableID);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
916
917
918
919
920
921
922
	}
    }

  streamptr->tsteps[tsID].nallrecs++;
  streamptr->nrecs++;

  if ( CDI_Debug )
923
    Message("varID = %d  param = %d  zaxistype = %d  gridID = %d  levelID = %d",
Uwe Schulzweida's avatar
Uwe Schulzweida committed
924
	    varID, param, zaxistype, gridID, levelID);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
925
926
927
}
#endif

928
int gribapiScanTimestep1(int streamID)
Uwe Schulzweida's avatar
Uwe Schulzweida committed
929
930
931
932
933
934
{
#if  defined  (HAVE_LIBGRIB_API)
  off_t recpos = 0;
  unsigned char *gribbuffer = NULL;
  long buffersize = 0;
  int iret = 0, ipunp = 0, iword = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
935
  int rstatus;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
936
937
  int status;
  int fileID;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
938
  int rtabnum = 0;
939
940
  int rcode = 0, level1 = 0, level2 = 0;
  int vdate = 0, vtime = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
941
  int param = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
942
943
944
945
946
  DateTime datetime, datetime0;
  int tsID;
  int varID;
  size_t readsize;
  int nrecords, nrecs, recID;
947
  int datatype;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
948
949
950
951
952
  long recsize = 0;
  int warn_time = TRUE;
  int warn_numavg = TRUE;
  int taxisID = -1;
  int rdate = 0, rtime = 0, tunit = 0, fcast = 0;
953
  taxis_t *taxis;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
954
  int vlistID;
955
  int comptype;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
956
957
958
959
960
  long unzipsize;
  compvar2_t compVar, compVar0;
  stream_t *streamptr;
  grib_handle *gh = NULL;
  int leveltype;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
961
  int pdis = 0, pcat = 0, pnum = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
962
963
  long editionNumber;
  long lpar;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
964
  int bitsPerValue;
965
  int lieee = FALSE;
966
  int lbounds;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
967
968
969

  streamptr = stream_to_pointer(streamID);

970
  stream_check_ptr(__func__, streamptr);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
971
972
973
974
975
976
977

  streamptr->curTsID = 0;

  tsID  = tstepsNewEntry(streamID);
  taxis = &streamptr->tsteps[tsID].taxis;

  if ( tsID != 0 )
978
    Error("Internal problem! tstepsNewEntry returns %d", tsID);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
979
980
981
982
983
984

  fileID = streamInqFileID(streamID);

  nrecs = 0;
  while ( TRUE )
    {
985
986
      level1 = 0;
      level2 = 0;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
1001
      recsize = gribGetSize(fileID);
      recpos  = fileGetPos(fileID);

      if ( recsize == 0 )
	{
	  streamptr->ntsteps = 1;
	  break;
	}
      if ( recsize > buffersize )
	{
	  buffersize = recsize;
	  gribbuffer = (unsigned char *) realloc(gribbuffer, buffersize);
	}

      readsize = recsize;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1002
1003
      rstatus = gribRead(fileID, gribbuffer, &readsize);
      if ( rstatus ) break;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1004

1005
1006
      lieee = FALSE;

1007
      comptype = COMPRESS_NONE;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1008
1009
      if ( gribGetZip(recsize, gribbuffer, &unzipsize) > 0 )
	{
1010
	  comptype = COMPRESS_SZIP;
1011
	  unzipsize += 100;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1012
1013
1014
1015
1016
1017
1018
1019
	  if ( (long) buffersize < unzipsize )
	    {
	      buffersize = unzipsize;
	      gribbuffer = (unsigned char *) realloc(gribbuffer, buffersize);
	    }
	}

      gh = grib_handle_new_from_message(NULL, (void *) gribbuffer, recsize);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1020
      GRIB_CHECK(grib_set_double(gh, "missingValue", GRIBAPI_MISSVAL), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1021
1022
1023

      GRIB_CHECK(grib_get_long(gh, "editionNumber", &editionNumber), 0);

1024
1025
      lbounds = gribapiGetZaxisHasBounds(gh);

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1026
1027
      if ( editionNumber <= 1 )
	{
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1028
1029
	  GRIB_CHECK(grib_get_long(gh, "table2Version", &lpar), 0);
	  rtabnum = (int) lpar;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1030
1031
	  GRIB_CHECK(grib_get_long(gh, "indicatorOfParameter", &lpar), 0);
	  rcode = (int) lpar;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1032

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1033
	  param = cdiEncodeParam(rcode, rtabnum, 255);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1034

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1035
1036
1037
1038
	  status = grib_get_long(gh, "indicatorOfTypeOfLevel", &lpar);
	  if ( status == 0 )
	    {
	      leveltype = (int) lpar;
1039
1040
	      grib1GetLevel(gh, leveltype, lbounds, &level1, &level2);
	      if ( leveltype == 99 ) leveltype = 100;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1041
	    }
1042
	  else
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1043
1044
1045
1046
1047
1048
	    {
	      leveltype = 0;
	    }
	}
      else
	{
1049
1050
	  size_t len = 256;
	  char typeOfPacking[256];
1051

1052
	  status = grib_get_string(gh, "packingType", typeOfPacking, &len);
1053
1054
	  if ( status == 0 )
	    {
1055
	      // fprintf(stderr, "packingType %d %s\n", len, typeOfPacking);
1056
1057
	      if      ( strncmp(typeOfPacking, "grid_jpeg", len) == 0 ) comptype = COMPRESS_JPEG;
	      else if ( strncmp(typeOfPacking, "grid_ieee", len) == 0 ) lieee = TRUE;
1058
	    }
1059

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1060
	  GRIB_CHECK(grib_get_long(gh, "discipline", &lpar), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1061
1062
	  pdis = (int) lpar;

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1063
	  GRIB_CHECK(grib_get_long(gh, "parameterCategory", &lpar), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1064
	  pcat = (int) lpar;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1065

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1066
	  GRIB_CHECK(grib_get_long(gh, "parameterNumber", &lpar), 0);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1067
1068
	  pnum = (int) lpar;

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1069
	  param = cdiEncodeParam(pnum, pcat, pdis);
1070

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1071
1072
1073
1074
	  status = grib_get_long(gh, "typeOfFirstFixedSurface", &lpar);
	  if ( status == 0 )
	    {
	      leveltype = (int) lpar;
1075
1076
	      grib2GetLevel(gh, leveltype, lbounds, &level1, &level2);
	      if ( leveltype == 99 ) leveltype = 100;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1077
	    }
1078
	  else
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1079
1080
1081
1082
1083
	    {
	      leveltype = 0;
	    }
	}

1084
1085
      gribapiGetValidityDateTime(gh, &vdate, &vtime);
      /*
1086
      printf("%d %d %d.%d.%d %d\n", vdate, vtime, pnum, pcat, pdis, leveltype);
1087
      */
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
      if ( lieee )
        {
          datatype = DATATYPE_FLT64;
          status = grib_get_long(gh, "precision", &lpar);
          if ( status == 0 && lpar == 1 ) datatype = DATATYPE_FLT32;
        }
      else
        {
          datatype = DATATYPE_PACK;
          status = grib_get_long(gh, "bitsPerValue", &lpar);
          if ( status == 0 )
            {
              bitsPerValue = (int) lpar;
              if ( bitsPerValue > 0 && bitsPerValue <= 32 )
                datatype = bitsPerValue;
            }
        }
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1105

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1106
1107
1108
1109
      if ( nrecs == 0 )
	{
	  datetime0.date = vdate;
	  datetime0.time = vtime;
1110
1111
1112
1113
1114
1115
	  GRIB_CHECK(grib_get_long(gh, "dataDate", &lpar), 0);
	  rdate = (int) lpar;
	  GRIB_CHECK(grib_get_long(gh, "dataTime", &lpar), 0);
	  rtime = (int) lpar*100;
	  fcast = gribapiTimeIsFC(gh);
	  if ( fcast ) tunit = gribapiGetTimeUnits(gh);
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1116
1117
1118
1119
1120
	}
      else
	{
	  datetime.date  = vdate;
	  datetime.time  = vtime;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1121
1122

	  compVar.param  = param;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1123
1124
          compVar.level1 = level1;
          compVar.level2 = level2;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1125
1126
	  compVar.ltype  = leveltype;

Uwe Schulzweida's avatar
Uwe Schulzweida committed
1127
1128
	  for ( recID = 0; recID < nrecs; recID++ )
	    {
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1129
	      compVar0.param  = streamptr->tsteps[0].records[recID].param;
Uwe Schulzweida's avatar
Uwe Schulzweida committed
1130
1131
1132
	      compVar0.level1 = streamptr->tsteps[0].records[recID].ilevel;
	      compVar0.level2 = streamptr->tsteps[0].records[recID].ilevel2;
	      compVar0.ltype  = streamptr->tsteps[0].records[recID].ltype;