Commit 48501709 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

stream_ieg: multiply pressure levels with 100

parent 565671ec
...@@ -2,6 +2,7 @@ ...@@ -2,6 +2,7 @@
* using GRIB library version 1.0.0 * using GRIB library version 1.0.0
* rename *New functions to *Create * rename *New functions to *Create
* stream_ieg: multiply pressure levels with 100 (report: Holger Goettel)
* add dummy functions gridDestroy, zaxisDestroy, taxisDestroy * add dummy functions gridDestroy, zaxisDestroy, taxisDestroy
* Version 1.0.0 released * Version 1.0.0 released
......
...@@ -605,6 +605,7 @@ void iegAddRecord(int streamID, int code, int *pdb, int *gdb, double *vct, ...@@ -605,6 +605,7 @@ void iegAddRecord(int streamID, int code, int *pdb, int *gdb, double *vct,
{ {
level1 = IEG_P_Level2(pdb); level1 = IEG_P_Level2(pdb);
level2 = 0; level2 = 0;
if ( IEG_P_LevelType(pdb) == 100 ) level1 *= 100;
} }
/* /*
fprintf(stderr, "%d %d %d %d\n", code, IEG_P_LevelType(pdb), IEG_P_Level1(pdb), IEG_P_Level2(pdb)); fprintf(stderr, "%d %d %d %d\n", code, IEG_P_LevelType(pdb), IEG_P_Level1(pdb), IEG_P_Level2(pdb));
...@@ -780,6 +781,8 @@ void iegScanTimestep1(int streamID) ...@@ -780,6 +781,8 @@ void iegScanTimestep1(int streamID)
else else
rlevel = IEG_P_Level2(iegp->ipdb); rlevel = IEG_P_Level2(iegp->ipdb);
if ( IEG_P_LevelType(iegp->ipdb) == 100 ) rlevel *= 100;
iegDateTime(iegp->ipdb, &vdate, &vtime); iegDateTime(iegp->ipdb, &vdate, &vtime);
if ( nrecs == 0 ) if ( nrecs == 0 )
...@@ -929,6 +932,8 @@ int iegScanTimestep2(int streamID) ...@@ -929,6 +932,8 @@ int iegScanTimestep2(int streamID)
else else
rlevel = IEG_P_Level2(iegp->ipdb); rlevel = IEG_P_Level2(iegp->ipdb);
if ( IEG_P_LevelType(iegp->ipdb) == 100 ) rlevel *= 100;
iegDateTime(iegp->ipdb, &vdate, &vtime); iegDateTime(iegp->ipdb, &vdate, &vtime);
if ( rindex == 0 ) if ( rindex == 0 )
...@@ -1103,6 +1108,8 @@ int iegScanTimestep(int streamID) ...@@ -1103,6 +1108,8 @@ int iegScanTimestep(int streamID)
else else
rlevel = IEG_P_Level2(iegp->ipdb); rlevel = IEG_P_Level2(iegp->ipdb);
if ( IEG_P_LevelType(iegp->ipdb) == 100 ) rlevel *= 100;
iegDateTime(iegp->ipdb, &vdate, &vtime); iegDateTime(iegp->ipdb, &vdate, &vtime);
if ( rindex == nrecs ) break; if ( rindex == nrecs ) break;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment