open TELEMAC-MASCARET The mathematically superior suite of solvers

Violeta
Could you post the listing file?
Maybe you should modify your fortran file to make the difference between the two lines WRITE(LU,*). Something like salinity initialization and temperature initialisation.

Regards

Thanks Christophe. I will do it and see if it works. Listing file - do you mean steering file?

Kind Regards!

Violeta

no the output file (-s option)

regards

Hi Christophe,

I have just got a new result file, please see attached. There is not value for tracers even after making changes in Fortran files: giving names of tracers in Write(LU,*) and swaping Salinity and Tracers. However other parameters :velocities seem fine.

hi
In your output file, you could see that condim is well treated as some information exists.
But the tracer initialization is missing.
You should investigate by adding a write before the line if (ntrac.NE.0) where you print a message and also the value of NTRAC

For me it looks like there was a problem with this value...

Hope this helps

A possible problem in your steering file is the names of tracers:
For each tracer, the name is introduced using 32 characters (16 for the name and 16 for the unit) thus for salinity, for e.g. you have to write "SALINITY

---

KG/M3

---

". and so for temperature.

Furthermore, in order to be to be more pragmatic, you can change the block
IF(TA%ADR(ITRAC)%P%NAME(1:4).EQ.'SALI') THEN
DO IPLAN=1,NPLAN
.
.
.
TA%ADR(ITRAC)%P%R(J)=0.D0

by

DO I=1,NPOIN2
DO IPLAN = 1 , NPLAN
J = I + ((IPLAN-1)*NPOIN2)
TA%ADR(1)%P%R(J) = value for temperature
TA%ADR(2)%P%R(J) = value for salinity

I'm not expert of T3D but it seems to me that there is a mandatory order of tracers when dealing with temperature and salinity. If I remember well, temperature must be the first one, then salinity and then the remaining tracers. To be confirmed (Cristophe or Chi-Tuân)

With my best regards

Riadh

Thank you very much Riadh, I will try it. You are right I have Salinity first and then temperature, I will change the way around and try again as per your suggestion.

Kind Regards!

Violeta

Hi Christophe/Riadh,

Thank you very much for your help. It worked. Here is the Fortran script that worked for me. So I can add more layers in similar way. Thanks again.

IF(NTRAC.NE.0) THEN
DO ITRAC=1,NTRAC
CALL OS('X=C ',X=TA%ADR(ITRAC)%P,C=TRAC0(ITRAC))
DO I=1,NPOIN2
DO IPLAN=1,NPLAN
J=I+((IPLAN-1)*NPOIN2)
IF(IPLAN.GT.18) THEN
TA%ADR(1)%P%R(J)=-0.29D0
TA%ADR(2)%P%R(J)=0.D0
ELSEIF(IPLAN.GT.15.AND.IPLAN.LE.28) THEN
TA%ADR(1)%P%R(J)=1.D0
TA%ADR(2)%P%R(J)=24.D0
ELSE
TA%ADR(1)%P%R(J)=4.3D0
TA%ADR(2)%P%R(J)=34.D0
ENDIF
ENDDO
ENDDO
ENDDO
ENDIF


Kind Regards!

Violeta

Looking at steering file again I have a question. Since I have three stratified layers then should I enter three values for salinity and temperature in Prescribed tracers vaules, as I have only one for each? I am not sure. Thanks.

Kind Regards!

Violet

No, you have 2 tracers so only 2 values in the prescribed tracer values (for boundary) and 2 values in the initial tracer values

About the Riadh remark on the order, I didn't check it but in my opinion, it's only important if you use the density law with the value 3. and this should never have an impact on the initialisation of the tracer ...

Regards