| 123456789101112131415161718192021222324252627282930313233343536373839 |
- #!/bin/sh
- # Test for t.rast3d.mapcalc
- # We need to set a specific region in the
- # @preprocess step of this test.
- # 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
- # Generate data
- r3.mapcalc --o expr="prec_1 = rand(0, 550)" -s
- r3.mapcalc --o expr="prec_2 = rand(0, 450)" -s
- r3.mapcalc --o expr="prec_3 = rand(0, 320)" -s
- r3.mapcalc --o expr="prec_4 = rand(0, 510)" -s
- r3.mapcalc --o expr="prec_5 = rand(0, 300)" -s
- r3.mapcalc --o expr="prec_6 = rand(0, 650)" -s
- t.create --o type=str3ds temporaltype=absolute output=precip_abs1 title="A test" descr="A test"
- t.create --o type=str3ds temporaltype=absolute output=precip_abs2 title="A test" descr="A test"
- t.register -i type=raster_3d input=precip_abs1 maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6 start="2001-01-01" increment="3 months"
- t.register --o type=raster_3d input=precip_abs2 maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6
- t.info type=str3ds input=precip_abs1
- t.info type=str3ds input=precip_abs2
- # The first @test
- t.rast3d.mapcalc --o --v -n inputs=precip_abs1,precip_abs2 output=precip_abs3 \
- expression=" precip_abs1 + precip_abs2" base=new_prec \
- method=equal nprocs=6
- t.info type=str3ds input=precip_abs3
- t.rast3d.mapcalc --o --v -s inputs=precip_abs1,precip_abs2,precip_abs3 output=precip_abs4 \
- expression=" (precip_abs1 + precip_abs2) / precip_abs2" base=new_prec \
- method=equal nprocs=6
- t.info type=str3ds input=precip_abs4
- t.unregister type=raster_3d maps=prec_1,prec_2,prec_3,prec_4,prec_5,prec_6
- t.unregister type=raster_3d maps=new_prec_1,new_prec_2,new_prec_3,new_prec_4,new_prec_5,new_prec_6
- t.remove type=str3ds input=precip_abs1,precip_abs2,precip_abs3,precip_abs4
|
