/[MITgcm]/MITgcm/pkg/exf/exf_check.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_check.F

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

revision 1.6 by jmc, Mon May 14 19:34:57 2007 UTC revision 1.31 by jmc, Sun Feb 12 00:50:58 2017 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "EXF_OPTIONS.h"  #include "EXF_OPTIONS.h"
5    
6        subroutine exf_check( mythid )  C--  File exf_check.F: Routines to check EXF settings
7    C--   Contents
8    C--   o EXF_CHECK
9    C--   o EXF_CHECK_INTERP
10    
11    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
12    CBOP
13    C     !ROUTINE: EXF_CHECK
14    C     !INTERFACE:
15    
16          SUBROUTINE EXF_CHECK( myThid )
17    
18    C     !DESCRIPTION: \bv
19    C     *==========================================================*
20    C     | S/R EXF_CHECK
21    C     | o Check parameters and other package dependences
22    C     *==========================================================*
23    C     \ev
24    
25  c     ==================================================================  C     !USES:
26  c     SUBROUTINE exf_check        IMPLICIT NONE
 c     ==================================================================  
 c  
       implicit none  
   
 c     == global variables ==  
27    
28    C     == Global variables ===
29  #include "EEPARAMS.h"  #include "EEPARAMS.h"
30  #include "SIZE.h"  #include "SIZE.h"
31  #include "PARAMS.h"  #include "PARAMS.h"
 c#include "FFIELDS.h"  
 c#include "GRID.h"  
32    
33  #include "EXF_PARAM.h"  #include "EXF_PARAM.h"
34  #include "EXF_CONSTANTS.h"  #include "EXF_CONSTANTS.h"
 c#include "EXF_FIELDS.h"  
 c     == routine arguments ==  
