/[MITgcm]/MITgcm/model/src/dynamics.F
ViewVC logotype

Diff of /MITgcm/model/src/dynamics.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.34 by cnh, Mon Sep 7 16:23:11 1998 UTC revision 1.37 by cnh, Tue Nov 3 15:28:04 1998 UTC
# Line 202  C     uninitialised but inert locations. Line 202  C     uninitialised but inert locations.
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    
# Line 232  C--     Set up work arrays that need val Line 233  C--     Set up work arrays that need val
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    
# Line 355  C--      Get temporary terms used by ten Line 357  C--      Get temporary terms used by ten
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,
# Line 384  C--      Calculate active tracer tendenc Line 380  C--      Calculate active tracer tendenc
380       I         xA,yA,uTrans,vTrans,rTrans,maskUp,maskC,       I         xA,yA,uTrans,vTrans,rTrans,maskUp,maskC,
381       I         K13,K23,KappaRT,KapGM,       I         K13,K23,KappaRT,KapGM,
382       U         aTerm,xTerm,fZon,fMer,fVerT,       U         aTerm,xTerm,fZon,fMer,fVerT,
383       I         myThid)       I         myTime, myThid)
384           ENDIF           ENDIF
385           IF ( saltStepping ) THEN           IF ( saltStepping ) THEN
386            CALL CALC_GS(            CALL CALC_GS(
# Line 392  C--      Calculate active tracer tendenc Line 388  C--      Calculate active tracer tendenc
388       I         xA,yA,uTrans,vTrans,rTrans,maskUp,maskC,       I         xA,yA,uTrans,vTrans,rTrans,maskUp,maskC,
389       I         K13,K23,KappaRS,KapGM,       I         K13,K23,KappaRS,KapGM,
390       U         aTerm,xTerm,fZon,fMer,fVerS,       U         aTerm,xTerm,fZon,fMer,fVerS,
391       I         myThid)       I         myTime, myThid)
392           ENDIF           ENDIF
393  C--      Prediction step (step forward all model variables)  C--      Prediction step (step forward all model variables)
394           CALL TIMESTEP(           CALL TIMESTEP(
# Line 454  C     write(0,*) 'dynamics: S  ',minval( Line 450  C     write(0,*) 'dynamics: S  ',minval(
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
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.37
  ViewVC Help
Powered by ViewVC 1.1.22