open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello

I want to simulate the "morphodynamics" in a 'city model with river and buildings' with Sisyphe.
But i will to define only the Bed of the river as erodable.

How can I do that?

Best regards
Rima

Hi
There is many options but I think one of the easiest is to define the level of the non erodible bed in BK as a new variable in the geometry file.
Then you modify the fortran (noerod.f) where this value is define by default to (bottom-100) to take in account the value you create in the geometry file

Dear Rima,

Expanding a little the previous answer, the one that I've tried was to create a new variable that states the erodible depth.

if you plan to follow this way:

1. Create a "NOER" variable in your geometry file
2. Map to each node the desired erodible depth (in your case, 0 to the urban area and a default value like 100 to anything else)
3. Place the keyword "FORTRAN FILE" in your steering case file and point the attached fortran file.

Regards,

José Diaz.

Thanks for both of you, Christophe and José

my model is so long. Therefore, the definition of the level of the non erodible bed is not representative of my situation.

I would try to follow the suggestion of José. That suits my situation.


Best wishes
Rima

Hi,

I addition to what José has suggested, I can mention that the procedure of reading a user variable stored in the geometry file can now be achieved by using the following keywords in your Sisyphe steering file (from version 7.1):

After doing that, you can access these private variables with the PRIVE structure, for example with PRIVE%ADR(1)%P%R for the private variable no 1.

If you do like above, defining your depth to erodible bed in a private variable called NOER, all you need to do in noerod.f is the following:

Add the following line just below USE BIEF:
Replace the default CALL OV that sets the erodible bed 100 m below the river bottom
By the following:
Pay attention to the X=Y-Z function, here it works when NOER (that is to say PRIVE%ADR(1)%P%R) is the thickness of the erodible bed, not the elevation of the non-erodible bed!

Hope it helps!

Best regards
PL

Hello
I need your support again.

1- I have already tried the suggestion of "José". In this case, I get an error 'Interface_hermes not found'.
Please see the attachment 'ERROR_interface_hermes'

2- the second suggestion of 'pilou1253' has not worked. Error 'Names of private variables not found'. Can it be, that there are three varaibles in my serafin file 'Zwoenitz'?
Please see the attachment 'Names of private variables'

Thanks for your support
Rima

Hi

You could find all the dictionnary in the source directory so you could have a look in it and try to understand the error messages.

Nonetheless, as you're using a previous version of telemac, this means that some advise are not necessarily relevant. In both case, Jose and Pilou gives you indication which will work only with the last version.
So you have to update telemac to use v7p1...

Regards

Hi Christophe

i have tried to update telemac. But it seems that v7p0 is the newest version.

where can I find the version v7p1?

Regards

Hi
I never used the installer so I don't know how it scan the versions...
The latest version is available on the svn.
You could brownse the repository to see the exact location of the latest version and then run a svn with this reference

Regards

It worked.
Thanks a lot