open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I assume you mean user_corfon.f rather than user_conf.f and CORFON rather than CORDON?

If you do not modify corfon subroutine (the one calling user_corfon),
"CORFON (TELEMAC2D): NO MODIFICATION OF BOTTOM" will always be written if BOTTOM SMOOTHINGS = 0 (default value). That does not mean what you have implemented in user_corfon subroutine is not applied.

You can check that your modifications in user_corfon are taken into account e.g.:
- add a print : WRITE(LU,*) 'BOTTOM MODIFIED BY GREGORY' in user_corfon,
- check that what you have implemented is reproduced in the RESULTS FILE for the variable B (bottom) after a 1 time step computation.

I suggest you to first implement something and run in serial/sequential mode (not parallel) and when you are OK, produce a source code working for parallel mode.
Moreover, it is always advised (when possible) not to produce source code with number of nodes but with geometry criteria such as operations for abscissae and/or ordinates. See many examples of user_corfon subroutine in the examples and the use of INPOLY subroutine in other FORTRAN FILE (to locate nodes inside a polygon).

Hope this helps,

Chi-Tuan

Hello again,

Perhaps not possible with v8p2, but perhaps possible with the latest release v8p4 (I am not a GAIA specialist). Many improvements and bug fixes have been done since v8p2 for GAIA.

Chi-Tuan

Hi Rebecca and Chi-Tuan. Thank you for your help!

Ok, i will try to install v8p4.

I think my user_corfon.f isn't considered because (as it write in manual user) i'm for COMPUTATION CONTINUED on TRUE. As i study breach ernergy, i need to computation continued on true.

I tested whithout "computation continued", it's ok. With, even if i see "bottom modife by greg", is not ok: the bottom isn't modified.

How can i call [user_corfon.f] even if i have COMPUTATION CONTINUED on TRUE?

still thank

Hello Gregory,

Indeed, if you activate COMPUTATION CONTINUED = YES, even if you add a user_corfon subroutine, with standard TELEMAC-2D, it will not be taken into account.

What you can do is to remove BOTTOM from your PREVIOUS COMPUTATION FILE.
To do so (if not following what is written in the TELEMAC-2D user manual), you can use run_telfile.py Python script with argument alter and option --vars (see the help run_telfile.py -h), e.g.:

run_telfile.py alter --vars 'VELOCITY U,VELOCITY V,FREE SURFACE' input_file output_file

Be careful as the velocity may not fit your new bottom elevation field.

Hope this helps,

Chi-Tuan

Great thank you Chi-Tuan, i try it and come back complete the post.

Hello,

Just for conclure on the subject "breach with gaia":

I try with v8p4r0,( that i succeeded install on linux but not on W10...): i did not notice ameliorations.

Without being a confirmed user, with my modest experience, to consider solid transport (erosion pit, fracture front viscosity,...) across a breach with gaia is not possible or difficult.

if anyone is using this kind of works successfully, i'm interested to see the *.cas files.

Thank you again for your help.

Hello,

For people who want to model breaches taking solid transport into account: if Gaia's module is having trouble, the Sisyphe module seems to work better with the breach module.

bye,
Grégory

Dear Gregory,

can you please give us more details? What does it mean that seems to work better?

Thanks,
Sara

Hello Sara,

It's more easy to simule solid transport across one breach with Sisyphe.

With Gaia (right on image): the breach does not form normally in spite of lot of differents setting.

With Sisyphe, there is no problem (left on image).





++
Grégory

Hello again,

thanks for these images.
Is it possible to have your case with GAIA's steering file, SISYPHE's steering file, as well as all the other input files please?
GAIA should not give worst results (results should be at least comparable and better for some cases) and I would like to understand if it's a problem of 'configuration' or if it comes from GAIA itself.
Can you please remember us with which version of GAIA and SISYPHE you obtained these results please?

Thanks and best regards,
Sara