/[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.14 by cnh, Sun Feb 4 14:38:45 2001 UTC revision 1.15 by jmc, Wed Feb 7 21:48:02 2001 UTC
# Line 6  C $Name$ Line 6  C $Name$
6  CStartOfInterFace  CStartOfInterFace
7        SUBROUTINE CALC_COMMON_FACTORS(        SUBROUTINE CALC_COMMON_FACTORS(
8       I        bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,       I        bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
9       O        xA,yA,uTrans,vTrans,rTrans,rVel,maskC,maskUp,       O        xA,yA,uTrans,vTrans,rTrans,maskC,maskUp,
10       I        myThid)       I        myThid)
11    
12  C     /==========================================================\  C     /==========================================================\
# Line 40  C     yA      - Tracer cell face area no Line 40  C     yA      - Tracer cell face area no
40  C     uTrans  - Zonal volume transport through cell face  C     uTrans  - Zonal volume transport through cell face
41  C     vTrans  - Meridional volume transport through cell face  C     vTrans  - Meridional volume transport through cell face
42  C     rTrans  - R-direction volume transport through cell face  C     rTrans  - R-direction volume transport through cell face
 C     rVel    - R-direction velocity at cell upper and lower faces  
43  C     maskC   - land/water mask for tracer points  C     maskC   - land/water mask for tracer points
44  C     maskUp  - land/water mask for Wvel points (above tracer level)  C     maskUp  - land/water mask for Wvel points (above tracer level)
45  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
# Line 51  C Line 50  C
50        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
52        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       _RL rVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)  
53        _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS maskC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
54        _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS maskUp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55  C  C
# Line 66  C                         into fluxUD. Line 64  C                         into fluxUD.
64        INTEGER i,j        INTEGER i,j
65        LOGICAL TOP_LAYER        LOGICAL TOP_LAYER
66    
 #ifdef ALLOW_AUTODIFF_TAMC  
 C--   rvel(:,:kDown) is still required  
       rVel(1,1,kDown) = rVel(1,1,kDown)  
 #endif  
   
67        TOP_LAYER = K .EQ. 1        TOP_LAYER = K .EQ. 1
68    
69  C--   Calculate mask for tracer cells  (0 => land, 1 => water)  C--   Calculate mask for tracer cells  (0 => land, 1 => water)
# Line 103  C--   tracer cell faces. Line 96  C--   tracer cell faces.
96         ENDDO         ENDDO
97        ENDDO        ENDDO
98    
 C--   Vertical velocity at upper face  
       DO j=jMin,jMax  
         DO i=iMin,iMax  
           rVel(i,j,kUp) = wVel(i,j,k,bi,bj)  
         ENDDO  
       ENDDO  
   
99  C--   Calculate vertical "volume transport" through  C--   Calculate vertical "volume transport" through
100  C--   tracer cell face *above* this level.  C--   tracer cell face *above* this level.
101        DO j=jMin,jMax        DO j=jMin,jMax
102          DO i=iMin,iMax          DO i=iMin,iMax
103            rTrans(i,j) = rVel(i,j,kUp)*rA(i,j,bi,bj)            rTrans(i,j) = wVel(i,j,k,bi,bj)*rA(i,j,bi,bj)
104          ENDDO          ENDDO
105        ENDDO        ENDDO
106    
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22