diff --git a/definitions/grib2/section.4.def b/definitions/grib2/section.4.def index fd65d076b..0c80854e4 100644 --- a/definitions/grib2/section.4.def +++ b/definitions/grib2/section.4.def @@ -49,16 +49,16 @@ transient genVertHeightCoords = 0; if (defined(typeOfFirstFixedSurface)) { if (typeOfFirstFixedSurface == 150) { transient genVertHeightCoords = 1; + transient PVPresent = 0; + ieeefloat nlev : dump ; + ieeefloat numberOfVGridUsed : dump; + byte[16] uuidOfVGrid : dump; } } if (genVertHeightCoords) { # Generalized vertical height coordinate case - transient PVPresent = 0; - ieeefloat nlev : dump ; alias numberOfVerticalCoordinateValues = nlev ; - ieeefloat numberOfVGridUsed : dump; - byte[16] uuidOfVGrid : dump; } else { if (NV == 0){ diff --git a/tests/filter.sh b/tests/filter.sh index 44307fee2..cedee9174 100755 --- a/tests/filter.sh +++ b/tests/filter.sh @@ -58,5 +58,20 @@ EOF result=`${tools_dir}/grib_filter ${data_dir}/formatint.rules $GRIB_SAMPLES_PATH/GRIB1.tmpl` [ "$result" = "edition=001, level= 500" ] +# Convert from grib1 to grib2 "Generalized vertical height coordinates" +cat >temp.filt <> local.log ${tools_dir}grib_get -p localDefinitionNumber loc1.grib2 >> local.log - + for localEps in 1 15 26 30 do ${tools_dir}grib_set -s localDefinitionNumber=$localEps,eps=1,perturbationNumber=2,numberOfForecastsInEnsemble=50 loc1.grib2 eps.grib2