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) |
|
CcnhDebugStarts |
|
|
C IF ( K .EQ. 1 ) THEN |
|
|
C CALL PLOT_FIELD_XYRL( rVel(1,1,1), 'K=1 Current rVel.1 ' , myIter, myThid ) |
|
|
C CALL PLOT_FIELD_XYRL( rVel(1,1,2), 'K=1 Current rVel.2 ' , myIter, myThid ) |
|
|
C ENDIF |
|
|
CcnhDebugEnds |
|
360 |
C-- Calculate the total vertical diffusivity |
C-- Calculate the total vertical diffusivity |
361 |
CALL CALC_DIFFUSIVITY( |
CALL CALC_DIFFUSIVITY( |
362 |
I bi,bj,iMin,iMax,jMin,jMax,K, |
I bi,bj,iMin,iMax,jMin,jMax,K, |
450 |
C & maxval(salt(1:sNx,1:sNy,:,:,:)) |
C & maxval(salt(1:sNx,1:sNy,:,:,:)) |
451 |
C write(0,*) 'dynamics: phiHyd ',minval(phiHyd/(Gravity*Rhonil),mask=phiHyd.NE.0.), |
C write(0,*) 'dynamics: phiHyd ',minval(phiHyd/(Gravity*Rhonil),mask=phiHyd.NE.0.), |
452 |
C & maxval(phiHyd/(Gravity*Rhonil)) |
C & maxval(phiHyd/(Gravity*Rhonil)) |
453 |
|
C CALL PLOT_FIELD_XYZRL( gU, ' GU exiting dyanmics ' , |
454 |
|
C &Nr, 1, myThid ) |
455 |
|
C CALL PLOT_FIELD_XYZRL( gV, ' GV exiting dyanmics ' , |
456 |
|
C &Nr, 1, myThid ) |
457 |
|
C CALL PLOT_FIELD_XYZRL( gS, ' GS exiting dyanmics ' , |
458 |
|
C &Nr, 1, myThid ) |
459 |
|
C CALL PLOT_FIELD_XYZRL( gT, ' GT exiting dyanmics ' , |
460 |
|
C &Nr, 1, myThid ) |
461 |
|
|
462 |
|
|
463 |
RETURN |
RETURN |
464 |
END |
END |