Dear ROM people,
I tried to calculate a passive tracer with define KPP option.
The tracer was calculated to be not 0, although I. C. and all B. C. of the tracer were set to 0.d0.
I read lines form 763 to 777 of pre_step3d.F as follows;
DO itrc=1,NT(ng)
ltrc=MIN(NAT,itrc)
DO k=1,N(ng)-1
DO i=Istr,Iend
cff=1.0_r8/(z_r(i,j,k+1)-z_r(i,j,k))
FC(i,k)=cff3*cff*Akt(i,j,k,ltrc)* &
& (t(i,j,k+1,nstp,itrc)- &
& t(i,j,k ,nstp,itrc))
# ifdef LMD_NONLOCAL
!
! Add in the nonlocal transport flux for unstable (convective)
! forcing conditions into matrix FC when using the Large et al.
! KPP scheme.
!
FC(i,k)=FC(i,k)-dt(ng)*Akt(i,j,k,ltrc)*ghats(i,j,k,ltrc)
# endif
END DO
END DO
.
If #define LMD_NONLOCAL, FC of the tracer is calculated to add the "ghats", related to temp. or sal. surface flux, calculated with lmd_skpp.F.
So, I added "IF(itrc.le.NAT) then" and "END IF" to before and after line-777 of pre_step3d.F and calculated the passive tracer. The calculated tracer became 0.
Is it O. K. ? I think it is not wrong.
Thanks a lot,
Tsubono