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

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

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

revision 1.31 by cnh, Tue Sep 8 01:37:49 1998 UTC revision 1.83 by adcroft, Thu Aug 15 17:25:31 2002 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6    CBOP
7    C     !ROUTINE: INI_PARMS
8    C     !INTERFACE:
9        SUBROUTINE INI_PARMS( myThid )        SUBROUTINE INI_PARMS( myThid )
10  C     /==========================================================\  C     !DESCRIPTION: \bv
11  C     | SUBROUTINE INI_PARMS                                     |  C     *==========================================================*
12  C     | o Routine to set model "parameters"                      |  C     | SUBROUTINE INI_PARMS                                      
13  C     |==========================================================|  C     | o Routine to set model "parameters"                      
14  C     | Notes:                                                   |  C     *==========================================================*
15  C     | ======                                                   |  C     | Notes:                                                    
16  C     | The present version of this routine is a place-holder.   |  C     | ======                                                    
17  C     | A production version needs to handle parameters from an  |  C     | The present version of this routine is a place-holder.    
18  C     | external file and possibly reading in some initial field |  C     | A production version needs to handle parameters from an  
19  C     | values.                                                  |  C     | external file and possibly reading in some initial field  
20  C     \==========================================================/  C     | values.                                                  
21    C     *==========================================================*
22    C     \ev
23    
24    C     !USES:
25          IMPLICIT NONE
26  C     === Global variables ===  C     === Global variables ===
27  #include "SIZE.h"  #include "SIZE.h"
28  #include "EEPARAMS.h"  #include "EEPARAMS.h"
29  #include "PARAMS.h"  #include "PARAMS.h"
30  #include "GRID.h"  #include "GRID.h"
31  #include "CG2D.h"  #include "EOS.h"
32    
33    C     !INPUT/OUTPUT PARAMETERS:
34  C     === Routine arguments ===  C     === Routine arguments ===
35  C     myThid - Number of this instance of INI_PARMS  C     myThid - Number of this instance of INI_PARMS
36        INTEGER myThid        INTEGER myThid
37    
38    C     !LOCAL VARIABLES:
39  C     === Local variables ===  C     === Local variables ===
40  C     dxSpacing, dySpacing - Default spacing in X and Y.  C     dxSpacing, dySpacing - Default spacing in X and Y.
41  C                            Units are that of coordinate system  C                            Units are that of coordinate system
# Line 41  C     K, I, J - Loop counters Line 51  C     K, I, J - Loop counters
51  C     xxxDefault - Default value for variable xxx  C     xxxDefault - Default value for variable xxx
52        _RL  dxSpacing        _RL  dxSpacing
53        _RL  dySpacing        _RL  dySpacing
54          CHARACTER*(MAX_LEN_FNAM) delXfile
55          CHARACTER*(MAX_LEN_FNAM) delYfile
56        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
57        CHARACTER*(MAX_LEN_PREC) record        CHARACTER*(MAX_LEN_PREC) record
58        INTEGER goptCount        INTEGER goptCount
59        INTEGER K, I, J, IL, iUnit        INTEGER K, I, IL, iUnit
60        INTEGER errIO        INTEGER errIO
61        INTEGER  IFNBLNK        INTEGER  IFNBLNK
62        EXTERNAL IFNBLNK        EXTERNAL IFNBLNK
# Line 56  C     dependency. Line 68  C     dependency.
68        _RL diffKrTDefault        _RL diffKrTDefault
69        _RL diffKrSDefault        _RL diffKrSDefault
70        _RL hFacMinDrDefault        _RL hFacMinDrDefault
71        _RL delRDefault        _RL delRDefault(Nr)
72  C     zCoordInputData - These are used to select between different coordinate systems.        _RS rkFacDefault
73  C     pCoordInputData   The vertical coordinate system in the rest of the model is  C     zCoordInputData :: Variables used to select between different coordinate systems.
74  C     rCoordInputData   written in terms of r. In the model "data" file input data can  C     pCoordInputData :: The vertical coordinate system in the rest of the model is
75  C     coordsSet         be interms of z, p or r.  C     rCoordInputData :: written in terms of r. In the model "data" file input data can
76  C                       e.g. delZ or delP or delR for the vertical grid spacing.  C     coordsSet       :: be interms of z, p or r.
77  C                       The following rules apply:  C                     :: e.g. delZ or delP or delR for the vertical grid spacing.
78  C                       All parameters must use the same vertical coordinate system.  C                     :: The following rules apply:
79  C                       e.g. delZ and viscAz is legal but  C                     :: All parameters must use the same vertical coordinate system.
80  C                            delZ and viscAr is an error.  C                     ::  e.g. delZ and viscAz is legal but
81  C                       Similarly specifyinh delZ and delP is an error.  C                     ::       delZ and viscAr is an error.
82  C                       zCoord..., pCoord..., rCoord... are used to flag when z, p or r are  C                     :: Similarly specifyinh delZ and delP is an error.
83  C                       used. coordsSet counts how many vertical coordinate systems have been  C                     :: zCoord..., pCoord..., rCoord... are used to flag when z, p or r are
84  C                       used to specify variables. coordsSet > 1 is an error.  C                     :: used. coordsSet counts how many vertical coordinate systems have been
85    C                        used to specify variables. coordsSet > 1 is an error.
86  C  C
87        LOGICAL zCoordInputData        LOGICAL zCoordInputData
88        LOGICAL pCoordInputData        LOGICAL pCoordInputData
89        LOGICAL rCoordInputData        LOGICAL rCoordInputData
90        INTEGER coordsSet        INTEGER coordsSet
91    
92    C     Retired main data file parameters. Kept here to trap use of old data files.
93    C     zonal_filt_lat  - Moved to package "zonal_filt"
94    C     nRetired :: Counter used to trap gracefully namelists containing "retired"
95    C              :: parameters. These are parameters that are either no-longer used
96    C                 or that have moved to a different input file and/or namelist.
97          _RL zonal_filt_lat
98          INTEGER nRetired
99    CEOP
100    
101  C--   Continuous equation parameters  C--   Continuous equation parameters
102        NAMELIST /PARM01/        NAMELIST /PARM01/
103       & gravity, gBaro, rhonil, tAlpha, sBeta, f0, beta,       & gravitySign,
104       &   viscAh,  viscAz,  viscA4,       & gravity, gBaro, rhonil, tAlpha, sBeta, f0, beta, omega,
105       &  diffKhT, diffKzT, diffK4T,       & viscAh,  viscAz,  viscA4, cosPower, viscAstrain, viscAtension,
106       &  diffKhS, diffKzS, diffK4S,       & diffKhT, diffKzT, diffK4T,
107       &  GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,GMmaxval,       & diffKhS, diffKzS, diffK4S,
108       &  tRef, sRef, eosType,       & tRef, sRef, eosType, Integr_GeoPot,
109         & no_slip_sides,no_slip_bottom,
110       & momViscosity,  momAdvection, momForcing, useCoriolis,       & momViscosity,  momAdvection, momForcing, useCoriolis,
111       & momPressureForcing, metricTerms,       & momPressureForcing, metricTerms, vectorInvariantMomentum,
112       & tempDiffusion, tempAdvection, tempForcing,       & tempDiffusion, tempAdvection, tempForcing,
113       & saltDiffusion, saltAdvection, saltForcing,       & saltDiffusion, saltAdvection, saltForcing,
114         & implicSurfPress, implicDiv2DFlow,
115       & implicitFreeSurface, rigidLid, freeSurfFac, hFacMin, hFacMinDz,       & implicitFreeSurface, rigidLid, freeSurfFac, hFacMin, hFacMinDz,
116       & tempStepping, saltStepping, momStepping, implicitDiffusion,       & exactConserv,uniformLin_PhiSurf,nonlinFreeSurf,hFacInf,hFacSup,
117         & staggerTimeStep,
118         & tempStepping, saltStepping, momStepping, tr1Stepping,
119         & implicitDiffusion, implicitViscosity,
120       & viscAr, diffKrT, diffKrS, hFacMinDr,       & viscAr, diffKrT, diffKrS, hFacMinDr,
121       & viscAp, diffKpT, diffKpS, hFacMinDp,       & viscAp, diffKpT, diffKpS, hFacMinDp,
122       & rhoConst, buoyancyRelation       & rhoConst, buoyancyRelation, HeatCapacity_Cp,
123         & writeBinaryPrec, readBinaryPrec, writeStatePrec,
124         & nonHydrostatic, globalFiles,
125         & allowFreezing, ivdc_kappa,
126         & bottomDragLinear,bottomDragQuadratic,
127         & usePickupBeforeC35, debugMode,
128         & readPickupWithTracer, writePickupWithTracer,
129         & tempAdvScheme, saltAdvScheme, tracerAdvScheme,
130         & multiDimAdvection, useEnergyConservingCoriolis,
131         & useJamartWetPoints,
132         & useRealFreshWaterFlux, convertFW2Salt,
133         & temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
134         & zonal_filt_lat
135    
136  C--   Elliptic solver parameters  C--   Elliptic solver parameters
137        NAMELIST /PARM02/        NAMELIST /PARM02/
138       & cg2dMaxIters, cg2dChkResFreq, cg2dTargetResidual, cg2dpcOffDFac       & cg2dMaxIters, cg2dChkResFreq, cg2dTargetResidual,
139         & cg2dTargetResWunit, cg2dpcOffDFac,
140         & cg3dMaxIters, cg3dChkResFreq, cg3dTargetResidual
141    
142  C--   Time stepping parammeters  C--   Time stepping parammeters
143        NAMELIST /PARM03/        NAMELIST /PARM03/
144       & nIter0, nTimeSteps, deltaT, deltaTmom, deltaTtracer, abEps, tauCD, rCD,       & nIter0, nTimeSteps, nEndIter,
145       & startTime, endTime, chkPtFreq, dumpFreq, taveFreq, deltaTClock,       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf,
146       & pChkPtFreq, cAdjFreq, tauThetaClimRelax, tauSaltClimRelax,       & forcing_In_AB, abEps, tauCD, rCD,
147         & startTime, endTime, chkPtFreq,
148         & dumpFreq, taveFreq, deltaTClock, diagFreq,
149         & monitorFreq, pChkPtFreq, cAdjFreq,
150         & tauThetaClimRelax, tauSaltClimRelax, tauTr1ClimRelax,
151       & periodicExternalForcing, externForcingPeriod, externForcingCycle       & periodicExternalForcing, externForcingPeriod, externForcingCycle
152    
153  C--   Gridding parameters  C--   Gridding parameters
154        NAMELIST /PARM04/        NAMELIST /PARM04/
155       & usingCartesianGrid, delZ, dxSpacing, dySpacing, delX, delY,       & usingCartesianGrid, dxSpacing, dySpacing, delX, delY, delZ,
156       & usingSphericalPolarGrid, phiMin, thetaMin, rSphere,       & usingSphericalPolarGrid, phiMin, thetaMin, rSphere,
157       & l, m, n, delP, delR, rkFac       & usingCurvilinearGrid, gravitySign,
158         & delP, delR, rkFac, Ro_SeaLevel, groundAtK1, delRc,
159         & delXfile, delYfile
160    
161  C--   Input files  C--   Input files
162        NAMELIST /PARM05/        NAMELIST /PARM05/
163       & bathyFile, hydrogThetaFile, hydrogSaltFile,       & bathyFile, topoFile, hydrogThetaFile, hydrogSaltFile,
164       & zonalWindFile, meridWindFile, thetaClimFile,       & zonalWindFile, meridWindFile,
165       & saltClimFile       & thetaClimFile, saltClimFile,
166         & surfQfile, EmPmRfile, surfQswfile,
167         & uVelInitFile, vVelInitFile, pSurfInitFile,
168         & dQdTFile, ploadFile
169    
170  C  C
171        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
172    
173    C     Defaults values for input parameters
174          CALL SET_DEFAULTS(
175         O   viscArDefault, diffKrTDefault, diffKrSDefault,
176         O   hFacMinDrDefault, delRdefault, rkFacDefault,
177         I   myThid )
178    
179  C--   Initialise "which vertical coordinate system used" flags.  C--   Initialise "which vertical coordinate system used" flags.
180        zCoordInputData = .FALSE.        zCoordInputData = .FALSE.
181        pCoordInputData = .FALSE.        pCoordInputData = .FALSE.
# Line 129  C--   Initialise "which vertical coordin Line 184  C--   Initialise "which vertical coordin
184        usingZCoords    = .FALSE.        usingZCoords    = .FALSE.
185        coordsSet       = 0        coordsSet       = 0
186    
187    C--   Iniialise retired parameters to unlikely value
188          nRetired = 0
189          zonal_filt_lat  = UNSET_RL
190    
191  C--   Open the parameter file  C--   Open the parameter file
192        OPEN(UNIT=scrUnit1,STATUS='SCRATCH')        OPEN(UNIT=scrUnit1,STATUS='SCRATCH')
193        OPEN(UNIT=scrUnit2,STATUS='SCRATCH')        OPEN(UNIT=scrUnit2,STATUS='SCRATCH')
194        OPEN(UNIT=modelDataUnit,FILE='data',STATUS='OLD',err=1,IOSTAT=errIO)        OPEN(UNIT=modelDataUnit,FILE='data',STATUS='OLD',
195        IF ( errIO .GE. 0 ) GOTO 2       &     IOSTAT=errIO)
196      1 CONTINUE        IF ( errIO .LT. 0 ) THEN
197         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
198       &  'S/R INI_PARMS'       &  'S/R INI_PARMS'
199         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
# Line 146  C--   Open the parameter file Line 205  C--   Open the parameter file
205         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
206         CALL MODELDATA_EXAMPLE( myThid )         CALL MODELDATA_EXAMPLE( myThid )
207         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
208      2 CONTINUE        ENDIF    
209    
210   1000 CONTINUE        DO WHILE ( .TRUE. )
211        READ(modelDataUnit,FMT='(A)',END=1001) RECORD         READ(modelDataUnit,FMT='(A)',END=1001) RECORD
212        IL = MAX(ILNBLNK(RECORD),1)         IL = MAX(ILNBLNK(RECORD),1)
213        IF ( RECORD(1:1) .NE. commentCharacter )         IF ( RECORD(1:1) .NE. commentCharacter )
214       &    WRITE(UNIT=scrUnit1,FMT='(A)') RECORD(:IL)       &     WRITE(UNIT=scrUnit1,FMT='(A)') RECORD(:IL)
215         WRITE(UNIT=scrUnit2,FMT='(A)') RECORD(:IL)          WRITE(UNIT=scrUnit2,FMT='(A)') RECORD(:IL)
216        GOTO 1000        ENDDO
217   1001 CONTINUE   1001 CONTINUE
218        CLOSE(modelDataUnit)        CLOSE(modelDataUnit)
219    
220  C--   Report contents of model parameter file  C--   Report contents of model parameter file
221        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
222       &'// ======================================================='       &'// ======================================================='
223        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
224         &                    SQUEEZE_RIGHT , 1)
225        WRITE(msgBuf,'(A)') '// Model parameter file "data"'        WRITE(msgBuf,'(A)') '// Model parameter file "data"'
226        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
227         &                    SQUEEZE_RIGHT , 1)
228        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
229       &'// ======================================================='       &'// ======================================================='
230        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
231       &  SQUEEZE_RIGHT , 1)       &  SQUEEZE_RIGHT , 1)
232        iUnit = scrUnit2        iUnit = scrUnit2
233        REWIND(iUnit)        REWIND(iUnit)
234   2000 CONTINUE        DO WHILE ( .TRUE. )
235         READ(UNIT=iUnit,FMT='(A)',END=2001) RECORD         READ(UNIT=iUnit,FMT='(A)',END=2001) RECORD
236         IL = MAX(ILNBLNK(RECORD),1)         IL = MAX(ILNBLNK(RECORD),1)
237         WRITE(msgBuf,'(A,A)') '>',RECORD(:IL)         WRITE(msgBuf,'(A,A)') '>',RECORD(:IL)
238         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, SQUEEZE_RIGHT , 1)         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
239        GOTO 2000       &                    SQUEEZE_RIGHT , 1)
240          ENDDO
241   2001 CONTINUE   2001 CONTINUE
242        CLOSE(iUnit)        CLOSE(iUnit)
243        WRITE(msgBuf,'(A)') ' '        WRITE(msgBuf,'(A)') ' '
# Line 188  C--   Read settings from model parameter Line 250  C--   Read settings from model parameter
250        REWIND(iUnit)        REWIND(iUnit)
251    
252  C--   Set default "physical" parameters  C--   Set default "physical" parameters
       DO K =1,Nr  
        tRef(K) = 30.D0 - FLOAT( K )  
       ENDDO  
       gravity  =   9.81D0  
       gBaro    = gravity  
       rhoNil   = 999.8D0  
       rhoConst = 999.8D0  
       f0       = 1.D-4  
       beta     = 1.D-11  
       viscAh   = 1.D3  
       diffKhT  = 1.D3  
       diffKhS  = 1.D3  
       viscArDefault = 1.D-3  
