1#!/bin/sh 2# This is a test to register and unregister raster3d maps in 3# space time raster3d datasets 4# The raster3d maps will be registered in different space time raster3d 5# datasets 6 7# We need to set a specific region in the 8# @preprocess step of this test. We generate 9# 3d raster with r3.mapcalc and create two space time raster3d datasets 10# with absolute time 11# The region setting should work for UTM and LL test locations 12g.region s=0 n=80 w=0 e=120 b=0 t=50 res=10 res3=10 -p3 13 14# Generate data 15r3.mapcalc --o expr="volume_1 = rand(0, 550)" -s 16r3.mapcalc --o expr="volume_2 = rand(0, 450)" -s 17r3.mapcalc --o expr="volume_3 = rand(0, 320)" -s 18r3.mapcalc --o expr="volume_4 = rand(0, 510)" -s 19r3.mapcalc --o expr="volume_5 = rand(0, 300)" -s 20r3.mapcalc --o expr="volume_6 = rand(0, 650)" -s 21 22# The first @test 23# We create the space time raster3d dataset and register the raster3d maps with absolute time interval 24 25t.create --v --o type=str3ds temporaltype=absolute output=volume_abs1 title="A test" descr="A test" 26t.create --v --o type=str3ds temporaltype=absolute output=volume_abs2 title="A test" descr="A test" 27t.create --v --o type=str3ds temporaltype=absolute output=volume_abs3 title="A test" descr="A test" 28t.create --v --o type=str3ds temporaltype=absolute output=volume_abs4 title="A test" descr="A test" 29t.create --v --o type=str3ds temporaltype=absolute output=volume_abs5 title="A test" descr="A test" 30t.create --v --o type=str3ds temporaltype=absolute output=volume_abs6 title="A test" descr="A test" 31t.create --v --o type=str3ds temporaltype=absolute output=volume_abs7 title="A test" descr="A test" 32 33t.register type=raster_3d --o --v -i input=volume_abs1 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="1 seconds" 34t.info type=str3ds input=volume_abs1 35t.unregister --v type=raster_3d input=volume_abs1 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 36t.info type=str3ds input=volume_abs1 37 38t.register type=raster_3d --o --v -i input=volume_abs2 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="20 seconds, 5 minutes" 39t.info type=str3ds input=volume_abs2 40r3.info volume_1 41r3.info volume_2 42r3.info volume_3 43r3.info volume_4 44r3.info volume_5 45r3.info volume_6 46 47t.register type=raster_3d --o --v -i input=volume_abs3 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="8 hours" 48t.info type=str3ds input=volume_abs3 49t.unregister --v type=raster_3d maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 50t.info type=str3ds input=volume_abs3 51 52t.register type=raster_3d --o input=volume_abs4 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="3 days" 53t.info type=str3ds input=volume_abs4 54 55t.register type=raster_3d --o input=volume_abs5 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="4 weeks" 56t.info type=str3ds input=volume_abs5 57 58t.register type=raster_3d --o input=volume_abs6 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-08-01" increment="2 months" 59t.info type=str3ds input=volume_abs6 60 61t.register type=raster_3d --o input=volume_abs7 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="20 years, 3 months, 1 days, 4 hours" 62t.info type=str3ds input=volume_abs7 63# Register with different valid time again 64t.register type=raster_3d --o input=volume_abs7 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="99 years, 9 months, 9 days, 9 hours" 65t.info type=str3ds input=volume_abs7 66# Register with different valid time again creating intervals 67t.register type=raster_3d -i --o input=volume_abs7 maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 start="2001-01-01" increment="99 years, 9 months, 9 days, 9 hours" 68t.info type=str3ds input=volume_abs7 69 70t.unregister --v type=raster_3d maps=volume_1,volume_2,volume_3,volume_4,volume_5,volume_6 71t.remove --v type=str3ds input=volume_abs1,volume_abs2,volume_abs3,volume_abs4,volume_abs5,volume_abs6,volume_abs7 72