o


r.sim.water -t elevin=sel9f dxin=sel9f.dx dyin=sel9f.dy rain=rain infil=infilt manin=man.05
r.sim.water [-mt] elevin=name dxin=name dyin=name rain=name
   infil=name [traps=name] manin=name [sites=name] [depth=name]
   [disch=name] [err=name] [outwalk=name] [nwalk=value] [niter=value]
   [outiter=value] [density=value] [diffc=value] [hmax=value]
   [halpha=value] [hbeta=value]

r.sim.sediment [-mt] elevin=name wdepth=name dxin=name dyin=name
   detin=name tranin=name tauin=name manin=name [sites=name] [tc=name]
   [et=name] [conc=name] [flux=name] [erdep=name] [nwalk=value]
   [niter=value] [outiter=value] [density=value] [diffc=value]

rain01=0.00001
man05=0.05
infil0=0.

r.sim.water -t elevin=elev.6 dxin=eldx.6 dyin=eldy.6 rain=rain01 infil=infil0 manin=man05 depth=hhtest.6
r.sim.water2 -t elevin=elev.6 dxin=eldx.6 dyin=eldy.6 rain=rain01 infil=infil0 manin=man05 depth=hhtestnul.6

 306  r.mapcalc tranin=0.001
  307  r.mapcalc detin=0.001
  308  r.mapcalc tauin=0.01

r.sim.sediment -t elevin=elev.6 dxin=eldx.6 dyin=eldy.6 rain=rain01 infil=infil0 manin=man05 depth=hhtest.6
r.sim.sediment2 -t elevin=elev.6 dxin=eldx.6 dyin=eldy.6 wdepth=hhtestnul.6.01199 detin=detin 
tranin=tranin tauin=tauin manin=man05 tc=tc.test2 et=et.test2 flux=flux.test2 erdep=erdep.test2

r.sim.sediment2 -t elevin=elev.6 dxin=eldx.6 dyin=eldy.6 wdepth=hhtestnul.6.01199 detin=detin
tranin=tranin tauin=tauin manin=man05 tc=tc.test2 et=et.test2 flux=flux.test2 erdep=erdep.test2

r.sim.sediment2 elevin=elev.6 dxin=eldx.6 dyin=eldy.6 wdepth=hhtestnul.6.01199 detin=detin
tranin=tranin tauin=tauin manin=man05 tc=tc.test3 et=et.test3 flux=flux.test3 erdep=erdep.test3

r.sim.sediment2 elevin=elev.6 dxin=eldx.6 dyin=eldy.6 wdepth=hhtestnul.6.01199 detin=detin
tranin=tranin tauin=tauin manin=man05 tc=tc.test4 et=et.test4 flux=flux.test4 erdep=erdep.test4

-t does not work, without it - it is OK.

r.mapcalc veloc="(1./man05)*exp(hhtestnul.6.01199,0.67)*sqrt(sin(elev6.slp))"
o======================================================================================

iblock=1 i=2399 miter=9599 nwalk=2120088 nwalka=7145 ERROR: cannot find file flux.test2
tries to write a color table for flux but flux was not written

NOTE deltap is HARD CODED for sediment

zmin,zmax 97.190437 131.692810
 simean,vmean,chmean,deltap 0.010124 0.202638 19.833569 0.100000
 sigmax,vmax,deltaw 3.849001 28.955183 0.007178
 MITER, timec 9599 1.250000
 maxwa, nblock 2000000 1
 barea,sarea,rwalk,sisum: 3.344509 0.371612 2000000.000000 2404.097944
 number of written walkers: 10600
 nwalk, maxw 2120088 7000000
 walkwe (walk weight),frac 2000000.000002 0.800000
 deldif,factor 0.252982 1.202049e-04

BUGS

r.sim.sediment
-t flag does not work (does not write the first file and then color table assignment cannot
find it and gives fatal error)

deltap is hardcoded at 0.1 - needs to be computed based on resolution and velocity
give warning if water depth is 

Jaro,

poslat  Lubosov material
sigma[k][l] = exp(-sigma[k][l] * deltap * slope[k][l]); /* not clear what's here :-\ */

v1[k][l] *= deltap; /* normalization */ toto je draha ktoru prejde prechadzka za jeden casovy krok,
otestuj ci nie je vacsia ako resolution - ak ano daj warning!


UNDEF is -9999, gives crazy values in erdep, sed. flow is OK
Change UNDEF to something else - theoretically it is OK, but somebody could use it
by mistake in the input parameters

Doxygen notes

pekny priklad dokumentacie (ale nie doxygen) je
libes/gis

doxygen je v libes/raster
pozri rastergraphicslib.dox - kombinacia html a latex