253        viscAz   = UNSET_RL            viscAz   = UNSET_RL    
254        viscAr   = UNSET_RL        viscAr   = UNSET_RL
255        viscAp   = UNSET_RL        viscAp   = UNSET_RL
       diffKrTDefault = 1.D-5  
256        diffKzT  = UNSET_RL        diffKzT  = UNSET_RL
257        diffKpT  = UNSET_RL        diffKpT  = UNSET_RL
258        diffKrT  = UNSET_RL        diffKrT  = UNSET_RL
       diffKrSDefault = 1.D-5  
259        diffKzS  = UNSET_RL        diffKzS  = UNSET_RL
260        diffKpS  = UNSET_RL        diffKpS  = UNSET_RL
261        diffKrS  = UNSET_RL        diffKrS  = UNSET_RL
262        viscA4   = 0.        gBaro    = UNSET_RL
263        diffK4T  = 0.        rhoConst = UNSET_RL
       diffK4S  = 0.  
       GMmaxslope   =   1.D-2  
       GMlength     = 200.D3  
       GMalpha      = 0.D0  
       GMdepth      = 1000.D0  
       GMkbackground= 0.D0  
       GMmaxval     = 2500.D0  
       tAlpha       = 2.D-4  
       sBeta        = 7.4D-4  
       eosType      = 'LINEAR'  
       buoyancyRelation    = 'OCEANIC'  
       implicitFreeSurface = .TRUE.  
       rigidLid            = .FALSE.  
       freeSurfFac         = 1.D0  
       hFacMin             = 0.D0  
       hFacMinDrDefault    = 0.D0  
