[GRASS-user] Usage of r.sim.water (error - no ouput grid created)

Wed Mar 5 09:18:33 PST 2014

Dear GRASS Users,

Being quite a newbie to GRASS, I have the task to model overlandflow via GIS-tools for simulating flash floods. This lead me to r.sim.water, which I tried to run for a ~5km² area with a 25m raster dgm, but with no success until now.

I have a grid with the DEM, two grids with the first partial derivative x/y, a grid with manning´s coeff, a grid with effective rainfall and that´s it - should work, shouldn´t it?

Do you have a hint for me, why this does not work in my case?
(see parameters in png, log in txt)

THANK YOU!!!
Regards,
Andreas

Mit freundlichen Grüßen
Mag. Andreas Reithofer
Akad. Geoinformatiker
email: Reithofer at ib-humer.at
Tel.: 07732-4146-35
________________________________________________________________________

Ingenieurbüro für Kulturtechnik und Wasserwirtschaft, Consulting Engineer
Feld 16, 4682 Geboltskirchen, Austria
Tel.: 07732-4146, Fax.: DW - 22

www.ib-humer.at
office at ib-humer.at
Fb.Nr.: FN 312673 p

g.proj -c proj4="+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\DEM_leonding_nosinks.tif" band=1 output=tmp13940388217816 --overwrite -o
r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\DEM_leonding_nosinks_1st_order_dx.tif" band=1 output=tmp13940388217817 --overwrite -o
r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\DEM_leonding_nosinks_1st_order_dy.tif" band=1 output=tmp13940388217818 --overwrite -o
r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\BNS_60min.tif" band=1 output=tmp13940388217819 --overwrite -o
r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\dummygrid_null.tif" band=1 output=tmp13940388217820 --overwrite -o
r.external input="K:\Proj_2012\12092_switchOn\work_temp\CLC2006\Raster_100m\m_coeff_25m.tif" band=1 output=tmp13940388217821 --overwrite -o
g.region n=2806700.0 s=2803000.0 e=4636212.5 w=4633100.0 res=25.0
r.sim.water elevin=tmp13940388217816 dxin=tmp13940388217817 dyin=tmp13940388217818 rain=tmp13940388217819 infil=tmp13940388217820 manin=tmp13940388217821 traps=tmp13940388217820 niter="15" outiter="15" diffc="0.8" hmax="0.3" halpha="4.0" hbeta="0.5" -t depth=depthd5667d4a1b7a4136a850a17a410fbc54 disch=dischd5667d4a1b7a4136a850a17a410fbc54 err=errd5667d4a1b7a4136a850a17a410fbc54 --overwrite
g.region rast=depthd5667d4a1b7a4136a850a17a410fbc54
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=depthd5667d4a1b7a4136a850a17a410fbc54 output="K:/Proj_2012/12092_switchOn/work_temp/SAGA_DEM_EU/_output.tif"
g.region rast=dischd5667d4a1b7a4136a850a17a410fbc54
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=dischd5667d4a1b7a4136a850a17a410fbc54 output="C:\Users\user_x\AppData\Local\Temp\processing\2219fb8bc6784cd0903ee0a0f3004e47\disch.tif"
g.region rast=errd5667d4a1b7a4136a850a17a410fbc54
r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=errd5667d4a1b7a4136a850a17a410fbc54 output="C:\Users\user_x\AppData\Local\Temp\processing\b7920c941552441f9b2cba64799bd835\err.tif"

C:\PROGRA~1\QGISDU~1\bin>set HOME=C:\Users\user_x 

C:\PROGRA~1\QGISDU~1\bin>set GISRC=C:\Users\user_x\.qgis2\\processing\processing.gisrc 

C:\PROGRA~1\QGISDU~1\bin>set GRASS_SH=C:/PROGRA~1/QGISDU~1/apps\msys\bin\sh.exe 

C:\PROGRA~1\QGISDU~1\bin>set PATH=C:/PROGRA~1/QGISDU~1/apps\msys\bin;C:/PROGRA~1/QGISDU~1/apps\msys\lib;C:/PROGRA~1/QGISDU~1/apps/grass/GRASS-~1.3/bin;C:/PROGRA~1/QGISDU~1/apps/grass/GRASS-~1.3/scripts;C:/PROGRA~1/QGISDU~1/apps/qgis/./grass/scripts/;C:/PROGRA~1/QGISDU~1/bin;C:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\bin;C:\PROGRA~1\QGISDU~1\apps\msys\bin;C:\PROGRA~1\QGISDU~1\apps\Python27\Scripts;C:\PROGRA~1\QGISDU~1\bin;C:\Windows\system32;C:\Windows;C:\Windows\WBem;C:\PROGRA~1\QGISDU~1\apps\msys\bin;C:\PROGRA~1\QGISDU~1\apps\qgis\bin;C:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\lib 