35    
36  c     mythid - thread number for this instance of the routine.  C     !INPUT/OUTPUT PARAMETERS:
37    C     myThid   :: my Thread Id number
38          INTEGER myThid
39    
40    C     !LOCAL VARIABLES:
41    C     msgBuf   :: Informational/error message buffer
42          CHARACTER*(MAX_LEN_MBUF) msgBuf
43          INTEGER errCount
44    CEOP
45    
46          _BEGIN_MASTER(myThid)
47          errCount = 0
48    
49          WRITE(msgBuf,'(A)') 'EXF_CHECK: #define ALLOW_EXF'
50          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
51         &                    SQUEEZE_RIGHT, myThid )
52    
53    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
54    C     check for consistency
55          IF (.NOT.
56         &     (exf_iprec.EQ.precFloat32 .OR. exf_iprec.EQ.precFloat64)
57         &     ) THEN
58           WRITE(msgBuf,'(A)')
59         &      'S/R EXF_CHECK: value of exf_iprec not allowed'
60           CALL PRINT_ERROR( msgBuf, myThid )
61           errCount = errCount + 1
62          ENDIF
63    
64          IF ( repeatPeriod.LT.0. ) THEN
65           WRITE(msgBuf,'(A)')
66         &      'S/R EXF_CHECK: repeatPeriod must be positive'
67           CALL PRINT_ERROR( msgBuf, myThid )
68           errCount = errCount + 1
69          ENDIF
70    
71        integer mythid        IF ( useExfYearlyFields ) THEN
72           IF ( .NOT.useCAL ) THEN
73            WRITE(msgBuf,'(2A)') 'S/R EXF_CHECK: ',
74         &       'useExfYearlyFields requires to use pkg/cal (useCAL=T)'
75            CALL PRINT_ERROR( msgBuf, myThid )
76            errCount = errCount + 1
77           ENDIF
78           IF ( repeatPeriod.NE.0. ) THEN
79            WRITE(msgBuf,'(2A)') 'S/R EXF_CHECK: The use of ',
80         &       'useExfYearlyFields AND repeatPeriod is not implemented'
81            CALL PRINT_ERROR( msgBuf, myThid )
82            errCount = errCount + 1
83           ENDIF
84          ENDIF
85          IF ( useOBCS .AND. useOBCSYearlyFields ) THEN
86           IF ( .NOT.useCAL ) THEN
87            WRITE(msgBuf,'(2A)') 'S/R EXF_CHECK: ',
88         &       'useOBCSYearlyFields requires to use pkg/cal (useCAL=T)'
89            CALL PRINT_ERROR( msgBuf, myThid )
90            errCount = errCount + 1
91           ENDIF
92           IF ( repeatPeriod.NE.0. ) THEN
93            WRITE(msgBuf,'(2A)') 'S/R EXF_CHECK: ',
94         &       'useOBCSYearlyFields not implemented for repeatPeriod <> 0'
95            CALL PRINT_ERROR( msgBuf, myThid )
96            errCount = errCount + 1
97           ENDIF
98          ENDIF
99    
100  c     == local variables ==        IF ( useAtmWind ) THEN
101           IF ( ustressfile .NE. ' ' .OR. vstressfile .NE. ' ' ) THEN
102            WRITE(msgBuf,'(A)')
103         &   'S/R EXF_CHECK: use u,v_wind components but not wind-stress'
104            CALL PRINT_ERROR( msgBuf, myThid )
105            errCount = errCount + 1
106           ENDIF
107          ENDIF
108    
109  c     == end of interface ==        IF ( .NOT.useAtmWind ) THEN
110           IF ( uwindfile .NE. ' ' .OR. vwindfile .NE. ' ' ) THEN
111  c     check for consistency          WRITE(msgBuf,'(A)')
       if (.NOT.  
      &     (exf_iprec .EQ. 32 .OR. exf_iprec .EQ. 64)  
      &     ) then  
          stop 'stop in exf_readparms: value of exf_iprec not allowed'  
       else if (.NOT.  
      &        (exf_yftype .EQ. 'RS' .OR.  
      &        exf_yftype .EQ. 'RL')  
      &        ) then  
          stop 'stop in exf_readparms: value of exf_yftype not allowed'  
       end if  
   
       if ( useCubedSphereExchange ) then  
 cph         if ( uvecfile .NE. ' ' .and. vvecfile .NE. ' ' ) then  
 c     some restrictions that can be relaxed later on  
 cph            if ( uvecstartdate .ne. vvecstartdate .or.  
 cph     &           uvecperiod    .ne. vvecperiod ) then  
 cph               print*,'For useCubedSphereExchange, S/R exf_set_uv.F'  
 cph               print*,'assumes that the u and v wind or wind stress'  
 cph               print*,'files have the same startdate and period.'  
 cph               stop  
 cph            endif  
 cph         endif  
       endif  
   
 #if ( defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_WIND) )  
       IF ( ustressfile .NE. ' ' .OR. ustressfile .NE. ' ' ) THEN  
         STOP  
      & 'S/R EXF_CHECK: use u,v_wind components but not wind-stress'  
       ENDIF  
 #endif  
 #ifndef ALLOW_ATM_WIND  
       IF ( uwindfile .NE. ' ' .OR. vwindfile .NE. ' ' ) THEN  
         STOP  