264        hFacMinDr           = UNSET_RL        hFacMinDr           = UNSET_RL
265        hFacMinDz           = UNSET_RL        hFacMinDz           = UNSET_RL
266        hFacMinDp           = UNSET_RL        hFacMinDp           = UNSET_RL
267        momViscosity        = .TRUE.        convertFW2Salt = UNSET_RL
268        momAdvection        = .TRUE.        tAlpha              = UNSET_RL
269        momForcing          = .TRUE.        sBeta               = UNSET_RL
270        useCoriolis         = .TRUE.        READ(UNIT=iUnit,NML=PARM01) !,IOSTAT=errIO)
271        momPressureForcing  = .TRUE.        IF ( errIO .LT. 0 ) THEN
       momStepping         = .TRUE.  
       tempStepping        = .TRUE.  
       saltStepping        = .TRUE.  
       metricTerms         = .TRUE.  
       implicitDiffusion   = .FALSE.  
       READ(UNIT=iUnit,NML=PARM01,IOSTAT=errIO,err=3)  
       IF ( errIO .GE. 0 ) GOTO 4  
     3 CONTINUE  
272         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
273       &  'S/R INI_PARMS'       &  'S/R INI_PARMS'
274         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
# Line 260  C--   Set default "physical" parameters Line 283  C--   Set default "physical" parameters
283         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
284         CALL MODELDATA_EXAMPLE( myThid )         CALL MODELDATA_EXAMPLE( myThid )
285         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
286     4  CONTINUE        ELSE
287           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM01 : OK'
288           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
289         &                    SQUEEZE_RIGHT , 1)
290          ENDIF
291        IF ( implicitFreeSurface ) freeSurfFac = 1.D0        IF ( implicitFreeSurface ) freeSurfFac = 1.D0
292        IF ( rigidLid            ) freeSurfFac = 0.D0        IF ( rigidLid            ) freeSurfFac = 0.D0
293          IF ( gBaro .EQ. UNSET_RL ) gBaro=gravity
294          IF ( rhoConst .EQ. UNSET_RL ) rhoConst=rhoNil
295  C--   Momentum viscosity on/off flag.  C--   Momentum viscosity on/off flag.
296        IF ( momViscosity        ) THEN        IF ( momViscosity        ) THEN
297         vfFacMom = 1.D0         vfFacMom = 1.D0
# Line 297  C--   Metric terms on/off flag. Line 326  C--   Metric terms on/off flag.
326        IF ( metricTerms ) THEN        IF ( metricTerms ) THEN
327         mTFacMom = 1.D0         mTFacMom = 1.D0
328        ELSE        ELSE
329         mTFacMom = 1.D0         mTFacMom = 0.D0
330        ENDIF        ENDIF
331    C--   Advection and Forcing for Temp and salt  on/off flags
332          tempAdvection = tempStepping .AND. tempAdvection
333          tempForcing   = tempStepping .AND. tempForcing
334          saltAdvection = saltStepping .AND. saltAdvection
335          saltForcing   = saltStepping .AND. saltForcing
336  C--   z,p,r coord input switching.  C--   z,p,r coord input switching.
337        IF ( viscAz .NE. UNSET_RL ) zCoordInputData = .TRUE.        IF ( viscAz .NE. UNSET_RL ) zCoordInputData = .TRUE.
338        IF ( viscAp .NE. UNSET_RL ) pCoordInputData = .TRUE.        IF ( viscAp .NE. UNSET_RL ) pCoordInputData = .TRUE.
# Line 324  C--   z,p,r coord input switching. Line 358  C--   z,p,r coord input switching.
358        IF ( hFacMinDz .NE. UNSET_RL ) zCoordInputData = .TRUE.        IF ( hFacMinDz .NE. UNSET_RL ) zCoordInputData = .TRUE.
359        IF ( hFacMinDp .NE. UNSET_RL ) pCoordInputData = .TRUE.        IF ( hFacMinDp .NE. UNSET_RL ) pCoordInputData = .TRUE.
360        IF ( hFacMinDr .NE. UNSET_RL ) rCoordInputData = .TRUE.        IF ( hFacMinDr .NE. UNSET_RL ) rCoordInputData = .TRUE.
361        IF ( hFacMinDz .EQ. UNSET_RL ) hFacMinDr       = hFacMinDz        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDz
362        IF ( hFacMinDp .EQ. UNSET_RL ) hFacMinDr       = hFacMinDp        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDp
363        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDrDefault        IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr       = hFacMinDrDefault
364    
365        IF ( implicitFreeSurface .AND.  rigidLid ) THEN        IF (convertFW2Salt.EQ.UNSET_RL) THEN
366         WRITE(msgBuf,'(A)')          convertFW2Salt = 35.
367       &  'S/R INI_PARMS: Cannot select both implicitFreeSurface and rigidLid.'          IF (useRealFreshWaterFlux) convertFW2Salt=-1
368          ENDIF
369    
370          IF ( ivdc_kappa .NE. 0. .AND. .NOT. implicitDiffusion ) THEN
371            WRITE(msgBuf,'(A,A)')
372         &  'S/R INI_PARMS: To use ivdc_kappa you must enable implicit',
373         &  ' vertical diffusion.'
374         CALL PRINT_ERROR( msgBuf , myThid)         CALL PRINT_ERROR( msgBuf , myThid)
375         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
376        ENDIF        ENDIF
377    
378        coordsSet = 0        coordsSet = 0
379        IF ( zCoordInputData ) coordsSet = coordsSet + 1        IF ( zCoordInputData ) coordsSet = coordsSet + 1
380        IF ( pCoordInputData ) coordsSet = coordsSet + 1        IF ( pCoordInputData ) coordsSet = coordsSet + 1
# Line 352  C--   z,p,r coord input switching. Line 393  C--   z,p,r coord input switching.
393        ELSE        ELSE
394         recip_rhoConst = 1.D0 / rhoConst         recip_rhoConst = 1.D0 / rhoConst
395        ENDIF        ENDIF
396          IF ( rhoNil .LE. 0. ) THEN
397           WRITE(msgBuf,'(A)')
398         &  'S/R INI_PARMS: rhoNil must be greater than 0.'
399           CALL PRINT_ERROR( msgBuf , myThid)
400           STOP 'ABNORMAL END: S/R INI_PARMS'
401          ELSE
402           recip_rhoNil = 1.D0 / rhoNil
403          ENDIF
404          IF ( HeatCapacity_Cp .LE. 0. ) THEN
405           WRITE(msgBuf,'(A)')
406         &  'S/R INI_PARMS: HeatCapacity_Cp must be greater than 0.'
407           CALL PRINT_ERROR( msgBuf , myThid)
408           STOP 'ABNORMAL END: S/R INI_PARMS'
409          ELSE
410           recip_Cp = 1.D0 / HeatCapacity_Cp
411          ENDIF
412          IF ( gravity .LE. 0. ) THEN
413           WRITE(msgBuf,'(A)')
414         &  'S/R INI_PARMS: gravity must be greater than 0.'
415           CALL PRINT_ERROR( msgBuf , myThid)
416           STOP 'ABNORMAL END: S/R INI_PARMS'
417          ELSE
418           recip_gravity = 1.D0 / gravity
419          ENDIF
420    C     Set globalFiles flag for READ_WRITE_FLD package
421          CALL SET_WRITE_GLOBAL_FLD( globalFiles )
422    C     Set globalFiles flag for READ_WRITE_REC package
423          CALL SET_WRITE_GLOBAL_REC( globalFiles )
424    C     Set globalFiles flag for READ_WRITE_REC package
425          CALL SET_WRITE_GLOBAL_PICKUP( globalFiles )
426    
427    C     Check for retired parameters still being used
428          nRetired = 0
429          IF ( zonal_filt_lat .NE. UNSET_RL ) THEN
430           nRetired = nRetired+1
431           WRITE(msgBuf,'(A,A)')
432         &  'S/R INI_PARMS: Paramater "zonal_filt_lat" is',
433         &  ' no longer allowed in file "data".'
434           CALL PRINT_ERROR( msgBuf , myThid)
435           WRITE(msgBuf,'(A,A)')
436         &  'S/R INI_PARMS: Paramater "zonal_filt_lat" is',
437         &  ' now read from file "data.zonfilt".'
438           CALL PRINT_ERROR( msgBuf , myThid)
439          ENDIF
440    
441  C--   Elliptic solver parameters  C--   Elliptic solver parameters
442        cg2dMaxIters       = 150        READ(UNIT=iUnit,NML=PARM02) !,IOSTAT=errIO)
443        cg2dTargetResidual = 1.D-7        IF ( errIO .LT. 0 ) THEN
       cg2dChkResFreq     = 1  
       cg2dpcOffDFac      = 0.51D0  
       READ(UNIT=iUnit,NML=PARM02,IOSTAT=errIO,err=5)  
       IF ( errIO .GE. 0 ) GOTO 6  
     5 CONTINUE  
