#!/bin/bash
#nastavi retezec rastry na prazdnou hodnotu, inicializace promenne
rastry="";
#spusti reklasifikaci vrstvy vyuzitin
#vznika tak treci povrch
#3 les absorbuje plyn
#1 vodni tok kvuli okolni vegeteaci take absorbuje castecne plyn
#0, 4, silnice a pole umoznuji dobre sireni plynu s malou absorbci vegetaci
#2 zastavene plochy - snadny prunik plynu dovnitr budov je velmi omezeny
#<<EOF umoznuje vstup dat do prikazu, tak jako by je zadaval uzivatel na klavesnici
r.reclass input=vyuzitin output=prostupnost --o<<EOF
3 = 50
1 = 10
0 4 = 5
2 = 100
end
EOF
#pomocna promena pom je vyuzivana pro omezeni sireni, tak aby v krocich mohla vzniknout animace
pom=1;
#cyklus provadi vypocet akumulace sireni plynu na zaklade treciho povrchu a to vzdy s omezenim
#ktere odpovida aktualni iteraci
for i in `seq 1 10`; do
	#definice omezeni vypoctu pro danou iteraci
	pom=$((100 * $i));
	echo $pom;
	#vypocet akumulace
	r.cost -v input=prostupnost output=sireni$i coordinate=-512356.92,-1112811.31 max_cost=$pom --o
	#sestavovani seznamu vrstev vzniklych vypoctem 
	if [ $i -lt 1 ]; then
		rastry="sireni$i";
		#zpomaleni animace tim, ze se jeden vystup vypoctu zapise jako 9 snimku animace
		for j in `seq 1 10`; do
			rastry="$rastry,sireni$i";
		done
	else
		for j in `seq 1 10`; do
			rastry="$rastry,sireni$i";
		done
	fi
	echo $i;
	
done
#vygenerovani animace na zaklade seznamu vrstev
#vyzaduje netpbm
r.out.mpeg view1=$rastry;
#spusteni animace s vyuzitim externiho programu xmms
mplayer gmovie.mpg