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

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

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

revision 1.1 by jmc, Thu Jan 31 14:52:28 2002 UTC revision 1.9 by jmc, Thu Apr 17 13:16:23 2003 UTC
# Line 39  CEOP Line 39  CEOP
39    
40  C-  check that CPP option is "defined" when running-flag parameter is on:  C-  check that CPP option is "defined" when running-flag parameter is on:
41    
42    #ifndef INCLUDE_CD_CODE
43          IF (useCDscheme) THEN
44            WRITE(msgBuf,'(A)')
45         &   'CONFIG_CHECK: #undef INCLUDE_CD_CODE and'
46            CALL PRINT_ERROR( msgBuf , myThid)
47            WRITE(msgBuf,'(A)')
48         &   'CONFIG_CHECK: useCDscheme is TRUE'
49            CALL PRINT_ERROR( msgBuf , myThid)
50            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
51          ENDIF
52    #endif
53    
54  #ifndef ALLOW_NONHYDROSTATIC  #ifndef ALLOW_NONHYDROSTATIC
55        IF (nonHydrostatic) THEN        IF (nonHydrostatic) THEN
56          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
57       &   'CONFIG_CHECK: #undef ALLOW_NONHYDROSTATIC and'       &   'CONFIG_CHECK: #undef ALLOW_NONHYDROSTATIC and'
58          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
59          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
60       &   'CONFIG_CHECK: nonHydrostatic is TRUE'       &   'CONFIG_CHECK: nonHydrostatic is TRUE'
61          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
62          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
63        ENDIF        ENDIF
64  #endif  #endif
# Line 55  C-  check that CPP option is "defined" w Line 67  C-  check that CPP option is "defined" w
67        IF (exactConserv) THEN        IF (exactConserv) THEN
68          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
69       &   'CONFIG_CHECK: #undef EXACT_CONSERV and'       &   'CONFIG_CHECK: #undef EXACT_CONSERV and'
70          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
71          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
72       &   'CONFIG_CHECK: exactConserv is TRUE'       &   'CONFIG_CHECK: exactConserv is TRUE'
73          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
74          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
75        ENDIF        ENDIF
76  #endif  #endif
# Line 67  C-  check that CPP option is "defined" w Line 79  C-  check that CPP option is "defined" w
79        IF (nonlinFreeSurf.NE.0) THEN        IF (nonlinFreeSurf.NE.0) THEN
80          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
81       &   'CONFIG_CHECK: #undef NONLIN_FRSURF and'       &   'CONFIG_CHECK: #undef NONLIN_FRSURF and'
82          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
83          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
84       &   'CONFIG_CHECK: nonlinFreeSurf is non-zero'       &   'CONFIG_CHECK: nonlinFreeSurf is non-zero'
85          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
86          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
87        ENDIF        ENDIF
88  #endif  #endif
89    
90  C-  check parameter consistency :  #ifndef NONLIN_FRSURF
91          IF (select_rStar .NE. 0) THEN
92            WRITE(msgBuf,'(A)')
93         &   'CONFIG_CHECK: rStar is part of NonLin-FS '
94            CALL PRINT_ERROR( msgBuf, myThid)                      
95            WRITE(msgBuf,'(A)')
96         &   'CONFIG_CHECK: ==> use #define NONLIN_FRSURF to use it'
97            CALL PRINT_ERROR( msgBuf, myThid)                      
98            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
99          ENDIF
100    #endif /* NONLIN_FRSURF */
101    
102  #ifdef USE_NATURAL_BCS  #ifdef USE_NATURAL_BCS
       IF (rigidLid) THEN  
