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

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

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

revision 1.7 by jmc, Mon Jan 26 01:17:58 2004 UTC revision 1.8 by jmc, Mon Jan 26 23:38:53 2004 UTC
# Line 26  C     === Global variables === Line 26  C     === Global variables ===
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
27  #include "PARAMS.h"  #include "PARAMS.h"
28  #include "GRID.h"  #include "GRID.h"
 #include "EESUPPORT.h"  
29  #ifdef USE_W2  #ifdef USE_W2
30  #include "W2_EXCH2_TOPOLOGY.h"  #include "W2_EXCH2_TOPOLOGY.h"
31  #include "W2_EXCH2_PARAMS.h"  #include "W2_EXCH2_PARAMS.h"
# Line 49  C     == Local variables == Line 48  C     == Local variables ==
48        INTEGER I,J        INTEGER I,J
49        CHARACTER*(15) fName        CHARACTER*(15) fName
50        _RL buf(sNx+1,sNy+1)        _RL buf(sNx+1,sNy+1)
51          INTEGER iG, iL
52          CHARACTER*(MAX_LEN_MBUF) msgBuf
53          INTEGER  ILNBLNK
54          EXTERNAL ILNBLNK
55  CEOP  CEOP
56    
57  C--   Set everything to zero everywhere  C--   Set everything to zero everywhere
# Line 245  c     write(10) DYG Line 248  c     write(10) DYG
248    
249  #else  #else
250    
251        DO bj = myByLo(myThid), myByHi(myThid)  C--   Only do I/O if I am the master thread
        DO bi = myBxLo(myThid), myBxHi(myThid)  
