/[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.6 by cnh, Mon Jun 8 21:43:00 1998 UTC revision 1.7 by cnh, Tue Aug 18 16:32:41 1998 UTC
# Line 5  C $Header$ Line 5  C $Header$
5  CStartOfInterFace  CStartOfInterFace
6        SUBROUTINE CALC_COMMON_FACTORS(        SUBROUTINE CALC_COMMON_FACTORS(
7       I        bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,       I        bi,bj,iMin,iMax,jMin,jMax,k,kM1,kUp,kDown,
8       O        xA,yA,uTrans,vTrans,wTrans,wVel,maskC,maskUp,       O        xA,yA,uTrans,vTrans,rTrans,rVel,maskC,maskUp,
9       I        myThid)       I        myThid)
10    
11  C     /==========================================================\  C     /==========================================================\
# Line 37  C     xA      - Tracer cell face area no Line 37  C     xA      - Tracer cell face area no
37  C     yA      - Tracer cell face area normal to X  C     yA      - Tracer cell face area normal to X
38  C     uTrans  - Zonal volume transport through cell face  C     uTrans  - Zonal volume transport through cell face
39  C     vTrans  - Meridional volume transport through cell face  C     vTrans  - Meridional volume transport through cell face
40  C     wTrans  - Vertical volume transport through cell face  C     rTrans  - R-direction volume transport through cell face
41  C     wVel    - Vertical velocity at cell upper and lower faces  C     rVel    - R-direction velocity at cell upper and lower faces
42  C     maskC   - land/water mask for tracer points  C     maskC   - land/water mask for tracer points
43  C     maskUp  - land/water mask for Wvel points (above tracer level)  C     maskUp  - land/water mask for Wvel points (above tracer level)
44  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
# Line 48  C Line 48  C
48        _RS yA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS yA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
49        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
50        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51        _RL wTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
52        _RL wVel  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _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 69  C                         into fluxUD. Line 69  C                         into fluxUD.
69  C--   Calculate tracer cell face open areas  C--   Calculate tracer cell face open areas
70        DO j=jMin,jMax        DO j=jMin,jMax
71         DO i=iMin,iMax         DO i=iMin,iMax
72          xA(i,j) = _dyG(i,j,bi,bj)*dzF(k)*_hFacW(i,j,k,bi,bj)          xA(i,j) = _dyG(i,j,bi,bj)*drF(k)*_hFacW(i,j,k,bi,bj)
73          yA(i,j) = _dxG(i,j,bi,bj)*dzF(k)*_hFacS(i,j,k,bi,bj)          yA(i,j) = _dxG(i,j,bi,bj)*drF(k)*_hFacS(i,j,k,bi,bj)
74         ENDDO         ENDDO
75        ENDDO        ENDDO
76    
# Line 87  C--   Calculate vertical "volume transpo Line 87  C--   Calculate vertical "volume transpo
87  C--   tracer cell face *above* this level.  C--   tracer cell face *above* this level.
88        DO j=jMin,jMax        DO j=jMin,jMax
89         DO i=iMin,iMax         DO i=iMin,iMax
90          wTrans(i,j) = uTrans(i,j)-uTrans(i+1,j)          rTrans(i,j) = uTrans(i,j)*recip_rkFac-uTrans(i+1,j)*_recip_rkFac
91       &               +vTrans(i,j)-vTrans(i,j+1)       &               +vTrans(i,j)*_recip_rkFac-vTrans(i,j+1)*_recip_rkFac
92       &               +wTrans(i,j)       &               +rTrans(i,j)
93         ENDDO         ENDDO
94        ENDDO        ENDDO
95    
96  C--   Vertical velocity at upper face  C--   Vertical velocity at upper face
97        DO j=jMin,jMax        DO j=jMin,jMax
98         DO i=iMin,iMax         DO i=iMin,iMax
99          wVel(i,j,kUp) = wTrans(i,j)/_zA(i,j,bi,bj)          rVel(i,j,kUp) = rTrans(i,j)/_rA(i,j,bi,bj)
100         ENDDO         ENDDO
101        ENDDO        ENDDO
102    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22