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
|
#!/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_multi_test"
tmpdata=temp.$label.grib
rm -f $tmpdata
parameterNumber=`${tools_dir}/grib_get -p parameterNumber -w parameterCategory=2,parameterNumber=3 ${data_dir}/multi.grib2`
if [ -z "$parameterNumber" ]; then
echo 'ERROR: grib_get parameterNumber failed'
exit 1
fi
${tools_dir}/grib_copy -w parameterCategory=2,parameterNumber=3 ${data_dir}/multi.grib2 $tmpdata.1
${tools_dir}/grib_compare ${data_dir}/v.grib2 $tmpdata.1
cat > $tmpdata.rules <<EOF
if ( parameterCategory==2 && parameterNumber==3) {
print "[parameterNumber]";
}
EOF
parameterNumber=`${tools_dir}/grib_filter $tmpdata.rules ${data_dir}/multi.grib2`
if [ -z "$parameterNumber" ]; then
echo 'ERROR: grib_filter failed'
exit 1
fi
# Convert to single-field
# -----------------------
infile=${data_dir}/multi.grib2
${tools_dir}/grib_copy $infile $tmpdata
c1=$(${tools_dir}/grib_count $infile)
c2=$(${tools_dir}/grib_count $tmpdata)
[ $c1 -eq 30 ]
[ $c2 -eq 56 ]
# Clean up
rm -f $tmpdata
rm -f $tmpdata.1 $tmpdata.rules
|
