/[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.33 by jmc, Wed Jul 26 22:01:57 2006 UTC revision 1.34 by jmc, Sun Aug 13 21:43:48 2006 UTC
# Line 48  C     tmpfld  :: Temporary array used to Line 48  C     tmpfld  :: Temporary array used to
48        INTEGER I, J, K        INTEGER I, J, K
49        _RL hFacCtmp        _RL hFacCtmp
50        _RL hFacMnSz        _RL hFacMnSz
51        _RL tileArea        _RL tileArea, threadArea
52  CEOP  CEOP
53    
54  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 186  C    &         'Model Depths K Index' , Line 186  C    &         'Model Depths K Index' ,
186       &         'Model Ro_surf (ini_masks_etc)', 1, myThid)       &         'Model Ro_surf (ini_masks_etc)', 1, myThid)
187    
188  C     Calculate quantities derived from XY depth map  C     Calculate quantities derived from XY depth map
189        globalArea = 0. _d 0        threadArea = 0. _d 0
190        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
191         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
192          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
# Line 209  C-      Compute the domain Area: Line 209  C-      Compute the domain Area:
209            tileArea = tileArea + rA(i,j,bi,bj)*maskH(i,j,bi,bj)            tileArea = tileArea + rA(i,j,bi,bj)*maskH(i,j,bi,bj)
210           ENDDO           ENDDO
211          ENDDO          ENDDO
212          globalArea = globalArea + tileArea          threadArea = threadArea + tileArea
213         ENDDO         ENDDO
214        ENDDO        ENDDO
215  C     _EXCH_XY_R4(   recip_Rcol, myThid )  C     _EXCH_XY_R4(   recip_Rcol, myThid )
216        _GLOBAL_SUM_R8( globalArea, myThid )  C_jmc: apply GLOBAL_SUM to thread-local variable (not in common block)
217          _GLOBAL_SUM_R8( threadArea, myThid )
218          _BEGIN_MASTER( myThid )
219          globalArea = threadArea
220          _END_MASTER( myThid )
221    
222  C     hFacW and hFacS (at U and V points)  C     hFacW and hFacS (at U and V points)
223        DO bj=myByLo(myThid), myByHi(myThid)        DO bj=myByLo(myThid), myByHi(myThid)
# Line 247  C     We should really supply a flag for Line 251  C     We should really supply a flag for
251          DO K=1, Nr          DO K=1, Nr
252           DO J=1-Oly,sNy+Oly           DO J=1-Oly,sNy+Oly
253            DO I=1-Olx,sNx+Olx            DO I=1-Olx,sNx+Olx
254             IF (DYG(I,J,bi,bj).EQ.0.) hFacW(I,J,K,bi,bj)=0.             IF (dyG(I,J,bi,bj).EQ.0.) hFacW(I,J,K,bi,bj)=0.
255             IF (DXG(I,J,bi,bj).EQ.0.) hFacS(I,J,K,bi,bj)=0.             IF (dxG(I,J,bi,bj).EQ.0.) hFacS(I,J,K,bi,bj)=0.
256            ENDDO            ENDDO
257           ENDDO           ENDDO
258          ENDDO          ENDDO
Legend:
Removed from v.1.33  
changed lines
  Added in v.1.34
  ViewVC Help
Powered by ViewVC 1.1.22