444         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
445       &  'S/R INI_PARMS'       &  'S/R INI_PARMS'
446         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
# Line 375  C--   Elliptic solver parameters Line 455  C--   Elliptic solver parameters
455         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
456         CALL MODELDATA_EXAMPLE( myThid )         CALL MODELDATA_EXAMPLE( myThid )
457         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
458     6  CONTINUE        ELSE
459           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM02 : OK'
460           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
461         &                    SQUEEZE_RIGHT , 1)
462          ENDIF    
463    
464  C--   Time stepping parameters  C--   Time stepping parameters
       startTime         = 0.  
       nTimeSteps        = 0  
       endTime           = 0.  
       nIter0            = 0  
       deltaT            = 0.  
       deltaTClock       = 0.  
       deltaTtracer      = 0.  
       deltaTMom         = 0.  
       abEps             = 0.01  
       pchkPtFreq        = 0.  
       chkPtFreq         = 3600.*25  
       dumpFreq          = 3600.*100  
       taveFreq          = 0.  
       writeStatePrec    = precFloat64  
       nCheckLev         = 1  
       checkPtSuff(1)    = 'ckptA'  
       checkPtSuff(2)    = 'ckptB'  
       cAdjFreq          = -1.D0  
465        rCD               = -1.D0        rCD               = -1.D0
466        tauCD             =  0.D0        READ(UNIT=iUnit,NML=PARM03) !,IOSTAT=errIO)
467        tauThetaClimRelax =  0.D0        IF ( errIO .LT. 0 ) THEN
       doThetaClimRelax  = .FALSE.  
       tauSaltClimRelax  =  0.D0  
       doSaltClimRelax   = .FALSE.  
       periodicExternalForcing = .FALSE.  
       externForcingPeriod     = 0.  
       externForcingCycle      = 0.  
       READ(UNIT=iUnit,NML=PARM03,IOSTAT=errIO,err=7)  
       IF ( errIO .GE. 0 ) GOTO 8  
     7 CONTINUE  