C:\PROGRA~1\QGISDU~1\bin>set WINGISBASE=C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3 

C:\PROGRA~1\QGISDU~1\bin>set GISBASE=C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3 

C:\PROGRA~1\QGISDU~1\bin>set GRASS_PROJSHARE=C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\share\proj 

C:\PROGRA~1\QGISDU~1\bin>set GRASS_MESSAGE_FORMAT=gui 

C:\PROGRA~1\QGISDU~1\bin>if "" == "" set PATH=C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\bin;C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\lib;C:/PROGRA~1/QGISDU~1/apps\msys\bin;C:/PROGRA~1/QGISDU~1/apps\msys\lib;C:/PROGRA~1/QGISDU~1/apps/grass/GRASS-~1.3/bin;C:/PROGRA~1/QGISDU~1/apps/grass/GRASS-~1.3/scripts;C:/PROGRA~1/QGISDU~1/apps/qgis/./grass/scripts/;C:/PROGRA~1/QGISDU~1/bin;C:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\bin;C:\PROGRA~1\QGISDU~1\apps\msys\bin;C:\PROGRA~1\QGISDU~1\apps\Python27\Scripts;C:\PROGRA~1\QGISDU~1\bin;C:\Windows\system32;C:\Windows;C:\Windows\WBem;C:\PROGRA~1\QGISDU~1\apps\msys\bin;C:\PROGRA~1\QGISDU~1\apps\qgis\bin;C:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\lib 

C:\PROGRA~1\QGISDU~1\bin>if not "" == "" set PATH=C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\bin;C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\lib;;C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\bin;C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\lib;C:/PROGRA~1/QGISDU~1/apps\msys\bin;C:/PROGRA~1/QGISDU~1/apps\msys\lib;C:/PROGRA~1/QGISDU~1/apps/grass/GRASS-~1.3/bin;C:/PROGRA~1/QGISDU~1/apps/grass/GRASS-~1.3/scripts;C:/PROGRA~1/QGISDU~1/apps/qgis/./grass/scripts/;C:/PROGRA~1/QGISDU~1/bin;C:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\bin;C:\PROGRA~1\QGISDU~1\apps\msys\bin;C:\PROGRA~1\QGISDU~1\apps\Python27\Scripts;C:\PROGRA~1\QGISDU~1\bin;C:\Windows\system32;C:\Windows;C:\Windows\WBem;C:\PROGRA~1\QGISDU~1\apps\msys\bin;C:\PROGRA~1\QGISDU~1\apps\qgis\bin;C:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\lib 

C:\PROGRA~1\QGISDU~1\bin>set GRASS_VERSION=6.4.0 

C:\PROGRA~1\QGISDU~1\bin>if not "" == "" goto langset 

C:\PROGRA~1\QGISDU~1\bin>FOR /F "usebackq delims==" %i IN (`"C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\etc\winlocale"`) DO @set LANG=%i 

C:\PROGRA~1\QGISDU~1\bin>set PATHEXT=.COM;.EXE;.BAT;.CMD;.VBS;.VBE;.JS;.JSE;.WSF;.WSH;.MSC;.PY 

C:\PROGRA~1\QGISDU~1\bin>set PYTHONPATH=C:/PROGRA~1/QGISDU~1/apps/grass/GRASS-~1.3/etc/python;;C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\etc\python;C:/PROGRA~1/QGISDU~1/apps\grass\grass-6.4.3\etc\wxpython\n 

C:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="MAPSET=PERMANENT" 

C:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="LOCATION=temp_location" 

C:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="LOCATION_NAME=temp_location" 

C:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="GISDBASE=C:\Users\user_x\AppData\Local\Temp\processing\grassdata" 

C:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="GRASS_GUI=text" 

C:\PROGRA~1\QGISDU~1\bin>g.proj -c proj4="+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs" 

GRASS_INFO_MESSAGE(4204,1): Die Standard Region wurde auf die neue Projektion aktualisiert. Wenn Sie aber mehrere Mapsets haben, sollten Sie `g.region -d` in jedem ausführen, um die Einstellungen von der Standardregion zu übernehmen. 
GRASS_INFO_END(4204,1) 

GRASS_INFO_MESSAGE(4204,2): Projektionsinformationen aktualisiert 
GRASS_INFO_END(4204,2) 