103          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
104       &   'CONFIG_CHECK: #define USE_NATURAL_BCS and'       &   'CONFIG_CHECK: USE_NATURAL_BCS option has been replaced'
105          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
106          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
107       &   'CONFIG_CHECK: rigidLid are not compatible'       &   'CONFIG_CHECK: by useRealFreshWaterFlux=TRUE in data file'
108          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
109            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
110    #endif
111    
112    C     o If pLoadFile is set, then we should make sure the corresponing
113    C       code is being compiled
114    #ifndef ATMOSPHERIC_LOADING
115          IF (pLoadFile.NE.' ') THEN
116            WRITE(msgBuf,'(A)')
117         &  'CONFIG_CHECK: pLoadFile is set but you have not'
118            CALL PRINT_ERROR( msgBuf , myThid)
119            WRITE(msgBuf,'(A)')
120         &  'compiled the model with the pressure loading code.'
121            CALL PRINT_ERROR( msgBuf , myThid)
122            WRITE(msgBuf,'(A,A)')
123         &  'Re-compile with:  #define ATMOSPHERIC_LOADING',
124         &  '              or  -DATMOSPHERIC_LOADING'
125            CALL PRINT_ERROR( msgBuf , myThid)
126          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
127        ENDIF        ENDIF
128        IF (exactConserv) THEN  #endif
129          WRITE(msgBuf,'(A)')  
130       &   'CONFIG_CHECK: #define USE_NATURAL_BCS with'  C     o If taveFreq is finite, then we must make sure the diagnostics
131    C       code is being compiled
132    #ifndef ALLOW_TIMEAVE
133          IF (taveFreq.NE.0.) THEN
134            WRITE(msgBuf,'(A)')
135         &  'CONFIG_CHECK: taveFreq <> 0  but you have'
136          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , 1)
137          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
138       &   'CONFIG_CHECK: exactConserv not yet implemented'       &  'not compiled the model with the diagnostics routines.'
139            CALL PRINT_ERROR( msgBuf , 1)
140            WRITE(msgBuf,'(A,A)')
141         &  'Re-compile with:  #define ALLOW_TIMEAVE',
142         &  '              or  -DALLOW_TIMEAVE'
143          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , 1)
144          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
145        ENDIF        ENDIF
146  #endif  #endif
147    
148    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
149    
150    C-  check parameter consistency :
151    
152          IF ( viscA4.NE.0. .AND. (Olx.LT.3 .OR. Oly.LT.3)) THEN
153            WRITE(msgBuf,'(A,A)')
154         &  'CONFIG_CHECK: cannot use Biharmonic Visc. (viscA4) with',
155         &  ' overlap (Olx,Oly) smaller than 3'
156            CALL PRINT_ERROR( msgBuf , myThid)
157            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
158          ENDIF                
159    
160          IF ( rigidLid .AND. implicitFreeSurface ) THEN
161            WRITE(msgBuf,'(A,A)')
162         &  'CONFIG_CHECK: Cannot select both implicitFreeSurface',
163         &  ' and rigidLid.'
164            CALL PRINT_ERROR( msgBuf , myThid)
165            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
166          ENDIF                
167    
168        IF (rigidLid .AND. exactConserv) THEN        IF (rigidLid .AND. exactConserv) THEN
169          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
170       &   'CONFIG_CHECK: exactConserv not compatible with'       &   'CONFIG_CHECK: exactConserv not compatible with'
171          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
172          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
173       &   'CONFIG_CHECK: rigidLid (meaningless in that case)'       &   'CONFIG_CHECK: rigidLid (meaningless in that case)'
174          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
175            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
176          ENDIF
177    
178          IF (rigidLid .AND. useRealFreshWaterFlux) THEN
179            WRITE(msgBuf,'(A)')
180         &   'CONFIG_CHECK: useRealFreshWaterFlux not compatible with'
181            CALL PRINT_ERROR( msgBuf , myThid)
182            WRITE(msgBuf,'(A)')
183         &   'CONFIG_CHECK: rigidLid (meaningless in that case)'
184            CALL PRINT_ERROR( msgBuf , myThid)
185            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
186          ENDIF
187    
188          IF ( (implicSurfPress.NE.1. .OR. implicDiv2DFlow.NE.1.)
189         &    .AND. nonHydrostatic ) THEN
190            WRITE(msgBuf,'(A,A)') 'CONFIG_CHECK: nonHydrostatic',
191         & ' NOT SAFE with non-fully implicit Barotropic solver'
192            CALL PRINT_ERROR( msgBuf , myThid)
193            WRITE(msgBuf,'(A,A)') 'CONFIG_CHECK: To by-pass this',
194         &    'STOP, comment this test and re-compile config_check'
195            CALL PRINT_ERROR( msgBuf , myThid)
196          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
197        ENDIF        ENDIF
198    
199        IF (nonlinFreeSurf.NE.0 .AND. .NOT.exactConserv) THEN        IF (nonlinFreeSurf.NE.0 .AND. .NOT.exactConserv) THEN
200          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
201       &   'CONFIG_CHECK: nonlinFreeSurf cannot be used'       &   'CONFIG_CHECK: nonlinFreeSurf cannot be used'
202          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
203          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
204       &   'CONFIG_CHECK: without exactConserv'       &   'CONFIG_CHECK: without exactConserv'
205          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
206          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
207        ENDIF        ENDIF
208    
209  C- note : not implemented in Release1_beta1 but it's done now (since 01-30-02)        IF (select_rStar.NE.0 .AND. .NOT.exactConserv) THEN
       IF (nonlinFreeSurf.NE.0 .AND. useOBCS ) THEN  
