mirror of https://github.com/ecmwf/eccodes.git
307 lines
9.4 KiB
Bash
Executable File
307 lines
9.4 KiB
Bash
Executable File
#!/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_ls_test"
|
|
tempLog=temp.$label.log
|
|
temp1=temp.$label.1.txt
|
|
temp2=temp.$label.2.txt
|
|
tempText=temp.$label.txt
|
|
rm -f $temp1 $temp2 $tempText $tempLog
|
|
|
|
sample_g1="$ECCODES_SAMPLES_PATH/GRIB1.tmpl"
|
|
sample_g2="$ECCODES_SAMPLES_PATH/GRIB2.tmpl"
|
|
|
|
cd ${data_dir}
|
|
infile=regular_gaussian_model_level.grib1
|
|
|
|
# Basic test of man page
|
|
set +e
|
|
${tools_dir}/grib_ls > $tempLog
|
|
status=$?
|
|
set -e
|
|
[ $status -eq 1 ]
|
|
grep -q "Full documentation and examples at" $tempLog
|
|
grep -q "https://confluence.ecmwf.int/display/ECC/grib_ls" $tempLog
|
|
|
|
set +e
|
|
${tools_dir}/grib_ls -? > $tempLog
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
|
|
set +e
|
|
${tools_dir}/grib_ls -h > $tempLog
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
|
|
set +e
|
|
DOXYGEN_USAGE=1 ${tools_dir}/grib_ls > $tempLog
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
|
|
|
|
${tools_dir}/grib_ls -P count $infile > $tempLog
|
|
${tools_dir}/grib_ls -p count,step $infile >> $tempLog
|
|
${tools_dir}/grib_ls $infile >> $tempLog
|
|
if [ $HAVE_GEOGRAPHY -eq 1 ]; then
|
|
${tools_dir}/grib_ls -l 0,0,1 $infile >> $tempLog
|
|
${tools_dir}/grib_get -l 0,0,1 $infile >> $tempLog
|
|
fi
|
|
${tools_dir}/grib_get -p count,step $infile >> $tempLog
|
|
${tools_dir}/grib_get -P count $infile >> $tempLog
|
|
|
|
# ECC-786 and ECC-791
|
|
result=$( ${tools_dir}/grib_get -p shortName -i 0 $infile )
|
|
[ "$result" = "t 199.078 " ]
|
|
result=$( ${tools_dir}/grib_get -i 8191 $infile )
|
|
[ "$result" = "160.852 " ]
|
|
|
|
|
|
files=" reduced_gaussian_lsm.grib1
|
|
reduced_gaussian_model_level.grib1
|
|
reduced_gaussian_model_level.grib2
|
|
reduced_gaussian_pressure_level.grib1
|
|
reduced_gaussian_pressure_level.grib2
|
|
reduced_gaussian_pressure_level_constant.grib1
|
|
reduced_gaussian_pressure_level_constant.grib2
|
|
reduced_gaussian_sub_area.grib1
|
|
reduced_gaussian_sub_area.grib2
|
|
reduced_gaussian_surface.grib1
|
|
reduced_gaussian_surface.grib2
|
|
reduced_latlon_surface.grib1
|
|
reduced_latlon_surface.grib2
|
|
regular_gaussian_model_level.grib1
|
|
regular_gaussian_model_level.grib2
|
|
regular_gaussian_pressure_level.grib1
|
|
regular_gaussian_pressure_level.grib2
|
|
regular_gaussian_pressure_level_constant.grib1
|
|
regular_gaussian_pressure_level_constant.grib2
|
|
regular_gaussian_surface.grib1
|
|
regular_gaussian_surface.grib2
|
|
regular_latlon_surface.grib1
|
|
regular_latlon_surface.grib2
|
|
"
|
|
if [ $HAVE_GEOGRAPHY -eq 1 ]; then
|
|
for file in $files; do
|
|
[ -f "$file" ]
|
|
echo $file >> $tempLog
|
|
${tools_dir}/grib_ls -l 40,28 $file | grep index | awk '{print $4;}' >> $tempLog
|
|
done
|
|
|
|
diff $tempLog ls.log
|
|
rm -f $tempLog
|
|
fi
|
|
|
|
echo "Test for bug GRIB-56..."
|
|
# ------------------------------
|
|
${tools_dir}/grib_set -s typeOfLevel=depthBelowLandLayer,topLevel=missing regular_latlon_surface.grib1 tmp_rlls.grib1
|
|
${tools_dir}/grib_ls -plevels tmp_rlls.grib1 | grep MISSING>/dev/null
|
|
|
|
${tools_dir}/grib_set -s typeOfLevel=depthBelowLandLayer,bottomLevel=missing regular_latlon_surface.grib1 tmp_rlls.grib1
|
|
${tools_dir}/grib_ls -plevels tmp_rlls.grib1 | grep MISSING>/dev/null
|
|
rm -f tmp_rlls.grib1
|
|
|
|
echo "GRIB-305. GRIB edition 1 file with one large message..."
|
|
# ----------------------------------------------------------
|
|
if [ -f "sst_globus0083.grib" ]; then
|
|
${tools_dir}/grib_ls sst_globus0083.grib > /dev/null
|
|
fi
|
|
|
|
echo "GRIB-387 printing key of type byte..."
|
|
# ----------------------------------------------------------
|
|
${tools_dir}/grib_ls -p uuidOfVGrid test_uuid.grib2 > /dev/null
|
|
|
|
type=`${tools_dir}/grib_get -wcount=1 -p typeOfLevel test_uuid.grib2`
|
|
[ "$type" = "generalVertical" ]
|
|
|
|
|
|
if [ $HAVE_GEOGRAPHY -eq 1 ]; then
|
|
echo "GRIB-213 nearest with land-sea mask..."
|
|
# ----------------------------------------------------------
|
|
${tools_dir}/grib_ls -l 85,13,1,reduced_gaussian_lsm.grib1 reduced_gaussian_surface.grib1 >$tempText
|
|
grep -q 'Point chosen #3 index=21 .* distance=11\.' $tempText
|
|
|
|
${tools_dir}/grib_ls -l 53,2,1,reduced_gaussian_lsm.grib1 reduced_gaussian_surface.grib1 >$tempText
|
|
grep -q 'Point chosen #2 index=749 .* distance=204\.' $tempText
|
|
|
|
${tools_dir}/grib_get -F%.2f -l 85,13,1,reduced_gaussian_lsm.grib1 reduced_gaussian_surface.grib1 >$tempText
|
|
grep -q '252.88' $tempText
|
|
fi
|
|
|
|
echo "ECC-278: grib_ls -n namespace..."
|
|
# ----------------------------------------------------------
|
|
${tools_dir}/grib_ls -n geography $ECCODES_SAMPLES_PATH/reduced_ll_sfc_grib2.tmpl
|
|
${tools_dir}/grib_ls -n data $sample_g1
|
|
|
|
set +e
|
|
${tools_dir}/grib_ls -n nosuchnamespace $sample_g1 > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "does not contain any key" $tempText
|
|
|
|
|
|
# Angle subdivisions
|
|
grib_check_key_equals $sample_g1 angleSubdivisions 1000
|
|
grib_check_key_equals $sample_g2 angleSubdivisions 1000000
|
|
|
|
# Print 'offset' key as string and integer
|
|
${tools_dir}/grib_ls -p offset:s tigge_cf_ecmwf.grib2 > $temp1
|
|
${tools_dir}/grib_ls -p offset:i tigge_cf_ecmwf.grib2 > $temp2
|
|
diff $temp1 $temp2
|
|
|
|
# Section pointers
|
|
grib_check_key_equals $sample_g2 'section0Pointer,section1Pointer,section3Pointer,section4Pointer' '0_16 16_21 37_72 109_34'
|
|
|
|
|
|
echo "Test constraints: -w option..."
|
|
# ----------------------------------------------------------
|
|
file=tigge_pf_ecmwf.grib2 # Has 38 messages
|
|
${tools_dir}/grib_ls -w count!=1 $file > $tempText
|
|
grep -q "37 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w count=1/2/38 $file > $tempText
|
|
grep -q "3 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w shortName=u $file > $tempText
|
|
grep -q "3 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w shortName=u/v $file > $tempText
|
|
grep -q "5 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w shortName!=u/v $file > $tempText
|
|
grep -q "33 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w shortName=u/v,level=925 $file > $tempText
|
|
grep -q "2 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w shortName=u/v/10u/10v,level=925/10 $file > $tempText
|
|
grep -q "4 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w packingType!=grid_simple $file > $tempText
|
|
grep -q "0 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w units!=K $file > $tempText
|
|
grep -q "30 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w scaleFactorOfSecondFixedSurface=missing $file > $tempText
|
|
grep -q "36 of 38 messages" $tempText
|
|
${tools_dir}/grib_ls -w scaleFactorOfSecondFixedSurface!=missing $file > $tempText
|
|
grep -q "2 of 38 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w referenceValue:d=0 $file > $tempText
|
|
grep -q "5 of 38 messages" $tempText
|
|
|
|
file=mixed.grib # Has 14 messages
|
|
${tools_dir}/grib_ls -w packingType=grid_simple,gridType=regular_ll/regular_gg $file > $tempText
|
|
grep -q "12 of 14 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w packingType=grid_simple/grid_simple_matrix,gridType=regular_ll/regular_gg $file > $tempText
|
|
grep -q "12 of 14 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w typeOfLevel=surface,centre=7 $file > $tempText
|
|
grep -q "3 of 14 messages" $tempText
|
|
|
|
${tools_dir}/grib_ls -w shortName=t/10u,gridType=regular_gg $file > $tempText
|
|
grep -q "5 of 14 messages" $tempText
|
|
|
|
|
|
# ECC-1528: Try some invalid values for the count constraint
|
|
for cval in -1 0 xx; do
|
|
set +e
|
|
rm -f $tempText
|
|
${tools_dir}/grib_ls -w count=$cval $file > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "Invalid value for key 'count'" $tempText
|
|
done
|
|
|
|
|
|
# ECC-1562: Segmentation fault: Invalid orderby directive
|
|
set +e
|
|
${tools_dir}/grib_ls -B'shortName: asc' tigge_af_ecmwf.grib2 > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "Invalid type for key=shortName" $tempText
|
|
|
|
# Do list after an offset
|
|
file=tigge_pf_ecmwf.grib2
|
|
${tools_dir}/grib_ls -X 62414 $file
|
|
set +e
|
|
${tools_dir}/grib_ls -X -1 $file > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "Invalid file offset" $tempText
|
|
|
|
# Decode an ascii key as double
|
|
file=$ECCODES_SAMPLES_PATH/reduced_gg_pl_32_grib2.tmpl
|
|
grib_check_key_equals $file 'expver:d' 1
|
|
grib_check_key_equals $file 'expver:s' '0001'
|
|
|
|
if [ $HAVE_GEOGRAPHY -eq 1 ]; then
|
|
# JSON and lat/lon
|
|
${tools_dir}/grib_ls -j -l0,0 -p referenceValue:d $data_dir/sample.grib2
|
|
${tools_dir}/grib_ls -j -l0,0 -p referenceValue:i $data_dir/sample.grib2
|
|
${tools_dir}/grib_ls -j -l0,0 -p bitmap $data_dir/simple_bitmap.grib > $tempText 2>&1
|
|
grep -q "invalid_type" $tempText
|
|
${tools_dir}/grib_ls -j -l0,0 -p nosuchkey $data_dir/sample.grib2 > $tempText 2>&1
|
|
grep -q "nosuchkey.* null" $tempText
|
|
fi
|
|
|
|
# -M and -g options
|
|
${tools_dir}/grib_ls -M -g $data_dir/gts.grib
|
|
|
|
|
|
if [ $HAVE_GEOGRAPHY -eq 1 ]; then
|
|
${tools_dir}/grib_get -l0,0,4 $data_dir/sample.grib2
|
|
|
|
set +e
|
|
${tools_dir}/grib_ls -l0,0,666 $data_dir/sample.grib2 > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "Wrong mode given" $tempText
|
|
|
|
set +e
|
|
${tools_dir}/grib_ls -l poo $data_dir/sample.grib2 > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "Wrong latitude value" $tempText
|
|
|
|
set +e
|
|
${tools_dir}/grib_ls -l0,0,1,nonexistingmask $data_dir/sample.grib2 > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "unable to open mask file" $tempText
|
|
fi
|
|
|
|
# ----------------------
|
|
# Printing array keys
|
|
# ----------------------
|
|
set +e
|
|
${tools_dir}/grib_ls -p bitmap $data_dir/reduced_latlon_surface.grib2 > $tempText 2>&1
|
|
status=$?
|
|
set -e
|
|
[ $status -ne 0 ]
|
|
grep -q "Hint: Tool grib_ls cannot print keys of array type" $tempText
|
|
|
|
|
|
# Clean up
|
|
rm -f $temp1 $temp2 $tempText $tempLog
|