/[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.2 by cnh, Fri Apr 24 02:05:41 1998 UTC revision 1.18 by heimbach, Wed Jun 30 23:45:35 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CBOP
7    C     !ROUTINE: INI_GRID
8    C     !INTERFACE:
9        SUBROUTINE INI_GRID( myThid )        SUBROUTINE INI_GRID( myThid )
10  C     /==========================================================\  C     !DESCRIPTION: \bv
11  C     | SUBROUTINE INI_GRID                                      |  C     *==========================================================*
12  C     | o Initialise model grid                                  |  C     | SUBROUTINE INI_GRID                                      
13  C     |==========================================================|  C     | o Initialise model grid                                  
14  C     | These arrays are used throughout the code in evaluating  |  C     *==========================================================*
15  C     | gradients, integrals and spatial avarages. This routine  |  C     | These arrays are used throughout the code in evaluating  
16  C     | is called separately by each thread and initialise only  |  C     | gradients, integrals and spatial avarages. This routine  
17  C     | the region of the domain it is "responsible" for.        |  C     | is called separately by each thread and initialise only  
18  C     | Notes:                                                   |  C     | the region of the domain it is "responsible" for.        
19  C     | Two examples are shown in this code. One illustrates the |  C     | Notes:                                                    
20  C     | initialisation of a cartesian grid. The other shows the  |  C     | Two examples are shown in this code. One illustrates the  
21  C     | inialisation of a spherical polar grid. Other orthonormal|  C     | initialisation of a cartesian grid. The other shows the  
22  C     | grids can be fitted into this design. In this case       |  C     | inialisation of a spherical polar grid. Other orthonormal
23  C     | custom metric terms also need adding to account for the  |  C     | grids can be fitted into this design. In this case        
24  C     | projections of velocity vectors onto these grids.        |  C     | custom metric terms also need adding to account for the  
25  C     | The structure used here also makes it possible to        |  C     | projections of velocity vectors onto these grids.        
26  C     | implement less regular grid mappings. In particular      |  C     | The structure used here also makes it possible to        
27  C     | o Schemes which leave out blocks of the domain that are  |  C     | implement less regular grid mappings. In particular      
28  C     |   all land could be supported.                           |  C     | o Schemes which leave out blocks of the domain that are  
29  C     | o Multi-level schemes such as icosohedral or cubic       |  C     |   all land could be supported.                            
30  C     |   grid projectedions onto a sphere can also be fitted    |  C     | o Multi-level schemes such as icosohedral or cubic        
31  C     |   within the strategy we use.                            |  C     |   grid projectedions onto a sphere can also be fitted    
32  C     |   Both of the above also require modifying the support   |  C     |   within the strategy we use.                            
33  C     |   routines that map computational blocks to simulation   |  C     |   Both of the above also require modifying the support    
34  C     |   domain blocks.                                         |  C     |   routines that map computational blocks to simulation    
35  C     \==========================================================/  C     |   domain blocks.                                          
36    C     *==========================================================*
37    C     \ev
38    
39    C     !USES:
40          IMPLICIT NONE
41  C     === Global variables ===  C     === Global variables ===
42  #include "SIZE.h"  #include "SIZE.h"
43  #include "EEPARAMS.h"  #include "EEPARAMS.h"
44  #include "PARAMS.h"  #include "PARAMS.h"
45  #include "GRID.h"  #include "GRID.h"
46    
47    C     !INPUT/OUTPUT PARAMETERS:
48  C     == Routine arguments ==  C     == Routine arguments ==
49  C     myThid -  Number of this instance of INI_GRID  C     myThid -  Number of this instance of INI_GRID
50        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
51    
52    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
56    CEOP
57    
58  C--   Set up grid coordinate system  C--   Set up vertical grid and coordinate system
59          CALL INI_VERTICAL_GRID( myThid )
60    
61    C--   Set up horizontal grid and coordinate system
62        IF ( usingCartesianGrid ) THEN        IF ( usingCartesianGrid ) THEN
63         CALL INI_CARTESIAN_GRID( myThid )         CALL INI_CARTESIAN_GRID( myThid )
64        ELSEIF ( usingSphericalPolarGrid ) THEN        ELSEIF ( usingSphericalPolarGrid ) THEN
65         CALL INI_SPHERICAL_POLAR_GRID( myThid )         CALL INI_SPHERICAL_POLAR_GRID( myThid )
66          ELSEIF ( usingCurvilinearGrid ) THEN
67           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 59  C--   Set up grid coordinate system Line 75  C--   Set up grid coordinate system
75           STOP 'ABNORMAL END: S/R INI_GRID'           STOP 'ABNORMAL END: S/R INI_GRID'
76         _END_MASTER(myThid)         _END_MASTER(myThid)
77        ENDIF        ENDIF
78  C  
79    C--   Write certain grid data to files (useful for creating netCDF
80    C     and general post-analysis)
81          CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid)
82          CALL WRITE_FLD_XY_RS( 'YC',' ',YC,0,myThid)
83          CALL WRITE_FLD_XY_RS( 'XG',' ',XG,0,myThid)
84          CALL WRITE_FLD_XY_RS( 'YG',' ',YG,0,myThid)
85          CALL WRITE_FLD_XY_RS( 'RAC',' ',rA,0,myThid)
86          CALL WRITE_FLD_XY_RS( 'RAW',' ',rAw,0,myThid)
87          CALL WRITE_FLD_XY_RS( 'RAS',' ',rAs,0,myThid)
88          CALL WRITE_FLD_XY_RS( 'DXG',' ',DXG,0,myThid)
89          CALL WRITE_FLD_XY_RS( 'DYG',' ',DYG,0,myThid)
90          CALL WRITE_FLD_XY_RS( 'DXC',' ',DXC,0,myThid)
91          CALL WRITE_FLD_XY_RS( 'DYC',' ',DYC,0,myThid)
92    
93    C--   Print out statistics of each horizontal grid array (helps when debugging)
94    
95    #ifdef ALLOW_MONITOR
96          CALL MON_PRINTSTATS_RS(1,XC,'XC',myThid)
97          CALL MON_PRINTSTATS_RS(1,XG,'XG',myThid)
98          CALL MON_PRINTSTATS_RS(1,DXC,'DXC',myThid)
99          CALL MON_PRINTSTATS_RS(1,DXF,'DXF',myThid)
100          CALL MON_PRINTSTATS_RS(1,DXG,'DXG',myThid)
101          CALL MON_PRINTSTATS_RS(1,DXV,'DXV',myThid)
102          CALL MON_PRINTSTATS_RS(1,YC,'YC',myThid)
103          CALL MON_PRINTSTATS_RS(1,YG,'YG',myThid)
104          CALL MON_PRINTSTATS_RS(1,DYC,'DYC',myThid)
105          CALL MON_PRINTSTATS_RS(1,DYF,'DYF',myThid)
106          CALL MON_PRINTSTATS_RS(1,DYG,'DYG',myThid)
107          CALL MON_PRINTSTATS_RS(1,DYU,'DYU',myThid)
108          CALL MON_PRINTSTATS_RS(1,RA,'RA',myThid)
109          CALL MON_PRINTSTATS_RS(1,RAW,'RAW',myThid)
110          CALL MON_PRINTSTATS_RS(1,RAS,'RAS',myThid)
111          CALL MON_PRINTSTATS_RS(1,RAZ,'RAZ',myThid)
112    #endif
113    
114        RETURN        RETURN
115        END        END
   
 C     $Id$  
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22