2023-06-16 10:24:21 +00:00
|
|
|
|
#!/bin/sh
|
|
|
|
|
|
# (C) Copyright 2005- ECMWF.
|
|
|
|
|
|
#
|
|
|
|
|
|
# This software is licensed under the terms of the Apache Licence Version 2.0
|
|
|
|
|
|
# which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
|
|
|
|
|
|
#
|
|
|
|
|
|
# In applying this licence, ECMWF does not waive the privileges and immunities granted to it by
|
|
|
|
|
|
# virtue of its status as an intergovernmental organisation nor does it submit to any jurisdiction.
|
|
|
|
|
|
#
|
|
|
|
|
|
|
|
|
|
|
|
. ./include.ctest.sh
|
|
|
|
|
|
|
|
|
|
|
|
|
2023-10-09 10:12:44 +00:00
|
|
|
|
label="grib_ecPoint_test"
|
2023-06-16 10:24:21 +00:00
|
|
|
|
tempGrib=temp.$label.grib
|
|
|
|
|
|
|
|
|
|
|
|
sample_grib2=$ECCODES_SAMPLES_PATH/GRIB2.tmpl
|
|
|
|
|
|
|
2023-10-09 10:33:48 +00:00
|
|
|
|
latest=`${tools_dir}/grib_get -p tablesVersionLatest $sample_grib2`
|
|
|
|
|
|
|
2023-06-16 10:24:21 +00:00
|
|
|
|
# 1. bias-corrected
|
|
|
|
|
|
## we check if we can set PDTN 73 and typeOfPostProcessing to ecPoint, ecPoint versioning and MARS type
|
2023-10-09 10:33:48 +00:00
|
|
|
|
${tools_dir}/grib_set -s tablesVersion=$latest,productDefinitionTemplateNumber=73,backgroundProcess=1,\
|
2023-06-16 10:24:21 +00:00
|
|
|
|
setLocalDefinition=1,grib2LocalSectionNumber=1,type=gbf,class=od,stream=enfo,\
|
|
|
|
|
|
generatingProcessIdentifier=2,typeOfPostProcessing=206 $sample_grib2 $tempGrib
|
|
|
|
|
|
grib_check_key_equals $tempGrib backgroundProcess,generatingProcessIdentifier,typeOfPostProcessing,type "1 2 206 gbf"
|
|
|
|
|
|
rm -f $tempGrib
|
|
|
|
|
|
|
|
|
|
|
|
# 2. Point values
|
2023-06-20 09:47:52 +00:00
|
|
|
|
## we check if we can set PDTN 90 and typeOfPostProcessing to ecPoint, ecPoint versioning, MARS type, quantile
|
2023-10-09 10:33:48 +00:00
|
|
|
|
${tools_dir}/grib_set -s tablesVersion=$latest,productDefinitionTemplateNumber=90,backgroundProcess=1,\
|
2023-06-16 10:24:21 +00:00
|
|
|
|
setLocalDefinition=1,grib2LocalSectionNumber=1,type=pfc,class=od,stream=enfo,\
|
|
|
|
|
|
totalNumberOfQuantiles=1000,quantileValue=995,\
|
|
|
|
|
|
generatingProcessIdentifier=2,typeOfPostProcessing=206 $sample_grib2 $tempGrib
|
|
|
|
|
|
grib_check_key_equals $tempGrib backgroundProcess,generatingProcessIdentifier,typeOfPostProcessing,type,quantile "1 2 206 pfc 995:1000"
|
|
|
|
|
|
rm -f $tempGrib
|
|
|
|
|
|
|
|
|
|
|
|
# 3. Point value metrics
|
2023-06-20 09:47:52 +00:00
|
|
|
|
## we check if we can set PDTN 90 and typeOfPostProcessing to ecPoint, ecPoint versioning, MARS type, quantile
|
2023-10-09 10:33:48 +00:00
|
|
|
|
${tools_dir}/grib_set -s tablesVersion=$latest,productDefinitionTemplateNumber=90,backgroundProcess=1,\
|
2023-06-16 10:24:21 +00:00
|
|
|
|
setLocalDefinition=1,grib2LocalSectionNumber=1,type=ppm,class=od,stream=enfo,\
|
|
|
|
|
|
totalNumberOfQuantiles=1000,quantileValue=995,\
|
|
|
|
|
|
generatingProcessIdentifier=2,typeOfPostProcessing=206 $sample_grib2 $tempGrib
|
|
|
|
|
|
grib_check_key_equals $tempGrib backgroundProcess,generatingProcessIdentifier,typeOfPostProcessing,type,quantile "1 2 206 ppm 995:1000"
|
|
|
|
|
|
rm -f $tempGrib
|
|
|
|
|
|
|
|
|
|
|
|
# 4. weather types
|
|
|
|
|
|
## we check if we can set PDTN 73 and typeOfPostProcessing to ecPoint, ecPoint versioning and MARS type
|
2023-10-09 10:33:48 +00:00
|
|
|
|
${tools_dir}/grib_set -s tablesVersion=$latest,productDefinitionTemplateNumber=73,backgroundProcess=1,\
|
2023-06-16 10:24:21 +00:00
|
|
|
|
setLocalDefinition=1,grib2LocalSectionNumber=1,type=gwt,class=od,stream=enfo,\
|
|
|
|
|
|
generatingProcessIdentifier=2,typeOfPostProcessing=206 $sample_grib2 $tempGrib
|
|
|
|
|
|
grib_check_key_equals $tempGrib backgroundProcess,generatingProcessIdentifier,typeOfPostProcessing,type "1 2 206 gwt"
|
2023-10-09 10:12:44 +00:00
|
|
|
|
|
|
|
|
|
|
# Clean up
|
2023-06-16 10:24:21 +00:00
|
|
|
|
rm -f $tempGrib
|
