open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,all

I'm trying to run tidal simulation with TELEMAC3D, and to provide tidal boundary conditions, I want to use formatted data file with tidal constans for each boundary nodes then use BORD3D.f to assign those values. I thought it would be similar to that in TELEMAC2D since this is the 2D boundary condition, but I got errors about the 'shape of the arry' in compilation of BORD.f, can someone give me a hint...
Thanks in advance!

Merry Christmas O(∩_∩)O~

With best regards.
Robinson

Hello Robinson,

Your problem may come from the declaration of the variable Z. In bord.f, it is declared as a double precision variable inside the subroutine whereas in bord3d.f it is not declared inside but in declarations_telemac3d.f (as an alias): Z is then a pointer on an array of double precision containing the z coordinates of the nodes. Just try to change the name of the variable Z in your subroutine and declare it as a double precision and tell me if it works or not.

Hope this helps,

Chi-Tuan

Hi Chi-Tuan,

Thank you very much for your help.

Yes, just like what you've said, if using another variable instead of Z or not using this intermediate variable at all, the compilation is ok, but when it comes to run the computation, it always stop with this kind of error message:
MORE PRESCRIBED ELEVATIONS ARE REQUIRED
IN THE PARAMETER FILE
AT LEAST 0 MUST BE GIVEN
OTHER POSSIBILITY
STAGE-DISCHARGE CURVES FILE MISSING
so I disable those sentences and the call of PLANTE(1),then it seems work fine(i think).

With best regard!
Robinson

Hi,

I think I have solve my problem, it's about the variable 'Z' that I used which may conflict with the system defined ones. If this intermediate variable is not introduced, and also make 'CALL PLANTE(1)' and LNG .EQ.1/2 disabled, it works fine if there is no other potential risk in it.

with best regard.
Robinson

Hello again,

I think that in your model, you have 2 liquid boundaries or more, do you? If yes, in your steering file, if N is the number of liquid boundaries, you should have N values (at least) for the keyword PRESCRIBED ELEVATION (e.g. PRESCRIBED ELEVATION = 0.;0. if there are two of them). These values are default values if you do not assign values for specific locations.
FYI, the error messages and calls of PLANTE subroutine exist in TELEMAC to prevent from errors that can occur not only in the subroutine where there are present but also in other subroutines. Be careful when you try to delete them.

Regards,

Chi-Tuan