2
# Test for t.rast.mapcalc
4
export GRASS_OVERWRITE=1
6
# We need to set a specific region in the
7
# @preprocess step of this test. We generate
8
# raster with r.mapcalc and create several space time raster datasets
9
# The region setting should work for UTM and LL test locations
10
g.region s=0 n=80 w=0 e=120 b=0 t=50 res=10 res3=10 -p3
13
r.mapcalc expr="prec_1 = rand(0, 550)" -s
14
r.mapcalc expr="prec_2 = rand(0, 450)" -s
15
r.mapcalc expr="prec_3 = rand(0, 320)" -s
16
r.mapcalc expr="prec_4 = rand(0, 510)" -s
17
r.mapcalc expr="prec_5 = rand(0, 300)" -s
18
r.mapcalc expr="prec_6 = rand(0, 650)" -s
20
t.create type=strds temporaltype=absolute output=precip_abs1 title="A test" descr="A test"
21
t.create type=strds temporaltype=absolute output=precip_abs2 title="A test" descr="A test"
22
t.register -i type=raster input=precip_abs1 maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6 start="2001-01-01" increment="3 months"
23
t.register type=raster input=precip_abs2 maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6
29
t.rast.mapcalc -n inputs=precip_abs1,precip_abs2 output=precip_abs3 \
30
expression=" precip_abs1 + precip_abs2" base=new_prec \
32
t.info type=strds input=precip_abs3
34
t.rast.mapcalc -s inputs=precip_abs1,precip_abs2,precip_abs3 output=precip_abs4 \
35
expression=" (precip_abs1 + precip_abs2) / precip_abs2" base=new_prec \
37
t.info type=strds input=precip_abs4
39
t.rast.mapcalc -s inputs=precip_abs1,precip_abs2 output=precip_abs4 \
40
expression=" (precip_abs1 + precip_abs2) * null()" base=new_prec \
42
t.info type=strds input=precip_abs4
44
t.rast.mapcalc -sn inputs=precip_abs1,precip_abs2 output=precip_abs4 \
45
expression=" (precip_abs1 + precip_abs2) * null()" base=new_prec \
47
t.info type=strds input=precip_abs4
50
g.remove -f type=raster name=prec_1
51
t.rast.mapcalc -sn inputs=precip_abs1,precip_abs2 output=precip_abs4 \
52
expression=" (precip_abs1 + precip_abs2) * null()" base=new_prec \
56
t.remove -rf type=strds input=precip_abs1,precip_abs2,precip_abs3,precip_abs4