Fortran example for codes_copy_key

This commit is contained in:
Shahram Najm 2016-11-10 17:05:02 +00:00
parent c5b2132413
commit 88ecf50464
3 changed files with 60 additions and 0 deletions

View File

@ -45,6 +45,7 @@ list( APPEND tests
bufr_read_temp
bufr_read_tropical_cyclone
bufr_set_keys
bufr_copy_keys
bufr_subset
get_product_kind
)

View File

@ -0,0 +1,34 @@
! Copyright 2005-2016 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.
!
! Description: how to copy a BUFR key from a message to another.
!
program bufr_copy_keys
use eccodes
implicit none
integer :: file1, file2, file3
integer :: ibufr1,ibufr2,ibufr3
call codes_open_file(file1, '../../data/bufr/PraticaTemp.bufr', 'r')
call codes_open_file(file2, '../../data/bufr/aaen_55.bufr', 'r')
call codes_open_file(file3, 'out.bufr_copy_keys_test_f.bufr', 'w')
call codes_bufr_new_from_file(file1, ibufr1)
call codes_bufr_new_from_file(file2, ibufr2)
call codes_clone(ibufr2, ibufr3)
call codes_copy_key(ibufr1, 'bufrHeaderCentre', ibufr3)
call codes_write(ibufr3, file3)
call codes_close_file(file1)
call codes_close_file(file2)
call codes_close_file(file3)
end program bufr_copy_keys

25
examples/F90/bufr_copy_keys.sh Executable file
View File

@ -0,0 +1,25 @@
#!/bin/sh
# Copyright 2005-2016 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.sh
#Define a common label for all the tmp files
label="bufr_copy_keys_test_f"
TEMP=$label.out.bufr
cd ${data_dir}/bufr
#The input and output BUFR files are hardcoded in the f90 example!!!
OUTPUT=out.$label.bufr
${examples_dir}eccodes_f_bufr_copy_keys
res=`${tools_dir}bufr_get -p bufrHeaderCentre $OUTPUT`
[ "$res" = "80" ]
rm -f $OUTPUT