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

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

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

revision 1.16 by adcroft, Thu May 13 15:40:53 2004 UTC revision 1.17 by afe, Thu Jun 24 20:25:44 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
 #include "PACKAGES_CONFIG.h"  
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CBOP  CBOP
# Line 54  C     !LOCAL VARIABLES: Line 53  C     !LOCAL VARIABLES:
53  C     == Local variables ==  C     == Local variables ==
54  C     msgBuf - Used for informational I/O.  C     msgBuf - Used for informational I/O.
55        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
 #ifdef ALLOW_EXF  
       INTEGER i, j, bi, bj  
 #endif /* ALLOW_EXF */  
56  CEOP  CEOP
57    
58  C--   Set up vertical grid and coordinate system  C--   Set up vertical grid and coordinate system
# Line 69  C--   Set up horizontal grid and coordin Line 65  C--   Set up horizontal grid and coordin
65         CALL INI_SPHERICAL_POLAR_GRID( myThid )         CALL INI_SPHERICAL_POLAR_GRID( myThid )
66        ELSEIF ( usingCurvilinearGrid ) THEN        ELSEIF ( usingCurvilinearGrid ) THEN
67         CALL INI_CURVILINEAR_GRID( myThid )         CALL INI_CURVILINEAR_GRID( myThid )
68          ELSEIF ( usingCylindricalGrid ) THEN
69           CALL INI_CYLINDER( myThid )
70        ELSE        ELSE
71         _BEGIN_MASTER(myThid)         _BEGIN_MASTER(myThid)
72          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
# Line 78  C--   Set up horizontal grid and coordin Line 76  C--   Set up horizontal grid and coordin
76         _END_MASTER(myThid)         _END_MASTER(myThid)
77        ENDIF        ENDIF
78    
 #ifdef ALLOW_EXF  
 C--   exf_interp assumes that 0 <= xG, xC <= 360  
 C     This is a quick fix until this assumption is relaxed  
 C     and the interpolation weights are pre-computed.  
       DO bj = myByLo(myThid), myByHi(myThid)  
        DO bi = myBxLo(myThid), myBxHi(myThid)  
         DO J=1-Oly,sNy+Oly  
          DO I=1-Olx,sNx+Olx  
           IF ( xG(I,J,bi,bj) .LT. 0 )  
      &           xG(I,J,bi,bj) = xG(I,J,bi,bj) + 360  
           IF ( xC(I,J,bi,bj) .LT. 0 )  
      &           xC(I,J,bi,bj) = xC(I,J,bi,bj) + 360  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
 #endif /* ALLOW_EXF */  
   
79  C--   Write certain grid data to files (useful for creating netCDF  C--   Write certain grid data to files (useful for creating netCDF
80  C     and general post-analysis)  C     and general post-analysis)
81  C     This I/O is now done in write_grid.F        CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)
82  c     CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)        CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)
83  c     CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)        CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)
84  c     CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)        CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)
85  c     CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)        CALL WRITE_FLD_XY_RS( 'RAC',' ',rA,0,myThid)
86  c     CALL WRITE_FLD_XY_RS( 'RAC',' ',rA,0,myThid)        CALL WRITE_FLD_XY_RS( 'RAW',' ',rAw,0,myThid)
87  c     CALL WRITE_FLD_XY_RS( 'RAW',' ',rAw,0,myThid)        CALL WRITE_FLD_XY_RS( 'RAS',' ',rAs,0,myThid)
88  c     CALL WRITE_FLD_XY_RS( 'RAS',' ',rAs,0,myThid)        CALL WRITE_FLD_XY_RS( 'DXG',' ',DXG,0,myThid)
89  c     CALL WRITE_FLD_XY_RS( 'DXG',' ',DXG,0,myThid)        CALL WRITE_FLD_XY_RS( 'DYG',' ',DYG,0,myThid)
90  c     CALL WRITE_FLD_XY_RS( 'DYG',' ',DYG,0,myThid)        CALL WRITE_FLD_XY_RS( 'DXC',' ',DXC,0,myThid)
91  c     CALL WRITE_FLD_XY_RS( 'DXC',' ',DXC,0,myThid)        CALL WRITE_FLD_XY_RS( 'DYC',' ',DYC,0,myThid)
 c     CALL WRITE_FLD_XY_RS( 'DYC',' ',DYC,0,myThid)  
92    
93  C--   Print out statistics of each horizontal grid array (helps when debugging)  C--   Print out statistics of each horizontal grid array (helps when debugging)
94    
95  #ifdef ALLOW_MONITOR  #ifndef EXCLUDE_MONITOR
96        CALL MON_PRINTSTATS_RS(1,XC,'XC',myThid)        CALL MON_PRINTSTATS_RS(1,XC,'XC',myThid)
97        CALL MON_PRINTSTATS_RS(1,XG,'XG',myThid)        CALL MON_PRINTSTATS_RS(1,XG,'XG',myThid)
98        CALL MON_PRINTSTATS_RS(1,DXC,'DXC',myThid)        CALL MON_PRINTSTATS_RS(1,DXC,'DXC',myThid)
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22