gribapiDefDateTimeRel: forecastTime [Bug #5435]

ChangeLog

@@ -4,6 +4,7 @@
 
 2015-02-02  Uwe Schulzweida
 
+	* gribapiDefDateTimeRel: forecastTime [Bug #5435]
 	* cdfCopyRecord: incorrect checksum for freed object [Bug #5461]
 
 2015-01-16  Uwe Schulzweida

src/mo_cdi.f90

@@ -1540,9 +1540,10 @@ module mo_cdi
         end function gridInqMask
       end interface
       interface
-        subroutine gridPrint(gridID,opt) bind(c,name='gridPrint')
+        subroutine gridPrint(gridID,index,opt) bind(c,name='gridPrint')
           import :: c_int
           integer(kind=c_int), value :: gridID
+          integer(kind=c_int), value :: index
           integer(kind=c_int), value :: opt
         end subroutine gridPrint
       end interface
@@ -2199,9 +2200,10 @@ module mo_cdi
         end subroutine zaxisResize
       end interface
       interface
-        subroutine zaxisPrint(zaxisID) bind(c,name='zaxisPrint')
+        subroutine zaxisPrint(zaxisID,index) bind(c,name='zaxisPrint')
           import :: c_int
           integer(kind=c_int), value :: zaxisID
+          integer(kind=c_int), value :: index
         end subroutine zaxisPrint
       end interface
       interface

src/stream_gribapi.c

@@ -2307,7 +2307,11 @@ int gribapiDefDateTimeRel(int editionNumber, grib_handle *gh, int rdate, int rti
 
       // printf(">>>>> tsteptype %d  startStep %ld  endStep %ld\n", tsteptype, startStep, endStep);
 
-      if ( proDefTempNum == 0 ) startStep = endStep;
+      // Product Definition Template Number: defined in GRIB_API file 4.0.table
+      // point in time products:
+      if ( (proDefTempNum >= 0 && proDefTempNum <=  7) || 
+           proDefTempNum == 55 || proDefTempNum == 40055 ) // Tile
+        startStep = endStep;
 
       if ( editionNumber > 1 ) GRIB_CHECK(my_grib_set_long(gh, "forecastTime", startStep), 0);
       GRIB_CHECK(my_grib_set_long(gh, "endStep", endStep), 0);