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

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

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

Revision 1.3 - (hide annotations) (download)
Sun Feb 10 00:35:30 2002 UTC (22 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint46l_post, checkpoint46g_pre, checkpoint46f_post, checkpoint44f_post, checkpoint46b_post, checkpoint46l_pre, chkpt44d_post, checkpoint44e_pre, checkpoint47a_post, checkpoint46d_pre, checkpoint45d_post, checkpoint46j_pre, checkpoint44h_pre, checkpoint46a_post, checkpoint46j_post, checkpoint46k_post, checkpoint45a_post, checkpoint44e_post, checkpoint44g_post, checkpoint46e_pre, checkpoint45b_post, checkpoint46b_pre, release1_final_v1, checkpoint46c_pre, checkpoint46, checkpoint47b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint45c_post, checkpoint44h_post, checkpoint46g_post, checkpoint46i_post, checkpoint46c_post, checkpoint46e_post, checkpoint47, checkpoint45, checkpoint46h_post, checkpoint44f_pre, checkpoint46d_post
Branch point for: release1_final, release1
Changes since 1.2: +84 -29 lines 
add parameters to use real fresh water flux (replace USE_NATURALBCS option)
  with non-linear free-surface.
1 jmc 1.3 C $Header: /u/gcmpack/models/MITgcmUV/model/src/config_check.F,v 1.2 2002/01/31 15:02:18 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "CPP_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: CONFIG_CHECK
8     C !INTERFACE:
9     SUBROUTINE CONFIG_CHECK( myThid )
10     C !DESCRIPTION: \bv
11     C *=========================================================*
12     C | SUBROUTINE CONFIG_CHECK
13     C | o Check model parameter settings.
14     C *=========================================================*
15     C | This routine help to prevent the use of parameters
16     C | that are not compatible with the model configuration.
17     C *=========================================================*
18     C \ev
19    
20     C !USES:
21     IMPLICIT NONE
22     C === Global variables ===
23     #include "SIZE.h"
24     #include "EEPARAMS.h"
25     #include "PARAMS.h"
26     c #include "GRID.h"
27    
28     C !INPUT/OUTPUT PARAMETERS:
29     C === Routine arguments ===
30     C myThid - Number of this instances of CONFIG_CHECK
31     INTEGER myThid
32     CEndOfInterface
33    
34     C !LOCAL VARIABLES:
35     C == Local variables ==
36     C msgBuf :: Informational/error meesage buffer
37     CHARACTER*(MAX_LEN_MBUF) msgBuf
38     CEOP
39    
40     C- check that CPP option is "defined" when running-flag parameter is on:
41    
42     #ifndef ALLOW_NONHYDROSTATIC
43     IF (nonHydrostatic) THEN
44     WRITE(msgBuf,'(A)')
45     & 'CONFIG_CHECK: #undef ALLOW_NONHYDROSTATIC and'
46 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
47 jmc 1.1 WRITE(msgBuf,'(A)')
48     & 'CONFIG_CHECK: nonHydrostatic is TRUE'
49 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
50 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
51     ENDIF
52     #endif
53    
54     #ifndef EXACT_CONSERV
55     IF (exactConserv) THEN
56     WRITE(msgBuf,'(A)')
57     & 'CONFIG_CHECK: #undef EXACT_CONSERV and'
58 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
59 jmc 1.1 WRITE(msgBuf,'(A)')
60     & 'CONFIG_CHECK: exactConserv is TRUE'
61 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
62 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
63     ENDIF
64     #endif
65    
66     #ifndef NONLIN_FRSURF
67     IF (nonlinFreeSurf.NE.0) THEN
68     WRITE(msgBuf,'(A)')
69     & 'CONFIG_CHECK: #undef NONLIN_FRSURF and'
70 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
71 jmc 1.1 WRITE(msgBuf,'(A)')
72     & 'CONFIG_CHECK: nonlinFreeSurf is non-zero'
73 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
74 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
75     ENDIF
76     #endif
77    
78     #ifdef USE_NATURAL_BCS
79     WRITE(msgBuf,'(A)')
80 jmc 1.3 & 'CONFIG_CHECK: USE_NATURAL_BCS option has been replaced'
81     CALL PRINT_ERROR( msgBuf , myThid)
82 jmc 1.1 WRITE(msgBuf,'(A)')
83 jmc 1.3 & 'CONFIG_CHECK: by useRealFreshWaterFlux=TRUE in data file'
84     CALL PRINT_ERROR( msgBuf , myThid)
85 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
86 jmc 1.3 #endif
87    
88     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
89    
90     C- check parameter consistency :
91    
92     IF ( rigidLid .AND. implicitFreeSurface ) THEN
93     WRITE(msgBuf,'(A,A)')
94     & 'CONFIG_CHECK: Cannot select both implicitFreeSurface',
95     & ' and rigidLid.'
96     CALL PRINT_ERROR( msgBuf , myThid)
97     STOP 'ABNORMAL END: S/R CONFIG_CHECK'
98     ENDIF
99    
100     IF (rigidLid .AND. exactConserv) THEN
101 jmc 1.1 WRITE(msgBuf,'(A)')
102 jmc 1.3 & 'CONFIG_CHECK: exactConserv not compatible with'
103     CALL PRINT_ERROR( msgBuf , myThid)
104 jmc 1.1 WRITE(msgBuf,'(A)')
105 jmc 1.3 & 'CONFIG_CHECK: rigidLid (meaningless in that case)'
106     CALL PRINT_ERROR( msgBuf , myThid)
107 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
108     ENDIF
109    
110 jmc 1.3 IF (rigidLid .AND. useRealFreshWaterFlux) THEN
111 jmc 1.1 WRITE(msgBuf,'(A)')
112 jmc 1.3 & 'CONFIG_CHECK: useRealFreshWaterFlux not compatible with'
113     CALL PRINT_ERROR( msgBuf , myThid)
114 jmc 1.1 WRITE(msgBuf,'(A)')
115     & 'CONFIG_CHECK: rigidLid (meaningless in that case)'
116 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
117     STOP 'ABNORMAL END: S/R CONFIG_CHECK'
118     ENDIF
119    
120     IF ( (implicSurfPress.NE.1. .OR. implicDiv2DFlow.NE.1.)
121     & .AND. nonHydrostatic ) THEN
122     WRITE(msgBuf,'(A,A)') 'CONFIG_CHECK: nonHydrostatic',
123     & ' NOT SAFE with non-fully implicit Barotropic solver'
124     CALL PRINT_ERROR( msgBuf , myThid)
125     WRITE(msgBuf,'(A,A)') 'CONFIG_CHECK: To by-pass this',
126     & 'STOP, comment this test and re-compile config_check'
127     CALL PRINT_ERROR( msgBuf , myThid)
128 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
129     ENDIF
130    
131     IF (nonlinFreeSurf.NE.0 .AND. .NOT.exactConserv) THEN
132     WRITE(msgBuf,'(A)')
133     & 'CONFIG_CHECK: nonlinFreeSurf cannot be used'
134 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
135 jmc 1.1 WRITE(msgBuf,'(A)')
136     & 'CONFIG_CHECK: without exactConserv'
137 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
138 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
139     ENDIF
140    
141     C- note : not implemented in Release1_beta1 but it's done now (since 01-30-02)
142 jmc 1.2 c IF (nonlinFreeSurf.NE.0 .AND. useOBCS ) THEN
143     c WRITE(msgBuf,'(A)')
144     c & 'CONFIG_CHECK: nonlinFreeSurf not yet implemented'
145     c CALL PRINT_ERROR( msgBuf , 1)
146     c WRITE(msgBuf,'(A)')
147     c & 'CONFIG_CHECK: in OBC package'
148     c CALL PRINT_ERROR( msgBuf , 1)
149     c STOP 'ABNORMAL END: S/R CONFIG_CHECK'
150     c ENDIF
151 jmc 1.1
152     IF (nonlinFreeSurf.NE.0 .AND. nonHydrostatic) THEN
153     WRITE(msgBuf,'(A)')
154     & 'CONFIG_CHECK: nonlinFreeSurf not yet implemented'
155 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
156 jmc 1.1 WRITE(msgBuf,'(A)')
157     & 'CONFIG_CHECK: in nonHydrostatic code'
158 jmc 1.3 CALL PRINT_ERROR( msgBuf , myThid)
159 jmc 1.1 STOP 'ABNORMAL END: S/R CONFIG_CHECK'
160     ENDIF
161 jmc 1.3
162     IF (nonlinFreeSurf.NE.0 .AND. deltaTmom.NE.deltaTtracer) THEN
163     WRITE(msgBuf,'(A)')
164     & 'CONFIG_CHECK: WARNING: nonlinFreeSurf might cause'
165     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
166     & SQUEEZE_RIGHT , 1)
167     WRITE(msgBuf,'(A)')
168     & 'CONFIG_CHECK: problems with different time-steps'
169     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
170     & SQUEEZE_RIGHT , 1)
171     ENDIF
172    
173     IF (useRealFreshWaterFlux .AND. exactConserv
174     & .AND.startTime.NE.0. .AND. implicSurfPress.EQ.0. _d 0) THEN
175     WRITE(msgBuf,'(A)')
176     & 'CONFIG_CHECK: RealFreshWaterFlux+implicSurfP=0+exactConserv:'
177     CALL PRINT_ERROR( msgBuf , myThid)
178     WRITE(msgBuf,'(A)')
179     & 'CONFIG_CHECK: restart not implemented in this config'
180     CALL PRINT_ERROR( msgBuf , myThid)
181     STOP 'ABNORMAL END: S/R CONFIG_CHECK'
182     ENDIF
183    
184     #ifndef NONLIN_FRSURF
185     IF (useRealFreshWaterFlux .AND. exactConserv
186     & .AND. implicSurfPress.NE.1. _d 0 ) THEN
187     WRITE(msgBuf,'(A)')
188     & 'CONFIG_CHECK: Pb with restart in this config'
189     CALL PRINT_ERROR( msgBuf , myThid)
190     WRITE(msgBuf,'(A)')
191     & 'CONFIG_CHECK: ==> use #define NONLIN_FRSURF to fix it'
192     CALL PRINT_ERROR( msgBuf , myThid)
193     STOP 'ABNORMAL END: S/R CONFIG_CHECK'
194     ENDIF
195     #endif
196 jmc 1.1
197     WRITE(msgBuf,'(A)') 'CONFIG_CHECK: OK'
198     CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,
199     & SQUEEZE_RIGHT,myThid)
200    
201     RETURN
202     END
  ViewVC Help
Powered by ViewVC 1.1.22