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

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

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

Revision 1.23 - (hide annotations) (download)
Wed Oct 13 04:37:37 2004 UTC (19 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint57g_pre, checkpoint57b_post, checkpoint57g_post, checkpoint56b_post, checkpoint57d_post, checkpoint57i_post, checkpoint57, checkpoint56, checkpoint55i_post, checkpoint57f_post, checkpoint57a_post, checkpoint57h_pre, checkpoint57h_post, checkpoint55g_post, checkpoint57c_post, checkpoint57c_pre, checkpoint55j_post, checkpoint55h_post, checkpoint57e_post, checkpoint55f_post, eckpoint57e_pre, checkpoint56a_post, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint56c_post, checkpoint57a_pre, checkpoint57k_post
Changes since 1.22: +4 -1 lines 
 o I am *sick* of moving these variables around.  But Jean-Michel has
   all but threatened a hissy fit if they aren't removed from PARAMS.h.
   So now here they are *back* in MNC_PARAMS.h where they were just a
   few days ago.
1 edhill 1.23 C $Header: /u/gcmpack/MITgcm/model/src/ini_grid.F,v 1.22 2004/10/10 06:08:47 edhill Exp $
2 adcroft 1.8 C $Name: $
3 cnh 1.1
4 jmc 1.19 #include "PACKAGES_CONFIG.h"
5 cnh 1.5 #include "CPP_OPTIONS.h"
6 cnh 1.1
7 cnh 1.10 CBOP
8     C !ROUTINE: INI_GRID
9 edhill 1.21
10 cnh 1.10 C !INTERFACE:
11 cnh 1.1 SUBROUTINE INI_GRID( myThid )
12 edhill 1.21 C !DESCRIPTION:
13     C These arrays are used throughout the code in evaluating gradients,
14     C integrals and spatial avarages. This routine is called separately
15     C by each thread and initializes only the region of the domain it is
16     C "responsible" for.
17 cnh 1.10
18     C !USES:
19 adcroft 1.6 IMPLICIT NONE
20 cnh 1.1 #include "SIZE.h"
21     #include "EEPARAMS.h"
22     #include "PARAMS.h"
23 edhill 1.23 #ifdef ALLOW_MNC
24     #include "MNC_PARAMS.h"
25     #endif
26 cnh 1.1 #include "GRID.h"
27 edhill 1.21 #ifdef ALLOW_MONITOR
28     #include "MONITOR.h"
29     #endif
30 cnh 1.1
31 cnh 1.10 C !INPUT/OUTPUT PARAMETERS:
32 cnh 1.1 INTEGER myThid
33 edhill 1.21 CEOP
34 cnh 1.1
35 cnh 1.10 C !LOCAL VARIABLES:
36 cnh 1.1 C msgBuf - Used for informational I/O.
37     CHARACTER*(MAX_LEN_MBUF) msgBuf
38 edhill 1.21 INTEGER i
39 cnh 1.1
40 edhill 1.21 C Set up vertical grid and coordinate system
41 adcroft 1.4 CALL INI_VERTICAL_GRID( myThid )
42    
43 edhill 1.21 C Two examples are shown in this code. One illustrates the
44     C initialization of a cartesian grid. The other shows the
45     C inialization of a spherical polar grid. Other orthonormal grids
46     C can be fitted into this design. In this case custom metric terms
47     C also need adding to account for the projections of velocity
48     C vectors onto these grids. The structure used here also makes it
49     C possible to implement less regular grid mappings. In particular:
50     C o Schemes which leave out blocks of the domain that are
51     C all land could be supported.
52     C o Multi-level schemes such as icosohedral or cubic
53     C grid projectedions onto a sphere can also be fitted
54     C within the strategy we use.
55     C Both of the above also require modifying the support
56     C routines that map computational blocks to simulation
57     C domain blocks.
58    
59     C Set up horizontal grid and coordinate system
60 cnh 1.1 IF ( usingCartesianGrid ) THEN
61 edhill 1.21 CALL INI_CARTESIAN_GRID( myThid )
62 cnh 1.1 ELSEIF ( usingSphericalPolarGrid ) THEN
63 edhill 1.21 CALL INI_SPHERICAL_POLAR_GRID( myThid )
64 adcroft 1.8 ELSEIF ( usingCurvilinearGrid ) THEN
65 edhill 1.21 CALL INI_CURVILINEAR_GRID( myThid )
66 afe 1.17 ELSEIF ( usingCylindricalGrid ) THEN
67 edhill 1.21 CALL INI_CYLINDER( myThid )
68 cnh 1.1 ELSE
69 edhill 1.21 _BEGIN_MASTER(myThid)
70     WRITE(msgBuf,'(2A)') 'S/R INI_GRID: ',
71     & 'No grid coordinate system has been selected'
72     CALL PRINT_ERROR( msgBuf , myThid)
73     STOP 'ABNORMAL END: S/R INI_GRID'
74     _END_MASTER(myThid)
75 cnh 1.1 ENDIF
76 dimitri 1.15
77 adcroft 1.20 #ifdef ALLOW_MONITOR
78 edhill 1.21 mon_write_stdout = .FALSE.
79     mon_write_mnc = .FALSE.
80     IF (monitor_stdio) THEN
81     mon_write_stdout = .TRUE.
82     ENDIF
83    
84     #ifdef ALLOW_MNC
85     IF (useMNC .AND. monitor_mnc) THEN
86     DO i = 1,MAX_LEN_MBUF
87     mon_fname(i:i) = ' '
88     ENDDO
89     mon_fname(1:12) = 'monitor_grid'
90     CALL MNC_CW_SET_UDIM(mon_fname, 1, myThid)
91     mon_write_mnc = .TRUE.
92     ENDIF
93     #endif /* ALLOW_MNC */
94    
95     C Print out statistics of each horizontal grid array (helps when
96     C debugging)
97 adcroft 1.9 CALL MON_PRINTSTATS_RS(1,XC,'XC',myThid)
98     CALL MON_PRINTSTATS_RS(1,XG,'XG',myThid)
99     CALL MON_PRINTSTATS_RS(1,DXC,'DXC',myThid)
100     CALL MON_PRINTSTATS_RS(1,DXF,'DXF',myThid)
101     CALL MON_PRINTSTATS_RS(1,DXG,'DXG',myThid)
102     CALL MON_PRINTSTATS_RS(1,DXV,'DXV',myThid)
103     CALL MON_PRINTSTATS_RS(1,YC,'YC',myThid)
104     CALL MON_PRINTSTATS_RS(1,YG,'YG',myThid)
105     CALL MON_PRINTSTATS_RS(1,DYC,'DYC',myThid)
106     CALL MON_PRINTSTATS_RS(1,DYF,'DYF',myThid)
107     CALL MON_PRINTSTATS_RS(1,DYG,'DYG',myThid)
108     CALL MON_PRINTSTATS_RS(1,DYU,'DYU',myThid)
109     CALL MON_PRINTSTATS_RS(1,RA,'RA',myThid)
110     CALL MON_PRINTSTATS_RS(1,RAW,'RAW',myThid)
111     CALL MON_PRINTSTATS_RS(1,RAS,'RAS',myThid)
112     CALL MON_PRINTSTATS_RS(1,RAZ,'RAZ',myThid)
113 edhill 1.21
114     mon_write_stdout = .FALSE.
115     mon_write_mnc = .FALSE.
116 heimbach 1.12 #endif
117 adcroft 1.8
118 cnh 1.1 RETURN
119     END
  ViewVC Help
Powered by ViewVC 1.1.22