eccodes/tests/grib_tigge_check.sh

#!/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

label="grib_tigge_check_test"
REDIRECT=/dev/null
TEMP=temp.$label.tigge

dir="${data_dir}/tigge/"

# Parameter 168: at 2m and 1.5m above ground
# -------------------------------------------
sample_g2=$ECCODES_SAMPLES_PATH/GRIB2.tmpl
${tools_dir}/grib_filter -o $TEMP - $sample_g2 << EOF
 set centre = "ammc";
 set productionStatusOfProcessedData=4; # TIGGE
 set discipline=0;
 set parameterCategory=0;
 set parameterNumber=6;
 set typeOfFirstFixedSurface=103;
 set scaledValueOfFirstFixedSurface=15;
 set scaleFactorOfFirstFixedSurface=1;
 write;
EOF
grib_check_key_equals $TEMP paramId,shortName '168 2d'

${tools_dir}/grib_filter -o $TEMP - $sample_g2 << EOF
 set centre = "ammc";
 set productionStatusOfProcessedData=4; # TIGGE
 set discipline=0;
 set parameterCategory=0;
 set parameterNumber=6;
 set typeOfFirstFixedSurface=103;
 set scaledValueOfFirstFixedSurface=2;
 set scaleFactorOfFirstFixedSurface=0;
 write;
EOF
grib_check_key_equals $TEMP paramId,shortName '168 2d'


# Check tigge global
# ------------------
for file in ${dir}tigge_*.grib; do
   ${tigge_dir}/tigge_check ${file} 2> $REDIRECT > $REDIRECT
done

# Check tigge-lam
# ------------------
for file in ${dir}tiggelam_*.grib; do
   ${tigge_dir}/tigge_check -l ${file} 2> $REDIRECT > $REDIRECT
done


# Test non-TIGGE files too. We now expect tigge_check to fail!
# All the GRIB files in the samples are non-TIGGE
# -------------------------------------------------------------
for file in ${ECCODES_SAMPLES_PATH}/regular_*.tmpl; do
   set +e
   ${tigge_dir}/tigge_check ${file} 2> $REDIRECT > $REDIRECT
   status=$?
   set -e
   if [ $status -eq 0 ]; then
      # should have failed and returned a non-zero exit code
      exit 1
   fi
done


# GRIB-531
# ---------
${tools_dir}/grib_get -nparameter ${data_dir}/tigge_pf_ecmwf.grib2 > $TEMP
diff ${data_dir}/tigge_pf_ecmwf.grib2.ref $TEMP

# GRIB-205. Changing productionStatusOfProcessedData
# should not change anything else
# ---------------------------------------------------
input=${dir}/tigge_ecmf_sfc_sd.grib
${tools_dir}/grib_set -s productionStatusOfProcessedData=5 $input $TEMP
${tools_dir}/grib_compare -bproductionStatusOfProcessedData $input $TEMP

rm -f $TEMP

# GRIB-757 validity date/time check fails for the following:
tigge_bad_validity="
  tigge_kwbc_sfc_sf.grib
  tigge_kwbc_sfc_slhf.grib
  tigge_kwbc_sfc_sshf.grib
  tigge_kwbc_sfc_ssr.grib
  tigge_kwbc_sfc_str.grib
  tigge_kwbc_sfc_ttr.grib
"
for file in $tigge_bad_validity; do
   set +e
   ${tigge_dir}/tigge_check -w ${dir}${file} > $TEMP
   status=$?
   set -e
   [ $status -eq 1 ]
   cat $TEMP
   grep -q "invalid validity Date/Time" $TEMP
done


# ECC-1428
# ----------
set +e
${tools_dir}/grib_set -s productionStatusOfProcessedData=5,marsType=xx $sample_g2 $TEMP
status=$?
set -e
[ $status -ne 0 ]

# ECC-1081
# ---------
${tools_dir}/grib_set -s level=2 ${data_dir}/tigge/tigge_ecmf_pv_pt.grib $TEMP
grib_check_key_equals $TEMP level 2
# ${tools_dir}/grib_dump -O -p section_4 $TEMP
${tools_dir}/grib_compare ${data_dir}/tigge/tigge_ecmf_pv_pt.grib $TEMP


rm -f $TEMP