~ubuntu-branches/ubuntu/vivid/grass/vivid-proposed

« back to all changes in this revision

Viewing changes to temporal/t.rast.mapcalc/test.t.rast.mapcalc.sh

  • Committer: Package Import Robot
  • Author(s): Bas Couwenberg
  • Date: 2015-02-20 23:12:08 UTC
  • mfrom: (8.2.6 experimental)
  • Revision ID: package-import@ubuntu.com-20150220231208-1u6qvqm84v430b10
Tags: 7.0.0-1~exp1
* New upstream release.
* Update python-ctypes-ternary.patch to use if/else instead of and/or.
* Drop check4dev patch, rely on upstream check.
* Add build dependency on libpq-dev to grass-dev for libpq-fe.h.
* Drop patches applied upstream, refresh remaining patches.
* Update symlinks for images switched from jpg to png.

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
#!/bin/sh
 
2
# Test for t.rast.mapcalc 
 
3
 
 
4
export GRASS_OVERWRITE=1
 
5
 
 
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
 
11
 
 
12
# Generate data
 
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
 
19
 
 
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
 
24
 
 
25
t.info precip_abs1
 
26
t.info precip_abs2
 
27
 
 
28
# The first @test
 
29
t.rast.mapcalc -n inputs=precip_abs1,precip_abs2 output=precip_abs3 \
 
30
           expression=" precip_abs1 + precip_abs2" base=new_prec \
 
31
           method=equal nprocs=5
 
32
t.info type=strds input=precip_abs3
 
33
 
 
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 \
 
36
           method=equal nprocs=5
 
37
t.info type=strds input=precip_abs4
 
38
 
 
39
t.rast.mapcalc -s inputs=precip_abs1,precip_abs2 output=precip_abs4 \
 
40
           expression=" (precip_abs1 + precip_abs2) * null()" base=new_prec \
 
41
           method=equal nprocs=5
 
42
t.info type=strds input=precip_abs4
 
43
 
 
44
t.rast.mapcalc -sn inputs=precip_abs1,precip_abs2 output=precip_abs4 \
 
45
           expression=" (precip_abs1 + precip_abs2) * null()" base=new_prec \
 
46
           method=equal nprocs=5
 
47
t.info type=strds input=precip_abs4
 
48
 
 
49
# Let the test fail
 
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 \
 
53
           method=equal nprocs=5
 
54
 
 
55
# @postprocess
 
56
t.remove -rf type=strds input=precip_abs1,precip_abs2,precip_abs3,precip_abs4