open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello Santiago,

Sorry, you are right, the source term of the MICROPOL module for the bed sediment (or Second Member of the advection-diffusion equation of this tracer) was not taken into account in my suggested fix for CVDFTR subroutine.

Please try the enclosed one (the 2 other subroutines still has to be used) and tell me if you are OK with the results.

Best regards,

Chi-Tuan

Hello Chi-Tuan,

I tried the new proposal for CVDFTR subroutine and it seems to be working just fine.

If I happen to run into any more issues I will reply to this (or start another) topic.

Best regards,

Santiago

Good news Santiago,

Do not hesitate if you see issues/bugs, it helps us.

FYI, all these fixes have been pushed in future tagged v8p2r0.

Best regards,

Chi-Tuan

Hello Chi-Tuan,

I have been thinking in the modification on the EROSION_FX subroutine.

Maybe using max is not the best solution, because then you are multiplying the erosion flux by the bed sediment concentration (in case it is > 0) but in the Krone and Partheniades equations the erosion flux is independent of it.

Best regards,

Santiago

Hello,

thanks for supplying all this information on waqtel/micropol.

For compatibility reasons, I want to use micropol on 7p3r1.
I updated the files according to your suggestions.
The gouttedo-case is runable. But the propagation is restricted to the x-direction. I made a test with other waqtel modules. Here the propagation is symmetrical.

Can some give a hint?

Kind regards,
Karl

Hello Karl,

No assistance for release before the latest one (v8p3r1) or the main.
There were some bugs fixed between v7p3r1 and v8p3r1 (4 years!).
Please upgrade your release and see if your issue still occurs.

Hope this heps,

Chi-Tuan

I have some question about the units concerning the parameters and variables of MICROPOL, specifically in calc2d_micropol.f the explicit calculation step of pollutant density [C]. As it goes, X=CY => T1 = Css*K-1. Then new T1 is calculated as T1 = T1 - k1* C * SS and finally the explicit part TEXP(IND_C) = TEXP(IND_C) + T1 is calculated. From the technical manual I find that k1 has [1/g/s] and k-1 has [1/s] unit, while C [kg/m3], SS [kg/m3] and Css [kg/kg]. My question is related to units at each of this step. What will be the unit for T1 each step of calculation ?

Regards,
Nithin

Hello Nithin,

For each step, the unit of T1 is [kg/m^3/s].
Be careful, for the 2nd step, what is done is T1 = T1 - K-1 * K_d * C * SS (you have forgotten K_d constant [m3/kg].

Hope this helps,

Chi-Tuan

Thanks Chi-Tuan. I mentioned K1 in T1 = T1 - k1* C * SS in my previous question since K-1*Kd = K1.
The units mentioned in the technical manual might be incorrect since it says k1[1/kg s], k-1[1/s] and kd [1/kg]. But what I understand from you is k1 [m3/kgs], k-1[1/s] and kd[m3/kg], am I correct ?
Thanks
Nithin

Hello Nithin,

You are almost all right except for the units in the Theory guide where it should be read l (L) for litter rather than 1 (one). I recognize typo is a little bit confusing but it would be written g^-1 rather than 1/g for the units.

Chi-Tuan