/[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.10 by cnh, Fri Nov 6 22:44:43 1998 UTC revision 1.11 by adcroft, Mon Mar 22 15:54:03 1999 UTC
# Line 27  C     == GLobal variables == Line 27  C     == GLobal variables ==
27  #include "EEPARAMS.h"  #include "EEPARAMS.h"
28  #include "PARAMS.h"  #include "PARAMS.h"
29  #include "GRID.h"  #include "GRID.h"
30    #ifdef ALLOW_NONHYDROSTATIC
31    #include "GW.h"
32    #endif
33    
34  C     == Routine arguments ==  C     == Routine arguments ==
35  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
# Line 64  C                         into fluxUD. Line 67  C                         into fluxUD.
67    
68        TOP_LAYER = K .EQ. 1        TOP_LAYER = K .EQ. 1
69    
70    C--   Calculate mask for tracer cells  (0 => land, 1 => water)
71          DO j=jMin,jMax
72           DO i=iMin,iMax
73            maskC(i,j) = 1.
74            IF (_hFacC(i,j,k,bi,bj).eq.0.) maskC(i,j)=0.
75            maskUp(i,j) = 1.
76            IF (_hFacC(i,j,k,bi,bj).eq.0. .OR. TOP_LAYER )
77         &    maskUp(i,j)=0.
78           ENDDO
79          ENDDO
80    
81  C--   Calculate tracer cell face open areas  C--   Calculate tracer cell face open areas
82        DO j=jMin,jMax        DO j=jMin,jMax
83         DO i=iMin,iMax         DO i=iMin,iMax
# Line 85  C--   tracer cell faces. Line 99  C--   tracer cell faces.
99    
100  C--   Calculate vertical "volume transport" through  C--   Calculate vertical "volume transport" through
101  C--   tracer cell face *above* this level.  C--   tracer cell face *above* this level.
102        DO j=jMin,jMax        IF (TOP_LAYER .AND. rigidLid) THEN
103         DO i=iMin,iMax         DO j=jMin,jMax
104          rTrans(i,j) =          DO i=iMin,iMax
105       &   uTrans(i,j)*recip_rkFac-uTrans(i+1,j)*recip_rkFac           rTrans(i,j) = 0.
106       &  +vTrans(i,j)*recip_rkFac-vTrans(i,j+1)*recip_rkFac          ENDDO
      &  +rTrans(i,j)  
107         ENDDO         ENDDO
108        ENDDO        ELSE
109           DO j=jMin,jMax
110            DO i=iMin,iMax
111             rTrans(i,j) =
112         &    uTrans(i,j)*recip_rkFac-uTrans(i+1,j)*recip_rkFac
113         &   +vTrans(i,j)*recip_rkFac-vTrans(i,j+1)*recip_rkFac
114         &   +rTrans(i,j)
115            ENDDO
116           ENDDO
117          ENDIF
118    
119  C--   Vertical velocity at upper face  C--   Vertical velocity at upper face
120        DO j=jMin,jMax        DO j=jMin,jMax
# Line 101  C--   Vertical velocity at upper face Line 123  C--   Vertical velocity at upper face
123         ENDDO         ENDDO
124        ENDDO        ENDDO
125    
126  C--    Calculate mask for tracer cells  (0 => land, 1 => water)  #ifdef ALLOW_NONHYDROSTATIC
127        DO j=jMin,jMax  C--   Vertical velocity at upper face
128         DO i=iMin,iMax  C     IF ( nonHydrostatic ) THEN
129          maskC(i,j) = 1.          DO j=jMin,jMax
130          if (_hFacC(i,j,k,bi,bj).eq.0.) maskC(i,j)=0.           DO i=iMin,iMax
131          maskUp(i,j) = 1.            wVel(i,j,k,bi,bj)=rVel(i,j,kUp)
132          if (_hFacC(i,j,k,bi,bj).eq.0. .or. TOP_LAYER )           ENDDO
133       &   maskUp(i,j)=0.          ENDDO
134         ENDDO  C     ENDIF
135        ENDDO  #endif
136    
137        RETURN        RETURN
138        END        END
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22