112       & 'S/R EXF_CHECK: read-in wind-stress but not u,v_wind components'       & 'S/R EXF_CHECK: read-in wind-stress but not u,v_wind components'
113            CALL PRINT_ERROR( msgBuf, myThid )
114            errCount = errCount + 1
115           ENDIF
116          ENDIF
117    
118    #ifdef ALLOW_SALTFLX
119          IF ( useSEAICE .OR. useThSIce )  THEN
120           IF ( saltflxfile .NE. ' ' ) THEN
121            WRITE(msgBuf,'(2A)') 'S/R EXF_CHECK: exf salt flux is not',
122         &       ' allowed when using either pkg/seaice or pkg/thsice'
123            CALL PRINT_ERROR( msgBuf, myThid )
124            errCount = errCount + 1
125           ENDIF
126        ENDIF        ENDIF
127  #endif  #endif /* ALLOW_SALTFLX */
128    
129    #ifdef ALLOW_ZENITHANGLE
130          IF ( ( useExfZenIncoming .OR. select_ZenAlbedo.NE.0 ) .AND.
131         &     ( usingCartesianGrid .OR. usingCylindricalGrid ) ) THEN
132           WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ZENITHANGLE code ',
133         &      'does not work for cartesian and cylindrical grids'
134           CALL PRINT_ERROR( msgBuf, myThid )
135           errCount = errCount + 1
136          ENDIF
137          IF ( select_ZenAlbedo.LT.0 .OR. select_ZenAlbedo.GT.3 ) THEN
138           WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: unsupported ',
139         &      'select_ZenAlbedo choice'
140           CALL PRINT_ERROR( msgBuf, myThid )
141           errCount = errCount + 1
142          ENDIF
143          IF ( select_ZenAlbedo.EQ.2 ) THEN
144           WRITE(msgBuf,'(A,A)')
145         &      'S/R EXF_CHECK: *** WARNING *** for daily mean albedo, ',
146         &      'it is advised to use select_ZenAlbedo.EQ.1 instead of 2'
147           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
148         &                     SQUEEZE_RIGHT, myThid )
149          ENDIF
150          IF ( select_ZenAlbedo.EQ.3 .AND. swdownperiod.GT.21600. ) THEN
151           WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: using diurnal albedo ',
152         &        'formula requires diurnal downward shortwave forcing'
153           CALL PRINT_ERROR( msgBuf, myThid )
154           errCount = errCount + 1
155          ENDIF
156          IF ( select_ZenAlbedo.EQ.3 .AND. swdownperiod.GT.3600. ) THEN
157           WRITE(msgBuf,'(A,A,A)')
158         &      'S/R EXF_CHECK: *** WARNING *** ',
159         &      'the diurnal albedo formula is likely not safe for such ',
160         &      'coarse temporal resolution downward shortwave forcing'
161           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
162         &                     SQUEEZE_RIGHT, myThid )
163          ENDIF
164    #else /* ALLOW_ZENITHANGLE */
165          IF ( useExfZenIncoming .OR. select_ZenAlbedo.NE.0 ) THEN
166            WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: unsupported option',
167         &        ' when ALLOW_ZENITHANGLE is not defined'
168            CALL PRINT_ERROR( msgBuf, myThid )
169            errCount = errCount + 1
170          ENDIF
171    #endif /* ALLOW_ZENITHANGLE */
172    
173  #ifdef USE_EXF_INTERPOLATION  #ifdef USE_EXF_INTERPOLATION
174        if ( climsst_nlat .GT. MAX_LAT_INC )        IF ( usingCartesianGrid ) THEN
175       & stop 'stopped in exf_readparms: climsst_nlat > MAX_LAT_INC'         WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ',
176        if ( climsss_nlat .GT. MAX_LAT_INC )       &      'USE_EXF_INTERPOLATION assumes latitude/longitude'
177       & stop 'stopped in exf_readparms: climsss_nlat > MAX_LAT_INC'         CALL PRINT_ERROR( msgBuf, myThid )
178        if ( usingCartesianGrid ) then         WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ',
179         print*,'USE_EXF_INTERPOLATION assumes latitude/longitude'       &      'input and output coordinates. Trivial to extend to'
180         print*,'input and output coordinates.  Trivial to extend to'         CALL PRINT_ERROR( msgBuf, myThid )
181         print*,'cartesian coordinates, but has not yet been done.'         WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ',
182         stop       &      'cartesian coordinates, but has not yet been done.'
183        endif         CALL PRINT_ERROR( msgBuf, myThid )
184        IF ( readStressOnAgrid.OR.readStressOnCgrid ) THEN         errCount = errCount + 1
185          STOP 'wind-stress position defined by Long/Lat'        ENDIF
186        ENDIF  
187  #else        CALL EXF_CHECK_INTERP('ustress',ustressfile,ustress_interpMethod,
188        IF ( ustressfile .NE. ' ' .OR. ustressfile .NE. ' ' ) THEN       &     ustress_nlat,ustress_lon_inc,errCount,myThid)
189          IF ( (readStressOnAgrid.AND.readStressOnCgrid) .OR.        CALL EXF_CHECK_INTERP('vstress',vstressfile,vstress_interpMethod,
190       &   .NOT.(readStressOnAgrid.OR.readStressOnCgrid) ) THEN       &     vstress_nlat,vstress_lon_inc,errCount,myThid)
191            STOP        CALL EXF_CHECK_INTERP('hflux',hfluxfile,hflux_interpMethod,
192       & 'S/R EXF_CHECK: Select 1 wind-stress position: A or C-grid'       &     hflux_nlat,hflux_lon_inc,errCount,myThid)
193          CALL EXF_CHECK_INTERP('sflux',sfluxfile,sflux_interpMethod,
194         &     sflux_nlat,sflux_lon_inc,errCount,myThid)
195          CALL EXF_CHECK_INTERP('swflux',swfluxfile,swflux_interpMethod,
196         &     swflux_nlat,swflux_lon_inc,errCount,myThid)
197          CALL EXF_CHECK_INTERP('runoff',runofffile,runoff_interpMethod,
198         &     runoff_nlat,runoff_lon_inc,errCount,myThid)
199          CALL EXF_CHECK_INTERP('saltflx',saltflxfile,saltflx_interpMethod,
200         &     saltflx_nlat,saltflx_lon_inc,errCount,myThid)
201          CALL EXF_CHECK_INTERP('atemp',atempfile,atemp_interpMethod,
202         &     atemp_nlat,atemp_lon_inc,errCount,myThid)
203          CALL EXF_CHECK_INTERP('aqh',aqhfile,aqh_interpMethod,
204         &     aqh_nlat,aqh_lon_inc,errCount,myThid)
205          CALL EXF_CHECK_INTERP( 'hs', hs_file, hs_interpMethod,
206         &     hs_nlat, hs_lon_inc, errCount, myThid )
207          CALL EXF_CHECK_INTERP( 'hl', hl_file, hl_interpMethod,
208         &     hl_nlat, hl_lon_inc, errCount, myThid )
209          CALL EXF_CHECK_INTERP('evap',evapfile,evap_interpMethod,
210         &     evap_nlat,evap_lon_inc,errCount,myThid)
211          CALL EXF_CHECK_INTERP('precip',precipfile,precip_interpMethod,
212         &     precip_nlat,precip_lon_inc,errCount,myThid)
213          CALL EXF_CHECK_INTERP('snowprecip',
214         &     snowprecipfile,snowprecip_interpMethod,
215         &     snowprecip_nlat,snowprecip_lon_inc,errCount,myThid)
216          CALL EXF_CHECK_INTERP('uwind',uwindfile,uwind_interpMethod,
217         &     uwind_nlat,uwind_lon_inc,errCount,myThid)
218          CALL EXF_CHECK_INTERP('vwind',vwindfile,vwind_interpMethod,
219         &     vwind_nlat,vwind_lon_inc,errCount,myThid)
220          CALL EXF_CHECK_INTERP('wspeed',wspeedfile,wspeed_interpMethod,
221         &     wspeed_nlat,wspeed_lon_inc,errCount,myThid)
222          CALL EXF_CHECK_INTERP('lwflux',lwfluxfile,lwflux_interpMethod,
223         &     lwflux_nlat,lwflux_lon_inc,errCount,myThid)
224          CALL EXF_CHECK_INTERP('swdown',swdownfile,swdown_interpMethod,
225         &     swdown_nlat,swdown_lon_inc,errCount,myThid)
226          CALL EXF_CHECK_INTERP('lwdown',lwdownfile,lwdown_interpMethod,
227         &     lwdown_nlat,lwdown_lon_inc,errCount,myThid)
228          CALL EXF_CHECK_INTERP('apressure',
229         &     apressurefile,apressure_interpMethod,
230         &     apressure_nlat,apressure_lon_inc,errCount,myThid)
231          CALL EXF_CHECK_INTERP('areamask',
232         &     areamaskfile,areamask_interpMethod,
233         &     areamask_nlat,areamask_lon_inc,errCount,myThid)
234          CALL EXF_CHECK_INTERP('climsst',climsstfile,climsst_interpMethod,
235         &     climsst_nlat,climsst_lon_inc,errCount,myThid)
236          CALL EXF_CHECK_INTERP('climsss',climsssfile,climsss_interpMethod,
237         &     climsss_nlat,climsss_lon_inc,errCount,myThid)
238          CALL EXF_CHECK_INTERP('climustr',
239         &     climustrfile,climustr_interpMethod,
240         &     climustr_nlat,climustr_lon_inc,errCount,myThid)
241          CALL EXF_CHECK_INTERP('climvstr',
242         &     climvstrfile,climvstr_interpMethod,
243         &     climvstr_nlat,climvstr_lon_inc,errCount,myThid)
244    
245    C-    some restrictions on 2-component vector field (might be relaxed later on)
246          IF ( ( uwind_interpMethod.GE.1 .AND. uwindfile.NE.' ' ) .OR.
247         &     ( vwind_interpMethod.GE.1 .AND. vwindfile.NE.' ' ) ) THEN
248           IF ( usingCurvilinearGrid .OR. rotateGrid ) THEN
249            IF ( uwind_interpMethod.EQ.0 .OR. uwindfile.EQ.' ' .OR.
250         &       vwind_interpMethod.EQ.0 .OR. vwindfile.EQ.' ' ) THEN
251    C-    stop if one expects interp+rotation (Curvilin-G) which will not happen
252             WRITE(msgBuf,'(A)')
253         &        'S/R EXF_CHECK: interp. needs 2 components (wind)'
254             CALL PRINT_ERROR( msgBuf, myThid )
255             errCount = errCount + 1
256          ENDIF          ENDIF
257        ELSE          IF ( uwindStartTime .NE. vwindStartTime .OR.
258          IF ( readStressOnAgrid .OR. readStressOnCgrid ) THEN       &       uwindperiod    .NE. vwindperiod   ) THEN
259            STOP 'S/R EXF_CHECK: wind-stress position irrelevant'           WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ',
260         &        'For CurvilinearGrid/RotatedGrid, the u and v wind '
261             CALL PRINT_ERROR( msgBuf, myThid )
262             WRITE(msgBuf,'(A,A,A)') 'S/R EXF_CHECK: ',
263         &        'files have to have the same StartTime and period, ',
264         &        'because S/R EXF_SET_UV assumes that.'
265             CALL PRINT_ERROR( msgBuf, myThid )
266             errCount = errCount + 1
267            ENDIF
268           ENDIF
269          ENDIF
270          IF ( (ustress_interpMethod.GE.1 .AND. ustressfile.NE.' ') .OR.
271         &     (vstress_interpMethod.GE.1 .AND. vstressfile.NE.' ') ) THEN
272           IF ( readStressOnCgrid ) THEN
273            WRITE(msgBuf,'(A,A)')
274         &       'S/R EXF_CHECK: readStressOnCgrid=.TRUE. ',
275         &       'and interp wind-stress (=A-grid) are not compatible'
276            CALL PRINT_ERROR( msgBuf, myThid )
277            errCount = errCount + 1
278           ENDIF
279           IF ( usingCurvilinearGrid .OR. rotateGrid ) THEN
280            IF ( ustress_interpMethod.EQ.0 .OR. ustressfile.EQ.' ' .OR.
281         &       vstress_interpMethod.EQ.0 .OR. vstressfile.EQ.' ' ) THEN
282    C-    stop if one expects interp+rotation (Curvilin-G) which will not happen
283             WRITE(msgBuf,'(A)')
284         &        'S/R EXF_CHECK: interp. needs 2 components (wind-stress)'
285             CALL PRINT_ERROR( msgBuf, myThid )
286             errCount = errCount + 1
287            ENDIF
288            IF ( ustressStartTime .NE. vstressStartTime .OR.
289         &       ustressperiod    .NE. vstressperiod   ) THEN
290             WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ',
291         &    'For CurvilinearGrid/RotatedGrid, the u and v wind stress '
292             CALL PRINT_ERROR( msgBuf, myThid )
293             WRITE(msgBuf,'(A,A,A)') 'S/R EXF_CHECK: ',
294         &        'files have to have the same StartTime and period, ',
295         &        'because S/R EXF_SET_UV assumes that.'
296             CALL PRINT_ERROR( msgBuf, myThid )
297             errCount = errCount + 1
298          ENDIF          ENDIF
299           ENDIF
300        ENDIF        ENDIF
 #endif  
