open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi

I'm running Telemac2D v7p3r0 and have discovered that my computation continued is not retaining the tracer conditions even though they are present in the previous computation file.

I have read previous threads on this problem for older versions (i.e. v7p1) where the problem was the tracer variable number was incorrect due to new variables being added before the tracers.

In telemac2d.f it gives the tracer variables as:

IF(NTRAC.GT.0) THEN
DO ITRAC=1,NTRAC
! SEE NOMVAR_TELEMAC2D
ALIRE(34+ITRAC) = 1
ENDDO
ENDIF


Which matches nomvar_telemac2d.f:

ILAST = 34
INEXT = ILAST+1

IF(NTRAC.GT.0) THEN
DO I=1,NTRAC
TEXTE(ILAST+I) = NAMETRAC(I)
TEXTPR(ILAST+I) = NAMETRAC(I)
WRITE(CHAR2,'(I2)') I
MNEMO(ILAST+I) = 'T'//ADJUSTL(CHAR2)//' '
ENDDO

Are there any other additional variables that have been added that have been missed?

I also notice though in rescue.f that the tracer reference in Trouve is effectively 36 for one tracer although ALIRE is 35, I don't know if ALIRE and TROUVE array sizes are the same or not?:

IF(NTRAC.GT.0) THEN
DO ITRAC=1,NTRAC
IF(TROUVE(34+1+ITRAC).EQ.0) THEN
IF(LNG.EQ.1) WRITE(LU,900)
IF(LNG.EQ.2) WRITE(LU,901)
900 FORMAT(1X,'RESCUE : CALCUL PRECEDENT SANS TRACEUR,',
& /,1X,' ON PREND TRAC0')
901 FORMAT(1X,'RESCUE : PREVIOUS CALCULATION WITHOUT TRACER',
& /,1X,' WE FIX IT TO TRAC0')
CALL OS( 'X=C ' , X=T%ADR(ITRAC)%P,C=TRAC0(ITRAC))
ENDIF
ENDDO
ENDIF

Kind regards

David

I have found the solution.

There is an error in the rescue.f routine.

If the line:

IF(TROUVE(34+1+ITRAC).EQ.0) THEN

is replaced with:

IF(TROUVE(34+ITRAC).EQ.0) THEN

then the TROUVE array matches the ALIRE array and the tracer values in the previous computation file are carried forward with COMPUTATION CONTINUED.

Regards

David

Hello David

I fixed the probem tentatively in the trunk version. The management of tracer index is being chenged deeply with the introduction of Water Quality and Ice new processes. Hopefully we will have much more efficient way to manage them.

with my kind regards

Riadh ATA