202 |
ENDDO |
ENDDO |
203 |
ENDDO |
ENDDO |
204 |
|
|
205 |
|
|
206 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
207 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
208 |
|
|
233 |
jMin = 1-OLy+1 |
jMin = 1-OLy+1 |
234 |
jMax = sNy+OLy |
jMax = sNy+OLy |
235 |
|
|
236 |
|
|
237 |
K = 1 |
K = 1 |
238 |
BOTTOM_LAYER = K .EQ. Nr |
BOTTOM_LAYER = K .EQ. Nr |
239 |
|
|
275 |
I myThid ) |
I myThid ) |
276 |
ENDIF |
ENDIF |
277 |
C-- Calculate buoyancy |
C-- Calculate buoyancy |
278 |
CALL CALC_BUOY( |
CALL CALC_BUOYANCY( |
279 |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1, |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1, |
280 |
O buoyKm1, |
O buoyKm1, |
281 |
I myThid ) |
I myThid ) |
315 |
I myThid ) |
I myThid ) |
316 |
ENDIF |
ENDIF |
317 |
C-- Calculate buoyancy |
C-- Calculate buoyancy |
318 |
CALL CALC_BUOY( |
CALL CALC_BUOYANCY( |
319 |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoK, |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoK, |
320 |
O buoyK, |
O buoyK, |
321 |
I myThid ) |
I myThid ) |
357 |
I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown, |
I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown, |
358 |
O xA,yA,uTrans,vTrans,rTrans,rVel,maskC,maskUp, |
O xA,yA,uTrans,vTrans,rTrans,rVel,maskC,maskUp, |
359 |
I myThid) |
I myThid) |
360 |
|
CcnhDebugStarts |
361 |
|
C IF ( K .EQ. 1 ) THEN |
362 |
|
C CALL PLOT_FIELD_XYRL( rVel(1,1,1), 'K=1 Current rVel.1 ' , myIter, myThid ) |
363 |
|
C CALL PLOT_FIELD_XYRL( rVel(1,1,2), 'K=1 Current rVel.2 ' , myIter, myThid ) |
364 |
|
C ENDIF |
365 |
|
CcnhDebugEnds |
366 |
C-- Calculate the total vertical diffusivity |
C-- Calculate the total vertical diffusivity |
367 |
CALL CALC_DIFFUSIVITY( |
CALL CALC_DIFFUSIVITY( |
368 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |