gmt-fork/test/blockmode/etest.sh

#!/usr/bin/env bash
#
# Test to make sure the -E option works as advertised.
# We have data that will fall in to 4 separate blocks
# in a -R0/2/0/2 -I2 -r situation (2x2 blocks)
# In all blocks mean = median = mode = 5.
# With the new mode calculation we are unable to get 5 for the SW block, it comes out as 5.5
# This test is only for blockmode
log=etest.log

cat << EOF > data.d
# Block NW (1 value)
0.4	1.35	5
# Block NE (2 values)
1.3	1.7	4
1.7	1.3	6
# Block SW (5 values)
0.2	0.2	3
0.4	0.4	4
0.6	0.6	5
0.5	0.5	6
0.2	0.7	7
# Block SE (6 values)
1.3	0.2	2
1.4	0.4	3
1.7	0.6	4
1.5	0.6	6
1.2	0.8	7
1.9	0.18	8
EOF
echo "All the mode values should be 5" > $log
echo "Plain modes" >> $log
gmt blockmode -R0/2/0/2 -I1 -r data.d >> $log
echo "Extended modes" >> $log
gmt blockmode -R0/2/0/2 -I1 -r -E data.d >> $log
$AWK '{if (NF == 6 && $3 != 5) print $0}' $log > fail