468         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
469       &  'S/R INI_PARMS'       &  'S/R INI_PARMS'
470         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
# Line 422  C--   Time stepping parameters Line 479  C--   Time stepping parameters
479         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
480         CALL MODELDATA_EXAMPLE( myThid )         CALL MODELDATA_EXAMPLE( myThid )
481         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
482     8  CONTINUE        ELSE
483           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM03 : OK'
484           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
485         &                    SQUEEZE_RIGHT , 1)
486          ENDIF  
487  C     Process "timestepping" params  C     Process "timestepping" params
488  C     o Time step size  C     o Time step size
489        IF ( deltaT       .EQ. 0. ) deltaT       = deltaTmom        IF ( deltaT       .EQ. 0. ) deltaT       = deltaTmom
# Line 430  C     o Time step size Line 491  C     o Time step size
491        IF ( deltaTmom    .EQ. 0. ) deltaTmom    = deltaT        IF ( deltaTmom    .EQ. 0. ) deltaTmom    = deltaT
492        IF ( deltaTtracer .EQ. 0. ) deltaTtracer = deltaT        IF ( deltaTtracer .EQ. 0. ) deltaTtracer = deltaT
493        IF ( deltaTClock  .EQ. 0. ) deltaTClock  = deltaT        IF ( deltaTClock  .EQ. 0. ) deltaTClock  = deltaT
494    C Note that this line should set deltaTreesurf=deltaTmom
495    C but this would change a lot of existing set-ups so we are
496    C obliged to set the default inappropriately.
497    C Be advised that when using asynchronous time stepping
498    C it is better to set deltaTreesurf=deltaTtracer
499          IF ( deltaTfreesurf .EQ. 0. ) deltaTfreesurf  = deltaTmom
500        IF ( periodicExternalForcing ) THEN        IF ( periodicExternalForcing ) THEN
501         IF ( externForcingCycle*externForcingPeriod .EQ. 0. ) THEN         IF ( externForcingCycle*externForcingPeriod .EQ. 0. ) THEN
502          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
# Line 461  C     o Convection frequency Line 528  C     o Convection frequency
528        IF ( cAdjFreq .LT. 0. ) THEN        IF ( cAdjFreq .LT. 0. ) THEN
529         cAdjFreq = deltaTClock         cAdjFreq = deltaTClock
530        ENDIF        ENDIF
531          IF ( ivdc_kappa .NE. 0. .AND. cAdjFreq .NE. 0. ) THEN
532           WRITE(msgBuf,'(A,A)')
533         &  'S/R INI_PARMS: You have enabled both ivdc_kappa and',
534         &  ' convective adjustment.'
535           CALL PRINT_ERROR( msgBuf , myThid)
536           STOP 'ABNORMAL END: S/R INI_PARMS'
537          ENDIF
538  C     o CD coupling  C     o CD coupling
539        IF ( tauCD .EQ. 0.D0 ) THEN        IF ( tauCD .EQ. 0.D0 ) THEN
540          tauCD = deltaTmom          tauCD = deltaTmom
# Line 484  C     o Salinity climatology relaxation Line 558  C     o Salinity climatology relaxation
558         doSaltClimRelax     = .TRUE.         doSaltClimRelax     = .TRUE.
559         lambdaSaltClimRelax = 1./tauSaltClimRelax         lambdaSaltClimRelax = 1./tauSaltClimRelax
560        ENDIF        ENDIF
561  C     o Time step count  C     o Tracer 1 climatology relaxation time scale
562        IF ( endTime .NE. 0 ) THEN        IF ( tauTr1ClimRelax .EQ. 0.D0 ) THEN
563         IF ( deltaTClock .NE. 0 ) nTimeSteps =         doTr1ClimRelax     = .FALSE.
564       &  INT((endTime-startTime)/deltaTClock)         lambdaTr1ClimRelax = 0.D0
565          ELSE
566           doTr1ClimRelax     = .TRUE.
567           lambdaTr1ClimRelax = 1./tauTr1ClimRelax
568          ENDIF
569    
570    C     o Start time
571          IF ( nIter0 .NE. 0 .AND. startTime .EQ. 0. )
572         &   startTime = deltaTClock*float(nIter0)
573    C     o nIter0
574          IF ( nIter0 .EQ. 0 .AND. startTime .NE. 0. )
575         &   nIter0 = INT( startTime/deltaTClock )
576    
577    C     o nTimeSteps 1
578          IF ( nTimeSteps .EQ. 0 .AND. nEndIter .NE. 0 )
579         &     nTimeSteps = nEndIter-nIter0
580    C     o nTimeSteps 2
581          IF ( nTimeSteps .EQ. 0 .AND. endTime .NE. 0. )
582         &     nTimeSteps = int(0.5+(endTime-startTime)/deltaTclock)
583    C     o nEndIter 1
584          IF ( nEndIter .EQ. 0 .AND. nTimeSteps .NE. 0 )
585         &     nEndIter = nIter0+nTimeSteps
586    C     o nEndIter 2
587          IF ( nEndIter .EQ. 0 .AND. endTime .NE. 0. )
588         &     nEndIter = int(0.5+endTime/deltaTclock)
589    C     o End Time 1
590          IF ( endTime .EQ. 0. .AND. nTimeSteps .NE. 0 )
591         &     endTime = startTime + deltaTClock*float(nTimeSteps)
592    C     o End Time 2
593          IF ( endTime .EQ. 0. .AND. nEndIter .NE. 0 )
594         &     endTime = deltaTClock*float(nEndIter)
595    
596    C     o Consistent?
597          IF ( nEndIter .NE. nIter0+nTimeSteps ) THEN
598           WRITE(msgBuf,'(A)')
599         & 'S/R INI_PARMS: nIter0, nTimeSteps and nEndIter are inconsistent'
600           CALL PRINT_ERROR( msgBuf , 1)
601           WRITE(msgBuf,'(A)')
602         & 'S/R INI_PARMS: Perhaps more than two were set at once'
603           CALL PRINT_ERROR( msgBuf , 1)
604           STOP 'ABNORMAL END: S/R INI_PARMS'
605          ENDIF
606          IF ( nTimeSteps .NE. int(0.5+(endTime-startTime)/deltaTClock) )
607         & THEN
608            WRITE(msgBuf,'(A)')
609         &  'S/R INI_PARMS: both endTime and nTimeSteps have been set'
610            CALL PRINT_ERROR( msgBuf , 1)
611            WRITE(msgBuf,'(A)')
612         &  'S/R INI_PARMS: but are inconsistent'
613            CALL PRINT_ERROR( msgBuf , 1)
614            STOP 'ABNORMAL END: S/R INI_PARMS'
615          ENDIF
616    
617    C     o Monitor (should also add CPP flag for monitor?)
618          IF (monitorFreq.LT.0.) THEN
619           monitorFreq=0.
620           IF (dumpFreq.NE.0.) monitorFreq=dumpFreq
621           IF (diagFreq.NE.0..AND.diagFreq.LT.monitorFreq)
622         &         monitorFreq=diagFreq
623           IF (taveFreq.NE.0..AND.taveFreq.LT.monitorFreq)
624         &         monitorFreq=taveFreq
625           IF (chkPtFreq.NE.0..AND.chkPtFreq.LT.monitorFreq)
626         &         monitorFreq=chkPtFreq
627           IF (pChkPtFreq.NE.0..AND.pChkPtFreq.LT.monitorFreq)
628         &         monitorFreq=pChkPtFreq
629           IF (monitorFreq.EQ.0.) monitorFreq=deltaTclock
630        ENDIF        ENDIF
 C     o Finish time  
       IF ( endTime .EQ. 0. ) endTime = FLOAT(nTimeSteps)*deltaTClock  
