/[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.34 by jmc, Sun Aug 13 21:43:48 2006 UTC revision 1.35 by jmc, Tue Sep 4 17:04:02 2007 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, threadArea        _RL tileArea(nSx,nSy), 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 203  C         Inverse of fluid column thickn Line 203  C         Inverse of fluid column thickn
203           ENDDO           ENDDO
204          ENDDO          ENDDO
205  C-      Compute the domain Area:  C-      Compute the domain Area:
206          tileArea = 0. _d 0          tileArea(bi,bj) = 0. _d 0
207          DO j=1,sNy          DO j=1,sNy
208           DO i=1,sNx           DO i=1,sNx
209            tileArea = tileArea + rA(i,j,bi,bj)*maskH(i,j,bi,bj)            tileArea(bi,bj) = tileArea(bi,bj)
210         &                    + rA(i,j,bi,bj)*maskH(i,j,bi,bj)
211           ENDDO           ENDDO
212          ENDDO          ENDDO
213          threadArea = threadArea + tileArea          threadArea = threadArea + tileArea(bi,bj)
214         ENDDO         ENDDO
215        ENDDO        ENDDO
216  C     _EXCH_XY_R4(   recip_Rcol, myThid )  C     _EXCH_XY_R4(   recip_Rcol, myThid )
217    #ifdef ALLOW_AUTODIFF_TAMC
218  C_jmc: apply GLOBAL_SUM to thread-local variable (not in common block)  C_jmc: apply GLOBAL_SUM to thread-local variable (not in common block)
219        _GLOBAL_SUM_R8( threadArea, myThid )        _GLOBAL_SUM_R8( threadArea, myThid )
220    #else
221          CALL GLOBAL_SUM_TILE_RL( tileArea, threadArea, myThid )
222    #endif
223        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
224        globalArea = threadArea        globalArea = threadArea
225        _END_MASTER( myThid )        _END_MASTER( myThid )
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.35
  ViewVC Help
Powered by ViewVC 1.1.22