open TELEMAC-MASCARET The mathematically superior suite of solvers

I am trying to adapt the micropol module of waqtel. I am fairly new to working with fortran files, therefore any guidance is welcome. In micropol there is one suspended sediment and one deposited sediment class. However for my application I want to add two more for each of these sediment classes. I am wondering what would be the most clean and efficient way to do this. Also which fortran files needs to be modified? Thank you. Looking forward to your answers.
Regards
Nithin

Hello,

Zied Amama did a great job this year for validating, debugging and improving the MICROPOL module, in particular for 3D. Do you work with TELEMAC-2D or TELEMAC-3D?

In particular, he has introduced a 2 step reversible model, maybe you can have a look at what he did (source code + manuals, in particular WAQTEL theory guide) and adapt his development of what you need if it is not enough.

In that case, you have to clone the main version (it will be available for next release v8p4, but is already available in the main) and compiling the documentation with command:
doc_telemac.py -m waqtel --theory

The main implementation is done in CALCS2D_MICROPOL (for 2D) or CALCS3D_MICROPOLV and CALCS3D_MICROPOLS (for 3D) subroutines located in $HOMETEL/sources/waqtel folder.

Hope this helps,

Chi-Tuan

thank you Chi-Tuan.
I managed to add the additional sediment classes. I am using telemac2d. The work of Zied is quite interesting and will come in handy for me sometime later.
However I have a question, are the masses of suspended and bed sediments in micropol conserved? I did some checks for my test model. The mass balance is quite off.
regards,
Nithin

Hello Nithin,

If you use dedicated options, mass balance should be conserved. E.g. use a distributive scheme for tracers (default = 1 in the 3rd item of TYPE OF ADVECTION or for SCHEME FOR ADVECTION OF TRACERS, i.e. characteristics which are not conservative), decrease accuracies (e.g. 1.E-8 or less), no use of masking...

Running the test case t2d_waq2d_micropol.cas seems OK to me with the (not so good) choices used to check this mass balance criterion.

As suggested in the rules, you can upload your steering files.

Hope this helps,

Chi-Tuan

Thank you Chi-Tuan. Indeed the mass balance that you showed is quite good. But is it from the test case in telemac folder? When I run the test case from the telemac folder, I don't get same results (I use v8p1r2, does it make a difference ?). The initial and final masses are different and the mass of suspended and deposited sediments are not balanced. is it possible to share the steering file that you used? Thank you
Attached are the case file of test case from telemac folder, its mass balance info and my steering file.

Hello Nithin,

Sorry, I ran the MICROPOL example from the main and forgot that we changed the geometry file to save time for this basin.

Anyway, please upgrade your release to the latest (currently v8p3r1 or main) before more questions because some bugs may have already been solved and as we only assist for this latest release.

If using v8p3 and modified steering file t2d_waq2d_micropol.cas with TYPE OF ADVECTION = 1;5;4 rather than default choice (with characteristics = 1 for tracers at the 3rd column), you can get rather good mass conservation for most of the tracer variables.

Chi-Tuan