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
|
#!/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_paramid_chemid_split_test"
tempGribA=temp1.$label.grib
tempGribB=temp2.$label.grib
tempFilt=temp.$label.filt
sample_grib2=$ECCODES_SAMPLES_PATH/GRIB2.tmpl
# Legacy using older tablesVersion
# We expect paramId 210121
${tools_dir}/grib_set -s productDefinitionTemplateNumber=40,discipline=0,parameterCategory=20,parameterNumber=2,constituentType=5 \
$sample_grib2 $tempGribA
grib_check_key_equals $tempGribA tablesVersion 4
grib_check_key_equals $tempGribA paramId,is_chemical "210121 1"
# The new keys should not be present
result=$( ${tools_dir}/grib_get -f -p chemName,chemId $tempGribA )
[ "$result" = "not_found not_found" ]
# Switch to the newer tablesVersion so now the chemId/paramId split is activated
# Now expect paramId 402000 and chemId 17
${tools_dir}/grib_set -s setLocalDefinition=1,class=a5,tablesVersion=34 $tempGribA $tempGribB
grib_check_key_equals $tempGribB paramId,chemId "402000 17"
grib_check_key_equals $tempGribB chemName "Nitrogen dioxide"
grib_check_key_equals $tempGribB name "Mass mixing ratio"
# The is_chemical key should not be present
result=$( ${tools_dir}/grib_get -f -p is_chemical $tempGribB )
[ "$result" = "1" ]
# Test non-chemical e.g. temperature. Now chem keys are present but unknown
${tools_dir}/grib_set -s tablesVersion=34,setLocalDefinition=1,class=a5,discipline=0,parameterCategory=0,parameterNumber=0 $sample_grib2 $tempGribA
grib_check_key_equals $tempGribA "chemId,chemName,chemShortName" "-1 unknown unknown"
# Conversion in memory using transient 'enableChemSplit'
cat >$tempFilt<<EOF
set tablesVersion = 35;
set productDefinitionTemplateNumber = 50;
set class = 'cr';
set enableChemSplit = 1;
assert( paramId == 402000 );
assert( chemId == 914 );
write;
EOF
${tools_dir}/grib_set -s \
setLocalDefinition=1,class=2,stream=1025,type=2,productDefinitionTemplateNumber=44,constituentType=62003,parameterCategory=20,parameterNumber=2 \
$sample_grib2 $tempGribA
grib_check_key_equals $tempGribA "paramId,shortName" "210249 aermr18"
${tools_dir}/grib_filter -o $tempGribB $tempFilt $tempGribA
grib_check_key_equals $tempGribB "paramId,shortName" "402000 mass_mixrat"
# Clean up
rm -f $tempGribA $tempGribB $tempFilt
|
