[GRASS-user] r.damflood question

[Codrina Maria Ilie]

Hi Maxi,

Yes, it is used and works well :)
You remember correctly, that was what I used.

Thank you and your colleague for the fine work at r.damflood.

Cheers,
Codrina


On 30/08/2017 11:56, Massimiliano Cannata wrote:
> Hi, glad to see the module is used :-D
> 
> If I correctly remember, you have to use a DEM with bathimetry and dam, 
> then the dam break is passed trough a raster with decrements of dam 
> height (it is use to calculate the new DEM after the dam breck by DEM - 
> BREAK).
> 
> Best,
> Maxi
> 
> 2017-08-30 10:44 GMT+02:00 Codrina Maria Ilie <codrina at geo-spatial.org 
> <mailto:codrina at geo-spatial.org>>:
> 
>     Hi Anna,
> 
>     You were perfectly right!
>     I changed from null to 0 and it works.
> 
>     Thank you!
>     Codrina
> 
> 
>     On 29/08/2017 00:05, Anna Petrášová wrote:
> 
>         Hi,
> 
>         On Mon, Aug 28, 2017 at 4:53 PM, Codrina Maria Ilie
>         <codrina at geo-spatial.org <mailto:codrina at geo-spatial.org>> wrote:
> 
>             Hi all,
> 
>             I have been trying to run GRASS7 module r.damflood using the
>             following
>             datasets:
>             [https://drive.google.com/drive/folders/0B2Yx_1shUSx3VThlY3BrV1JyWjg?usp=sharing
>             <https://drive.google.com/drive/folders/0B2Yx_1shUSx3VThlY3BrV1JyWjg?usp=sharing>],
>             with this command:
> 
>             r.damflood --overwrite --verbose elev=dem lake=water_depth
>             dambreak=dam_breach manning=manning tstop=3 deltat=1
>             h=b_depth vel=b_vel
>             hmax=b_mwd vmax=b_mwv imax=b_mi wavefront=b_wf
> 
>             To make it easier, the lake has a constant depth.
> 
>             As far as I can understand, the rasters should be correctly
>             built, but the
>             outputs of r.damflood are obviously wrong (rasters with
>             nodata or one
>             category = 0; at second 1 the depth raster is half the
>             lake's size (s00.png,
>             s01.png)and at second 2, it completely disappears).
> 
>             During the processing, I get the Courant-Friedrich-Lewy
>             stability condition
>             warning message.
> 
>             If I try to chose as computational method for initial
>             velocity estimation,
>             uniform drop in of lake or small dam breach, I get the
>             following error
>             "ERROR: Don't find the dambreak - Please select a correct
>             map or adjust the
>             computational region".
> 
>             Does anyone have an idea what could be the reason of these
>             results, please?
>             Any suggestion would be highly appreciated.
> 
> 
>         Your water_depth should have 0 instead of nulls, I think that's the
>         problem in your case. The stability condition warning probably means
>         you should reduce your time step, on the other hand, it will
>         take more
>         time.
> 
>         Anna
> 
> 
>             Cheers,
>             Codrina
> 
> 
> 
> 
> 
> 
>             _______________________________________________
>             grass-user mailing list
>             grass-user at lists.osgeo.org <mailto:grass-user at lists.osgeo.org>
>             https://lists.osgeo.org/mailman/listinfo/grass-user
>             <https://lists.osgeo.org/mailman/listinfo/grass-user>
> 
>     _______________________________________________
>     grass-user mailing list
>     grass-user at lists.osgeo.org <mailto:grass-user at lists.osgeo.org>
>     https://lists.osgeo.org/mailman/listinfo/grass-user
>     <https://lists.osgeo.org/mailman/listinfo/grass-user>
> 
> 
> 
> 
> -- 
> *Massimiliano Cannata*
> 
> Professore SUPSI in ingegneria Geomatica
> 
> Responsabile settore Geomatica
> 
> 
> Istituto scienze della Terra
> 
> Dipartimento ambiente costruzione e design
> 
> Scuola universitaria professionale della Svizzera italiana
> 
> Campus Trevano, CH - 6952 Canobbio
> 
> Tel. +41 (0)58 666 62 14____
> 
> Fax +41 (0)58 666 62 09____
> 
> massimiliano.cannata at supsi.ch <mailto:massimiliano.cannata at supsi.ch>
> 
> _www.supsi.ch/ist <http://www.supsi.ch/ist>_
>