631    
632  C     o If taveFreq is finite, then we must make sure the diagnostics  C     o If taveFreq is finite, then we must make sure the diagnostics
633  C       code is being compiled  C       code is being compiled
634  #ifndef ALLOW_DIAGNOSTICS  #ifndef ALLOW_TIMEAVE
635        IF (taveFreq.NE.0.) THEN        IF (taveFreq.NE.0.) THEN
636          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
637       &  'S/R INI_PARMS: taveFreq <> 0  but you have'       &  'S/R INI_PARMS: taveFreq <> 0  but you have'
# Line 502  C       code is being compiled Line 639  C       code is being compiled
639          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
640       &  'not compiled the model with the diagnostics routines.'       &  'not compiled the model with the diagnostics routines.'
641          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , 1)
642          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A,A)')
643       &  'Re-compile with:  #define ALLOW_DIAGNOSTICS  or  -DALLOW_DIAGNOSTICS'       &  'Re-compile with:  #define ALLOW_TIMEAVE',
644         &  '              or  -DALLOW_TIMEAVE'
645          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , 1)
646          STOP 'ABNORMAL END: S/R INI_PARMS'          STOP 'ABNORMAL END: S/R INI_PARMS'
647        ENDIF        ENDIF
# Line 511  C       code is being compiled Line 649  C       code is being compiled
649    
650  C--   Grid parameters  C--   Grid parameters
651  C     In cartesian coords distances are in metres  C     In cartesian coords distances are in metres
652        rkFac = UNSET_I        rkFac = UNSET_RS
       usingCartesianGrid = .TRUE.  
       delRDefault =  1.D2  
653        DO K =1,Nr        DO K =1,Nr
654         delZ(K) = UNSET_RL         delZ(K) = UNSET_RL
655         delP(K) = UNSET_RL         delP(K) = UNSET_RL
656         delR(K) = UNSET_RL         delR(K) = UNSET_RL
657        ENDDO        ENDDO
       dxSpacing = 20.D0 * 1000.D0  
       dySpacing = 20.D0 * 1000.D0  
       DO i=1,Nx  
        delX(i) = dxSpacing  
       ENDDO  
       DO j=1,Ny  
        delY(j) = dySpacing  
       ENDDO  
658  C     In spherical polar distances are in degrees  C     In spherical polar distances are in degrees
       usingSphericalPolarGrid = .FALSE.  
       phiMin    = -5.0  
       thetaMin  = 0.  
       rSphere   = 6370. * 1.D3  
659        recip_rSphere  = 1.D0/rSphere        recip_rSphere  = 1.D0/rSphere
660        IF ( usingSphericalPolarGrid ) THEN        dxSpacing = UNSET_RL
661         dxSpacing = 1.        dySpacing = UNSET_RL
662         dySpacing = 1.        delXfile = ' '
663         DO I=1,Nx        delYfile = ' '
664          delX(I) = dxSpacing        READ(UNIT=iUnit,NML=PARM04,IOSTAT=errIO)
665         ENDDO        IF ( errIO .LT. 0 ) THEN
        DO J=1,Ny  
         delY(J) = dySpacing  
        ENDDO  
       ENDIF  
   
       READ(UNIT=iUnit,NML=PARM04,IOSTAT=errIO,err=9)  
       IF ( errIO .GE. 0 ) GOTO 10  
     9 CONTINUE  
666         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
667       &  'S/R INI_PARMS'       &  'S/R INI_PARMS'
668         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
# Line 561  C     In spherical polar distances are i Line 677  C     In spherical polar distances are i
677         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
678         CALL MODELDATA_EXAMPLE( myThid )         CALL MODELDATA_EXAMPLE( myThid )
679         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
   10  CONTINUE  
 C  
       IF ( rSphere .NE. 0 ) THEN  
        recip_rSphere = 1.D0/rSphere  
