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

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

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

revision 1.3 by jmc, Fri Apr 23 20:21:06 2010 UTC revision 1.4 by jmc, Sat Jul 9 21:53:35 2011 UTC
# Line 101  C-    Only 1rst dim is set: assuming a r Line 101  C-    Only 1rst dim is set: assuming a r
101          nBl = nRd          nBl = nRd
102        ELSEIF ( nRd+nGr+nBl.EQ.0 ) THEN        ELSEIF ( nRd+nGr+nBl.EQ.0 ) THEN
103  C-    try to get cube size from number of tiles, assuming a regular Cube  C-    try to get cube size from number of tiles, assuming a regular Cube
104          nGr = nTiles*sNx*sNy          nGr = exch2_nTiles*sNx*sNy
105          tmpVar = FLOAT(nGr)/6.          tmpVar = FLOAT(nGr)/6.
106          tmpVar = SQRT(tmpVar)          tmpVar = SQRT(tmpVar)
107          nRd = NINT(tmpVar)          nRd = NINT(tmpVar)
# Line 113  C-    try to get cube size from number o Line 113  C-    try to get cube size from number o
113            CALL PRINT_MESSAGE( msgBuf, W2_oUnit,SQUEEZE_RIGHT,myThid )            CALL PRINT_MESSAGE( msgBuf, W2_oUnit,SQUEEZE_RIGHT,myThid )
114          ELSE          ELSE
115            WRITE(msgBuf,'(3(A,I4),A,I10,A,I6,A)')            WRITE(msgBuf,'(3(A,I4),A,I10,A,I6,A)')
116       &     ' nTiles*sNx*sNy=', nTiles, ' x',sNx,' x',sNy,' =',nGr       &     ' nTiles*sNx*sNy=', exch2_nTiles,' x',sNx,' x',sNy,' =',nGr
117            CALL PRINT_MESSAGE( msgBuf, W2_oUnit,SQUEEZE_RIGHT,myThid )            CALL PRINT_MESSAGE( msgBuf, W2_oUnit,SQUEEZE_RIGHT,myThid )
118            WRITE(msgBuf,'(A,I6,A,I10)')            WRITE(msgBuf,'(A,I6,A,I10)')
119       &     '       not equal to: 6 x',nRd,'^2 =', nRd*nRd*6       &     '       not equal to: 6 x',nRd,'^2 =', nRd*nRd*6
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22