/[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.9 by adcroft, Mon Jun 4 19:46:22 2001 UTC revision 1.15 by dimitri, Wed Dec 10 19:37:25 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CStartOfInterface  CBOP
8    C     !ROUTINE: INI_GRID
9    C     !INTERFACE:
10        SUBROUTINE INI_GRID( myThid )        SUBROUTINE INI_GRID( myThid )
11  C     /==========================================================\  C     !DESCRIPTION: \bv
12  C     | SUBROUTINE INI_GRID                                      |  C     *==========================================================*
13  C     | o Initialise model grid                                  |  C     | SUBROUTINE INI_GRID                                      
14  C     |==========================================================|  C     | o Initialise model grid                                  
15  C     | These arrays are used throughout the code in evaluating  |  C     *==========================================================*
16  C     | gradients, integrals and spatial avarages. This routine  |  C     | These arrays are used throughout the code in evaluating  
17  C     | is called separately by each thread and initialise only  |  C     | gradients, integrals and spatial avarages. This routine  
18  C     | the region of the domain it is "responsible" for.        |  C     | is called separately by each thread and initialise only  
19  C     | Notes:                                                   |  C     | the region of the domain it is "responsible" for.        
20  C     | Two examples are shown in this code. One illustrates the |  C     | Notes:                                                    
21  C     | initialisation of a cartesian grid. The other shows the  |  C     | Two examples are shown in this code. One illustrates the  
22  C     | inialisation of a spherical polar grid. Other orthonormal|  C     | initialisation of a cartesian grid. The other shows the  
23  C     | grids can be fitted into this design. In this case       |  C     | inialisation of a spherical polar grid. Other orthonormal
24  C     | custom metric terms also need adding to account for the  |  C     | grids can be fitted into this design. In this case        
25  C     | projections of velocity vectors onto these grids.        |  C     | custom metric terms also need adding to account for the  
26  C     | The structure used here also makes it possible to        |  C     | projections of velocity vectors onto these grids.        
27  C     | implement less regular grid mappings. In particular      |  C     | The structure used here also makes it possible to        
28  C     | o Schemes which leave out blocks of the domain that are  |  C     | implement less regular grid mappings. In particular      
29  C     |   all land could be supported.                           |  C     | o Schemes which leave out blocks of the domain that are  
30  C     | o Multi-level schemes such as icosohedral or cubic       |  C     |   all land could be supported.                            
31  C     |   grid projectedions onto a sphere can also be fitted    |  C     | o Multi-level schemes such as icosohedral or cubic        
32  C     |   within the strategy we use.                            |  C     |   grid projectedions onto a sphere can also be fitted    
33  C     |   Both of the above also require modifying the support   |  C     |   within the strategy we use.                            
34  C     |   routines that map computational blocks to simulation   |  C     |   Both of the above also require modifying the support    
35  C     |   domain blocks.                                         |  C     |   routines that map computational blocks to simulation    
36  C     \==========================================================/  C     |   domain blocks.                                          
37        IMPLICIT NONE  C     *==========================================================*
38    C     \ev
39    
40    C     !USES:
41          IMPLICIT NONE
42  C     === Global variables ===  C     === Global variables ===
43  #include "SIZE.h"  #include "SIZE.h"
44  #include "EEPARAMS.h"  #include "EEPARAMS.h"
45  #include "PARAMS.h"  #include "PARAMS.h"
46  #include "GRID.h"  #include "GRID.h"
47    
48    C     !INPUT/OUTPUT PARAMETERS:
49  C     == Routine arguments ==  C     == Routine arguments ==
50  C     myThid -  Number of this instance of INI_GRID  C     myThid -  Number of this instance of INI_GRID
51        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
52    
53    C     !LOCAL VARIABLES:
54  C     == Local variables ==  C     == Local variables ==
55  C     msgBuf - Used for informational I/O.  C     msgBuf - Used for informational I/O.
56        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
57    #ifdef ALLOW_EXF
58          INTEGER i, j, bi, bj
59    #endif /* ALLOW_EXF */
60    CEOP
61    
62  C--   Set up vertical grid and coordinate system  C--   Set up vertical grid and coordinate system
63        CALL INI_VERTICAL_GRID( myThid )        CALL INI_VERTICAL_GRID( myThid )
# Line 67  C--   Set up horizontal grid and coordin Line 78  C--   Set up horizontal grid and coordin
78         _END_MASTER(myThid)         _END_MASTER(myThid)
79        ENDIF        ENDIF
80    
81    #ifdef ALLOW_EXF
82    C--   exf_interp assumes that 0 <= xG, xC <= 360
83    C     This is a quick fix until this assumption is relaxed
84    C     and the interpolation weights are pre-computed.
85          DO bj = myByLo(myThid), myByHi(myThid)
86           DO bi = myBxLo(myThid), myBxHi(myThid)
87            DO J=1-Oly,sNy+Oly
88             DO I=1-Olx,sNx+Olx
89              IF ( xG(I,J,bi,bj) .LT. 0 )
90         &           xG(I,J,bi,bj) = xG(I,J,bi,bj) + 360
91              IF ( xC(I,J,bi,bj) .LT. 0 )
92         &           xC(I,J,bi,bj) = xC(I,J,bi,bj) + 360
93             ENDDO
94            ENDDO
95           ENDDO
96          ENDDO
97    #endif /* ALLOW_EXF */
98    
99  C--   Write certain grid data to files (useful for creating netCDF  C--   Write certain grid data to files (useful for creating netCDF
100  C     and general post-analysis)  C     and general post-analysis)
101        CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)        CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)
102        CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)        CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)
103        CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)        CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)
104        CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)        CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)
105        CALL WRITE_FLD_XY_RS( 'AC',' ',rA,0,myThid)        CALL WRITE_FLD_XY_RS( 'RAC',' ',rA,0,myThid)
106          CALL WRITE_FLD_XY_RS( 'RAW',' ',rAw,0,myThid)
107          CALL WRITE_FLD_XY_RS( 'RAS',' ',rAs,0,myThid)
108          CALL WRITE_FLD_XY_RS( 'DXG',' ',DXG,0,myThid)
109          CALL WRITE_FLD_XY_RS( 'DYG',' ',DYG,0,myThid)
110          CALL WRITE_FLD_XY_RS( 'DXC',' ',DXC,0,myThid)
111          CALL WRITE_FLD_XY_RS( 'DYC',' ',DYC,0,myThid)
112    
113  C--   Print out statistics of each horizontal grid array (helps when debugging)  C--   Print out statistics of each horizontal grid array (helps when debugging)
114    
115    #ifdef ALLOW_MONITOR
116        CALL MON_PRINTSTATS_RS(1,XC,'XC',myThid)        CALL MON_PRINTSTATS_RS(1,XC,'XC',myThid)
117        CALL MON_PRINTSTATS_RS(1,XG,'XG',myThid)        CALL MON_PRINTSTATS_RS(1,XG,'XG',myThid)
118        CALL MON_PRINTSTATS_RS(1,DXC,'DXC',myThid)        CALL MON_PRINTSTATS_RS(1,DXC,'DXC',myThid)
# Line 92  C--   Print out statistics of each horiz Line 129  C--   Print out statistics of each horiz
129        CALL MON_PRINTSTATS_RS(1,RAW,'RAW',myThid)        CALL MON_PRINTSTATS_RS(1,RAW,'RAW',myThid)
130        CALL MON_PRINTSTATS_RS(1,RAS,'RAS',myThid)        CALL MON_PRINTSTATS_RS(1,RAS,'RAS',myThid)
131        CALL MON_PRINTSTATS_RS(1,RAZ,'RAZ',myThid)        CALL MON_PRINTSTATS_RS(1,RAZ,'RAZ',myThid)
132    #endif
133    
134        RETURN        RETURN
135        END        END
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22