open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I want to simulate the 3D flow in square domain(4 open boundaries) by imposing TEMPERATURE & SALINITY(Space & time dependent) conditions along the open boundaries. I wish to know the subroutines that deal with such type of boundary condition.

Thanks
Sathya

Hello,

Once you boundary conditions types are given in the boundary conditions file, the only subroutine to be changed is bord3d.f. You can look in this subroutine what is done on tracers for boundaries, namely from line 396 to 467. The only task done here is to fill the prescribed values at the boundaries, which are:

TABORL%ADR(ITRAC)%P%R(IBORD)

where ITRAC is the tracer number and IBORD the boundary numbering. You can replace all this block or rather put after a new block that gives these values in your case, depending on space and time. The time is variable TIME and the coordinates are MESH3D%X%R(I) and MESH3D%Y%R(I) for a global point I.

With best regards,

Jean-Michel Hervouet

Hello,

1) Would you tell about the VERSION of telemac where the lines 396 to 467 belong to?

2) What will be the BOUNDARY TYPE CODE (for tracers) to be written in the BOUNDARY CONDITION FILE to impose TIME & SPACE DEPENDENT Temperature and salinity boundary condition, ?

2) Will the bord3d subroutine take care if i consider the following as tracers?

NUMBER OF TRACERS : 2
NAMES of TRACERS : 'TEMPERATURE';'SALINITY'

i.e. itrac=1 for TEMPERATURE
and itrac=2 for SALINITY.



Thanks
Sathya

Hello,

1) Version 7.1, but line 400 in 7.0, namely all the block:

IF(NTRAC.GT.0.AND.NTRACER.GT.0) THEN
...
ENDIF

2) no specific type, except that it is prescribed values, so 5 or KENT. It is never looked after bord3d.f if there is a variation in time and/or space or not.

3) OK, temperature and salinity are recognised by their name, but the names must be given in 32 characters, 16 for the name and 16 for the unit.

Regards,

JMH