open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Costas,

You are right, I attached the wrong file. Here is the correct one:





Sorry for that.

Antoine

Dear Antoine,

The version of streamline you modified is from v6p3. Can these be directly applied to version v7p0r0? Or if it is possible, can you make the modifications in the v7p0 version so that I can compare directly?

Regarding the interpolation in derive.f, I notice that this is done only once for the first algae drop, line 286:
Is this limit related to the issues you discussed earlier? Because this certainly stops the algae transport module for dealing with algae drops at multiple times.

Regards,
Costas

Hello Antoine,

I modified line 286 in derive.f so that the interpolation condition coincides with the times I release algae in my domain and it seems that the problem, no longer occurs. The computation with ALGAE TRANSPORT seems to carry on like with plain floats. I will check more over the weekend.

Can you confirm this as the culprit of the error?

Regards,
Costas

Hello Costas,

I think you should be able to apply the modification to streamline directly to the v7p0r0, as I don't think it has evolved greatly.

Modifying line 286 in derive.f will work numerically, but this means that you are reinterpolating the values of variables ending with "_0" for all the particles. Therefore you will introduce a small error as you are erasing the historical values kept for those values (the "*_0" variables are linked to the previous time step). This was why it was only done for the first algae drop.

Furthermore you will notice in the subroutine DISP_ALGAE of algae_transp.f that ALGAE_START is used as well as there are other historical terms that are calculated. Multiple releases will not have these terms properly taken into account. The way to go around that is to add a variable to algae module so that each particle has its own ALGAE_START variable, but it is quite a lot of work to do this. In the report describing the algae transport module it is explained how to add a variable:

www.opentelemac.org/index.php/manuals/su...looms-and-oil-spills

All the best,

Antoine

Hello Antoine,

Thank you for your explanations, it is clearer now. I believe this limitation should be mentioned within flot.f, because ordinary float transport works with multiple drops but algae transport doesn't, and first time users (like me) could get confused.

I find that multiple drop (in time) functionality would greatly increase the applicability of the algae transport module, however it seems beyond my programming skills, at least for now.

Regarding your modification in streamline.f, how necessary is it? The problem with algae transport module was occurring at the time of the second drop, which was due to the discussed limitations in derive.f.

Regards,
Costas

Hello Costas,

If you only release particles once, then the modification is unnecessary.

Furthermore, I agree with the fact that multiple release do not work with the algae module should be mentioned in flot.f and the user manual. I will try and get it in the next version.

If I have time I might do the necessary developpements for the algae transport module.

Kind regards,

Antoine

Hi Dear,

Is possible user your code for modelling of transport of wood in river.


Regards,