200 |
jp1 = j + 1 |
jp1 = j + 1 |
201 |
DO i = imin, imax |
DO i = imin, imax |
202 |
ip1 = i+1 |
ip1 = i+1 |
|
|
|
203 |
bo(I,J) = - gravity * |
bo(I,J) = - gravity * |
204 |
& ( TTALPHA(I,J,1) * (surfForcT(i,j,bi,bj)+ |
& ( TTALPHA(I,J,1) * (surfForcT(i,j,bi,bj)+ |
205 |
& surfForcTice(i,j,bi,bj)) + |
& surfForcTice(i,j,bi,bj)) + |
206 |
& SSBETA(I,J,1) * surfForcS(i,j,bi,bj) ) |
& SSBETA(I,J,1) * surfForcS(i,j,bi,bj) ) |
207 |
& / rhoSurf(I,J) |
& / rhoSurf(I,J) |
|
|
|
208 |
bosol(I,J) = gravity * TTALPHA(I,J,1) * Qsw(i,j,bi,bj) * |
bosol(I,J) = gravity * TTALPHA(I,J,1) * Qsw(i,j,bi,bj) * |
209 |
& recip_Cp*recip_rhoConst |
& recip_Cp*recip_rhoConst |
210 |
& / rhoSurf(I,J) |
& / rhoSurf(I,J) |
|
#ifdef ALLOW_SALT_PLUME |
|
|
boplume(I,J) = - gravity * SSBETA(I,J,1) |
|
|
& * saltPlumeFlux(i,j,bi,bj) |
|
|
& * recip_rhoConst / rhoSurf(I,J) |
|
|
#endif /* ALLOW_SALT_PLUME */ |
|
211 |
END DO |
END DO |
212 |
END DO |
END DO |
213 |
|
|
214 |
|
#ifdef ALLOW_SALT_PLUME |
215 |
|
IF ( useSALT_PLUME ) THEN |
216 |
|
DO j = jmin, jmax |
217 |
|
jp1 = j + 1 |
218 |
|
DO i = imin, imax |
219 |
|
ip1 = i+1 |
220 |
|
boplume(I,J) = - gravity * SSBETA(I,J,1) |
221 |
|
& * saltPlumeFlux(i,j,bi,bj) |
222 |
|
& * recip_rhoConst / rhoSurf(I,J) |
223 |
|
END DO |
224 |
|
END DO |
225 |
|
ENDIF |
226 |
|
#endif /* ALLOW_SALT_PLUME */ |
227 |
|
|
228 |
cph( |
cph( |
229 |
CADJ store ustar = comlev1_kpp, key = ikppkey |
CADJ store ustar = comlev1_kpp, key = ikppkey |
230 |
cph) |
cph) |
231 |
|
|
232 |
|
#ifdef ALLOW_DIAGNOSTICS |
233 |
|
IF ( useDiagnostics ) THEN |
234 |
|
CALL DIAGNOSTICS_FILL(bo ,'KPPbo ',0,1,1,bi,bj,myThid) |
235 |
|
CALL DIAGNOSTICS_FILL(bosol ,'KPPbosol',0,1,1,bi,bj,myThid) |
236 |
|
#ifdef ALLOW_SALT_PLUME |
237 |
|
CALL DIAGNOSTICS_FILL(boplume,'KPPboplm',0,1,1,bi,bj,myThid) |
238 |
|
#endif /* ALLOW_SALT_PLUME */ |
239 |
|
ENDIF |
240 |
|
#endif /* ALLOW_DIAGNOSTICS */ |
241 |
|
|
242 |
c------------------------------------------------------------------------ |
c------------------------------------------------------------------------ |
243 |
c velocity shear |
c velocity shear |
244 |
c -------------- |
c -------------- |