C:\PROGRA~1\QGISDU~1\bin>r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\DEM_leonding_nosinks.tif" band=1 output=tmp13940388217816 --overwrite -o 

GRASS_INFO_WARNING(2564,1): Übersteuere die Überprüfung der Projektion. 
GRASS_INFO_END(2564,1) 

GRASS_INFO_MESSAGE(2564,2): Importing band 1 of 1... 
GRASS_INFO_END(2564,2) 

GRASS_INFO_MESSAGE(2564,3): Link zur Rasterkarte erzeugt 
GRASS_INFO_END(2564,3) 

C:\PROGRA~1\QGISDU~1\bin>r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\DEM_leonding_nosinks_1st_order_dx.tif" band=1 output=tmp13940388217817 --overwrite -o 

GRASS_INFO_WARNING(4232,1): Übersteuere die Überprüfung der Projektion. 
GRASS_INFO_END(4232,1) 

GRASS_INFO_MESSAGE(4232,2): Importing band 1 of 1... 
GRASS_INFO_END(4232,2) 

GRASS_INFO_MESSAGE(4232,3): Link zur Rasterkarte erzeugt 
GRASS_INFO_END(4232,3) 

C:\PROGRA~1\QGISDU~1\bin>r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\DEM_leonding_nosinks_1st_order_dy.tif" band=1 output=tmp13940388217818 --overwrite -o 

GRASS_INFO_WARNING(7000,1): Übersteuere die Überprüfung der Projektion. 
GRASS_INFO_END(7000,1) 

GRASS_INFO_MESSAGE(7000,2): Importing band 1 of 1... 
GRASS_INFO_END(7000,2) 

GRASS_INFO_MESSAGE(7000,3): Link zur Rasterkarte erzeugt 
GRASS_INFO_END(7000,3) 

C:\PROGRA~1\QGISDU~1\bin>r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\BNS_60min.tif" band=1 output=tmp13940388217819 --overwrite -o 

GRASS_INFO_WARNING(6372,1): Übersteuere die Überprüfung der Projektion. 
GRASS_INFO_END(6372,1) 

GRASS_INFO_MESSAGE(6372,2): Importing band 1 of 1... 
GRASS_INFO_END(6372,2) 

GRASS_INFO_MESSAGE(6372,3): Link zur Rasterkarte erzeugt 
GRASS_INFO_END(6372,3) 

C:\PROGRA~1\QGISDU~1\bin>r.external input="K:\Proj_2012\12092_switchOn\work_temp\SAGA_DEM_EU\dummygrid_null.tif" band=1 output=tmp13940388217820 --overwrite -o 

GRASS_INFO_WARNING(5760,1): Übersteuere die Überprüfung der Projektion. 
GRASS_INFO_END(5760,1) 

GRASS_INFO_MESSAGE(5760,2): Importing band 1 of 1... 
GRASS_INFO_END(5760,2) 

GRASS_INFO_MESSAGE(5760,3): Link zur Rasterkarte erzeugt 
GRASS_INFO_END(5760,3) 

C:\PROGRA~1\QGISDU~1\bin>r.external input="K:\Proj_2012\12092_switchOn\work_temp\CLC2006\Raster_100m\m_coeff_25m.tif" band=1 output=tmp13940388217821 --overwrite -o 

GRASS_INFO_WARNING(6596,1): Übersteuere die Überprüfung der Projektion. 
GRASS_INFO_END(6596,1) 

GRASS_INFO_MESSAGE(6596,2): Importing band 1 of 1... 
GRASS_INFO_END(6596,2) 

GRASS_INFO_MESSAGE(6596,3): Link zur Rasterkarte erzeugt 
GRASS_INFO_END(6596,3) 

C:\PROGRA~1\QGISDU~1\bin>g.region n=2806700.0 s=2803000.0 e=4636212.5 w=4633100.0 res=25.0 

C:\PROGRA~1\QGISDU~1\bin>r.sim.water elevin=tmp13940388217816 dxin=tmp13940388217817 dyin=tmp13940388217818 rain=tmp13940388217819 infil=tmp13940388217820 manin=tmp13940388217821 traps=tmp13940388217820 niter="15" outiter="15" diffc="0.8" hmax="0.3" halpha="4.0" hbeta="0.5" -t depth=depthd5667d4a1b7a4136a850a17a410fbc54 disch=dischd5667d4a1b7a4136a850a17a410fbc54 err=errd5667d4a1b7a4136a850a17a410fbc54 --overwrite 

