Hello everybody,
I am trying to include a new variable as a coefficient for tracer degradation, K=2.303/T90*24 by using the private array that varies with concentrations of salinity and temperature in preres_telemac2d.f:
IF((LEO.AND.SORLEO(23)).OR.(IMP.AND.SORIMP(23))) THEN
DO N=1,NPOIN
PRIVE%ADR(1)%P%R(N) = (0.8+0.02*TN%ADR(1)%P%R(N))*1.07**
& (TN%ADR(2)%P%R(N)-20)
ENDDO
ENDIF
With K = PRIVE%ADR(1)%P%R(N), salinity = TN%ADR(1)%P%R(N), and temperature = TN%ADR(2)%P%R(N), non-conservative tracer = TN%ADR(3)%P%R(N).
In difsou.f, looks like to multiply directly the calculated decay rate, K with the non-conservative tracer = TN%ADR(3)%P%R(N) as follow, is not producing the decay result.
ELSEIF(LOITRAC(ITRAC).EQ.2) THEN
YASMI(ITRAC)=.TRUE.
CALL OS('X=CY ',X=TIMP%ADR(3)%P,
& Y=TN%ADR(3)%P,
& C=-PRIVE%ADR(1)%P/24.D0/3600.D0)
Can anybody do suggestions on this matter, either to replace the private array PRIVE%ADR(1)%P%R(N) with a new declared variable?
Thank you, Amyrhul
