#!/bin/bash
#	$Id: io.sh 12526 2013-11-16 20:38:52Z pwessel $
#
# Test that the -i, -o, -bi and bo works OK.

# 1. Prepare input files with 8 columns, both ascii and binary
cat << EOF > ascii_i.txt
0	4.1	0	0	4.1	0	41 	1
3	4.3	9	3	1.3	90	13 	91
4.67	4	6	4.67	4	60	4 	61
0.69	3.33	4	0.69	3.33	40	33	41
2	3.3	5	2	3.3	50	33 	51
5	3.35	8	5	3.35	80	35	81
0.25	2.25	3	0.25	2.25	30	25	31
3	2	7	3	2	70	2	71
1.33	1.8	2	1.33	1.8	20	18 	21
0.65	0.7	1	0.65	0.7	10	7 	11
EOF
gmt gmtconvert ascii_i.txt -bod > bin_i.b

# 2. do basic gmt blockmean ascii/bin i/o with no -i/-o
gmt blockmean -R0/5/0/5 -I1 -r ascii_i.txt > ascii_o.txt
gmt blockmean -R0/5/0/5 -I1 -r bin_i.b -bi8d -bod | gmt gmtconvert -bi3d > bin_o.txt
gmt gmtmath -T -Sl ascii_o.txt bin_o.txt SUB SUM = io_answer.txt

# 3. Same as 2, but with selecting cols 3-5 via -i
gmt blockmean -R0/5/0/5 -I1 -r ascii_i.txt -i3-5 > ascii_o.txt
gmt blockmean -R0/5/0/5 -I1 -r bin_i.b -bi8d -i3-5 -bod | gmt gmtconvert -bi3d > bin_o.txt
gmt gmtmath -T -Sl ascii_o.txt bin_o.txt SUB SUM = >> io_answer.txt

# 4. Same 2-3, but just output cols 2,0 via -o
gmt blockmean -R0/5/0/5 -I1 -r ascii_i.txt -o2,0 > ascii_o.txt
gmt blockmean -R0/5/0/5 -I1 -r bin_i.b -bi8d -o2,0 -bo2d | gmt gmtconvert -bi2d > bin_o.txt
gmt gmtmath -T -Sl ascii_o.txt bin_o.txt SUB SUM = >> io_answer.txt

# 5. Same 5, but with selecting cols 3-5 via -i and output cols 2,0 via -o
gmt blockmean -R0/5/0/5 -I1 -r ascii_i.txt -i3-5 -o2,0 > ascii_o.txt
gmt blockmean -R0/5/0/5 -I1 -r bin_i.b -bi8d -i3-5 -o2,0 -bo2d | gmt gmtconvert -bi2d > bin_o.txt
gmt gmtmath -T -Sl ascii_o.txt bin_o.txt SUB SUM = >> io_answer.txt

diff io_answer.txt  "${src:-.}"/io_answer.txt --strip-trailing-cr > fail
#rm -f ascii_i.txt ascii_o.txt bin_o.txt bin_i.b