/[MITgcm]/MITgcm/model/src/ini_masks_etc.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_masks_etc.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.29 by adcroft, Thu May 13 15:40:53 2004 UTC revision 1.30 by jmc, Wed Jul 20 22:31:47 2005 UTC
# Line 52  C     I,J,K Line 52  C     I,J,K
52  #endif  #endif
53        _RL hFacCtmp        _RL hFacCtmp
54        _RL hFacMnSz        _RL hFacMnSz
55          _RL tileArea
56  CEOP  CEOP
57    
58  C- Calculate lopping factor hFacC : over-estimate the part inside of the domain  C- Calculate lopping factor hFacC : over-estimate the part inside of the domain
# Line 160  C    &         'Model Depths K Index' , Line 161  C    &         'Model Depths K Index' ,
161       &         'Model Ro_surf (ini_masks_etc)', 1, myThid)       &         'Model Ro_surf (ini_masks_etc)', 1, myThid)
162    
163  C     Calculate quantities derived from XY depth map  C     Calculate quantities derived from XY depth map
164          globalArea = 0. _d 0
165        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
166         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
167          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
# Line 175  C         Inverse of fluid column thickn Line 177  C         Inverse of fluid column thickn
177            ENDIF            ENDIF
178           ENDDO           ENDDO
179          ENDDO          ENDDO
180    C-      Compute the domain Area:
181            tileArea = 0. _d 0
182            DO j=1,sNy
183             DO i=1,sNx
184              tileArea = tileArea + rA(i,j,bi,bj)*maskH(i,j,bi,bj)
185             ENDDO
186            ENDDO
187            globalArea = globalArea + tileArea
188         ENDDO         ENDDO
189        ENDDO        ENDDO
190  C     _EXCH_XY_R4(   recip_Rcol, myThid )  C     _EXCH_XY_R4(   recip_Rcol, myThid )
191          _GLOBAL_SUM_R8( globalArea, myThid )
192    
193  C     hFacW and hFacS (at U and V points)  C     hFacW and hFacS (at U and V points)
194        DO bj=myByLo(myThid), myByHi(myThid)        DO bj=myByLo(myThid), myByHi(myThid)
Legend:
Removed from v.1.29  
changed lines
  Added in v.1.30
  ViewVC Help
Powered by ViewVC 1.1.22