open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello dears,

I am trying to couple my own program with Telemac2d. I followed all the steps required in the programming manual and same files that used by other programs such as sisyphe and telemac2d. I named my program as ESTEL3D because I couldnt compile my program under different name eventhough I added the name to the compileTelemac.py file and cmd file too. In the telemac2d, for COUPLING WITH ='ESTEL3D', and when run the telemac2d, in the home_telemac2d, reading my coupling and said coupling=ESTEL3D. However, in the Telemac2d.f when running, doesnt show coupling for my program and when I have if condition, doesnt recognize my program. In lecdon_telemac2d, it reads my program.

Here I attached a log of running telemac2d coupled with my program.

Another program is when I call my program from telemac2d without if condition. Is working, but when I am writing the mesh for the initiol value, for example, the elevation of points such as we have in TELEMAC3D, it gives error in bluekenue.

Best Regards,

AMANJ

Hello dears,

I was able to fix the first part of my previous problems. The problem was I had have duplicated coupling.

Stil I have the second one which is related to writing the mesh into serafin file and then read by bluekenue. I looked at Telemac3D and I did exactly the same and just for test I used z as elevation of the mesh points which we have same in Telemac3d but I am getting this error message.

Hello Amanj

It depends on how your 3d code builds the 3D mesh, is it exactly as done by Telemac-3d ? In this case, just try to view the 3D mesh at its initial conditions.
If it is OK (initial conditions), the problem could come from the computation which can cause layer overlapping.

with my best regards

Riadh ATA

Hello dear Riadh,

I tried to create 3D mesh by taking benifit from the same way as Telemac3D doing using equally space and NPLAN. I have two problems:

first, If look at the pictures I attached, you can see the mesh is refined and overlayed at the toes of the dam (upstream and downstream). the question is how to avoid this, I mean how to correct the mesh and clip the layers to justify for example saying the min, depth of the layer as 0.01m. Second problem is, this is a case of coupling Telemac2D and sisyphe ( sediment erosion), so when the depth decreases, the mesh should update. This updtae either by creating new mesh at each time step which is not practice and layer depth will reduce and at the end all overlay each others, or have to remove the crashed layer,I mean starting from top layer, which is bathy mesh for T2D when this layer is washout totally.

Any solution?

Best Regards,
AMANJ

Hello,

I would like to add some suggestions. What if I mask those elements whose height lessthan the threshold value .Or mask the elements or nodes close to the toes.

Question: do the layers overlapped in the place where the depths are zero? for example, where the surface reach to the bottom, layers depth are zero. but in calcote we still consider as a layer right? how in Telemac3d clipped or masked those elements whose height are zero?

If you look at the pics attached. this is Telemac3D mesh at initial condition. I believe the 3D mesh layers extend from IPOIN=1, NPOIN. It means covers all points regardless of 3D element height.

Kind Regards,

AMANJ

Hello Amanj

When I read your posts, I'm mixed up between Telemac 2D and 3D :)

Yes you can mask elements which causes you troubles. You can do that using subroutine maskob.f. This will be efficient for you case, I think.

Other option, is to use the clipping (CLIPPING=YES), which allows you to remove water in elements when it is less than a threshol given be MINIMUM VALUE OF DEPTH.

good luck

with my kind regards

Riadh