open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello forum!

I am currently trying to model DO/BOD processes and after hustling with manual edits in SOURCE_TRAC I decided to give WAQTEL a shot.

However, the simulations won't start, no error messages are shown. Attached are my .cas files.

Any ideas?

/Schmirre

Hello

Can you send a copture of the listing you obtain?
Please add the keyword DEBUGGER=1

with kind regards

Riadh

Hi again,

Thanks for a fast answer. I solved the problem, I had some problems with my gfortran version after compiling the lastest release version!

Now I've had some time to look into the code and started test simulations of the O2-processes. I have a couple of questions, and I hope someone here could help me to clarify?

1. I did some simulation runs without any oxygen-demandning sources (no organic load, no NH4) just to make sure that the benthic oxygen demand and the reareation are corrently represented. In my model, the oxygen increases at the bottom and decreases at the surface. Why is this reversed? Can it have anything to do with how i specify my vertical layers? Files are attached.

2. In several fortran subroutines it is possible to account for the actual temperature in the model. There is a variable called:

"YATEMP" - if yes temperature is available.
"IND_T" - Index of temperature in tracer table

And this IND_T variable are called in several subroutines (for example calcs3d_o2.f, reaer.f, satur_o2.f). But how do I define temperature? When I've done it in the main .cas file - the tracers order is overridden by the subroutines. Also here my files are attached.

Thanks in advance!

/Schmirre

And here is an update:

When i run the simulations without any sources, WAQTEL automatically generates three tracers: DISSOLVED O2, ORGANIC LOAD, NH4 LOAD.

When I add a source, the simulation stops because it says I haven't defined the initial tracer values for all tracers, which I have..

I also tried to use in total 4 tracers, where i define the first one as 'TEMPERATURE', this doesn't change anything though..

Thanks in advance for your ideas!

Hi Schmirre,

I got the same problem which has been bothering me for a while, it seems that there's a conflict between prescribing Tracer TEMPERATURE for T3D and other Tracers for WAQTEL in the steering file at the same time, especially when taking density flow into account. I'm wondering if you've solved this problem, hope to hear from you soon!

With best regards,
Yunhao

Hi Yunhao!

I am still stuck with the same problem!

Do you have any experience in utilizing a varying temperature for calculations of temperature dependent constants in these subroutines?

Best regards, Schmirre

Hi Schmirre,

It seems that both of us get stuck with this problem , and I don't have much experience about what you said but there's a temperature difference in my domain aiming to study the density flow. Hope the problem can be solved ASAP.

Cheers,
Yunhao

Hi Schmirre,

I solved my problem, as Waqtel increases automatically the number of tracers when coupled with T3D, there's no need to add the number of tracers inside Waqtel (PO4, NO3, etc.) to 'number of tracers', if you use temperature, just set 'number of tracers=1'. For initial tracer values you need to give all of them a value, I just add 'zero' to the unused one.

Besides, there is a bug in the error message in bord3d.f in line 563, it should be 'NUMLIQ%I(K)*NTRAC' instread of 'NUMLIQ%I(K)'.

Hope this helps!
Yunhao