#!/bin/sh
# This is a test to list raster maps of a space time raster dataset

# We need to set a specific region in the
# @preprocess step of this test. We generate
# raster maps with r.mapcalc and create a space time raster datasets
# The region setting should work for UTM and LL test locations
g.region s=0 n=80 w=0 e=120 b=0 t=50 res=10 res3=10 -p3

export GRASS_OVERWRITE=1

# Generate data
r.mapcalc expr="prec_1 = rand(0, 550)" -s
r.mapcalc expr="prec_2 = rand(0, 450)" -s
r.mapcalc expr="prec_3 = rand(0, 320)" -s
r.mapcalc expr="prec_4 = rand(0, 510)" -s
r.mapcalc expr="prec_5 = rand(0, 300)" -s
r.mapcalc expr="prec_6 = rand(0, 650)" -s

n1=`g.tempfile pid=1 -d`
n2=`g.tempfile pid=2 -d`
n3=`g.tempfile pid=3 -d`
n4=`g.tempfile pid=4 -d`
n5=`g.tempfile pid=5 -d`

cat > "${n1}" << EOF
prec_1
prec_2
prec_3
prec_4
prec_5
prec_6
EOF

cat > "${n2}" << EOF
prec_1|2001-01-01
prec_2|2001-03-01
prec_3|2001-04-01
prec_4|2001-05-01
prec_5|2001-08-01
prec_6|2001-09-01
EOF

cat > "${n3}" << EOF
prec_1|2001-01-01|2001-04-01
prec_2|2001-05-01|2001-07-01
prec_3|2001-08-01|2001-10-01
prec_4|2001-11-01|2002-01-01
prec_5|2002-02-01|2002-04-01
prec_6|2002-05-01|2002-07-01
EOF

cat > "${n4}" << EOF
prec_1|2001-01-01|2001-07-01
prec_2|2001-02-01|2001-04-01
prec_3|2001-03-01|2001-04-01
prec_4|2001-04-01|2001-06-01
prec_5|2001-05-01|2001-06-01
prec_6|2001-06-01|2001-07-01
EOF

cat > "${n5}" << EOF
prec_1|2001-01-01|2001-03-11
prec_2|2001-02-01|2001-04-01
prec_3|2001-03-01|2001-06-02
prec_4|2001-04-01|2001-04-01
prec_5|2001-05-01|2001-05-01
prec_6|2001-06-01|2001-07-01
EOF

t.create type=strds temporaltype=absolute output=precip_abs0 title="A test with input files" descr="A test with input files"

# The @test
t.register type=raster -i input=precip_abs0 file="${n1}" start="2001-01-01" increment="1 month"
t.rast.list  separator=" | " method=comma     input=precip_abs0
t.rast.list  input=precip_abs0
t.rast.list  separator=" | " method=cols      input=precip_abs0
t.rast.list  separator=" | " method=delta     input=precip_abs0
t.rast.list  separator=" | " method=deltagaps input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="2 months"
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="1 day"

t.register type=raster input=precip_abs0 file="${n1}" start="2001-01-01" increment="1 month"
t.rast.list  separator=" | " method=comma     input=precip_abs0
t.rast.list  input=precip_abs0
t.rast.list  separator=" | " method=cols      input=precip_abs0
t.rast.list  separator=" | " method=delta     input=precip_abs0
t.rast.list  separator=" | " method=deltagaps input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="2 months"
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="6 days"

t.register type=raster input=precip_abs0 file="${n2}"
t.rast.list  separator=" | " method=comma     input=precip_abs0
t.rast.list  input=precip_abs0
t.rast.list  separator=" | " method=cols      input=precip_abs0
t.rast.list  separator=" | " method=delta     input=precip_abs0
t.rast.list  separator=" | " method=deltagaps input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="2 months"
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="6 days"

t.register type=raster input=precip_abs0 file="${n3}"
t.rast.list  separator=" | " method=comma     input=precip_abs0
t.rast.list  separator=" | " method=delta     input=precip_abs0
t.rast.list  separator=" | " method=deltagaps input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="2 months"
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="6 days"

t.register type=raster input=precip_abs0 file="${n4}"
t.rast.list  separator=" | " method=comma     input=precip_abs0
t.rast.list  separator=" | " method=delta     input=precip_abs0
t.rast.list  separator=" | " method=deltagaps input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="2 months"
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="6 days"

t.register type=raster input=precip_abs0 file="${n5}"
t.rast.list  separator=" | " method=comma     input=precip_abs0
t.rast.list  input=precip_abs0
t.rast.list  separator=" | " method=cols      input=precip_abs0
t.rast.list  separator=" | " method=delta     input=precip_abs0
t.rast.list  separator=" | " method=deltagaps input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="8 months"
t.rast.list  separator=" | " method=gran      input=precip_abs0 gran="13 days"

t.unregister type=raster maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6
t.remove type=strds input=precip_abs0
