open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

Yes indeed, the rain was not implemented with finite volumes in 7.2 but changes are being made by Riadh in the trunk.

You can get rain working with FV in 7.2.0 if you do the following (sorry not to have mentionned that earlier):

- use the attached prosou subroutine (the version 7.2 is not optimized for FV and the CALL to the subroutine RUNOFF_SCS_CN is within a IF(BANDEC) test (tidal flats which is not used in FV)). The attached routine is a quick fix. I just made the changes but did not actually try to compile it, check it though but I got something similar working.

- get the subroutine majzz.f from the trunk and add it without modifications to your fortran file (it has been updated in the trunk only)

I have a question for Riadh: I saw that there is a new prosou_fv.in the trunk. I guess it is called if using FV. In this version too there is the test IF(BANDEC)... in the rain section, so this should maybe be changed to be 100% compatible with FV.

Good luck,
PL

Good morning everyone,
I need to implement splash erosion model, E3 found in page 19 of the document found in this link;
www.opentelemac.org/downloads/Papers%20a...2015-proceedings.pdf
Question:
1. Which specific subroutine in Sisyphe should I modify?
2. The model depends on rainfall intensity, thus prosou.f in Telemac has to be engaged. How can it be done?


Koskei

Hi,

I cannot help with the Sisyphe part, sorry.
Otherwise, you need to run a coupled simulation between T2D and Sisyphe so prosou (and the runoff model if you are using it) will be called as if it was a single T2D run. The thing you need to now about coupling is that all your fortran subroutines (modified) should be added in T2D's FORTRAN FILE.

I hope I did not misunderstood your question.

Good luck,
PL

Hello

You can start by browsing the subroutine sisyphe.F and more particularly the subroutine suspension_computation.f in folder sources/sisyphe.
The code is well commented and if you follow carefully the program line per line you should be able to understand how things are implemented.
You should also read sisyphe manuals in order to catch the theory behind.

with my best regards

Riadh ATA

Dear everyone,

Based on recommendations by Riadh, we modified the following lines in suspension_computation.f (attached); line 539, 540, 563, 588, 597. we implemented equation 4 after introducing sediment density to ensure it's dimensions are correct (i.e. m/s). I have tried running the model (trunk version) with it but it doesn't compile. I get the error message shown in error.png (attached).

Please assist me in checking if modification is correct and remedy to the error msg.

Koskei

Hello Koskei

The error message is explicit. It is said that at line 806 one argument (flulimebe) is passed as an integer whether it is expected to be a bief_obj type.
try to verifiy each argument in the call of the new subroutine.

with my best regards

Riadh