open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,
I would like to change the bathymetry of a number of points.
I have modified CORFON this way :

IF(LISFON.GT.0) THEN
!
! !MAS=.TRUE.
! !CALL FILTER(ZF,MAS,T1,T2,AM1,'MATMAS ',
! & 1.D0,T1,T1,T1,T1,T1,T1,MESH,MSK,MASKEL,LISFON)
DO I1=43400,43500
ZF%R(I1)=10000
!WRITE(LU,*) 'I1 : ',I1,'Zf :',ZF%R(I1)
ENDDO
!
ENDIF

I added the key word LISSAGES DU FOND=1 /si différent de 0 : batardeaux

The computation runs, it seems to work but I have no change of bathymetry in the result file.
Could anyone help ?
Thank you,

Pauline

Hello Pauline,

If you use the keyword "BOTTOM SMOOTHINGS = 1" (LISSAGES DE FOND), Telemac will run the routine filter once, which will smooth the bathymetry. I saw in your code that you've commented these lines, which is not recommended.
The corfon subroutine is given to permit changing bathymetry at some (or all) mesh points. To introduce these changes, it is better to start coding just after the the 'CALL FILTER ...'

The question why you don't see any change in the bathymetry, is probably linked to the global/local numbering of nodes. I guess that you are working in parallel mode. If it is the case, you have to proceed by coordinates and not by node numbers. for Telemac, node number 43400 is the local number at the working subdomain. To switch to global number you need to use the table KNOGL (see guide for programming).
A much more easier way to proceed is to run your case 1 time step in serial mode, and to use the result file as a restart.

I hope that it was helpful

With my best regards
Riadh ATA