open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

I am trying to run a simple SISYPHE simulation coupled with TELEMAC2D.
But the case has failed to run.



The geometry and the boundary files are the same as mentioned in my previous post,
www.opentelemac.org/media/kunena/attachments/2535/testcase.zip

The tides at the open boundaries nodes are read from a file by modifying the bord.f subroutine. The simulation also has one discharge boundary.

Below are the modified, case file and the fortran file.







Do I need to change any subroutine, or is there a problem with the case file.

With Kind regards,
Gourish

Hi Gourish,
I'm not able (at this moment) to take a look to your files, but what about your "hydrodynamics-only" simulation? Did you try to run a steady or quasi-steady hydrodynamics solution that can be used as initial condition of your coupled problem ?

I hope it helps,

Pablo

Hi Pablo,

I have not yet tried running a steady or a quasi-steady hydrodynamics solution,
I shall try running it. Is there any example in the validation case, where a steady or a quasi-steady hydrodynamics solution is used as an initial condition for coupled problem.

I am attaching the error snapshot of the case.



Thanks and Best Regards,
Gourish

Hi,

I have changed the simulation by giving tide input at the open boundaries using a liquid file. Now the tide at open boundaries is not given from the formatted file. In this condition the simulation works properly when coupled with SISYPHE.

In order to check as where the problem may possible be, I included the following statement in the BORD subroutine
.

The outcome of the inquire command is that "Unit 26 is not open "
and program crashes once it tries to read the line from file unit number 26,
on the read statement


On using the inquire statement to get the file name i.e inquire (26,NAME=name), gives the file name as a binary output. Same problem is found even when I use FORMATTED DATA FILE 2, for imposing my tidal data.

Does this indicates that when TELEMAC is coupled with SISYPHE, the unit no. 26 file is allotted to some other file. Is there any solution for this problem.

I have attached the complete simulation file.



Regards,
Gourish

Hello,

Sure, in coupling there is a conflict between logical unit numbers, so they are distributed differently and you cannot rely on hardcoded values such as 26.

Basically (I take the example of FORMATTED DATA FILE 1 in Telemac-2D. In the dictionnary you can see the line SUBMIT=... corresponding to this file, and you see that its coded name in the temporary folder will be T2DFO1.

To access this file within a subroutine do the following :

Add at the beginning the command :

USE DECLARATIONS_TELEMAC2D, ONLY : T2D_FILES,T2DFO1

Declare :

INTEGER NFO1

Then the logical unit will be :

NFO1=T2D_FILES(T2DFO1)%LU

You can check that the integer T2DFO1 is 26, that without coupling NFO1 will also be 26, and that in coupling the value of NFO1 is changed.

T2D_FILES is an array of a structure called BIEF_FILE, which contains a number of data on files (names, logical unit, etc, see BIEF_DEF in library bief.

With best regards,

Jean-Michel Hervouet