Ticket #660: test_gts_headers

File test_gts_headers, 994 bytes (added by warrick, 4 years ago)

Script used to compare ECBUFR way of setting GTS headers with new way using perl script

Line 
1#!/bin/sh
2export ECCODES_BUFR_SET_TO_MISSING_IF_OUT_OF_RANGE=1
3
4# 1. run ropp2bufr -gi on test data
5old_dir=/data/users/fwarrick/ROPP/ropp_src/branches/dev/Share/ROPP10_prototype/ropp_io/tools
6$old_dir/ropp2bufr -gi ./input_data.nc -o old_result.bufr
7
8# 1. run ropp2bufr -g on test data
9old_dir=/data/users/fwarrick/ROPP/ropp_src/branches/dev/Share/ROPP10_prototype/ropp_io/tools
10$old_dir/ropp2bufr -g ./input_data.nc -o old_result_without_iph.bufr
11
12# 2a. run ropp2bufr_eccodes on test data
13new_dir=/data/users/fwarrick/ROPP/ropp_src/branches/dev/Share/ROPP10_ecCodes/ropp_io/tools
14$new_dir/ropp2bufr_eccodes ./input_data.nc -o new_result.bufr
15
16# 2b. run gtsheaders_bufr.pl on result of 2a
17perl ./gtsheaders_bufr.pl new_result.bufr new_result_plus_perl.bufr IUTB14 EKMI
18#190117
19
20# 3 compare octal dumps of each BUFR file
21od -c old_result.bufr > old_result_od.txt
22od -c old_result_without_iph.bufr > old_result_without_iph_od.txt
23od -c new_result_plus_perl.bufr > new_result_plus_perl_od.txt