open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

Trying to run a simulation in v7p2r0, which has already run successfully using previous versions, I got an error while the code was trying to read the FORMATTED DATA FILE 1 (invoked in the modified bord3d subroutine). After some debugging I found that:
Each of the logic unit variables (T3DFO1, T3DHYD, etc.) is given a value in subroutine lecdon_telemac.f. Then in bief_open_files.f, the unit numbers of the data files, defined also by variables T3D_FILES%LU, change due to the implementation of the function get_free_id.f. The T3Dxxx variables are never updated (to get the same values as T3D_FILES%LU) and this leads to an error in case a data file needs to be read (e.g. when a boundary conditions file has to be read in bord3d.f).
In my case, I solved it updating T3DFO1 variable in bord3d.f.

Cheers,
Makis

Hello,

The principle is that the values like T3DFO1 are never changed and that logical units T3D_FILES%ADR(T3DFO1)%LU may change depending on coupling, so T3DFO1 must never be updated, but to read the file we must always use T3D_FILES%ADR(T3DFO1)%LU, not T3DFO1. In view of this could you tell us if you still see a mistake? You should never open a file with logical unit T3DF01 but with T3D_FILES%ADR(T3DFO1)%LU.

With best regards,

Jean-Michel Hervouet

Hello Jean-Michel,

Actually the update of T3DFO1 that I mentioned in my first post is the same with what you explain here. It works fine!

The modified bord3d.f which I use, had been developed a long time before I started using it. It worked fine with the previous Telemac3d versions so I thought using T3DFO1 variable to open the data file wouldn't cause any problem.

Thanks,
Makis