[GRASS-dev] r.sim.water working?

Helena Mitasova hmitaso at unity.ncsu.edu
Sat Oct 17 11:04:13 EDT 2009


Please try to increase number of iterations to several hundred, run  
it without -t
and run it on something smaller (1000x1000) - it looks like it did  
not get
too far when it finished the computation, but it should have given  
you at least
a warning or error message (it might have run out of memory).
I need to try it out, I still run it in older versions of GRASS64,

Helena


Helena Mitasova
Associate Professor
Department of Marine, Earth
and Atmospheric Sciences
1125 Jordan Hall, Campus Box 8208
North Carolina State University
Raleigh NC 27695-8208
http://skagit.meas.ncsu.edu/~helena/



On Oct 17, 2009, at 3:49 AM, cgnicholas at alamedanet.net wrote:

>> From a previous post I realize r.terraflow isnt working in 6.4  
>> yet.., so
> thought I would try a dumbed-down-DEM with r.sim.water ... DEM, mask,
> slope, aspect, etc  look good; everything seem to run, but...no  
> outputs
> are generated.... hmmm...?
> ( I can take this to users' list if this is the wrong place, but  
> seems a
> lot of things are in motion here...)
>
> thanks in advance!
> Chris
> ----------------------------------------------------------------
> [Raster MASK present]
> GRASS 6.4.0RC5 (Haiti):~/tmp > r.sim.water -t elevin=DEM100m at PERMANENT
> dxin=dx at PERMANENT dyin=dy at PERMANENT rain_val=50 infil_val=0.0
> manin_val=0.1 depth=water_depth at PERMANENT  
> disch=water_discharge at PERMANENT
> err=error_map at PERMANENT niter=10 outiter=2 diffc=0.8 hmax=0.3  
> halpha=4.0
> hbeta=0.5 --overwrite
> default nwalk=68487324, rwalk=68487324.000000
>
> Min elevation = 1.00 m
> Max elevation = 3109.00 m
> Mean Source Rate (rainf. excess or sediment) = 0.000002 m/s or kg/m2s
> Mean flow velocity = 0.628966 m/s
> Mean Mannings = 0.606325
> Number of iterations = 3 cells
> Time step = 39.75 s
> [Raster MASK present]
> GRASS 6.4.0RC5 (Haiti):~/tmp > g.gisenv
> LOCATION_NAME=Haiti
> MAPSET=PERMANENT
> DIGITIZER=none
> GISDBASE=/home/cgn/GRASS
> MONITOR=x7
> GRASS_GUI=tcltk
> [Raster MASK present]
> GRASS 6.4.0RC5 (Haiti):~/tmp > g.region -p
> projection: 1 (UTM)
> zone:       0
> datum:      wgs84
> ellipsoid:  a=6378137 es=0.006694379990141316
> north:      2338174.90627377
> south:      1879574.90627377
> west:       499985.21620627
> east:       1246685.21620627
> nsres:      100
> ewres:      100
> rows:       4586
> cols:       7467
> cells:      34243662
> [Raster MASK present]
> GRASS 6.4.0RC5 (Haiti):~/tmp > g.list type=rast
> ----------------------------------------------
> raster files available in mapset <PERMANENT>:
> DEM100m         aspect          dx              dyy
> slope_steepness
> MASK            basinid         dxx             half_basins      
> streams
> accum           curvature       dxy             slope            
> tcurvature
> asin_visual     drainage_dir    dy              slope_length
>
> ----------------------------------------------
> [Raster MASK present]
> GRASS 6.4.0RC5 (Haiti):~/tmp >
>
>
> _______________________________________________
> grass-dev mailing list
> grass-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-dev



More information about the grass-dev mailing list