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

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

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

revision 1.12 by heimbach, Fri Jun 9 14:26:30 2000 UTC revision 1.13 by adcroft, Fri Feb 2 21:04:47 2001 UTC
# Line 102  C--   tracer cell faces. Line 102  C--   tracer cell faces.
102         ENDDO         ENDDO
103        ENDDO        ENDDO
104    
 C--   Calculate vertical "volume transport" through  
 C--   tracer cell face *above* this level.  
       IF (TOP_LAYER .AND. rigidLid) THEN  
        DO j=jMin,jMax  
         DO i=iMin,iMax  
          rTrans(i,j) = 0.  
         ENDDO  
        ENDDO  
       ELSE  
        DO j=jMin,jMax  
         DO i=iMin,iMax  
          rTrans(i,j) =  
      &    uTrans(i,j)*recip_rkFac-uTrans(i+1,j)*recip_rkFac  
      &   +vTrans(i,j)*recip_rkFac-vTrans(i,j+1)*recip_rkFac  
      &   +rTrans(i,j)  
         ENDDO  
        ENDDO  
       ENDIF  
   
105  C--   Vertical velocity at upper face  C--   Vertical velocity at upper face
106        DO j=jMin,jMax        DO j=jMin,jMax
107         DO i=iMin,iMax          DO i=iMin,iMax
108          rVel(i,j,kUp) = rTrans(i,j)/_rA(i,j,bi,bj)            rVel(i,j,kUp) = wVel(i,j,k,bi,bj)
109         ENDDO          ENDDO
110        ENDDO        ENDDO
111    
112  #ifdef ALLOW_NONHYDROSTATIC  C--   Calculate vertical "volume transport" through
113  C--   Vertical velocity at upper face  C--   tracer cell face *above* this level.
114  C     IF ( nonHydrostatic ) THEN        DO j=jMin,jMax
115          DO j=jMin,jMax          DO i=iMin,iMax
116           DO i=iMin,iMax            rTrans(i,j) = rVel(i,j,kUp)*rA(i,j,bi,bj)
           wVel(i,j,k,bi,bj)=rVel(i,j,kUp)  
          ENDDO  
117          ENDDO          ENDDO
118  C     ENDIF        ENDDO
 #endif  
119    
120        RETURN        RETURN
121        END        END
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22