69 |
df(i,j) = - _rA(i,j,bi,bj) |
df(i,j) = - _rA(i,j,bi,bj) |
70 |
& *( KPPdiffKzS(i,j,k,bi,bj) + tmpFac*Kwz(i,j,k,bi,bj) ) |
& *( KPPdiffKzS(i,j,k,bi,bj) + tmpFac*Kwz(i,j,k,bi,bj) ) |
71 |
& *KPPghat(i,j,km1,bi,bj) |
& *KPPghat(i,j,km1,bi,bj) |
72 |
& *surfaceForcingS(i,j,bi,bj) |
& *( surfaceForcingS(i,j,bi,bj) |
73 |
|
#ifdef ALLOW_SALT_PLUME |
74 |
|
& + saltPlumeFlux(i,j,bi,bj)*recip_rhoConst |
75 |
|
& *(1.-KPPplumefrac(i,j,bi,bj)) |
76 |
|
#endif /* ALLOW_SALT_PLUME */ |
77 |
|
& ) |
78 |
ENDDO |
ENDDO |
79 |
ENDDO |
ENDDO |
80 |
|
|
85 |
df(i,j) = - _rA(i,j,bi,bj) |
df(i,j) = - _rA(i,j,bi,bj) |
86 |
& *KPPdiffKzS(i,j,k,bi,bj) |
& *KPPdiffKzS(i,j,k,bi,bj) |
87 |
& *KPPghat(i,j,km1,bi,bj) |
& *KPPghat(i,j,km1,bi,bj) |
88 |
& *surfaceForcingS(i,j,bi,bj) |
& *( surfaceForcingS(i,j,bi,bj) |
89 |
|
#ifdef ALLOW_SALT_PLUME |
90 |
|
& + saltPlumeFlux(i,j,bi,bj)*recip_rhoConst |
91 |
|
& *(1.-KPPplumefrac(i,j,bi,bj)) |
92 |
|
#endif /* ALLOW_SALT_PLUME */ |
93 |
|
& ) |
94 |
ENDDO |
ENDDO |
95 |
ENDDO |
ENDDO |
96 |
|
|