GRASS_INFO_MESSAGE(6308,1): Standard nwalk=37000, rwalk=37000.000000 
GRASS_INFO_END(6308,1) 

GRASS_INFO_MESSAGE(6308,2): Min elevation = 273.88 m 

GRASS_INFO_MESSAGE(6308,2): Max elevation = 506.77 m 

GRASS_INFO_END(6308,2) 

GRASS_INFO_MESSAGE(6308,3): Mean Source Rate (rainf. excess or sediment) = 0.000009 m/s or kg/m2s 

GRASS_INFO_END(6308,3) 

GRASS_INFO_MESSAGE(6308,4): Mean flow velocity = 2.819856 m/s 

GRASS_INFO_END(6308,4) 

GRASS_INFO_MESSAGE(6308,5): Mean Mannings = 0.081520 

GRASS_INFO_END(6308,5) 

GRASS_INFO_MESSAGE(6308,6): Number of iterations = 101 cells 

GRASS_INFO_END(6308,6) 

GRASS_INFO_MESSAGE(6308,7): Time step = 2.21 s 

GRASS_INFO_END(6308,7) 

C:\PROGRA~1\QGISDU~1\bin>g.region rast=depthd5667d4a1b7a4136a850a17a410fbc54 

GRASS_INFO_ERROR(6300,1): Rasterkarte konnte nicht gefunden werden. 
GRASS_INFO_END(6300,1) 

C:\PROGRA~1\QGISDU~1\bin>r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=depthd5667d4a1b7a4136a850a17a410fbc54 output="K:/Proj_2012/12092_switchOn/work_temp/SAGA_DEM_EU/_output.tif" 

GRASS_INFO_ERROR(7036,1): Rasterkarte oder Gruppe nicht gefunden. 
GRASS_INFO_END(7036,1) 

C:\PROGRA~1\QGISDU~1\bin>g.region rast=dischd5667d4a1b7a4136a850a17a410fbc54 

GRASS_INFO_ERROR(6508,1): Rasterkarte konnte nicht gefunden werden. 
GRASS_INFO_END(6508,1) 

C:\PROGRA~1\QGISDU~1\bin>r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=dischd5667d4a1b7a4136a850a17a410fbc54 output="C:\Users\user_x\AppData\Local\Temp\processing\2219fb8bc6784cd0903ee0a0f3004e47\disch.tif" 

GRASS_INFO_ERROR(2612,1): Rasterkarte oder Gruppe nicht gefunden. 
GRASS_INFO_END(2612,1) 

C:\PROGRA~1\QGISDU~1\bin>g.region rast=errd5667d4a1b7a4136a850a17a410fbc54 

GRASS_INFO_ERROR(1436,1): Rasterkarte konnte nicht gefunden werden. 
GRASS_INFO_END(1436,1) 

C:\PROGRA~1\QGISDU~1\bin>r.out.gdal -c createopt="TFW=YES,COMPRESS=LZW" input=errd5667d4a1b7a4136a850a17a410fbc54 output="C:\Users\user_x\AppData\Local\Temp\processing\b7920c941552441f9b2cba64799bd835\err.tif" 

GRASS_INFO_ERROR(2244,1): Rasterkarte oder Gruppe nicht gefunden. 
GRASS_INFO_END(2244,1) 

C:\PROGRA~1\QGISDU~1\bin>exit 
Converting outputs
Loading resulting layers
Algorithm r.sim.water - Overland flow hydrologic simulation using path sampling method (SIMWE). finished
-------------- n?chster Teil --------------
Ein Dateianhang mit HTML-Daten wurde abgetrennt...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20140305/5537b43f/attachment-0001.html>
-------------- n?chster Teil --------------
Ein Dateianhang mit Bin?rdaten wurde abgetrennt...
Dateiname   : image001.gif
Dateityp    : image/gif
Dateigr??e  : 5601 bytes
Beschreibung: nicht verf?gbar
URL         : <http://lists.osgeo.org/pipermail/grass-user/attachments/20140305/5537b43f/attachment-0001.gif>
-------------- n?chster Teil --------------
Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt.
Name: logfile.txt
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20140305/5537b43f/attachment-0001.txt>
-------------- n?chster Teil --------------
Ein Dateianhang mit Bin?rdaten wurde abgetrennt...
Dateiname   : parameters.png
Dateityp    : image/png
Dateigr??e  : 73470 bytes
Beschreibung: nicht verf?gbar
URL         : <http://lists.osgeo.org/pipermail/grass-user/attachments/20140305/5537b43f/attachment-0001.png>