--- MITgcm/model/src/dynamics.F 1998/09/06 17:35:20 1.33 +++ MITgcm/model/src/dynamics.F 1998/09/29 18:50:57 1.35 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/dynamics.F,v 1.33 1998/09/06 17:35:20 cnh Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/dynamics.F,v 1.35 1998/09/29 18:50:57 cnh Exp $ #include "CPP_OPTIONS.h" @@ -202,6 +202,7 @@ ENDDO ENDDO + DO bj=myByLo(myThid),myByHi(myThid) DO bi=myBxLo(myThid),myBxHi(myThid) @@ -232,6 +233,7 @@ jMin = 1-OLy+1 jMax = sNy+OLy + K = 1 BOTTOM_LAYER = K .EQ. Nr @@ -356,10 +358,10 @@ O xA,yA,uTrans,vTrans,rTrans,rVel,maskC,maskUp, I myThid) CcnhDebugStarts - IF ( K .EQ. 1 ) THEN - CALL PLOT_FIELD_XYRL( rVel(1,1,1), 'K=1 Current rVel.1 ' , myIter, myThid ) - CALL PLOT_FIELD_XYRL( rVel(1,1,2), 'K=1 Current rVel.2 ' , myIter, myThid ) - ENDIF +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 C-- Calculate the total vertical diffusivity CALL CALC_DIFFUSIVITY(