/[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.8 by cnh, Sat Aug 22 17:51:07 1998 UTC revision 1.9 by cnh, Wed Oct 28 03:11:36 1998 UTC
# Line 18  C     | that are used at various points Line 18  C     | that are used at various points
18  C     | This reduces the amount of total work, total memory      |  C     | This reduces the amount of total work, total memory      |
19  C     | and therefore execution time and is generally a good     |  C     | and therefore execution time and is generally a good     |
20  C     | idea.                                                    |  C     | idea.                                                    |
 C     | We also think lower taxes are a good idea but we doubt   |  
 C     | whether we'll ever get them.                             |  
21  C     \==========================================================/  C     \==========================================================/
22        IMPLICIT NONE        IMPLICIT NONE
23    
# Line 69  C                         into fluxUD. Line 67  C                         into fluxUD.
67  C--   Calculate tracer cell face open areas  C--   Calculate tracer cell face open areas
68        DO j=jMin,jMax        DO j=jMin,jMax
69         DO i=iMin,iMax         DO i=iMin,iMax
70          xA(i,j) = _dyG(i,j,bi,bj)*drF(k)*_hFacW(i,j,k,bi,bj)          xA(i,j) = _dyG(i,j,bi,bj)
71          yA(i,j) = _dxG(i,j,bi,bj)*drF(k)*_hFacS(i,j,k,bi,bj)       &   *drF(k)*_hFacW(i,j,k,bi,bj)
72            yA(i,j) = _dxG(i,j,bi,bj)
73         &   *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          rTrans(i,j) = uTrans(i,j)*recip_rkFac-uTrans(i+1,j)*recip_rkFac          rTrans(i,j) =
91       &               +vTrans(i,j)*recip_rkFac-vTrans(i,j+1)*recip_rkFac       &   uTrans(i,j)*recip_rkFac-uTrans(i+1,j)*recip_rkFac
92       &               +rTrans(i,j)       &  +vTrans(i,j)*recip_rkFac-vTrans(i,j+1)*recip_rkFac
93         &  +rTrans(i,j)
94         ENDDO         ENDDO
95        ENDDO        ENDDO
96    
# Line 106  C--    Calculate mask for tracer cells Line 107  C--    Calculate mask for tracer cells
107          maskC(i,j) = 1.          maskC(i,j) = 1.
108          if (_hFacC(i,j,k,bi,bj).eq.0.) maskC(i,j)=0.          if (_hFacC(i,j,k,bi,bj).eq.0.) maskC(i,j)=0.
109          maskUp(i,j) = 1.          maskUp(i,j) = 1.
110          if (_hFacC(i,j,k,bi,bj).eq.0. .or. TOP_LAYER ) maskUp(i,j)=0.          if (_hFacC(i,j,k,bi,bj).eq.0. .or. TOP_LAYER )
111         &   maskUp(i,j)=0.
112         ENDDO         ENDDO
113        ENDDO        ENDDO
114    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22