/[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.7 by cnh, Tue Aug 18 16:32:41 1998 UTC revision 1.10 by cnh, Fri Nov 6 22:44:43 1998 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  CStartOfInterFace  CStartOfInterFace
6        SUBROUTINE CALC_COMMON_FACTORS(        SUBROUTINE CALC_COMMON_FACTORS(
# 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.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22