301    
302  #ifdef ALLOW_CLIMTEMP_RELAXATION        IF ( (ustress_interpMethod.EQ.0 .AND. ustressfile.NE.' ') .OR.
303        STOP 'ALLOW_CLIMTEMP_RELAXATION no longer supported. Use pkg/rbcs'       &     (vstress_interpMethod.EQ.0 .AND. vstressfile.NE.' ') ) THEN
304  #endif  #else /* ndef USE_EXF_INTERPOLATION */
305          IF     ( ustressfile .NE. ' ' .OR. vstressfile .NE. ' ' ) THEN
306    #endif /* USE_EXF_INTERPOLATION */
307           IF (     (readStressOnAgrid.AND.readStressOnCgrid) .OR.
308         &      .NOT.(readStressOnAgrid.OR.readStressOnCgrid) ) THEN
309            WRITE(msgBuf,'(A)')
310         & 'S/R EXF_CHECK: Select 1 wind-stress position: A or C-grid'
311            CALL PRINT_ERROR( msgBuf, myThid )
312            errCount = errCount + 1
313           ENDIF
314           IF (rotateStressOnAgrid.AND..NOT.readStressOnAgrid) THEN
315            WRITE(msgBuf,'(2A)') 'S/R EXF_CHECK: rotateStressOnAgrid ',
316         & 'only applies to cases readStressOnAgrid is true'
317            CALL PRINT_ERROR( msgBuf, myThid )
318            errCount = errCount + 1
319           ENDIF
320    
321          ELSE
322           IF ( readStressOnAgrid .OR. readStressOnCgrid .OR.
323         &      rotateStressOnAgrid) THEN
324            WRITE(msgBuf,'(A)')
325         &       'S/R EXF_CHECK: wind-stress position irrelevant'
326            CALL PRINT_ERROR( msgBuf, myThid )
327            errCount = errCount + 1
328           ENDIF
329          ENDIF
330    
331    #ifdef USE_NO_INTERP_RUNOFF
332          WRITE(msgBuf,'(A)')
333         &     'S/R EXF_CHECK: USE_NO_INTERP_RUNOFF code has been removed;'
334          CALL PRINT_ERROR( msgBuf, myThid )
335          WRITE(msgBuf,'(A,A)')
336         &     'S/R EXF_CHECK: use instead "runoff_interpMethod=0"',
337         &      ' in "data.exf" (EXF_NML_04)'
338          CALL PRINT_ERROR( msgBuf, myThid )
339          errCount = errCount + 1
340    #endif /* USE_NO_INTERP_RUNOFF */
341    
342    #ifdef ALLOW_CLIMTEMP_RELAXATION
343          WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ',
344         &     'ALLOW_CLIMTEMP_RELAXATION no longer supported. Use pkg/rbcs'
345          CALL PRINT_ERROR( msgBuf, myThid )
346          errCount = errCount + 1
347    #endif /* ALLOW_CLIMTEMP_RELAXATION */
348    
349  #ifdef ALLOW_CLIMSALT_RELAXATION  #ifdef ALLOW_CLIMSALT_RELAXATION
350        STOP 'ALLOW_CLIMSALT_RELAXATION no longer supported. Use pkg/rbcs'        WRITE(msgBuf,'(A,A)') 'S/R EXF_CHECK: ',
351  #endif       &     'ALLOW_CLIMSALT_RELAXATION no longer supported. Use pkg/rbcs'
352          CALL PRINT_ERROR( msgBuf, myThid )
353          errCount = errCount + 1
354    #endif /* ALLOW_CLIMSALT_RELAXATION */
355    
356          IF ( climsstTauRelax.NE.0. ) THEN
357    #ifndef ALLOW_CLIMSST_RELAXATION
358           WRITE(msgBuf,'(A)') 'S/R EXF_CHECK: climsstTauRelax > 0'
359           CALL PRINT_ERROR( msgBuf, myThid )
360           WRITE(msgBuf,'(A)')
361         &      'S/R EXF_CHECK: but ALLOW_CLIMSST_RELAXATION is not defined'
362           CALL PRINT_ERROR( msgBuf, myThid )
363           errCount = errCount + 1
364    #endif /* ndef ALLOW_CLIMSST_RELAXATION */
365           IF ( climsstfile.EQ.' ' ) THEN
366            WRITE(msgBuf,'(A)') 'S/R EXF_CHECK: climsstTauRelax > 0 but'
367            CALL PRINT_ERROR( msgBuf, myThid )
368            WRITE(msgBuf,'(A)') 'S/R EXF_CHECK: climsstfile is not set'
369            CALL PRINT_ERROR( msgBuf, myThid )
370            errCount = errCount + 1
371           ENDIf
372          ENDIf
373    
374          IF ( climsssTauRelax.NE.0. ) THEN
375    #ifndef ALLOW_CLIMSSS_RELAXATION
376           WRITE(msgBuf,'(A)') 'S/R EXF_CHECK: climsssTauRelax > 0'
377           CALL PRINT_ERROR( msgBuf, myThid )
378           WRITE(msgBuf,'(A)')
379         &      'S/R EXF_CHECK: but ALLOW_CLIMSSS_RELAXATION is not defined'
380           CALL PRINT_ERROR( msgBuf, myThid )
381           errCount = errCount + 1
382    #endif /* ALLOW_CLIMSSS_RELAXATION */
383           IF ( climsssfile.EQ.' ' ) THEN
384            WRITE(msgBuf,'(A)') 'S/R EXF_CHECK: climsssTauRelax > 0 but'
385            CALL PRINT_ERROR( msgBuf, myThid )
386            WRITE(msgBuf,'(A)') 'S/R EXF_CHECK: climsssfile is not set'
387            CALL PRINT_ERROR( msgBuf, myThid )
388            errCount = errCount + 1
389           ENDIF
390          ENDIF
391    
392    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
393    
394          IF ( errCount.GE.1 ) THEN
395            WRITE(msgBuf,'(A,I3,A)')
396         &       'EXF_CHECK: detected', errCount,' fatal error(s)'
397            CALL PRINT_ERROR( msgBuf, myThid )
398            CALL ALL_PROC_DIE( 0 )
399            STOP 'ABNORMAL END: S/R EXF_CHECK'
400          ENDIF
401    
402          _END_MASTER(myThid)
403    
404          RETURN
405          END
406    
407    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
408    
409    CBOP
410    C     !ROUTINE: EXF_CHECK_INTERP
411    C     !INTERFACE:
412    
413          SUBROUTINE EXF_CHECK_INTERP(
414         I               loc_name, loc_file,
415         I               loc_interpMethod, loc_nlat, loc_lon_inc,
416         U               errCount,
417         I               myThid )
418    
419    C     !DESCRIPTION: \bv
420    C     *==========================================================*
421    C     | S/R EXF_CHECK_INTERP
422    C     | o Check parameters for one of the pkg/exf variable
423    C     *==========================================================*
424    C     \ev
425    
426    C     !USES:
427          IMPLICIT NONE
428    
429    C     == Global variables ===
430    #include "EEPARAMS.h"
431    #include "EXF_PARAM.h"
432    
433    C     !INPUT/OUTPUT PARAMETERS:
434    C     myThid   :: my Thread Id number
435          CHARACTER*(*) loc_name
436          CHARACTER*(*) loc_file
437          INTEGER loc_interpMethod
438          INTEGER loc_nlat
439          _RL loc_lon_inc
440          INTEGER errCount
441          INTEGER myThid
442    
443    C     !LOCAL VARIABLES:
444    C     msgBuf   :: Informational/error message buffer
445          CHARACTER*(MAX_LEN_MBUF) msgBuf
446    CEOP
447    
448          IF ( loc_interpMethod.GE.1 .AND. loc_file.NE.' ' ) THEN
449           IF ( loc_nlat .GT. MAX_LAT_INC ) THEN
450            WRITE(msgBuf,'(3A)') 'S/R EXF_CHECK_INTERP: ',loc_name,
451         &                      '_nlat > MAX_LAT_INC'
452            CALL PRINT_ERROR( msgBuf, myThid )
453            errCount = errCount + 1
454           ENDIF
455           IF ( loc_lon_inc.GT.500. ) THEN
456            WRITE(msgBuf,'(4A,1PE16.8)') 'S/R EXF_CHECK_INTERP: ',
457         &    'Invalid value for: ',loc_name,'_lon_inc =', loc_lon_inc
458            CALL PRINT_ERROR( msgBuf, myThid )
459            WRITE(msgBuf,'(4A)') 'S/R EXF_CHECK_INTERP: Fix it ',
460         &    'or Turn off ',loc_name,'-interp (interpMethod=0)'
461            CALL PRINT_ERROR( msgBuf, myThid )
462            errCount = errCount + 1
463           ENDIF
464          ENDIF
465    
466        RETURN        RETURN
467        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22