210          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
211       &   'CONFIG_CHECK: nonlinFreeSurf not yet implemented'       &   'CONFIG_CHECK: r* Coordinate cannot be used'
212          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
213          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
214       &   'CONFIG_CHECK: in OBC package'       &   'CONFIG_CHECK: without exactConserv'
215          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
216          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
217        ENDIF        ENDIF
218    
219    C- note : not implemented in checkpoint48b but it's done now (since 01-28-03)
220    c     IF (select_rStar.GT.0 .AND. useOBCS ) THEN
221    c       WRITE(msgBuf,'(A)')
222    c    &   'CONFIG_CHECK: r* Coordinate not yet implemented'
223    c       CALL PRINT_ERROR( msgBuf , 1)
224    c       WRITE(msgBuf,'(A)')
225    c    &   'CONFIG_CHECK: in OBC package'
226    c       CALL PRINT_ERROR( msgBuf , 1)
227    c       STOP 'ABNORMAL END: S/R CONFIG_CHECK'
228    c     ENDIF
229    
230        IF (nonlinFreeSurf.NE.0 .AND. nonHydrostatic) THEN        IF (nonlinFreeSurf.NE.0 .AND. nonHydrostatic) THEN
231          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
232       &   'CONFIG_CHECK: nonlinFreeSurf not yet implemented'       &   'CONFIG_CHECK: nonlinFreeSurf not yet implemented'
233          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
234          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
235       &   'CONFIG_CHECK: in nonHydrostatic code'       &   'CONFIG_CHECK: in nonHydrostatic code'
236          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , myThid)
237            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
238          ENDIF
239    
240          IF (nonlinFreeSurf.NE.0.AND.deltaTfreesurf.NE.deltaTtracer) THEN
241            WRITE(msgBuf,'(A)')
242         &   'CONFIG_CHECK: WARNING: nonlinFreeSurf might cause problems'
243            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
244         &                    SQUEEZE_RIGHT , myThid)                      
245            WRITE(msgBuf,'(A)')
246         &   'CONFIG_CHECK: with different FreeSurf & Tracer time-steps'
247            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
248         &                    SQUEEZE_RIGHT , myThid)                      
249          ENDIF
250    
251          IF (useRealFreshWaterFlux .AND. exactConserv
252         &    .AND.startTime.NE.0. .AND. implicSurfPress.EQ.0. _d 0) THEN
253            WRITE(msgBuf,'(A)')
254         &   'CONFIG_CHECK: RealFreshWaterFlux+implicSurfP=0+exactConserv:'
255            CALL PRINT_ERROR( msgBuf , myThid)
256            WRITE(msgBuf,'(A)')
257         &   'CONFIG_CHECK: restart not implemented in this config'
258            CALL PRINT_ERROR( msgBuf , myThid)
259            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
260          ENDIF
261    
262    #ifdef NONLIN_FRSURF
263          IF (useRealFreshWaterFlux .AND. .NOT.exactConserv
264         &            .AND. buoyancyRelation.EQ.'OCEANICP' ) THEN
265            WRITE(msgBuf,'(A)')
266         &   'CONFIG_CHECK: RealFreshWaterFlux with OCEANICP'
267            CALL PRINT_ERROR( msgBuf , myThid)
268            WRITE(msgBuf,'(A)')
269         &   'CONFIG_CHECK: requires exactConserv=T'
270            CALL PRINT_ERROR( msgBuf , myThid)
271            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
272          ENDIF
273    #else
274          IF (useRealFreshWaterFlux .AND. exactConserv
275         &            .AND. implicSurfPress.NE.1. _d 0 ) THEN
276            WRITE(msgBuf,'(A)')
277         &   'CONFIG_CHECK: Pb with restart in this config'
278            CALL PRINT_ERROR( msgBuf , myThid)
279            WRITE(msgBuf,'(A)')
280         &   'CONFIG_CHECK: ==> use #define NONLIN_FRSURF to fix it'
281            CALL PRINT_ERROR( msgBuf , myThid)
282            STOP 'ABNORMAL END: S/R CONFIG_CHECK'
283          ENDIF
284    
285          IF (useRealFreshWaterFlux
286         &            .AND. buoyancyRelation.EQ.'OCEANICP' ) THEN
287            WRITE(msgBuf,'(A)')
288         &   'CONFIG_CHECK: E-P effects on wVel are not included'
289            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
290         &                    SQUEEZE_RIGHT , myThid)                      
291            WRITE(msgBuf,'(A)')
292         &   'CONFIG_CHECK: ==> use #define NONLIN_FRSURF to fix it'
293            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
294         &                    SQUEEZE_RIGHT , myThid)                      
295          ENDIF
296    #endif /* NONLIN_FRSURF */
297    
298          IF (.NOT.useCDscheme .AND. (tauCD.NE.0. .OR. rCD.NE.-1.) ) THEN
299    C- jmc: since useCDscheme is a new [04-13-03] flag (default=F),
300    C       put this WARNING to stress that even if CD-scheme parameters
301    C       (tauCD,rCD) are set, CD-scheme is not used without useCDscheme=T
302    C-    and STOP if using mom_fluxform (following Chris advise).
303    C- jmc: but ultimately, this block can/will be removed.
304           IF (.NOT.vectorInvariantMomentum.AND.momStepping) THEN
305            WRITE(msgBuf,'(A)')
306         &   'CONFIG_CHECK: CD-scheme is OFF but params(tauCD,rCD) are set'
307            CALL PRINT_ERROR( msgBuf , myThid)
308            WRITE(msgBuf,'(2A)')
309         &   'CONFIG_CHECK: to turn ON CD-scheme: => "useCDscheme=.TRUE."',
310         &   ' in "data", namelist PARM01'
311            CALL PRINT_ERROR( msgBuf , myThid)
312          STOP 'ABNORMAL END: S/R CONFIG_CHECK'          STOP 'ABNORMAL END: S/R CONFIG_CHECK'
313           ENDIF
314            WRITE(msgBuf,'(2A)') '**WARNNING** ',
315         &   'CONFIG_CHECK: CD-scheme is OFF but params(tauCD,rCD) are set'
316            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
317         &                    SQUEEZE_RIGHT , myThid)                      
318            WRITE(msgBuf,'(2A)')
319         &   'CONFIG_CHECK: to turn ON CD-scheme: => "useCDscheme=.TRUE."',
320         &   ' in "data", namelist PARM01'
321            CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
322         &                    SQUEEZE_RIGHT , myThid)                      
323        ENDIF        ENDIF
324    
325        WRITE(msgBuf,'(A)') 'CONFIG_CHECK: OK'        WRITE(msgBuf,'(A)') 'CONFIG_CHECK: OK'
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22