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 |
|
|
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 |
CcnhDebugStarts |
361 |
IF ( K .EQ. 1 ) THEN |
C IF ( K .EQ. 1 ) THEN |
362 |
CALL PLOT_FIELD_XYRL( rVel(1,1,1), 'K=1 Current rVel.1 ' , myIter, myThid ) |
C CALL PLOT_FIELD_XYRL( rVel(1,1,1), 'K=1 Current rVel.1 ' , myIter, myThid ) |
363 |
CALL PLOT_FIELD_XYRL( rVel(1,1,2), 'K=1 Current rVel.2 ' , myIter, myThid ) |
C CALL PLOT_FIELD_XYRL( rVel(1,1,2), 'K=1 Current rVel.2 ' , myIter, myThid ) |
364 |
ENDIF |
C ENDIF |
365 |
CcnhDebugEnds |
CcnhDebugEnds |
366 |
C-- Calculate the total vertical diffusivity |
C-- Calculate the total vertical diffusivity |
367 |
CALL CALC_DIFFUSIVITY( |
CALL CALC_DIFFUSIVITY( |