open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear all,
I want to correct the simulated 3D flow filed by Telemac-3D based on the velocity data measured by ADCP and data assimilation method. I want to do this with TelApy module. The flow velocity data of ADCP are different velocities values at different water depths on a vertical line.

I want to extract the flow velocity simulation values of different plans when the model run to the time where the ADCP data are available. However, I didn’t find any guidance in TelApy User Manual and previous posts.

There are two ways I can do this. My 2D mesh has a total of 14151 nodes and the number of plan is 6.

The first way is that the nodes at the same point of each layer = [4934, 4934+14151, 4934+2*14151, 4934+3*14151, 4934+4*14151, 4934+5*14151], then the simulated flow velocity of each point is obtained with function “myCase.get ('MODEL.VELOCITYU', i = node)”. But, the velocity value obtained in this way does not accord with the actual situation.

The second way is that the node number of the same point in each layer is the same, [4934, 4934, 4934, 4934, 4934, 4934]. I want to get the simulated flow velocity of the same point of each layer with the “For Loop” in each plan. However, I can’t specify a fixed layer.

So, my question is whether there is a function like “myCase.get ('MODEL.VELOCITYU', i = node)” that specifies a particular plan. For example, fisrt the plan 6 is extracted/specified , then the velocity value of the node on plan 6 is extracted. If there is no such function, what other way can it be implemented?

Looking forward to your reply. Thanks!


Yixiang Sun

Hi Yixiang Sun,

I am also interested in this topic. However, I was leaning towards trying to edit the Fortran Api instead of using the Telapy module based on my programming expertise, but it seems that using Telapy is much easier the more I read about it. Have you had any progress with this?

Cheers,

Nayef