Testing: Function to check grib_set and grib_filter

This commit is contained in:
shahramn 2024-08-19 11:25:37 +01:00
parent 1b532e9260
commit 22d4bff321
1 changed files with 33 additions and 0 deletions

View File

@ -49,3 +49,36 @@ grib_check_key_is_scalar()
result=$(echo "print '|[$a_key]|';" | ${tools_dir}/grib_filter - $a_file)
[ "$result" = "|$a_val|" ]
}
# This ensures setting the keys via grib_filter (similar to the API, one key at a time)
# works the same way as grib_set. The resultant GRIB outputs are compared
# E.g.,
# grib_encode_via_set_and_filter $label key1=val1,key2=val2,key3=val3 $ECCODES_SAMPLES_PATH/GRIB2.tmpl
#
grib_encode_via_set_and_filter()
{
_label=$1 # For generating temp files
_kvals=$2 # key1=val1,key2=val2
_input=$3 # Must be a GRIB file
temp_encode_filt=temp.encode.$_label.filt
temp_encode_gribA=temp.encode.$_label.A.grib
temp_encode_gribB=temp.encode.$_label.B.grib
for kv in $(echo $_kvals | tr ',' ' '); do
k=$(echo $kv | awk -F= '{print $1}')
v=$(echo $kv | awk -F= '{print $2}')
if [[ "$v" =~ ^[A-z] ]]; then
echo "set $k = \"$v\" ;" >> $temp_encode_filt
else
echo "set $k = $v ;" >> $temp_encode_filt
fi
done
echo "write;" >> $temp_encode_filt
# cat $temp_encode_filt
${tools_dir}/grib_filter -o $temp_encode_gribA $temp_encode_filt $_input
${tools_dir}/grib_set -s $_kvals $_input $temp_encode_gribB
${tools_dir}/grib_compare $temp_encode_gribA $temp_encode_gribB
rm -f $temp_encode_filt $temp_encode_gribB $temp_encode_gribA
}