252        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
253  #ifdef ALLOW_USE_MPI  
254        write(fName(1:15),'("tile",I3.3,".mitgrid")') myPid+1        DO bj = 1,nSy
255  #else         DO bi = 1,nSx
256        write(fName(1:15),'("tile",I3.3,".mitgrid")') bi          iG=bi+(myXGlobalLo-1)/sNx
257  #endif          WRITE(fName(1:15),'("tile",I3.3,".mitgrid")') iG
258            WRITE(msgBuf,'(A,I4)') 'tile:',iG
259  #ifdef USE_W2  #ifdef USE_W2
260        myTile = W2_myTileList(bi)        myTile = W2_myTileList(bi)
261        write(fName(1:15),'("tile",I3.3,".mitgrid")')        write(fName(1:15),'("tile",I3.3,".mitgrid")')
262       &  exch2_myface(myTile)       &  exch2_myface(myTile)
263          WRITE(msgBuf,'(A,I4)') 'tile:',myTile
264  #endif  #endif
265            iL = ILNBLNK(msgBuf)
266            WRITE(msgBuf,'(3A)') msgBuf(1:iL),
267         &                   ' ; Read from file ',fName(1:15)
268            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
269         &                      SQUEEZE_RIGHT , myThid)
270            WRITE(msgBuf,'(A)') '  =>'
271    
272            CALL READSYMTILE_RS(fName,1,XC,bi,bj,buf,myThid)
273            iL = ILNBLNK(msgBuf)
274            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'XC'
275            CALL READSYMTILE_RS(fName,2,YC,bi,bj,buf,myThid)
276            iL = ILNBLNK(msgBuf)
277            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'YC'
278            CALL READSYMTILE_RS(fName,3,DXF,bi,bj,buf,myThid)
279            iL = ILNBLNK(msgBuf)
280            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DXF'
281            CALL READSYMTILE_RS(fName,4,DYF,bi,bj,buf,myThid)
282            iL = ILNBLNK(msgBuf)
283            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DYF'
284            CALL READSYMTILE_RS(fName,5,RA,bi,bj,buf,myThid)
285            iL = ILNBLNK(msgBuf)
286            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'RA'
287            CALL READSYMTILE_RS(fName,6,XG,bi,bj,buf,myThid)
288            iL = ILNBLNK(msgBuf)
289            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'XG'
290            CALL READSYMTILE_RS(fName,7,YG,bi,bj,buf,myThid)
291            iL = ILNBLNK(msgBuf)
292            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'YG'
293            CALL READSYMTILE_RS(fName,8,DXV,bi,bj,buf,myThid)
294            iL = ILNBLNK(msgBuf)
295            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DXV'
296            CALL READSYMTILE_RS(fName,9,DYU,bi,bj,buf,myThid)
297            iL = ILNBLNK(msgBuf)
298            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DYU'
299            CALL READSYMTILE_RS(fName,10,RAZ,bi,bj,buf,myThid)
300            iL = ILNBLNK(msgBuf)
301            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'RAZ'
302            CALL READSYMTILE_RS(fName,11,DXC,bi,bj,buf,myThid)
303            iL = ILNBLNK(msgBuf)
304            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DXC'
305            CALL READSYMTILE_RS(fName,12,DYC,bi,bj,buf,myThid)
306            iL = ILNBLNK(msgBuf)
307            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DYC'
308            CALL READSYMTILE_RS(fName,13,RAW,bi,bj,buf,myThid)
309            iL = ILNBLNK(msgBuf)
310            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'RAW'
311            CALL READSYMTILE_RS(fName,14,RAS,bi,bj,buf,myThid)
312            iL = ILNBLNK(msgBuf)
313            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'RAS'
314            CALL READSYMTILE_RS(fName,15,DXG,bi,bj,buf,myThid)
315            iL = ILNBLNK(msgBuf)
316            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DXG'
317            CALL READSYMTILE_RS(fName,16,DYG,bi,bj,buf,myThid)
318            iL = ILNBLNK(msgBuf)
319            WRITE(msgBuf,'(A,1X,A)') msgBuf(1:iL),'DYG'
320    
321            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
322         &                      SQUEEZE_RIGHT , myThid)
323    
       CALL READSYMTILE_RS(fName,1,XC,bi,bj,buf,myThid)  
       write(0,*) 'Read XC'  
       CALL READSYMTILE_RS(fName,2,YC,bi,bj,buf,myThid)  
       write(0,*) 'Read YC'  
       CALL READSYMTILE_RS(fName,3,DXF,bi,bj,buf,myThid)  
       write(0,*) 'Read DXF'  
       CALL READSYMTILE_RS(fName,4,DYF,bi,bj,buf,myThid)  
       write(0,*) 'Read DYF'  
       CALL READSYMTILE_RS(fName,5,RA,bi,bj,buf,myThid)  
       write(0,*) 'Read RA'  
       CALL READSYMTILE_RS(fName,6,XG,bi,bj,buf,myThid)  
       write(0,*) 'Read XG'  
       CALL READSYMTILE_RS(fName,7,YG,bi,bj,buf,myThid)  
       write(0,*) 'Read YG'  
       CALL READSYMTILE_RS(fName,8,DXV,bi,bj,buf,myThid)  
       write(0,*) 'Read DXV'  
       CALL READSYMTILE_RS(fName,9,DYU,bi,bj,buf,myThid)  
       write(0,*) 'Read DYU'  
       CALL READSYMTILE_RS(fName,10,RAZ,bi,bj,buf,myThid)  
       write(0,*) 'Read RAZ'  
       CALL READSYMTILE_RS(fName,11,DXC,bi,bj,buf,myThid)  
       write(0,*) 'Read DXC'  
       CALL READSYMTILE_RS(fName,12,DYC,bi,bj,buf,myThid)  
       write(0,*) 'Read DYC'  
       CALL READSYMTILE_RS(fName,13,RAW,bi,bj,buf,myThid)  
       write(0,*) 'Read RAW'  
       CALL READSYMTILE_RS(fName,14,RAS,bi,bj,buf,myThid)  
       write(0,*) 'Read RAS'  
       CALL READSYMTILE_RS(fName,15,DXG,bi,bj,buf,myThid)  
       write(0,*) 'Read DXG'  
       CALL READSYMTILE_RS(fName,16,DYG,bi,bj,buf,myThid)  
       write(0,*) 'Read DYG'  
       _END_MASTER(myThid)  
324         ENDDO         ENDDO
325        ENDDO        ENDDO
326          _END_MASTER(myThid)
327    
328        CALL EXCH_XY_RS(XC,myThid)        CALL EXCH_XY_RS(XC,myThid)
329        CALL EXCH_XY_RS(YC,myThid)        CALL EXCH_XY_RS(YC,myThid)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22