/[MITgcm]/MITgcm/pkg/exch2/w2_set_map_tiles.F
ViewVC logotype

Diff of /MITgcm/pkg/exch2/w2_set_map_tiles.F

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

revision 1.5 by jmc, Sat Jul 9 21:53:35 2011 UTC revision 1.6 by jmc, Mon Jul 16 20:25:10 2012 UTC
# Line 105  C--   Check that domain size and (SIZE.h Line 105  C--   Check that domain size and (SIZE.h
105        IF ( W2_mapIO.EQ.1 ) THEN        IF ( W2_mapIO.EQ.1 ) THEN
106  C--   Compact IO map (mostly in Y dir): search for Greatest Common Divisor  C--   Compact IO map (mostly in Y dir): search for Greatest Common Divisor
107  C     of all x-size (faster to apply GCD to Nb of Tiles in X):  C     of all x-size (faster to apply GCD to Nb of Tiles in X):
108            k = 0
109            nnx(1) = 0
110          DO j=1,nFacets          DO j=1,nFacets
111            nnx(j) = facet_dims(2*j-1)/tNx  C     skip empty facet
112              IF ( facet_dims(2*j-1).GT.0 ) THEN
113                k = k + 1
114                nnx(k) = facet_dims(2*j-1)/tNx
115              ENDIF
116          ENDDO          ENDDO
117          divide = FIND_GCD_N( nnx, nFacets)          divide = FIND_GCD_N( nnx, k )
118          W2_mapIO = divide*tNx          W2_mapIO = divide*tNx
119          WRITE(msgBuf,'(A,2(I5,A))') ' W2_mapIO =', W2_mapIO,          WRITE(msgBuf,'(A,2(I5,A))') ' W2_mapIO =', W2_mapIO,
120       &                              ' (=', divide, '*sNx)'       &                              ' (=', divide, '*sNx)'
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22