273 |
I myThid ) |
I myThid ) |
274 |
ENDIF |
ENDIF |
275 |
C-- Calculate buoyancy |
C-- Calculate buoyancy |
276 |
CALL CALC_BUOY( |
CALL CALC_BUOYANCY( |
277 |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1, |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1, |
278 |
O buoyKm1, |
O buoyKm1, |
279 |
I myThid ) |
I myThid ) |
313 |
I myThid ) |
I myThid ) |
314 |
ENDIF |
ENDIF |
315 |
C-- Calculate buoyancy |
C-- Calculate buoyancy |
316 |
CALL CALC_BUOY( |
CALL CALC_BUOYANCY( |
317 |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoK, |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoK, |
318 |
O buoyK, |
O buoyK, |
319 |
I myThid ) |
I myThid ) |
355 |
I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown, |
I bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown, |
356 |
O xA,yA,uTrans,vTrans,rTrans,rVel,maskC,maskUp, |
O xA,yA,uTrans,vTrans,rTrans,rVel,maskC,maskUp, |
357 |
I myThid) |
I myThid) |
358 |
|
CcnhDebugStarts |
359 |
|
C IF ( K .EQ. 1 ) THEN |
360 |
|
C CALL PLOT_FIELD_XYRL( rVel(1,1,1), 'K=1 Current rVel.1 ' , myIter, myThid ) |
361 |
|
C CALL PLOT_FIELD_XYRL( rVel(1,1,2), 'K=1 Current rVel.2 ' , myIter, myThid ) |
362 |
|
C ENDIF |
363 |
|
CcnhDebugEnds |
364 |
C-- Calculate the total vertical diffusivity |
C-- Calculate the total vertical diffusivity |
365 |
CALL CALC_DIFFUSIVITY( |
CALL CALC_DIFFUSIVITY( |
366 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |