43 |
DO j=jMin,jMax |
DO j=jMin,jMax |
44 |
DO i=iMin,iMax |
DO i=iMin,iMax |
45 |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj) |
46 |
& + foFacMom*fu(i,j,bi,bj)*_maskW(i,j,kLev,bi,bj) |
& +foFacMom*fu(i,j,bi,bj) |
47 |
|
& *horiVertRatio*recip_rhoNil*recip_dRf(kLev) |
48 |
|
& *_maskW(i,j,kLev,bi,bj) |
49 |
ENDDO |
ENDDO |
50 |
ENDDO |
ENDDO |
51 |
ENDIF |
ENDIF |
93 |
DO j=jMin,jMax |
DO j=jMin,jMax |
94 |
DO i=iMin,iMax |
DO i=iMin,iMax |
95 |
gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj) |
gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj) |
96 |
& + foFacMom*fv(i,j,bi,bj)*_maskS(i,j,kLev,bi,bj) |
& +foFacMom*fv(i,j,bi,bj) |
97 |
|
& *horiVertRatio*recip_rhoNil*recip_dRf(kLev) |
98 |
|
& *_maskS(i,j,kLev,bi,bj) |
99 |
ENDDO |
ENDDO |
100 |
ENDDO |
ENDDO |
101 |
ENDIF |
ENDIF |
147 |
gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj) |
gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj) |
148 |
& +maskC(i,j)*( |
& +maskC(i,j)*( |
149 |
& -lambdaThetaClimRelax*(theta(i,j,kLev,bi,bj)-SST(i,j,bi,bj)) |
& -lambdaThetaClimRelax*(theta(i,j,kLev,bi,bj)-SST(i,j,bi,bj)) |
150 |
& -Qnet(i,j,bi,bj) ) |
& -Qnet(i,j,bi,bj)*recip_Cp*recip_rhoNil*recip_dRf(kLev) ) |
151 |
ENDDO |
ENDDO |
152 |
ENDDO |
ENDDO |
153 |
ENDIF |
ENDIF |
199 |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj) |
200 |
& +maskC(i,j)*( |
& +maskC(i,j)*( |
201 |
& -lambdaSaltClimRelax*(salt(i,j,kLev,bi,bj)-SSS(i,j,bi,bj)) |
& -lambdaSaltClimRelax*(salt(i,j,kLev,bi,bj)-SSS(i,j,bi,bj)) |
202 |
& +EmPmR(i,j,bi,bj) ) |
#ifndef USE_NATURAL_BCS |
203 |
|
& +EmPmR(i,j,bi,bj)*recip_dRf(1)*35. |
204 |
|
#endif |
205 |
|
& ) |
206 |
ENDDO |
ENDDO |
207 |
ENDDO |
ENDDO |
208 |
ENDIF |
ENDIF |