/[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.17 by jmc, Thu Sep 23 21:27:04 2004 UTC revision 1.18 by edhill, Mon Feb 7 18:37:19 2005 UTC
# Line 31  C     === Global variables === Line 31  C     === Global variables ===
31  #include "W2_EXCH2_TOPOLOGY.h"  #include "W2_EXCH2_TOPOLOGY.h"
32  #include "W2_EXCH2_PARAMS.h"  #include "W2_EXCH2_PARAMS.h"
33  #endif  #endif
34    #ifdef ALLOW_MNC
35    #include "MNC_PARAMS.h"
36    #endif
37    
38  #ifndef ALLOW_EXCH2  #ifndef ALLOW_EXCH2
39  C- note: default is to use "new" grid files (OLD_GRID_IO undef) with EXCH2  C- note: default is to use "new" grid files (OLD_GRID_IO undef) with EXCH2
# Line 48  C     == Local variables == Line 51  C     == Local variables ==
51        INTEGER bi,bj, myTile, myiter        INTEGER bi,bj, myTile, myiter
52        INTEGER I,J        INTEGER I,J
53        CHARACTER*(15) fName        CHARACTER*(15) fName
54    #ifdef ALLOW_MNC
55          CHARACTER*(80) mncFn
56    #endif
57        _RL buf(sNx+1,sNy+1)        _RL buf(sNx+1,sNy+1)
58        INTEGER iG, iL        INTEGER iG, iL
59        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
# Line 89  C--   Set everything to zero everywhere Line 95  C--   Set everything to zero everywhere
95         ENDDO         ENDDO
96        ENDDO        ENDDO
97    
98    
99    #ifdef ALLOW_MNC
100          IF (useMNC .AND. readgrid_mnc) THEN
101    
102            _BEGIN_MASTER(myThid)
103            DO i = 1,80
104              mncFn(i:i) = ' '
105            ENDDO
106            write(mncFn,'(a)') 'mitgrid'
107            DO i = 1,MAX_LEN_MBUF
108              msgBuf(i:i) = ' '
109            ENDDO
110            WRITE(msgBuf,'(2A)') msgBuf,' ; Reading grid info using MNC'
111            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
112         &       SQUEEZE_RIGHT , myThid)
113            CALL MNC_FILE_CLOSE_ALL_MATCHING(mncFn, myThid)
114            CALL MNC_CW_SET_UDIM(mncFn, 1, myThid)
115            CALL MNC_CW_RS_R('R',mncFn,0,0,'XC', XC,  myThid)
116            CALL MNC_CW_RS_R('R',mncFn,0,0,'XG', XG,  myThid)
117            CALL MNC_CW_RS_R('R',mncFn,0,0,'YC', YC,  myThid)
118            CALL MNC_CW_RS_R('R',mncFn,0,0,'YG', YG,  myThid)
119            CALL MNC_CW_RS_R('R',mncFn,0,0,'dxC',DXC, myThid)
120            CALL MNC_CW_RS_R('R',mncFn,0,0,'dyC',DYC, myThid)
121            CALL MNC_CW_RS_R('R',mncFn,0,0,'dxF',DXF, myThid)
122            CALL MNC_CW_RS_R('R',mncFn,0,0,'dyF',DYF, myThid)
123            CALL MNC_CW_RS_R('R',mncFn,0,0,'dxG',DXG, myThid)
124            CALL MNC_CW_RS_R('R',mncFn,0,0,'dyG',DYG, myThid)
125            CALL MNC_CW_RS_R('R',mncFn,0,0,'dxV',DXV, myThid)
126            CALL MNC_CW_RS_R('R',mncFn,0,0,'dyU',DYU, myThid)
127            CALL MNC_CW_RS_R('R',mncFn,0,0,'rA', RA,  myThid)
128            CALL MNC_CW_RS_R('R',mncFn,0,0,'rAz',RAZ, myThid)
129            CALL MNC_CW_RS_R('R',mncFn,0,0,'rAw',RAW, myThid)
130            CALL MNC_CW_RS_R('R',mncFn,0,0,'rAs',RAS, myThid)
131    
132            _END_MASTER(myThid)
133    
134            CALL EXCH_XY_RS(XC,myThid)
135            CALL EXCH_XY_RS(YC,myThid)
136    #ifdef HRCUBE
137            CALL EXCH_XY_RS(DXF,myThid)
138            CALL EXCH_XY_RS(DYF,myThid)
139    #endif
140            CALL EXCH_XY_RS(RA,myThid )
141            CALL EXCH_Z_XY_RS(XG,myThid)
142            CALL EXCH_Z_XY_RS(YG,myThid)
143            CALL EXCH_Z_XY_RS(RAZ,myThid)
144            CALL EXCH_UV_XY_RS(DXC,DYC,.FALSE.,myThid)
145            CALL EXCH_UV_XY_RS(RAW,RAS,.FALSE.,myThid)
146            CALL EXCH_UV_XY_RS(DYG,DXG,.FALSE.,myThid)
147    
148          ELSE
149    #endif
150    
151  C     Here we make no assumptions about grid symmetry and simply  C     Here we make no assumptions about grid symmetry and simply
152  C     read the raw grid data from files  C     read the raw grid data from files
153    
# Line 337  C--   Only do I/O if I am the master thr Line 396  C--   Only do I/O if I am the master thr
396    
397         ENDDO         ENDDO
398        ENDDO        ENDDO
399    
400        _END_MASTER(myThid)        _END_MASTER(myThid)
401    
402        CALL EXCH_XY_RS(XC,myThid)        CALL EXCH_XY_RS(XC,myThid)
# Line 359  c     CALL EXCH_Z_XY_RS(DYU,myThid) Line 419  c     CALL EXCH_Z_XY_RS(DYU,myThid)
419    
420  #endif /* OLD_GRID_IO */  #endif /* OLD_GRID_IO */
421    
422    #ifdef ALLOW_MNC
423          ENDIF
424    #endif /* ALLOW_MNC */
425    
426  c     CALL WRITE_FULLARRAY_RL('DXV',DXV,1,0,0,0,myThid)  c     CALL WRITE_FULLARRAY_RL('DXV',DXV,1,0,0,0,myThid)
427  c     CALL WRITE_FULLARRAY_RL('DYU',DYU,1,0,0,0,myThid)  c     CALL WRITE_FULLARRAY_RL('DYU',DYU,1,0,0,0,myThid)
428  c     CALL WRITE_FULLARRAY_RL('RAZ',RAZ,1,0,0,0,myThid)  c     CALL WRITE_FULLARRAY_RL('RAZ',RAZ,1,0,0,0,myThid)
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22