/[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.20 by cnh, Mon Jun 15 05:17:42 1998 UTC revision 1.21 by adcroft, Mon Jun 22 15:26:25 1998 UTC
# Line 223  C--     Calculate gradient of surface pr Line 223  C--     Calculate gradient of surface pr
223    
224  C--     Update fields in top level according to tendency terms  C--     Update fields in top level according to tendency terms
225          CALL CORRECTION_STEP(          CALL CORRECTION_STEP(
226       I       bi,bj,iMin,iMax,jMin,jMax,K,pSurfX,pSurfY,myThid)       I       bi,bj,iMin,iMax,jMin,jMax,K,pSurfX,pSurfY,myTime,myThid)
227    
228            IF ( .NOT. BOTTOM_LAYER ) THEN
229    C--      Update fields in layer below according to tendency terms
230             CALL CORRECTION_STEP(
231         I        bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)
232            ENDIF
233    
234  C--     Density of 1st level (below W(1)) reference to level 1  C--     Density of 1st level (below W(1)) reference to level 1
235          CALL FIND_RHO(          CALL FIND_RHO(
# Line 258  C--     Integrate hydrostatic balance fo Line 264  C--     Integrate hydrostatic balance fo
264    
265           BOTTOM_LAYER = K .EQ. Nz           BOTTOM_LAYER = K .EQ. Nz
266    
267  C--      Update fields in Kth level according to tendency terms           IF ( .NOT. BOTTOM_LAYER ) THEN
268           CALL CORRECTION_STEP(  C--       Update fields in layer below according to tendency terms
269       I        bi,bj,iMin,iMax,jMin,jMax,K,pSurfX,pSurfY,myThid)            CALL CORRECTION_STEP(
270         I         bi,bj,iMin,iMax,jMin,jMax,K+1,pSurfX,pSurfY,myTime,myThid)
271             ENDIF
272    C--      Update fields in layer below according to tendency terms
273    C        CALL CORRECTION_STEP(
274    C    I        bi,bj,iMin,iMax,jMin,jMax,K,pSurfX,pSurfY,myTime,myThid)
275    
276  C--      Density of K level (below W(K)) reference to K level  C--      Density of K level (below W(K)) reference to K level
277           CALL FIND_RHO(           CALL FIND_RHO(
278       I      bi, bj, iMin, iMax, jMin, jMax,  K, K, eosType,       I      bi, bj, iMin, iMax, jMin, jMax,  K, K, eosType,
# Line 343  C--      Calculate active tracer tendenc Line 355  C--      Calculate active tracer tendenc
355           IF ( tempStepping ) THEN           IF ( tempStepping ) THEN
356            CALL CALC_GT(            CALL CALC_GT(
357       I         bi,bj,iMin,iMax,jMin,jMax, k,kM1,kUp,kDown,       I         bi,bj,iMin,iMax,jMin,jMax, k,kM1,kUp,kDown,
358       I         xA,yA,uTrans,vTrans,wTrans,maskUp,       I         xA,yA,uTrans,vTrans,wTrans,maskUp,maskC,
359       I         K13,K23,KappaZT,KapGM,       I         K13,K23,KappaZT,KapGM,
360       U         aTerm,xTerm,fZon,fMer,fVerT,       U         aTerm,xTerm,fZon,fMer,fVerT,
361       I         myThid)       I         myThid)
# Line 351  C--      Calculate active tracer tendenc Line 363  C--      Calculate active tracer tendenc
363           IF ( saltStepping ) THEN           IF ( saltStepping ) THEN
364            CALL CALC_GS(            CALL CALC_GS(
365       I         bi,bj,iMin,iMax,jMin,jMax, k,kM1,kUp,kDown,       I         bi,bj,iMin,iMax,jMin,jMax, k,kM1,kUp,kDown,
366       I         xA,yA,uTrans,vTrans,wTrans,maskUp,       I         xA,yA,uTrans,vTrans,wTrans,maskUp,maskC,
367       I         K13,K23,KappaZS,KapGM,       I         K13,K23,KappaZS,KapGM,
368       U         aTerm,xTerm,fZon,fMer,fVerS,       U         aTerm,xTerm,fZon,fMer,fVerS,
369       I         myThid)       I         myThid)
# Line 383  C--     Implicit diffusion Line 395  C--     Implicit diffusion
395  C     write(0,*) 'dynamics: pS ',minval(cg2d_x(1:sNx,1:sNy,:,:)),  C     write(0,*) 'dynamics: pS ',minval(cg2d_x(1:sNx,1:sNy,:,:)),
396  C    &                           maxval(cg2d_x(1:sNx,1:sNy,:,:))  C    &                           maxval(cg2d_x(1:sNx,1:sNy,:,:))
397  C     write(0,*) 'dynamics: U  ',minval(uVel(1:sNx,1:sNy,1,:,:),mask=uVel(1:sNx,1:sNy,1,:,:).NE.0.),  C     write(0,*) 'dynamics: U  ',minval(uVel(1:sNx,1:sNy,1,:,:),mask=uVel(1:sNx,1:sNy,1,:,:).NE.0.),
398  C    &                           maxval(uVel(1:sNx,1:sNy,1,:,:))  C    &                           maxval(uVel(1:sNx,1:sNy,1,:,:),mask=uVel(1:sNx,1:sNy,1,:,:).NE.0.)
399  C     write(0,*) 'dynamics: V  ',minval(vVel(1:sNx,1:sNy,1,:,:),mask=vVel(1:sNx,1:sNy,1,:,:).NE.0.),  C     write(0,*) 'dynamics: V  ',minval(vVel(1:sNx,1:sNy,1,:,:),mask=vVel(1:sNx,1:sNy,1,:,:).NE.0.),
400  C    &                           maxval(vVel(1:sNx,1:sNy,1,:,:))  C    &                           maxval(vVel(1:sNx,1:sNy,1,:,:),mask=vVel(1:sNx,1:sNy,1,:,:).NE.0.)
401  C     write(0,*) 'dynamics: wVel(1) ',  C     write(0,*) 'dynamics: wVel(1) ',
402  C    &            minval(wVel(1:sNx,1:sNy,1),mask=wVel(1:sNx,1:sNy,1).NE.0.),  C    &            minval(wVel(1:sNx,1:sNy,1),mask=wVel(1:sNx,1:sNy,1).NE.0.),
403  C    &            maxval(wVel(1:sNx,1:sNy,1))  C    &            maxval(wVel(1:sNx,1:sNy,1),mask=wVel(1:sNx,1:sNy,1).NE.0.)
404  C     write(0,*) 'dynamics: wVel(2) ',  C     write(0,*) 'dynamics: wVel(2) ',
405  C    &            minval(wVel(1:sNx,1:sNy,2),mask=wVel(1:sNx,1:sNy,2).NE.0.),  C    &            minval(wVel(1:sNx,1:sNy,2),mask=wVel(1:sNx,1:sNy,2).NE.0.),
406  C    &            maxval(wVel(1:sNx,1:sNy,2))  C    &            maxval(wVel(1:sNx,1:sNy,2),mask=wVel(1:sNx,1:sNy,2).NE.0.)
407  cblk  write(0,*) 'dynamics: K13',minval(K13(1:sNx,1:sNy,:)),  cblk  write(0,*) 'dynamics: K13',minval(K13(1:sNx,1:sNy,:)),
408  cblk &                           maxval(K13(1:sNx,1:sNy,:))  cblk &                           maxval(K13(1:sNx,1:sNy,:))
409  cblk  write(0,*) 'dynamics: K23',minval(K23(1:sNx,1:sNy,:)),  cblk  write(0,*) 'dynamics: K23',minval(K23(1:sNx,1:sNy,:)),
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22