680        ELSE        ELSE
681         recip_rSphere = 0.         WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM04 : OK'
682        ENDIF         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
683  C--   Initialize EOS coefficients (3rd order polynomial)       &                    SQUEEZE_RIGHT , 1)
684        IF (eostype.eq.'POLY3') THEN        ENDIF    
685         OPEN(37,FILE='POLY3.COEFFS',STATUS='OLD',FORM='FORMATTED')  
686         READ(37,*) I  C     X coordinate
687         IF (I.NE.Nr) THEN        IF ( delXfile .NE. ' ' ) THEN
688          WRITE(msgBuf,'(A)')         IF ( delX(1) .NE. UNSET_RL .OR. dxSpacing .NE. UNSET_RL ) THEN
689       &  'ini_parms: attempt to read POLY3.COEFFS failed'           WRITE(msgBuf,'(A,A)') 'Too many specifications for delX:',
690         &   'Specify only one of delX, dxSpacing or delXfile'
691          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , 1)
692          WRITE(msgBuf,'(A)')          STOP 'ABNORMAL END: S/R INI_PARMS'
693       &  '           because bad # of levels in data'         ELSE
694            _BEGIN_MASTER( myThid )
695            IF (readBinaryPrec.EQ.precFloat32) THEN
696             OPEN(37,FILE=delXfile,STATUS='OLD',FORM='UNFORMATTED',
697         &        ACCESS='DIRECT',RECL=WORDLENGTH*Nx)
698             READ(37,rec=1) delX
699    #ifdef _BYTESWAPIO
700                call MDS_BYTESWAPR4( Nx, delX )
701    #endif
702             CLOSE(37)
703            ELSEIF (readBinaryPrec.EQ.precFloat64) THEN
704             OPEN(37,FILE=delXfile,STATUS='OLD',FORM='UNFORMATTED',
705         &        ACCESS='DIRECT',RECL=WORDLENGTH*2*Nx)
706             READ(37,rec=1) delX
707    #ifdef _BYTESWAPIO
708                call MDS_BYTESWAPR8( Nx, delX )
709    #endif
710             CLOSE(37)
711            ENDIF
712            _END_MASTER(myThid)
713           ENDIF
714          ENDIF
715          IF ( dxSpacing .NE. UNSET_RL ) THEN
716           DO i=1,Nx
717            delX(i) = dxSpacing
718           ENDDO
719          ENDIF
720    C     Y coordinate
721          IF ( delYfile .NE. ' ' ) THEN
722           IF ( delY(1) .NE. UNSET_RL .OR. dySpacing .NE. UNSET_RL ) THEN
723             WRITE(msgBuf,'(A,A)') 'Too many specifications for delY:',
724         &   'Specify only one of delY, dySpacing or delYfile'
725          CALL PRINT_ERROR( msgBuf , 1)          CALL PRINT_ERROR( msgBuf , 1)
726          STOP 'Bad data in POLY3.COEFFS'          STOP 'ABNORMAL END: S/R INI_PARMS'
727           ELSE
728            _BEGIN_MASTER( myThid )
729            IF (readBinaryPrec.EQ.precFloat32) THEN
730             OPEN(37,FILE=delYfile,STATUS='OLD',FORM='UNFORMATTED',
731         &        ACCESS='DIRECT',RECL=WORDLENGTH*Ny)
732             READ(37,rec=1) delY
733    #ifdef _BYTESWAPIO
734                call MDS_BYTESWAPR4( Ny, delY )
735    #endif
736             CLOSE(37)
737            ELSEIF (readBinaryPrec.EQ.precFloat64) THEN
738             OPEN(37,FILE=delYfile,STATUS='OLD',FORM='UNFORMATTED',
739         &        ACCESS='DIRECT',RECL=WORDLENGTH*2*Ny)
740             READ(37,rec=1) delY
741    #ifdef _BYTESWAPIO
742                call MDS_BYTESWAPR8( Ny, delY )
743    #endif
744             CLOSE(37)
745            ENDIF
746            _END_MASTER(myThid)
747         ENDIF         ENDIF
748         READ(37,*) (eosRefT(K),eosRefS(K),eosSig0(K),K=1,Nr)        ENDIF
749         DO K=1,Nr        IF ( dySpacing .NE. UNSET_RL ) THEN
750          READ(37,*) (eosC(I,K),I=1,9)         DO i=1,Ny
751            delY(i) = dySpacing
752         ENDDO         ENDDO
753         CLOSE(37)        ENDIF
754    C
755          IF ( rSphere .NE. 0 ) THEN
756           recip_rSphere = 1.D0/rSphere
757          ELSE
758           recip_rSphere = 0.
759        ENDIF        ENDIF
760  C--   Check for conflicting grid definitions.  C--   Check for conflicting grid definitions.
761        goptCount = 0        goptCount = 0
762        IF ( usingCartesianGrid )      goptCount = goptCount+1        IF ( usingCartesianGrid )      goptCount = goptCount+1
763        IF ( usingSphericalPolarGrid ) goptCount = goptCount+1        IF ( usingSphericalPolarGrid ) goptCount = goptCount+1
764        IF ( goptCount .NE. 1 ) THEN        IF ( usingCurvilinearGrid )    goptCount = goptCount+1
765          IF ( goptCount .GT. 1 ) THEN
766         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
767       &  'S/R INI_PARMS: More than one coordinate system requested'       &  'S/R INI_PARMS: More than one coordinate system requested'
768         CALL PRINT_ERROR( msgBuf , myThid)         CALL PRINT_ERROR( msgBuf , myThid)
769         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
770        ENDIF        ENDIF
771          IF ( goptCount .LT. 1 ) THEN
772           WRITE(msgBuf,'(A)')
773         &  'S/R INI_PARMS: No coordinate system requested'
774           CALL PRINT_ERROR( msgBuf , myThid)
775           STOP 'ABNORMAL END: S/R INI_PARMS'
776          ENDIF
777  C--   Make metric term settings consistent with underlying grid.  C--   Make metric term settings consistent with underlying grid.
778        IF ( usingCartesianGrid ) THEN        IF ( usingCartesianGrid ) THEN
779         usingSphericalPolarMterms = .FALSE.         usingSphericalPolarMterms = .FALSE.
780         metricTerms = .FALSE.         metricTerms = .FALSE.
781         mTFacMom = 0         mTFacMom = 0.
782         useBetaPlaneF = .TRUE.         useBetaPlaneF = .TRUE.
783        ENDIF        ENDIF
784        IF ( usingSphericalPolarGrid ) THEN        IF ( usingSphericalPolarGrid ) THEN
785         useConstantF  = .FALSE.         useConstantF  = .FALSE.
786         useBetaPlaneF = .FALSE.         useBetaPlaneF = .FALSE.
787         useSphereF    = .TRUE.         useSphereF    = .TRUE.
788         omega         = 2.D0 * PI / ( 3600.D0 * 24.D0 )         usingSphericalPolarMterms = metricTerms
        usingSphericalPolarMterms = .TRUE.  
        metricTerms = .TRUE.  
        mTFacMom = 1  
