/[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.31 by cnh, Sat Aug 22 17:51:08 1998 UTC revision 1.35 by cnh, Tue Sep 29 18:50:57 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 273  C--      Recompute density after mixing Line 275  C--      Recompute density after mixing
275       I      myThid )       I      myThid )
276          ENDIF          ENDIF
277  C--     Calculate buoyancy  C--     Calculate buoyancy
278          CALL CALC_BUOY(          CALL CALC_BUOYANCY(
279       I      bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1,       I      bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1,
280       O      buoyKm1,       O      buoyKm1,
281       I      myThid )       I      myThid )
# Line 313  C--       Recompute density after mixing Line 315  C--       Recompute density after mixing
315       I       myThid )       I       myThid )
316           ENDIF           ENDIF
317  C--      Calculate buoyancy  C--      Calculate buoyancy
318           CALL CALC_BUOY(           CALL CALC_BUOYANCY(
319       I       bi,bj,iMin,iMax,jMin,jMax,K,rhoK,       I       bi,bj,iMin,iMax,jMin,jMax,K,rhoK,
320       O       buoyK,       O       buoyK,
321       I       myThid )       I       myThid )
# 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)
360    CcnhDebugStarts
361    C        IF ( K .EQ. 1 ) THEN
362    C          CALL PLOT_FIELD_XYRL( rVel(1,1,1), 'K=1 Current rVel.1 ' , myIter, myThid )
363    C          CALL PLOT_FIELD_XYRL( rVel(1,1,2), 'K=1 Current rVel.2 ' , myIter, myThid )
364    C        ENDIF
365    CcnhDebugEnds
366  C--      Calculate the total vertical diffusivity  C--      Calculate the total vertical diffusivity
367           CALL CALC_DIFFUSIVITY(           CALL CALC_DIFFUSIVITY(
368       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
Legend:
Removed from v.1.31  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22