60 |
C Here, we'll copy the exact forcing used for salinity to allow us |
C Here, we'll copy the exact forcing used for salinity to allow us |
61 |
C to test everything is working. |
C to test everything is working. |
62 |
gPtracer(i,j,k,bi,bj)=gPtracer(i,j,k,bi,bj) |
gPtracer(i,j,k,bi,bj)=gPtracer(i,j,k,bi,bj) |
63 |
& + EmPmR(i,j,bi,bj)*convertFW2Salt |
& +maskC(i,j,k,bi,bj)*surfaceTendencyS(i,j,bi,bj) |
64 |
& *recip_drF(1)*recip_hFacC(i,j,1,bi,bj) |
|
65 |
& - lambdaSaltClimRelax* |
cph & + EmPmR(i,j,bi,bj)*convertFW2Salt |
66 |
& (pTracer(i,j,k,bi,bj,iTracer)-SSS(i,j,bi,bj)) |
cph & *recip_drF(1)*recip_hFacC(i,j,1,bi,bj) |
67 |
|
cph & - lambdaSaltClimRelax* |
68 |
|
cph & (pTracer(i,j,k,bi,bj,iTracer)-SSS(i,j,bi,bj)) |
69 |
|
|
70 |
ENDDO |
ENDDO |
71 |
ENDDO |
ENDDO |