open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello all!

In the framework of my Msc. thesis I have implemented a new flocculation model in Télémac-3D. I am currently trying to validate the model with data for the Belgian coast.

I have fitted a cosine function for the water depth during a tidal cycle with the data (along with a function for bottom shear) and have set up a test case. I don't have experience with TOMAWAC, so my idea is to overwrite values of the water depth in the fortran file.

Unfortunately this gives me some issues. The values of the water depth seem to vary according to up- and downstream boundary conditions (imposed flow 4 5 5 or free 4 4 4) and do not follow my function at all:

I would like to know what would be the optimal way to incorporate this varying water depth in a correct fashion.

Thanks in advance!

Shandro

Hello,

It is strange that you apply your formula only to the point NPOIN3 which is the last of the 3D points:

HN%R(NPOIN3)=...

HN%R is a 2-dimensional array so you should give values to indices that are below or equal to NPOIN2, corresponding to the boundary where you want to apply this formula.

With best regards,

Jean-Michel Hervouet