789        ENDIF        ENDIF
790          IF ( usingCurvilinearGrid ) THEN
791           useSphereF    = .TRUE.
792          ENDIF
793    C--   set cell Center depth and put Interface at the middle between 2 C
794          setCenterDr = .FALSE.
795          IF (delRc(1).NE.UNSET_RL) setCenterDr=.TRUE.
796          DO K=1,Nr+1
797            IF (delRc(K).EQ.UNSET_RL) setCenterDr = .FALSE.
798          ENDDO
799  C--   p, z, r coord parameters  C--   p, z, r coord parameters
800        DO K = 1, Nr        DO K = 1, Nr
801         IF ( delZ(K) .NE. UNSET_RL ) zCoordInputData = .TRUE.         IF ( delZ(K) .NE. UNSET_RL ) zCoordInputData = .TRUE.
# Line 620  C--   p, z, r coord parameters Line 803  C--   p, z, r coord parameters
803         IF ( delR(K) .NE. UNSET_RL ) rCoordInputData = .TRUE.         IF ( delR(K) .NE. UNSET_RL ) rCoordInputData = .TRUE.
804         IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delZ(K)         IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delZ(K)
805         IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delP(K)         IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delP(K)
806         IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delRDefault         IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delRDefault(K)
807           IF (.NOT.setCenterDr .AND. delR(K).EQ.delRDefault(K) ) THEN
808             WRITE(msgBuf,'(A,I4)')
809         &  'S/R INI_PARMS: No value for delZ/delP/delR at K = ',K
810             CALL PRINT_ERROR( msgBuf , 1)
811             STOP 'ABNORMAL END: S/R INI_PARMS'
812           ELSEIF ( setCenterDr .AND. delR(K).NE.delRDefault(K) ) THEN
813             WRITE(msgBuf,'(2A,I4)') 'S/R INI_PARMS:',
814         &    ' Cannot specify both delRc and delZ/delP/delR at K=', K
815             CALL PRINT_ERROR( msgBuf , 1)
816             STOP 'ABNORMAL END: S/R INI_PARMS'
817           ENDIF
818        ENDDO        ENDDO
819  C     Check for multiple coordinate systems  C     Check for multiple coordinate systems
820        coordsSet = 0        CoordsSet = 0
821        IF ( zCoordInputData ) coordsSet = coordsSet + 1        IF ( zCoordInputData ) coordsSet = coordsSet + 1
822        IF ( pCoordInputData ) coordsSet = coordsSet + 1        IF ( pCoordInputData ) coordsSet = coordsSet + 1
823        IF ( rCoordInputData ) coordsSet = coordsSet + 1        IF ( rCoordInputData ) coordsSet = coordsSet + 1
# Line 635  C     Check for multiple coordinate syst Line 829  C     Check for multiple coordinate syst
829        ENDIF        ENDIF
830    
831  C--   Input files  C--   Input files
832        bathyFile       = ' '        READ(UNIT=iUnit,NML=PARM05) !,IOSTAT=errIO)
833        hydrogSaltFile  = ' '        IF ( errIO .LT. 0 ) THEN    
       hydrogThetaFile = ' '  
       zonalWindFile   = ' '  
       meridWindFile   = ' '  
       thetaClimFile   = ' '  
       saltClimFile    = ' '  
       READ(UNIT=iUnit,NML=PARM05,IOSTAT=errIO,err=11)  
       IF ( errIO .GE. 0 ) GOTO 12  
    11 CONTINUE  
        WRITE(msgBuf,'(A)')  
      &  'S/R INI_PARMS'  
        CALL PRINT_ERROR( msgBuf , 1)  
834         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
835       &  'Error reading numerical model '       &  'Error reading numerical model '
836         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
# Line 659  C--   Input files Line 842  C--   Input files
842         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
843         CALL MODELDATA_EXAMPLE( myThid )         CALL MODELDATA_EXAMPLE( myThid )
844         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
845    12  CONTINUE        ELSE
846           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM05 : OK'
847           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
848         &                    SQUEEZE_RIGHT , 1)
849          ENDIF      
850    C     o If pLoadFile is set, then we should make sure the corresponing
851    C       code is being compiled
852    #ifndef ATMOSPHERIC_LOADING
853          IF (pLoadFile.NE.' ') THEN
854            WRITE(msgBuf,'(A)')
855         &  'S/R INI_PARMS: pLoadFile is set but you have not'
856            CALL PRINT_ERROR( msgBuf , 1)
857            WRITE(msgBuf,'(A)')
858         &  'compiled the model with the pressure loading code.'
859            CALL PRINT_ERROR( msgBuf , 1)
860            WRITE(msgBuf,'(A,A)')
861         &  'Re-compile with:  #define ATMOSPHERIC_LOADING',
862         &  '              or  -DATMOSPHERIC_LOADING'
863            CALL PRINT_ERROR( msgBuf , 1)
864            STOP 'ABNORMAL END: S/R INI_PARMS'
865          ENDIF
866    #endif
867    
868  C  C
869  C--   Set factors required for mixing pressure and meters as vertical coordinate.  C--   Set factors required for mixing pressure and meters as vertical coordinate.
# Line 670  C     rkFac = -1 applies when K and the Line 874  C     rkFac = -1 applies when K and the
874  C     horiVertRatio is a parameter that maps horizontal units to vertical units.  C     horiVertRatio is a parameter that maps horizontal units to vertical units.
875  C     It is used in certain special cases where lateral and vertical terms are  C     It is used in certain special cases where lateral and vertical terms are
876  C     being combined and a single frame of reference is needed.  C     being combined and a single frame of reference is needed.
877        IF ( zCoordInputData .AND. rkFac .EQ. UNSET_I ) THEN        IF ( zCoordInputData .AND. rkFac .EQ. UNSET_RS ) THEN
878         rkFac       = 1.D0         rkFac       = 1.D0
879         horiVertRatio = 1.D0         horiVertRatio = 1.D0
880        ENDIF        ENDIF
881        IF ( pCoordInputData .AND. rkFac .EQ. UNSET_I ) THEN        IF ( pCoordInputData .AND. rkFac .EQ. UNSET_RS ) THEN
882         rkFac = -1.D0         rkFac = -1.D0
883         horiVertRatio = Gravity * rhoConst         horiVertRatio = Gravity * rhoConst
884        ENDIF        ENDIF
885        IF ( rCoordInputData .AND. rkFac .EQ. UNSET_I ) THEN        IF ( rCoordInputData .AND. rkFac .EQ. UNSET_RS ) THEN
886         rkFac =  1.D0         rkFac =  1.D0
887         horiVertRatio = 1.D0         horiVertRatio = 1.D0
888        ENDIF        ENDIF
889          IF (buoyancyRelation.EQ.'ATMOSPHERIC')
890         &   horiVertRatio = Gravity * rhoConst                                    
891          IF (buoyancyRelation.EQ.'OCEANICP')
892         &   horiVertRatio = Gravity * rhoConst                                    
893          IF ( rkFac .EQ. UNSET_RS ) rkFac=rkFacDefault
894        recip_rkFac = 1.D0 / rkFac        recip_rkFac = 1.D0 / rkFac
895          recip_horiVertRatio = 1./horiVertRatio
896        IF ( zCoordInputData ) usingZCoords = .TRUE.        IF ( zCoordInputData ) usingZCoords = .TRUE.
897        IF ( pCoordInputData ) usingPCoords = .TRUE.        IF ( pCoordInputData ) usingPCoords = .TRUE.
898    
899  C  C
900        CLOSE(iUnit)        CLOSE(iUnit)
901    
902    C--   Check whether any retired parameters were found.
903    C--   Stop if they were
904          IF ( nRetired .GT. 0 ) THEN    
905           WRITE(msgBuf,'(A)')
906         &  'Error reading '
907           CALL PRINT_ERROR( msgBuf , 1)
908           WRITE(msgBuf,'(A)')
909         &  'parameter file "data"'
910           CALL PRINT_ERROR( msgBuf , 1)
911           WRITE(msgBuf,'(A)')
912         &  'some out of date parameters were found in the namelist'
913           CALL PRINT_ERROR( msgBuf , 1)
914           STOP 'ABNORMAL END: S/R INI_PARMS'
915          ENDIF
916    
917        _END_MASTER(myThid)        _END_MASTER(myThid)
918    
919  C--   Everyone else must wait for the parameters to be loaded  C--   Everyone else must wait for the parameters to be loaded
920        _BARRIER        _BARRIER
921  C  C
   
922        RETURN        RETURN
923        END        END
924    
Legend:
Removed from v.1.31  
changed lines
  Added in v.1.83
  ViewVC Help
Powered by ViewVC 1.1.22