1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
#!/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
set -u
label="grib_levtype_test"
if [ ! -d "$ECCODES_DEFINITION_PATH" ]; then
echo "Test $0 disabled. No definitions directory"
exit 0
fi
sample2=$ECCODES_SAMPLES_PATH/GRIB2.tmpl
tempGrib=temp.${label}.grib
tempText=temp.${label}.txt
latestAvailable=`${tools_dir}/grib_get -p tablesVersionLatest $sample2`
latestOfficial=`${tools_dir}/grib_get -p tablesVersionLatestOfficial $sample2`
# ECC-1333: levtype should be read-only in GRIB2
set +e
${tools_dir}/grib_set -s mars.levtype=o2d $sample2 $tempGrib 2>$tempText
status=$?
set -e
[ $status -ne 0 ]
grep -q "Value is read only" $tempText
# These level types are S2S ocean parameters and are dealt with differently (See products_s2s.def)
exclude="20 160 169"
paramId_file_wmo="$ECCODES_DEFINITION_PATH/grib2/paramId.def"
paramId_file_ecm="$ECCODES_DEFINITION_PATH/grib2/localConcepts/ecmf/paramId.def"
levtypes=`grep typeOfFirstFixedSurface $paramId_file_wmo $paramId_file_ecm |\
awk -F'=' '{print $2}'|tr -d ' '|tr -d ';'|sort -un`
for lt in $levtypes; do
process_type=1
for ex in $exclude; do
if [ "$lt" = "$ex" ]; then process_type=0; break; fi
done
if [ $process_type = 1 ]; then
${tools_dir}/grib_set -s tablesVersion=$latestAvailable,typeOfFirstFixedSurface=$lt $sample2 $tempGrib
result=`${tools_dir}/grib_get -p mars.levtype $tempGrib`
if [ "$result" = "$lt" ]; then
echo "ERROR: typeOfFirstFixedSurface of |$lt| not mapped to a string!"
exit 1
fi
${tools_dir}/grib_set -s tablesVersion=$latestOfficial,typeOfFirstFixedSurface=$lt $sample2 $tempGrib
result=`${tools_dir}/grib_get -p mars.levtype $tempGrib`
if [ "$result" = "$lt" ]; then
echo "ERROR: typeOfFirstFixedSurface of |$lt| not mapped to a string!"
exit 1
fi
fi
done
# ECC-1328
params='228007 228011'
for p in $params; do
${tools_dir}/grib_set -s paramId=$p $sample2 $tempGrib
grib_check_key_equals $tempGrib 'mars.levtype,typeOfLevel' 'sfc entireLake'
done
rm -f $tempGrib
rm -f $tempText
|
