/[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.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 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 378  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 386  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 448  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.31  
changed lines
  Added in v.1.37
  ViewVC Help
Powered by ViewVC 1.1.22