stream_cgribex: check range of pfactor

9b159681 · Uwe Schulzweida · ba57b8d9 · 9b159681 · 9b159681
Commit 9b159681 authored Feb 15, 2010 by Uwe Schulzweida
Hide whitespace changes
Inline Side-by-side

Showing with 9 additions and 1 deletion

ChangeLog ChangeLog +4 -0

src/stream_cgribex.c src/stream_cgribex.c +5 -1

No files found.
--- a/ChangeLog
+++ b/ChangeLog
+2010-02-15  Uwe Schulzweida  <Uwe.Schulzweida@zmaw.de>
+
+	* cgribexDefTime: bug fix for timerange = -1
+
 2010-02-10  Uwe Schulzweida  <Uwe.Schulzweida@zmaw.de>

 	* added support for complex packing of spectral data in GRIB

--- a/src/stream_cgribex.c
+++ b/src/stream_cgribex.c
@@ -2151,7 +2151,11 @@ size_t cgribexEncode(int varID, int levelID, int vlistID, int gridID, int zaxisI
    }

  if ( isec4[2] == 128 && isec4[3] == 64 )
-    isec4[16] = (int) (1000*calculate_pfactor(data, ISEC2_PentaJ, isec4[17]));
+    {
+      isec4[16] = (int) (1000*calculate_pfactor(data, ISEC2_PentaJ, isec4[17]));
+      if ( isec4[16] < -10000 ) isec4[16] = -10000;
+      if ( isec4[16] >  10000 ) isec4[16] =  10000;
+    }
  //printf("isec4[16] %d\n", isec4[16]);

  gribExDP(isec0, isec1, isec2, fsec2, isec3, fsec3, isec4, (double*) data,