File indexing completed on 2024-02-01 06:10:08 UTC
view on githubraw file Latest commit 427e24e1 on 2024-01-31 16:50:14 UTC
b9b591469d Jean*0001 #include "PACKAGES_CONFIG.h"
b9377e56be Alis*0002 #include "CPP_OPTIONS.h"
b9b591469d Jean*0003 #ifdef ALLOW_EXCH2
0004 # include "W2_OPTIONS.h"
0005 #endif /* ALLOW_EXCH2 */
924557e60a Chri*0006
a30418b6b9 Ed H*0007
9366854e02 Chri*0008
427e24e121 Jean*0009
0010
924557e60a Chri*0011 SUBROUTINE INI_PARMS( myThid )
0012
427e24e121 Jean*0013
809bdccbfc Jean*0014
0015
9366854e02 Chri*0016
0017 IMPLICIT NONE
924557e60a Chri*0018 #include "SIZE.h"
0019 #include "EEPARAMS.h"
0020 #include "PARAMS.h"
b9b591469d Jean*0021 #ifdef ALLOW_EXCH2
0022 # include "W2_EXCH2_SIZE.h"
0023 # include "W2_EXCH2_TOPOLOGY.h"
0024 #endif /* ALLOW_EXCH2 */
a37a13034c Mart*0025 #include "EOS.h"
97168926d1 Jean*0026
2ad79bdf32 Jean*0027 #include "GRID.h"
97168926d1 Jean*0028 #include "SET_GRID.h"
924557e60a Chri*0029
427e24e121 Jean*0030
0031
924557e60a Chri*0032 INTEGER myThid
0033
427e24e121 Jean*0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044
0045
0046
0047
910f05e765 Chri*0048 _RL dxSpacing
0049 _RL dySpacing
062a876ce5 Jean*0050 _RL deltaTtracer
924557e60a Chri*0051 CHARACTER*(MAX_LEN_MBUF) msgBuf
c2b6ed6bfd Jean*0052 LOGICAL forcing_In_AB
924557e60a Chri*0053 INTEGER goptCount
b9b591469d Jean*0054 INTEGER gridNx, gridNy
aa8f58c33d Jean*0055 INTEGER k, i, j, iUnit
90aed06079 Jean*0056 INTEGER errIO, errCount
427e24e121 Jean*0057
0058
910f05e765 Chri*0059 _RL viscArDefault
0060 _RL diffKrTDefault
0061 _RL diffKrSDefault
0062 _RL hFacMinDrDefault
0127add478 Alis*0063 _RL delRDefault(Nr)
427e24e121 Jean*0064
0065
0066
0067
0068
0069
0070
0071
0072
0073
0074
0075
0076
0077
0078
0079
0080
53f30fd0a1 Jean*0081
910f05e765 Chri*0082 LOGICAL zCoordInputData
0083 LOGICAL pCoordInputData
0084 LOGICAL rCoordInputData
0085 INTEGER coordsSet
7cd7754949 Jean*0086 INTEGER vertSetCount
13b4d13c98 Jean*0087 LOGICAL specifiedDiffKrT, specifiedDiffKrS
ede2cfe569 Chri*0088
427e24e121 Jean*0089
0090
0091
0092
0093
0094
0095
0096
0097
53f30fd0a1 Jean*0098 _RL delZ(Nr)
0099 _RL delP(Nr)
0100 _RL viscAz
0101 _RL viscAp
7cd7754949 Jean*0102 _RL viscAr
53f30fd0a1 Jean*0103 _RL diffKzT
0104 _RL diffKpT
8fc5af62b8 Jean*0105 _RL diffKrT
53f30fd0a1 Jean*0106 _RL diffKzS
0107 _RL diffKpS
8fc5af62b8 Jean*0108 _RL diffKrS
53f30fd0a1 Jean*0109
427e24e121 Jean*0110
0111
0112
0113
742cf4499c Jean*0114
427e24e121 Jean*0115
0116
0117
0118
0119
0120
0121
0122
0123
0124
0125
0126
0127
0128
0129
0130
0131
0132
0133
0134
0135
0136
0137
0138
0139
0140
0141
0142
aecc8b0f47 Mart*0143
427e24e121 Jean*0144
0145
0146
0147
742cf4499c Jean*0148
427e24e121 Jean*0149
0150
0151
742cf4499c Jean*0152
427e24e121 Jean*0153
0154
0155
0156
0157
0158
8194bc4a99 Mart*0159
427e24e121 Jean*0160
0161
0162
0163
8194bc4a99 Mart*0164
9780090eaa Jean*0165 INTEGER nRetired
c07cd3bfa8 Jean*0166 LOGICAL useConstantF, useBetaPlaneF, useSphereF
427e24e121 Jean*0167 LOGICAL useJamartWetPoints, useEnergyConservingCoriolis
7e00d7e8f9 Jean*0168 LOGICAL SadournyCoriolis
427e24e121 Jean*0169 LOGICAL use3dCoriolis, metricTerms
7a4ec66481 Jean*0170 LOGICAL tempDiffusion, saltDiffusion
0171 INTEGER tracerAdvScheme
0172 _RL trac_EvPrRn
2ad79bdf32 Jean*0173 _RL zonal_filt_lat
0174
8f934af0f4 Jean*0175 _RL viscAstrain, viscAtension
6d61b52b09 Jean*0176 LOGICAL useAnisotropicViscAgridMax
0177 LOGICAL usePickupBeforeC35
83513ae28b Jean*0178 LOGICAL saveDebugMode
e793395f17 Jean*0179 LOGICAL allowInteriorFreezing, useOldFreezing
7e00d7e8f9 Jean*0180 LOGICAL balanceEmPmR
427e24e121 Jean*0181 INTEGER writeStatePrec
742cf4499c Jean*0182
aecc8b0f47 Mart*0183 INTEGER cg2dChkResFreq, cg3dChkResFreq
0184
46979d29da Jean*0185 _RL tauThetaClimRelax3Dim, tauSaltClimRelax3Dim
742cf4499c Jean*0186 LOGICAL calendarDumps
0187
0188 LOGICAL groundAtK1
2ad79bdf32 Jean*0189 _RL Ro_SeaLevel
742cf4499c Jean*0190 _RL rkFac
9780090eaa Jean*0191 _RL thetaMin, phiMin
7e00d7e8f9 Jean*0192
8194bc4a99 Mart*0193 CHARACTER*(MAX_LEN_FNAM) shelfIceFile
275ad78b78 Jean*0194 CHARACTER*(MAX_LEN_FNAM) dQdTfile
924557e60a Chri*0195
0196
0197 NAMELIST /PARM01/
fe4f35d666 Alis*0198 & gravitySign, nh_Am2,
7fbd8d8c1c Jean*0199 & gravity, gBaro, gravityFile, rhonil, tAlpha, sBeta,
c07cd3bfa8 Jean*0200 & selectCoriMap, f0, beta, fPrime, omega, rotationPeriod,
8f106aecfd Bayl*0201 & viscAh, viscAhW, viscAhMax,
0202 & viscAhGrid, viscAhGridMax, viscAhGridMin,
dc3adfb09b Jean*0203 & viscC2leith, viscC4leith, smag3D_coeff, useSmag3D,
6d61b52b09 Jean*0204 & useFullLeith, useAnisotropicViscAgridMax, useStrainTensionVisc,
f59d76b0dd Ed D*0205 & useAreaViscLength, viscC2leithD, viscC4leithD, viscC2LeithQG,
0206 & viscC2smag, viscC4smag, viscAhD, viscAhZ, viscA4D, viscA4Z,
0b6cbae535 Mart*0207 & viscA4, viscA4W,
0208 & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,
1e9e8678dd Bayl*0209 & viscA4ReMax, viscAhReMax,
7cd7754949 Jean*0210 & cosPower, viscAstrain, viscAtension,
427e24e121 Jean*0211 & viscAz, diffKzT, diffKzS, viscAp, diffKpT, diffKpS,
0212 & viscAr, diffKrT, diffKrS, viscArNr, diffKrNrT, diffKrNrS,
2d5bb917cc Jean*0213 & diffKhT, diffK4T, diffKhS, diffK4S, smag3D_diffCoeff,
427e24e121 Jean*0214 & diffKr4T, diffKr4S, BL79LatVary,
0215 & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
0216 & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,
9330bd8273 Jean*0217 & surf_pRef, tRef, sRef, tRefFile, sRefFile, rhoRefFile,
7fbd8d8c1c Jean*0218 & eosType, selectP_inEOS_Zc, integr_GeoPot, selectFindRoSurf,
aad87ebb4c Jean*0219 & HeatCapacity_Cp, celsius2K, atm_Cp, atm_Rd, atm_Rq, atm_Po,
59b35dd864 Jean*0220 & no_slip_sides, sideDragFactor, no_slip_bottom, bottomVisc_pCell,
ab47de63dc Mart*0221 & bottomDragLinear, bottomDragQuadratic, zRoughBot,
0222 & selectBotDragQuadr,
427e24e121 Jean*0223 & momPressureForcing, momForcing, momTidalForcing,
0224 & momViscosity, momAdvection, vectorInvariantMomentum,
0225 & useConstantF, useBetaPlaneF, useSphereF, useCoriolis,
0226 & use3dCoriolis, select3dCoriScheme, selectCoriScheme,
0227 & useCDscheme, useJamartWetPoints, useEnergyConservingCoriolis,
0228 & useAbsVorticity, selectVortScheme, SadournyCoriolis,
0229 & useJamartMomAdv, upwindVorticity, highOrderVorticity,
0230 & upwindShear, selectKEscheme,
0231 & selectMetricTerms, metricTerms, useNHMTerms, addFrictionHeating,
46918f1b26 Jean*0232 & tempDiffusion, tempAdvection, tempForcing, temp_stayPositive,
0233 & saltDiffusion, saltAdvection, saltForcing, salt_stayPositive,
9e44938df7 Jean*0234 & implicSurfPress, implicDiv2DFlow, implicitNHPress,
7cd7754949 Jean*0235 & implicitFreeSurface, rigidLid, freeSurfFac,
0236 & hFacMin, hFacMinDz, hFacMinDp, hFacMinDr,
1e6181f584 Davi*0237 & exactConserv, linFSConserveTr, uniformLin_PhiSurf,
0238 & nonlinFreeSurf, hFacInf, hFacSup, select_rStar,
3fcd8a21e5 Jean*0239 & nonHydrostatic, selectNHfreeSurf, quasiHydrostatic,
627a21a418 Jean*0240 & implicitIntGravWave, staggerTimeStep, doResetHFactors,
8a28092f34 Patr*0241 & tempStepping, saltStepping, momStepping,
9e44938df7 Jean*0242 & implicitDiffusion, implicitViscosity, selectImplicitDrag,
b9d351b225 Jean*0243 & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
aad87ebb4c Jean*0244 & rhoConst, thetaConst, rhoConstFresh, buoyancyRelation,
7c642b2a2b Dimi*0245 & allowFreezing, allowInteriorFreezing, useOldFreezing, ivdc_kappa,
d1b0368d70 Davi*0246 & hMixCriteria, dRhoSmall, hMixSmooth,
427e24e121 Jean*0247 & tempAdvScheme, tempVertAdvScheme, tracerAdvScheme,
0248 & saltAdvScheme, saltVertAdvScheme, multiDimAdvection,
d18df35fee Jean*0249 & selectAddFluid, useRealFreshWaterFlux, convertFW2Salt,
745a0098ab Jean*0250 & temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
80d98e0151 Dimi*0251 & temp_addMass, salt_addMass, zonal_filt_lat,
0320e25227 Mart*0252 & smoothAbsFuncRange, sIceLoadFac,
7e00d7e8f9 Jean*0253 & selectBalanceEmPmR, balanceEmPmR, balanceQnet, balancePrintMean,
427e24e121 Jean*0254 & balanceThetaClimRelax, balanceSaltClimRelax,
0255 & readBinaryPrec, writeBinaryPrec, writeStatePrec, globalFiles,
0256 & useSingleCpuIO, useSingleCpuInput, usePickupBeforeC54,
0257 & usePickupBeforeC35, debugMode, debugLevel, plotLevel
924557e60a Chri*0258
0259
0260 NAMELIST /PARM02/
aecc8b0f47 Mart*0261 & cg2dMaxIters, cg2dMinItersNSA, cg2dChkResFreq, cg2dUseMinResSol,
b46f9da188 Jean*0262 & cg2dTargetResidual, cg2dTargetResWunit,
0263 & cg2dpcOffDFac, cg2dPreCondFreq,
23bce0bbb8 Mart*0264 & cg3dMaxIters, cg3dChkResFreq, cg3dTargetResidual,
aecc8b0f47 Mart*0265 & useNSACGSolver, useSRCGSolver, printResidualFreq
924557e60a Chri*0266
0267
0268 NAMELIST /PARM03/
6e6f314aa7 Patr*0269 & nIter0, nTimeSteps, nTimeSteps_l2, nEndIter,
f804abbd25 Jean*0270 & baseTime, startTime, endTime,
db8d49beca Jean*0271 & deltaT, deltaTClock, deltaTMom,
0272 & deltaTtracer, dTtracerLev, deltaTFreeSurf,
b9306711eb Jean*0273 & forcing_In_AB, momForcingOutAB, tracForcingOutAB,
c2b6ed6bfd Jean*0274 & momDissip_In_AB, doAB_onGtGs,
2a618ca8c5 Jean*0275 & abEps, alph_AB, beta_AB, startFromPickupAB2, applyExchUV_early,
8039e9b985 Jean*0276 & tauCD, rCD, epsAB_CD, cAdjFreq,
f804abbd25 Jean*0277 & chkPtFreq, pChkPtFreq, pickupSuff, pickupStrictlyMatch,
8c73a5b228 Mart*0278 & writePickupAtEnd,
82d0948361 Jean*0279 & dumpFreq, dumpInitAndLast, adjDumpFreq, taveFreq, tave_lastIter,
f804abbd25 Jean*0280 & diagFreq, monitorFreq, adjMonitorFreq, monitorSelect,
2e3729af6b Jean*0281 & outputTypesInclusive, rwSuffixType,
b9306711eb Jean*0282 & tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax,
a31e157718 Jean*0283 & tauThetaClimRelax3Dim, tauSaltClimRelax3Dim,
fc64655a31 Dimi*0284 & periodicExternalForcing, externForcingPeriod, externForcingCycle,
82d0948361 Jean*0285 & calendarDumps
924557e60a Chri*0286
0287
0288 NAMELIST /PARM04/
b9306711eb Jean*0289 & usingCartesianGrid, usingCylindricalGrid,
9780090eaa Jean*0290 & usingSphericalPolarGrid, usingCurvilinearGrid,
0291 & xgOrigin, ygOrigin, dxSpacing, dySpacing,
0292 & delX, delY, delXFile, delYFile, horizGridFile,
0293 & phiEuler, thetaEuler, psiEuler,
2ad79bdf32 Jean*0294 & rSphere, radius_fromHorizGrid, deepAtmosphere, seaLev_Z,
0295 & top_Pres, delZ, delP, delR, delRc, delRFile, delRcFile,
bb2fd3f1ad Jean*0296 & useMin4hFacEdges, interViscAr_pCell, interDiffKr_pCell,
0297 & pCellMix_select, pCellMix_maxFac, pCellMix_delR,
0298 & pCellMix_viscAr, pCellMix_diffKr,
0299 & selectSigmaCoord, rSigmaBnd, hybSigmFile,
2ad79bdf32 Jean*0300 & Ro_SeaLevel, rkFac, groundAtK1, thetaMin, phiMin
924557e60a Chri*0301
81bc00c2f0 Chri*0302
0303 NAMELIST /PARM05/
26565fca80 Jean*0304 & bathyFile, topoFile, addWwallFile, addSwallFile, shelfIceFile,
0305 & diffKrFile, viscAhDfile, viscAhZfile, viscA4Dfile, viscA4Zfile,
4c5bb1c88e Jean*0306 & hydrogThetaFile, hydrogSaltFile,
0307 & maskIniTemp, maskIniSalt, checkIniTemp, checkIniSalt,
88830be691 Alis*0308 & zonalWindFile, meridWindFile,
0309 & thetaClimFile, saltClimFile,
1e273d1bf5 Jean*0310 & surfQfile, surfQnetFile, surfQswFile, EmPmRfile, saltFluxFile,
2dcaa8b9a5 Patr*0311 & uVelInitFile, vVelInitFile, pSurfInitFile,
0320e25227 Mart*0312 & dQdTFile, ploadFile, geoPotAnomFile, addMassFile,
0313 & tCylIn, tCylOut,
90929f8806 Patr*0314 & eddyPsiXFile, eddyPsiYFile, geothermalFile,
7e00d7e8f9 Jean*0315 & lambdaThetaFile, lambdaSaltFile, wghtBalanceFile,
cbc417d429 Mart*0316 & mdsioLocalDir, adTapeDir,
9a263a84a8 Ed H*0317 & the_run_name
a30418b6b9 Ed H*0318
81bc00c2f0 Chri*0319
b9b591469d Jean*0320 #ifdef ALLOW_EXCH2
0321 gridNx = exch2_mydNx(1)
0322 gridNy = exch2_mydNy(1)
0323 #else /* ALLOW_EXCH2 */
0324 gridNx = Nx
0325 gridNy = Ny
0326 #endif /* ALLOW_EXCH2 */
0327
924557e60a Chri*0328 _BEGIN_MASTER(myThid)
0329
0127add478 Alis*0330
0331 CALL SET_DEFAULTS(
0332 O viscArDefault, diffKrTDefault, diffKrSDefault,
f15994caab Jean*0333 O hFacMinDrDefault, delRDefault,
0127add478 Alis*0334 I myThid )
7843dde2de jm-c 0335
427e24e121 Jean*0336 useJamartWetPoints = .FALSE.
7843dde2de jm-c 0337 useEnergyConservingCoriolis = .FALSE.
e35a804907 Jean*0338 SadournyCoriolis = .FALSE.
427e24e121 Jean*0339 use3dCoriolis = .TRUE.
0340 metricTerms = .TRUE.
7e00d7e8f9 Jean*0341 balanceEmPmR = .FALSE.
0127add478 Alis*0342
910f05e765 Chri*0343
0344 zCoordInputData = .FALSE.
0345 pCoordInputData = .FALSE.
0346 rCoordInputData = .FALSE.
0347 coordsSet = 0
0348
8fc5af62b8 Jean*0349
ede2cfe569 Chri*0350 nRetired = 0
c07cd3bfa8 Jean*0351 useConstantF = .FALSE.
0352 useBetaPlaneF = .FALSE.
0353 useSphereF = .TRUE.
7a4ec66481 Jean*0354 tempDiffusion = .TRUE.
0355 saltDiffusion = .TRUE.
0356 tracerAdvScheme = UNSET_I
0357 trac_EvPrRn = UNSET_RL
ede2cfe569 Chri*0358 zonal_filt_lat = UNSET_RL
36b12bb7ff Jean*0359 gravitySign = UNSET_RL
8f934af0f4 Jean*0360 viscAstrain = UNSET_RL
0361 viscAtension = UNSET_RL
742cf4499c Jean*0362 useAnisotropicViscAgridMax=.TRUE.
6d61b52b09 Jean*0363 usePickupBeforeC35 = .FALSE.
83513ae28b Jean*0364 saveDebugMode = debugMode
e0b3e1bdd8 Dimi*0365 allowInteriorFreezing = .FALSE.
e793395f17 Jean*0366 useOldFreezing = .FALSE.
427e24e121 Jean*0367 writeStatePrec = UNSET_I
7e00d7e8f9 Jean*0368 cg2dChkResFreq = UNSET_I
aecc8b0f47 Mart*0369 cg3dChkResFreq = UNSET_I
46979d29da Jean*0370 tauThetaClimRelax3Dim = UNSET_RL
0371 tauSaltClimRelax3Dim = UNSET_RL
742cf4499c Jean*0372 calendarDumps = .FALSE.
2ad79bdf32 Jean*0373 Ro_SeaLevel = UNSET_RL
742cf4499c Jean*0374 rkFac = UNSET_RL
0375 groundAtK1 = .FALSE.
9780090eaa Jean*0376 thetaMin = UNSET_RL
0377 phiMin = UNSET_RL
8194bc4a99 Mart*0378 shelfIceFile = ' '
275ad78b78 Jean*0379 dQdTFile = ' '
ede2cfe569 Chri*0380
924557e60a Chri*0381
b9306711eb Jean*0382 WRITE(msgBuf,'(A)')
aa8f58c33d Jean*0383 & ' INI_PARMS: opening model parameter file "data"'
924557e60a Chri*0384 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0385 & SQUEEZE_RIGHT, myThid )
924557e60a Chri*0386
aa8f58c33d Jean*0387 CALL OPEN_COPY_DATA_FILE( 'data', 'INI_PARMS',
0388 O iUnit, myThid )
924557e60a Chri*0389
aa8f58c33d Jean*0390
0391 errIO = 0
90aed06079 Jean*0392 errCount = 0
924557e60a Chri*0393
0394
6deb815953 Mart*0395 viscAhW = UNSET_RL
0b6cbae535 Mart*0396 viscA4W = UNSET_RL
915a82d1f5 Jean*0397 viscAhD = UNSET_RL
0398 viscAhZ = UNSET_RL
0399 viscA4D = UNSET_RL
0400 viscA4Z = UNSET_RL
809bdccbfc Jean*0401 viscAz = UNSET_RL
910f05e765 Chri*0402 viscAp = UNSET_RL
7cd7754949 Jean*0403 viscAr = UNSET_RL
910f05e765 Chri*0404 diffKzT = UNSET_RL
0405 diffKpT = UNSET_RL
0406 diffKrT = UNSET_RL
0407 diffKzS = UNSET_RL
0408 diffKpS = UNSET_RL
0409 diffKrS = UNSET_RL
aad87ebb4c Jean*0410 hFacMinDr = UNSET_RL
0411 hFacMinDz = UNSET_RL
0412 hFacMinDp = UNSET_RL
0413 tAlpha = UNSET_RL
0414 sBeta = UNSET_RL
0415 implicitNHPress = UNSET_RL
49efc6c1e0 Jean*0416 tempVertAdvScheme = 0
0417 saltVertAdvScheme = 0
337b46d524 Jean*0418 plotLevel = UNSET_I
49efc6c1e0 Jean*0419
aa8f58c33d Jean*0420 WRITE(msgBuf,'(A)') ' INI_PARMS ; starts to read PARM01'
6b6f69d6b0 Jean*0421 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0422 & SQUEEZE_RIGHT, myThid )
88830be691 Alis*0423 READ(UNIT=iUnit,NML=PARM01)
c206b07c1d Chri*0424 IF ( errIO .LT. 0 ) THEN
924557e60a Chri*0425 WRITE(msgBuf,'(A)')
aa8f58c33d Jean*0426 & 'S/R INI_PARMS: Error reading model parameter file "data"'
b9306711eb Jean*0427 CALL PRINT_ERROR( msgBuf, myThid )
aa8f58c33d Jean*0428 WRITE(msgBuf,'(A)') 'S/R INI_PARMS: Problem in namelist PARM01'
b9306711eb Jean*0429 CALL PRINT_ERROR( msgBuf, myThid )
924557e60a Chri*0430 STOP 'ABNORMAL END: S/R INI_PARMS'
41f1ac6de1 Jean*0431 ELSE
aa8f58c33d Jean*0432 WRITE(msgBuf,'(A)') ' INI_PARMS ; read PARM01 : OK'
41f1ac6de1 Jean*0433 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0434 & SQUEEZE_RIGHT, myThid )
c206b07c1d Chri*0435 ENDIF
463053c692 Jean*0436
599be48d9a Jean*0437
e99e7fe4e6 Jean*0438
0439 usingPCoords = .FALSE.
0440 usingZCoords = .FALSE.
0441 fluidIsAir = .FALSE.
0442 fluidIsWater = .FALSE.
0443 IF ( buoyancyRelation.EQ.'ATMOSPHERIC' ) THEN
0444 usingPCoords = .TRUE.
0445 fluidIsAir = .TRUE.
0446 ELSEIF ( buoyancyRelation.EQ.'OCEANICP') THEN
0447 usingPCoords = .TRUE.
0448 fluidIsWater = .TRUE.
0449 ELSEIF ( buoyancyRelation.EQ.'OCEANIC' ) THEN
0450 usingZCoords = .TRUE.
0451 fluidIsWater = .TRUE.
0452 ELSE
0453 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS:',
0454 & ' Bad value of buoyancyRelation '
b9306711eb Jean*0455 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0456 errCount = errCount + 1
e99e7fe4e6 Jean*0457 ENDIF
0458
f615684c22 Jean*0459 IF ( .NOT.rigidLid .AND.
0460 & .NOT.implicitFreeSurface ) THEN
0461
0462 WRITE(msgBuf,'(A)')
f76021331f Jean*0463 & 'S/R INI_PARMS: No request for barotropic solver'
f615684c22 Jean*0464 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0465 & SQUEEZE_RIGHT, myThid )
f615684c22 Jean*0466 WRITE(msgBuf,'(A)')
0467 & 'S/R INI_PARMS: => Use implicitFreeSurface as default'
0468 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0469 & SQUEEZE_RIGHT, myThid )
f615684c22 Jean*0470 implicitFreeSurface = .TRUE.
0471 ENDIF
8039e9b985 Jean*0472 IF ( implicitFreeSurface ) freeSurfFac = 1. _d 0
0473 IF ( rigidLid ) freeSurfFac = 0. _d 0
0127add478 Alis*0474 IF ( gBaro .EQ. UNSET_RL ) gBaro=gravity
0475 IF ( rhoConst .EQ. UNSET_RL ) rhoConst=rhoNil
01d019192b Jean*0476 IF ( rhoConstFresh .EQ. UNSET_RL ) rhoConstFresh=rhoConst
16f5093311 Jean*0477 IF ( implicitNHPress.EQ.UNSET_RL )
0478 & implicitNHPress = implicSurfPress
81a2602dac Jean*0479 IF ( omega .EQ. UNSET_RL ) THEN
f0bc3639a4 Jean*0480 omega = 0. _d 0
b9306711eb Jean*0481 IF ( rotationPeriod .NE. 0. _d 0 )
8039e9b985 Jean*0482 & omega = 2. _d 0 * PI / rotationPeriod
f0bc3639a4 Jean*0483 ELSEIF ( omega .EQ. 0. _d 0 ) THEN
0484 rotationPeriod = 0. _d 0
81a2602dac Jean*0485 ELSE
8039e9b985 Jean*0486 rotationPeriod = 2. _d 0 * PI / omega
81a2602dac Jean*0487 ENDIF
aad87ebb4c Jean*0488 IF ( atm_Rd .EQ. UNSET_RL ) THEN
463053c692 Jean*0489 atm_Rd = atm_Cp * atm_kappa
0490 ELSE
0491 atm_kappa = atm_Rd / atm_Cp
0492 ENDIF
775ad00e06 Alis*0493
0494 IF (nonHydrostatic.AND.quasiHydrostatic) THEN
0495 WRITE(msgBuf,'(A)')
0496 & 'Illegal: both nonHydrostatic = quasiHydrostatic = TRUE'
b9306711eb Jean*0497 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0498 errCount = errCount + 1
775ad00e06 Alis*0499 ENDIF
809bdccbfc Jean*0500
49efc6c1e0 Jean*0501 IF (tempVertAdvScheme.EQ.0) tempVertAdvScheme = tempAdvScheme
0502 IF (saltVertAdvScheme.EQ.0) saltVertAdvScheme = saltAdvScheme
915a82d1f5 Jean*0503
0504 IF ( viscAhD .EQ. UNSET_RL ) viscAhD = viscAh
0505 IF ( viscAhZ .EQ. UNSET_RL ) viscAhZ = viscAh
0506 IF ( viscA4D .EQ. UNSET_RL ) viscA4D = viscA4
0507 IF ( viscA4Z .EQ. UNSET_RL ) viscA4Z = viscA4
2d5bb917cc Jean*0508
92da82d27f Jean*0509 IF ( viscAhW .EQ. UNSET_RL ) viscAhW = viscAhD
0510 IF ( viscA4W .EQ. UNSET_RL ) viscA4W = viscA4D
910f05e765 Chri*0511
0512 IF ( viscAz .NE. UNSET_RL ) zCoordInputData = .TRUE.
0513 IF ( viscAp .NE. UNSET_RL ) pCoordInputData = .TRUE.
0514 IF ( viscAr .NE. UNSET_RL ) rCoordInputData = .TRUE.
0515 IF ( viscAr .EQ. UNSET_RL ) viscAr = viscAz
0516 IF ( viscAr .EQ. UNSET_RL ) viscAr = viscAp
7cd7754949 Jean*0517 vertSetCount = 0
0518 DO k=1,Nr
13b4d13c98 Jean*0519 IF ( viscArNr(k).NE.UNSET_RL ) vertSetCount = vertSetCount + 1
7cd7754949 Jean*0520 ENDDO
13b4d13c98 Jean*0521 IF ( vertSetCount.GT.0 .AND. vertSetCount.LT.Nr ) THEN
0522 WRITE(msgBuf,'(A,2(I5,A))') 'S/R INI_PARMS: Partial setting (',
0523 & vertSetCount, ' /', Nr, ') of viscArNr is not allowed'
0524 CALL PRINT_ERROR( msgBuf, myThid )
0525 errCount = errCount + 1
0526 ENDIF
0527 IF ( viscAr .EQ. UNSET_RL ) THEN
0528 viscAr = viscArDefault
0529 ELSEIF ( vertSetCount.GT.0 ) THEN
0530 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS: Cannot set both ',
0531 & 'viscArNr and viscAr (or Ap,Az) in param file data'
0532 CALL PRINT_ERROR( msgBuf, myThid )
0533 errCount = errCount + 1
0534 ENDIF
0535 IF ( vertSetCount.EQ.0 ) THEN
0536 DO k=1,Nr
0537 viscArNr(k) = viscAr
0538 ENDDO
7cd7754949 Jean*0539 ENDIF
f9ef55fa0c Jean*0540 #ifdef ALLOW_MOM_COMMON
0541
0542 IF ( selectBotDragQuadr.EQ.-1 .AND. bottomDragQuadratic.NE.0. )
0543 & selectBotDragQuadr = 0
ab47de63dc Mart*0544 IF ( selectBotDragQuadr.EQ.-1 .AND. zRoughBot.NE.0. )
0545 & selectBotDragQuadr = 0
f9ef55fa0c Jean*0546 #endif /* ALLOW_MOM_COMMON */
910f05e765 Chri*0547
2d5bb917cc Jean*0548 IF ( smag3D_diffCoeff.NE.zeroRL .AND. .NOT.useSmag3D ) THEN
0549 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: ',
0550 & 'will not use "smag3D_diffCoeff" without useSmag3D'
0551 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
0552 & SQUEEZE_RIGHT, myThid )
0553 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: ',
0554 & '==> reset "smag3D_diffCoeff" to zero'
0555 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
0556 & SQUEEZE_RIGHT, myThid )
0557 smag3D_diffCoeff = zeroRL
0558 ENDIF
910f05e765 Chri*0559 IF ( diffKzT .NE. UNSET_RL ) zCoordInputData = .TRUE.
0560 IF ( diffKpT .NE. UNSET_RL ) pCoordInputData = .TRUE.
0561 IF ( diffKrT .NE. UNSET_RL ) rCoordInputData = .TRUE.
0562 IF ( diffKrT .EQ. UNSET_RL ) diffKrT = diffKzT
0563 IF ( diffKrT .EQ. UNSET_RL ) diffKrT = diffKpT
7cd7754949 Jean*0564 vertSetCount = 0
8fc5af62b8 Jean*0565 DO k=1,Nr
13b4d13c98 Jean*0566 IF ( diffKrNrT(k).NE.UNSET_RL ) vertSetCount = vertSetCount + 1
8fc5af62b8 Jean*0567 ENDDO
13b4d13c98 Jean*0568 IF ( vertSetCount.GT.0 .AND. vertSetCount.LT.Nr ) THEN
0569 WRITE(msgBuf,'(A,2(I5,A))') 'S/R INI_PARMS: Partial setting (',
0570 & vertSetCount, ' /', Nr, ') of diffKrNrT is not allowed'
0571 CALL PRINT_ERROR( msgBuf, myThid )
0572 errCount = errCount + 1
0573 ENDIF
0574 specifiedDiffKrT = vertSetCount.EQ.Nr
0575 IF ( diffKrT .EQ. UNSET_RL ) THEN
0576 diffKrT = diffKrTDefault
0577 ELSEIF ( vertSetCount.GT.0 ) THEN
0578 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS: Cannot set both ',
0579 & 'diffKrNrT and diffKrT (or Kp,Kz) in param file data'
0580 CALL PRINT_ERROR( msgBuf, myThid )
0581 errCount = errCount + 1
0582 ELSE
0583 specifiedDiffKrT = .TRUE.
0584 ENDIF
0585 IF ( vertSetCount.EQ.0 ) THEN
0586 DO k=1,Nr
0587 diffKrNrT(k) = diffKrT
0588 ENDDO
8fc5af62b8 Jean*0589 ENDIF
910f05e765 Chri*0590
0591 IF ( diffKzS .NE. UNSET_RL ) zCoordInputData = .TRUE.
0592 IF ( diffKpS .NE. UNSET_RL ) pCoordInputData = .TRUE.
0593 IF ( diffKrS .NE. UNSET_RL ) rCoordInputData = .TRUE.
0594 IF ( diffKrS .EQ. UNSET_RL ) diffKrS = diffKzS
0595 IF ( diffKrS .EQ. UNSET_RL ) diffKrS = diffKpS
7cd7754949 Jean*0596 vertSetCount = 0
8fc5af62b8 Jean*0597 DO k=1,Nr
13b4d13c98 Jean*0598 IF ( diffKrNrS(k).NE.UNSET_RL ) vertSetCount = vertSetCount + 1
8fc5af62b8 Jean*0599 ENDDO
13b4d13c98 Jean*0600 IF ( vertSetCount.GT.0 .AND. vertSetCount.LT.Nr ) THEN
0601 WRITE(msgBuf,'(A,2(I5,A))') 'S/R INI_PARMS: Partial setting (',
0602 & vertSetCount, ' /', Nr, ') of diffKrNrS is not allowed'
0603 CALL PRINT_ERROR( msgBuf, myThid )
0604 errCount = errCount + 1
0605 ENDIF
a631d4754d Jean*0606 IF ( vertSetCount.EQ.Nr ) THEN
0607 specifiedDiffKrS = .TRUE.
0608 IF ( diffKrS.NE.UNSET_RL ) THEN
0609 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS: Cannot set both ',
13b4d13c98 Jean*0610 & 'diffKrNrS and diffKrS (or Kp,Kz) in param file data'
a631d4754d Jean*0611 CALL PRINT_ERROR( msgBuf, myThid )
0612 errCount = errCount + 1
0613 ENDIF
0614 ELSEIF ( diffKrS.NE.UNSET_RL ) THEN
13b4d13c98 Jean*0615 specifiedDiffKrS = .TRUE.
0616 DO k=1,Nr
0617 diffKrNrS(k) = diffKrS
0618 ENDDO
a631d4754d Jean*0619 ELSE
0620 specifiedDiffKrS = .FALSE.
0621 diffKrS = diffKrSDefault
0622
0623 DO k=1,Nr
0624 diffKrNrS(k) = diffKrNrT(k)
0625 ENDDO
8fc5af62b8 Jean*0626 ENDIF
910f05e765 Chri*0627
e40c34e398 Dimi*0628 IF (diffKrBLEQsurf .EQ. UNSET_RL) diffKrBLEQsurf = diffKrBL79surf
0629 IF (diffKrBLEQdeep .EQ. UNSET_RL) diffKrBLEQdeep = diffKrBL79deep
0630 IF (diffKrBLEQscl .EQ. UNSET_RL) diffKrBLEQscl = diffKrBL79scl
0631 IF (diffKrBLEQHo .EQ. UNSET_RL) diffKrBLEQHo = diffKrBL79Ho
0632
910f05e765 Chri*0633 IF ( hFacMinDz .NE. UNSET_RL ) zCoordInputData = .TRUE.
0634 IF ( hFacMinDp .NE. UNSET_RL ) pCoordInputData = .TRUE.
0635 IF ( hFacMinDr .NE. UNSET_RL ) rCoordInputData = .TRUE.
fe0a6d6615 Alis*0636 IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr = hFacMinDz
0637 IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr = hFacMinDp
910f05e765 Chri*0638 IF ( hFacMinDr .EQ. UNSET_RL ) hFacMinDr = hFacMinDrDefault
c0a4efc370 Chri*0639
745a0098ab Jean*0640 IF (convertFW2Salt.EQ.UNSET_RL) THEN
0641 convertFW2Salt = 35.
0642 IF (useRealFreshWaterFlux) convertFW2Salt=-1
d18df35fee Jean*0643 IF ( selectAddFluid.GE.1 ) convertFW2Salt=-1
745a0098ab Jean*0644 ENDIF
0645
7843dde2de jm-c 0646
0647 IF ( selectCoriScheme.EQ.UNSET_I ) THEN
0648 selectCoriScheme = 0
0649 IF ( useJamartWetPoints ) selectCoriScheme = 1
0650 IF ( useEnergyConservingCoriolis .AND.
0651 & .NOT.vectorInvariantMomentum )
0652 & selectCoriScheme = selectCoriScheme + 2
0653 ENDIF
0654 IF ( vectorInvariantMomentum ) THEN
0655 IF ( useJamartWetPoints .AND. selectCoriScheme.NE.1 ) THEN
0656 WRITE(msgBuf,'(A,I5,A)')
0657 & 'S/R INI_PARMS: selectCoriScheme=', selectCoriScheme,
0658 & ' conflicts with "useJamartWetPoints"'
0659 CALL PRINT_ERROR( msgBuf, myThid )
0660 errCount = errCount + 1
0661 ENDIF
0662 ELSE
0663 IF ( useEnergyConservingCoriolis
0664 & .AND. selectCoriScheme.LT.2 ) THEN
0665 WRITE(msgBuf,'(A,I5,A)')
0666 & 'S/R INI_PARMS: selectCoriScheme=', selectCoriScheme,
0667 & ' conflicts with "useEnergyConservingCoriolis"'
0668 CALL PRINT_ERROR( msgBuf, myThid )
0669 errCount = errCount + 1
0670 ENDIF
0671 IF ( useJamartWetPoints .AND. selectCoriScheme.NE.1
0672 & .AND. selectCoriScheme.NE.3 ) THEN
0673 WRITE(msgBuf,'(A,I5,A)')
0674 & 'S/R INI_PARMS: selectCoriScheme=', selectCoriScheme,
0675 & ' conflicts with "useJamartWetPoints"'
0676 CALL PRINT_ERROR( msgBuf, myThid )
0677 errCount = errCount + 1
0678 ENDIF
0679 ENDIF
e35a804907 Jean*0680 IF ( SadournyCoriolis ) THEN
0681
0682 IF ( selectVortScheme.EQ.UNSET_I ) selectVortScheme = 2
0683 IF ( selectVortScheme.NE.2 ) THEN
7fe6343684 Jean*0684 WRITE(msgBuf,'(A,I5,A)')
0685 & 'S/R INI_PARMS: selectVortScheme=', selectVortScheme,
0686 & ' conflicts with "SadournyCoriolis"'
e35a804907 Jean*0687 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0688 errCount = errCount + 1
e35a804907 Jean*0689 ENDIF
0690 ENDIF
427e24e121 Jean*0691 IF ( select3dCoriScheme.EQ.UNSET_I ) THEN
0692
0693 select3dCoriScheme = 0
0694 IF ( use3dCoriolis ) select3dCoriScheme = 1
0695 ELSEIF ( select3dCoriScheme.NE.0 .AND. .NOT.use3dCoriolis ) THEN
0696 WRITE(msgBuf,'(A,I5,A)')
0697 & 'S/R INI_PARMS: select3dCoriScheme=', select3dCoriScheme,
0698 & ' conflicts with "use3dCoriolis=F"'
0699 CALL PRINT_ERROR( msgBuf, myThid )
0700 errCount = errCount + 1
0701 ENDIF
0702 IF ( selectMetricTerms.EQ.UNSET_I ) THEN
0703
0704 selectMetricTerms = 0
0705 IF ( metricTerms ) selectMetricTerms = 1
0706 ELSEIF ( selectMetricTerms.NE.0 .AND. .NOT.metricTerms ) THEN
0707 WRITE(msgBuf,'(A,I5,A)')
0708 & 'S/R INI_PARMS: selectMetricTerms=', selectMetricTerms,
0709 & ' conflicts with "metricTerms=F"'
0710 CALL PRINT_ERROR( msgBuf, myThid )
0711 errCount = errCount + 1
0712 ENDIF
7e00d7e8f9 Jean*0713 IF ( selectBalanceEmPmR.EQ.UNSET_I ) THEN
0714 selectBalanceEmPmR = 0
0715 IF ( balanceEmPmR ) selectBalanceEmPmR = 1
0716 ELSEIF ( selectBalanceEmPmR.NE.1 .AND. balanceEmPmR ) THEN
0717 WRITE(msgBuf,'(A,I5,A)')
0718 & 'S/R INI_PARMS: selectBalanceEmPmR=', selectBalanceEmPmR,
0719 & ' conflicts with "balanceEmPmR"'
0720 CALL PRINT_ERROR( msgBuf, myThid )
0721 errCount = errCount + 1
0722 ENDIF
e35a804907 Jean*0723
aad87ebb4c Jean*0724 IF ( ivdc_kappa.NE.zeroRL .AND. .NOT.implicitDiffusion ) THEN
e1c6dcc4ea Jean*0725 WRITE(msgBuf,'(A,A)')
d4701cb6da Alis*0726 & 'S/R INI_PARMS: To use ivdc_kappa you must enable implicit',
0727 & ' vertical diffusion.'
b9306711eb Jean*0728 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0729 errCount = errCount + 1
d4701cb6da Alis*0730 ENDIF
0731
910f05e765 Chri*0732 coordsSet = 0
0733 IF ( zCoordInputData ) coordsSet = coordsSet + 1
0734 IF ( pCoordInputData ) coordsSet = coordsSet + 1
0735 IF ( rCoordInputData ) coordsSet = coordsSet + 1
0736 IF ( coordsSet .GT. 1 ) THEN
0737 WRITE(msgBuf,'(A)')
0738 & 'S/R INI_PARMS: Cannot mix z, p and r in the input data.'
b9306711eb Jean*0739 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0740 errCount = errCount + 1
c0a4efc370 Chri*0741 ENDIF
910f05e765 Chri*0742 IF ( rhoConst .LE. 0. ) THEN
0743 WRITE(msgBuf,'(A)')
0744 & 'S/R INI_PARMS: rhoConst must be greater than 0.'
b9306711eb Jean*0745 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0746 errCount = errCount + 1
0747 recip_rhoConst = 1. _d 0
910f05e765 Chri*0748 ELSE
8039e9b985 Jean*0749 recip_rhoConst = 1. _d 0 / rhoConst
910f05e765 Chri*0750 ENDIF
53a1bebe0a Jean*0751 IF ( eosType.EQ.'LINEAR' .AND. rhoNil.LE.0. ) THEN
9bef66796c Alis*0752 WRITE(msgBuf,'(A)')
0753 & 'S/R INI_PARMS: rhoNil must be greater than 0.'
b9306711eb Jean*0754 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0755 errCount = errCount + 1
9bef66796c Alis*0756 ENDIF
0127add478 Alis*0757 IF ( HeatCapacity_Cp .LE. 0. ) THEN
0758 WRITE(msgBuf,'(A)')
0759 & 'S/R INI_PARMS: HeatCapacity_Cp must be greater than 0.'
b9306711eb Jean*0760 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0761 errCount = errCount + 1
0127add478 Alis*0762 ENDIF
b0f9f620c0 Chri*0763 IF ( gravity .LE. 0. ) THEN
0764 WRITE(msgBuf,'(A)')
0765 & 'S/R INI_PARMS: gravity must be greater than 0.'
b9306711eb Jean*0766 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*0767 errCount = errCount + 1
0768 recip_gravity = 1. _d 0
b0f9f620c0 Chri*0769 ELSE
8039e9b985 Jean*0770 recip_gravity = 1. _d 0 / gravity
b0f9f620c0 Chri*0771 ENDIF
631fe75038 Dimi*0772
83513ae28b Jean*0773
0774 printResidualFreq = 0
0775 IF ( debugLevel.GE.debLevE ) printResidualFreq = 1
337b46d524 Jean*0776 IF ( plotLevel.EQ.UNSET_I ) plotLevel = debugLevel
924557e60a Chri*0777
631fe75038 Dimi*0778
0779 IF ( useSingleCpuIO ) useSingleCpuInput=.TRUE.
0780
ede2cfe569 Chri*0781
0782 nRetired = 0
c07cd3bfa8 Jean*0783 IF ( useConstantF ) THEN
0784 nRetired = nRetired+1
0785 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "useConstantF" ',
0786 & 'is no longer allowed in file "data"'
0787 CALL PRINT_ERROR( msgBuf, myThid )
0788 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: set "selectCoriMap"',
0789 & ' [0,1,2,3] to impose a setting over grid default'
0790 CALL PRINT_ERROR( msgBuf, myThid )
0791 ENDIF
0792 IF ( useBetaPlaneF ) THEN
0793 nRetired = nRetired+1
0794 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "useBetaPlaneF" ',
0795 & 'is no longer allowed in file "data"'
0796 CALL PRINT_ERROR( msgBuf, myThid )
0797 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: set "selectCoriMap"',
0798 & ' [0,1,2,3] to impose a setting over grid default'
0799 CALL PRINT_ERROR( msgBuf, myThid )
0800 ENDIF
0801 IF ( .NOT. useSphereF ) THEN
0802 nRetired = nRetired+1
0803 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "useSphereF" ',
0804 & 'is no longer allowed in file "data"'
0805 CALL PRINT_ERROR( msgBuf, myThid )
0806 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: set "selectCoriMap"',
0807 & ' [0,1,2,3] to impose a setting over grid default'
0808 CALL PRINT_ERROR( msgBuf, myThid )
0809 ENDIF
ede2cfe569 Chri*0810 IF ( zonal_filt_lat .NE. UNSET_RL ) THEN
0811 nRetired = nRetired+1
0812 WRITE(msgBuf,'(A,A)')
0813 & 'S/R INI_PARMS: Paramater "zonal_filt_lat" is',
0814 & ' no longer allowed in file "data".'
b9306711eb Jean*0815 CALL PRINT_ERROR( msgBuf, myThid )
ede2cfe569 Chri*0816 WRITE(msgBuf,'(A,A)')
0817 & 'S/R INI_PARMS: Paramater "zonal_filt_lat" is',
0818 & ' now read from file "data.zonfilt".'
b9306711eb Jean*0819 CALL PRINT_ERROR( msgBuf, myThid )
ede2cfe569 Chri*0820 ENDIF
36b12bb7ff Jean*0821 IF ( gravitySign .NE. UNSET_RL ) THEN
0822 nRetired = nRetired+1
0823 WRITE(msgBuf,'(A,A)')
0824 & 'S/R INI_PARMS: "gravitySign" is set according to vertical ',
0825 & ' coordinate and is no longer allowed in file "data".'
b9306711eb Jean*0826 CALL PRINT_ERROR( msgBuf, myThid )
36b12bb7ff Jean*0827 ENDIF
7a4ec66481 Jean*0828 IF ( tracerAdvScheme .NE. UNSET_I ) THEN
0829 nRetired = nRetired+1
0830 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "tracerAdvScheme" ',
0831 & '(old passive tracer code) is no longer allowed in file "data"'
b9306711eb Jean*0832 CALL PRINT_ERROR( msgBuf, myThid )
7a4ec66481 Jean*0833 ENDIF
0834 IF ( trac_EvPrRn .NE. UNSET_RL ) THEN
0835 nRetired = nRetired+1
0836 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "trac_EvPrRn" ',
0837 & '(old passive tracer code) is no longer allowed in file "data"'
b9306711eb Jean*0838 CALL PRINT_ERROR( msgBuf, myThid )
7a4ec66481 Jean*0839 ENDIF
0840 IF ( .NOT. tempDiffusion ) THEN
0841 nRetired = nRetired+1
0842 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "tempDiffusion" ',
0843 & 'is no longer allowed in file "data"'
b9306711eb Jean*0844 CALL PRINT_ERROR( msgBuf, myThid )
7a4ec66481 Jean*0845 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: to turn off diffusion',
0846 & ' => set diffusivity to zero'
b9306711eb Jean*0847 CALL PRINT_ERROR( msgBuf, myThid )
7a4ec66481 Jean*0848 ENDIF
0849 IF ( .NOT. saltDiffusion ) THEN
0850 nRetired = nRetired+1
0851 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "saltDiffusion" ',
0852 & 'is no longer allowed in file "data"'
b9306711eb Jean*0853 CALL PRINT_ERROR( msgBuf, myThid )
7a4ec66481 Jean*0854 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: to turn off diffusion',
0855 & ' => set diffusivity to zero'
b9306711eb Jean*0856 CALL PRINT_ERROR( msgBuf, myThid )
7a4ec66481 Jean*0857 ENDIF
8f934af0f4 Jean*0858 IF ( viscAstrain .NE. UNSET_RL ) THEN
0859 nRetired = nRetired+1
0860 WRITE(msgBuf,'(A,A)')
0861 & 'S/R INI_PARMS: "viscAstrain" ',
0862 & 'is no longer allowed in file "data"'
b9306711eb Jean*0863 CALL PRINT_ERROR( msgBuf, myThid )
8f934af0f4 Jean*0864 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: to use Strain & Tension',
0865 & ' formulation => set useStrainTensionVisc to TRUE'
b9306711eb Jean*0866 CALL PRINT_ERROR( msgBuf, myThid )
8f934af0f4 Jean*0867 ENDIF
0868 IF ( viscAtension .NE. UNSET_RL ) THEN
0869 nRetired = nRetired+1
0870 WRITE(msgBuf,'(A,A)')
0871 & 'S/R INI_PARMS: "viscAtension" ',
0872 & 'is no longer allowed in file "data"'
b9306711eb Jean*0873 CALL PRINT_ERROR( msgBuf, myThid )
8f934af0f4 Jean*0874 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: to use Strain & Tension',
0875 & ' formulation => set useStrainTensionVisc to TRUE'
b9306711eb Jean*0876 CALL PRINT_ERROR( msgBuf, myThid )
8f934af0f4 Jean*0877 ENDIF
6d61b52b09 Jean*0878 IF ( .NOT.useAnisotropicViscAgridMax ) THEN
cef83a8f90 Bayl*0879 nRetired = nRetired+1
0880 WRITE(msgBuf,'(A,A)')
0881 & 'S/R INI_PARMS: "useAnisotropicViscAgridMax" ',
0882 & 'is not allowed in "data" substitute useAreaViscLength=true'
b9306711eb Jean*0883 CALL PRINT_ERROR( msgBuf, myThid )
cef83a8f90 Bayl*0884 ENDIF
6d61b52b09 Jean*0885 IF ( usePickupBeforeC35 ) THEN
0886 nRetired = nRetired+1
0887 WRITE(msgBuf,'(A,A)')
0888 & 'S/R INI_PARMS: "usePickupBeforeC35" ',
0889 & 'is no longer supported & not longer allowed in file "data"'
0890 CALL PRINT_ERROR( msgBuf, myThid )
0891 ENDIF
83513ae28b Jean*0892 IF ( debugMode.NEQV.saveDebugMode ) THEN
0893 nRetired = nRetired+1
0894 WRITE(msgBuf,'(A,A)')
0895 & 'S/R INI_PARMS: "debugMode" has been moved to "eedata"',
0896 & ' and is no longer allowed in file "data"'
0897 CALL PRINT_ERROR( msgBuf, myThid )
0898 ENDIF
e0b3e1bdd8 Dimi*0899 IF ( allowInteriorFreezing ) THEN
0900 nRetired = nRetired+1
0901 WRITE(msgBuf,'(A,A)')
0902 & 'S/R INI_PARMS: "allowInteriorFreezing" has been replaced',
0903 & ' by pkg/frazil and is no longer allowed in file "data"'
0904 CALL PRINT_ERROR( msgBuf, myThid )
0905 ENDIF
e793395f17 Jean*0906 IF ( useOldFreezing ) THEN
0907 nRetired = nRetired+1
0908 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "useOldFreezing" ',
0909 & 'is no longer supported & not longer allowed in file "data"'
0910 CALL PRINT_ERROR( msgBuf, myThid )
0911 ENDIF
427e24e121 Jean*0912 IF ( writeStatePrec .NE. UNSET_I ) THEN
0913 nRetired = nRetired+1
0914 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "writeStatePrec" ',
0915 & '(un-used) is no longer allowed in file "data"'
0916 CALL PRINT_ERROR( msgBuf, myThid )
0917 ENDIF
ede2cfe569 Chri*0918
924557e60a Chri*0919
aa8f58c33d Jean*0920 WRITE(msgBuf,'(A)') ' INI_PARMS ; starts to read PARM02'
6b6f69d6b0 Jean*0921 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0922 & SQUEEZE_RIGHT, myThid )
88830be691 Alis*0923 READ(UNIT=iUnit,NML=PARM02)
c206b07c1d Chri*0924 IF ( errIO .LT. 0 ) THEN
924557e60a Chri*0925 WRITE(msgBuf,'(A)')
aa8f58c33d Jean*0926 & 'S/R INI_PARMS: Error reading model parameter file "data"'
b9306711eb Jean*0927 CALL PRINT_ERROR( msgBuf, myThid )
aa8f58c33d Jean*0928 WRITE(msgBuf,'(A)') 'S/R INI_PARMS: Problem in namelist PARM02'
b9306711eb Jean*0929 CALL PRINT_ERROR( msgBuf, myThid )
924557e60a Chri*0930 STOP 'ABNORMAL END: S/R INI_PARMS'
41f1ac6de1 Jean*0931 ELSE
aa8f58c33d Jean*0932 WRITE(msgBuf,'(A)') ' INI_PARMS ; read PARM02 : OK'
41f1ac6de1 Jean*0933 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0934 & SQUEEZE_RIGHT, myThid )
809bdccbfc Jean*0935 ENDIF
aecc8b0f47 Mart*0936
0937 IF ( cg2dChkResFreq .NE. UNSET_I ) THEN
0938 nRetired = nRetired+1
0939 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: unused "cg2dChkResFreq"',
0940 & ' is no longer allowed in file "data"'
0941 CALL PRINT_ERROR( msgBuf, myThid )
0942 ENDIF
0943 IF ( cg3dChkResFreq .NE. UNSET_I ) THEN
0944 nRetired = nRetired+1
0945 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: unused "cg3dChkResFreq"',
0946 & ' is no longer allowed in file "data"'
0947 CALL PRINT_ERROR( msgBuf, myThid )
0948 ENDIF
924557e60a Chri*0949
0950
8039e9b985 Jean*0951 rCD = -1. _d 0
0952 epsAB_CD = UNSET_RL
36b12bb7ff Jean*0953 latBandClimRelax = UNSET_RL
062a876ce5 Jean*0954 deltaTtracer = 0. _d 0
c2b6ed6bfd Jean*0955 forcing_In_AB = .TRUE.
aa8f58c33d Jean*0956 WRITE(msgBuf,'(A)') ' INI_PARMS ; starts to read PARM03'
6b6f69d6b0 Jean*0957 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0958 & SQUEEZE_RIGHT, myThid )
88830be691 Alis*0959 READ(UNIT=iUnit,NML=PARM03)
c206b07c1d Chri*0960 IF ( errIO .LT. 0 ) THEN
924557e60a Chri*0961 WRITE(msgBuf,'(A)')
aa8f58c33d Jean*0962 & 'S/R INI_PARMS: Error reading model parameter file "data"'
b9306711eb Jean*0963 CALL PRINT_ERROR( msgBuf, myThid )
aa8f58c33d Jean*0964 WRITE(msgBuf,'(A)') 'S/R INI_PARMS: Problem in namelist PARM03'
b9306711eb Jean*0965 CALL PRINT_ERROR( msgBuf, myThid )
924557e60a Chri*0966 STOP 'ABNORMAL END: S/R INI_PARMS'
41f1ac6de1 Jean*0967 ELSE
aa8f58c33d Jean*0968 WRITE(msgBuf,'(A)') ' INI_PARMS ; read PARM03 : OK'
41f1ac6de1 Jean*0969 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*0970 & SQUEEZE_RIGHT, myThid )
809bdccbfc Jean*0971 ENDIF
46979d29da Jean*0972
0973 IF ( tauThetaClimRelax3Dim .NE. UNSET_RL ) THEN
0974 nRetired = nRetired+1
0975 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "tauThetaClimRelax3Dim" ',
0976 & 'is no longer allowed in file "data"'
b9306711eb Jean*0977 CALL PRINT_ERROR( msgBuf, myThid )
46979d29da Jean*0978 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: 3-dim. relaxation code',
0979 & ' has moved to separate pkg/rbcs.'
b9306711eb Jean*0980 CALL PRINT_ERROR( msgBuf, myThid )
46979d29da Jean*0981 ENDIF
0982 IF ( tauSaltClimRelax3Dim .NE. UNSET_RL ) THEN
0983 nRetired = nRetired+1
0984 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "tauSaltClimRelax3Dim" ',
0985 & 'is no longer allowed in file "data"'
b9306711eb Jean*0986 CALL PRINT_ERROR( msgBuf, myThid )
46979d29da Jean*0987 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: 3-dim. relaxation code',
0988 & ' has moved to separate pkg/rbcs.'
b9306711eb Jean*0989 CALL PRINT_ERROR( msgBuf, myThid )
46979d29da Jean*0990 ENDIF
742cf4499c Jean*0991 IF ( calendarDumps ) THEN
0992 nRetired = nRetired+1
0993 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: "calendarDumps" ',
0994 & 'is no longer allowed in file "data"'
b9306711eb Jean*0995 CALL PRINT_ERROR( msgBuf, myThid )
742cf4499c Jean*0996 WRITE(msgBuf,'(A,A)') 'S/R INI_PARMS: calendarDumps',
0997 & ' has moved to "data.cal"'
b9306711eb Jean*0998 CALL PRINT_ERROR( msgBuf, myThid )
742cf4499c Jean*0999 ENDIF
46979d29da Jean*1000
7a7a4899b4 Chri*1001
1002
062a876ce5 Jean*1003 IF ( deltaTtracer .NE. dTtracerLev(1) .AND.
1004 & deltaTtracer .NE. 0. .AND. dTtracerLev(1) .NE. 0. ) THEN
1005 WRITE(msgBuf,'(A)')
1006 & 'S/R INI_PARMS: deltaTtracer & dTtracerLev(1) not equal'
b9306711eb Jean*1007 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1008 errCount = errCount + 1
062a876ce5 Jean*1009 ELSEIF ( dTtracerLev(1) .NE. 0. ) THEN
1010 deltaTtracer = dTtracerLev(1)
1011 ENDIF
ad330e45ef Jean*1012 IF ( deltaT .EQ. 0. ) deltaT = deltaTClock
7a7a4899b4 Chri*1013 IF ( deltaT .EQ. 0. ) deltaT = deltaTtracer
db8d49beca Jean*1014 IF ( deltaT .EQ. 0. ) deltaT = deltaTMom
1015 IF ( deltaT .EQ. 0. ) deltaT = deltaTFreeSurf
90aed06079 Jean*1016 IF ( deltaT .EQ. 0. ) THEN
1017 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS: ',
1018 & 'need to specify in file "data", namelist "PARM03"'
1019 CALL PRINT_ERROR( msgBuf, myThid )
1020 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS: ',
1021 & ' a model timestep (in s) deltaT or deltaTClock= ?'
1022 CALL PRINT_ERROR( msgBuf, myThid )
1023 errCount = errCount + 1
1024 deltaT = 1.
1025 ENDIF
db8d49beca Jean*1026 IF ( deltaTMom .EQ. 0. ) deltaTMom = deltaT
7a7a4899b4 Chri*1027 IF ( deltaTtracer .EQ. 0. ) deltaTtracer = deltaT
66dc79a095 Chri*1028 IF ( deltaTClock .EQ. 0. ) deltaTClock = deltaT
062a876ce5 Jean*1029 DO k=1,Nr
1030 IF (dTtracerLev(k).EQ.0.) dTtracerLev(k)= deltaTtracer
1031 ENDDO
ad330e45ef Jean*1032
d48ed2b20a Alis*1033
1034
1035
1036
db8d49beca Jean*1037 IF ( deltaTFreeSurf .EQ. 0. ) deltaTFreeSurf = deltaTMom
4edf45584c Alis*1038 IF ( periodicExternalForcing ) THEN
1039 IF ( externForcingCycle*externForcingPeriod .EQ. 0. ) THEN
1040 WRITE(msgBuf,'(A)')
1041 & 'S/R INI_PARMS: externForcingCycle,externForcingPeriod =0'
b9306711eb Jean*1042 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1043 errCount = errCount + 1
1044 ELSEIF ( INT(externForcingCycle/externForcingPeriod) .NE.
1045 & externForcingCycle/externForcingPeriod ) THEN
4edf45584c Alis*1046 WRITE(msgBuf,'(A)')
1047 & 'S/R INI_PARMS: externForcingCycle <> N*externForcingPeriod'
b9306711eb Jean*1048 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1049 errCount = errCount + 1
1050 ELSEIF ( externForcingCycle.LT.externForcingPeriod ) THEN
4edf45584c Alis*1051 WRITE(msgBuf,'(A)')
1052 & 'S/R INI_PARMS: externForcingCycle < externForcingPeriod'
b9306711eb Jean*1053 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1054 errCount = errCount + 1
4edf45584c Alis*1055 ENDIF
db8d49beca Jean*1056 IF ( externForcingPeriod.LT.deltaTClock ) THEN
4edf45584c Alis*1057 WRITE(msgBuf,'(A)')
db8d49beca Jean*1058 & 'S/R INI_PARMS: externForcingPeriod < deltaTClock'
b9306711eb Jean*1059 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1060 errCount = errCount + 1
4edf45584c Alis*1061 ENDIF
1062 ENDIF
c2b6ed6bfd Jean*1063
b9306711eb Jean*1064 IF ( momForcingOutAB .EQ. UNSET_I ) THEN
c2b6ed6bfd Jean*1065 momForcingOutAB = 1
1066 IF ( forcing_In_AB ) momForcingOutAB = 0
1067 ENDIF
b9306711eb Jean*1068 IF ( tracForcingOutAB .EQ. UNSET_I ) THEN
c2b6ed6bfd Jean*1069 tracForcingOutAB = 1
1070 IF ( forcing_In_AB ) tracForcingOutAB = 0
1071 ENDIF
cf8488c0fd Chri*1072
1073 IF ( cAdjFreq .LT. 0. ) THEN
1074 cAdjFreq = deltaTClock
1075 ENDIF
d4701cb6da Alis*1076 IF ( ivdc_kappa .NE. 0. .AND. cAdjFreq .NE. 0. ) THEN
1077 WRITE(msgBuf,'(A,A)')
1078 & 'S/R INI_PARMS: You have enabled both ivdc_kappa and',
1079 & ' convective adjustment.'
b9306711eb Jean*1080 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1081 errCount = errCount + 1
d4701cb6da Alis*1082 ENDIF
a97467b673 Jean*1083 IF (useCDscheme) THEN
1084
db8d49beca Jean*1085 IF ( tauCD .EQ. 0. _d 0 ) tauCD = deltaTMom
a97467b673 Jean*1086 IF ( rCD .LT. 0. ) rCD = 1. _d 0 - deltaTMom/tauCD
8039e9b985 Jean*1087 IF ( epsAB_CD .EQ. UNSET_RL ) epsAB_CD = abEps
b05b067368 Chri*1088 ENDIF
88830be691 Alis*1089
19eda8b250 Jean*1090 IF ( startTime.EQ.UNSET_RL .AND. nIter0.EQ.-1 ) THEN
1091
1092 startTime = baseTime
1093 nIter0 = 0
1094 ELSEIF ( startTime.EQ.UNSET_RL ) THEN
1095
d987e9fc35 Jean*1096 startTime = baseTime + deltaTClock*DFLOAT(nIter0)
19eda8b250 Jean*1097 ELSEIF ( nIter0.EQ.-1 ) THEN
1098
1099 nIter0 = NINT( (startTime-baseTime)/deltaTClock )
1100 ELSEIF ( baseTime.EQ.0. ) THEN
1101
d987e9fc35 Jean*1102 baseTime = startTime - deltaTClock*DFLOAT(nIter0)
19eda8b250 Jean*1103 ENDIF
d4701cb6da Alis*1104
6e6f314aa7 Patr*1105 nTimeSteps_l2 = 4
d4701cb6da Alis*1106
1107 IF ( nTimeSteps .EQ. 0 .AND. nEndIter .NE. 0 )
1108 & nTimeSteps = nEndIter-nIter0
1109
88830be691 Alis*1110 IF ( nTimeSteps .EQ. 0 .AND. endTime .NE. 0. )
db8d49beca Jean*1111 & nTimeSteps = NINT((endTime-startTime)/deltaTClock)
d4701cb6da Alis*1112
1113 IF ( nEndIter .EQ. 0 .AND. nTimeSteps .NE. 0 )
1114 & nEndIter = nIter0+nTimeSteps
1115
1116 IF ( nEndIter .EQ. 0 .AND. endTime .NE. 0. )
db8d49beca Jean*1117 & nEndIter = NINT((endTime-baseTime)/deltaTClock)
d4701cb6da Alis*1118
1119 IF ( endTime .EQ. 0. .AND. nTimeSteps .NE. 0 )
d987e9fc35 Jean*1120 & endTime = startTime + deltaTClock*DFLOAT(nTimeSteps)
d4701cb6da Alis*1121
1122 IF ( endTime .EQ. 0. .AND. nEndIter .NE. 0 )
d987e9fc35 Jean*1123 & endTime = baseTime + deltaTClock*DFLOAT(nEndIter)
d4701cb6da Alis*1124
88830be691 Alis*1125
d987e9fc35 Jean*1126 IF ( startTime .NE. baseTime+deltaTClock*DFLOAT(nIter0) ) THEN
16c445e700 Jean*1127 WRITE(msgBuf,'(A)')
1128 & 'S/R INI_PARMS: startTime, baseTime and nIter0 are inconsistent'
b9306711eb Jean*1129 CALL PRINT_ERROR( msgBuf, myThid )
16c445e700 Jean*1130 WRITE(msgBuf,'(A)')
1131 & 'S/R INI_PARMS: Perhaps more than two were set at once'
b9306711eb Jean*1132 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1133 errCount = errCount + 1
16c445e700 Jean*1134 ENDIF
d4701cb6da Alis*1135 IF ( nEndIter .NE. nIter0+nTimeSteps ) THEN
1136 WRITE(msgBuf,'(A)')
1137 & 'S/R INI_PARMS: nIter0, nTimeSteps and nEndIter are inconsistent'
b9306711eb Jean*1138 CALL PRINT_ERROR( msgBuf, myThid )
d4701cb6da Alis*1139 WRITE(msgBuf,'(A)')
1140 & 'S/R INI_PARMS: Perhaps more than two were set at once'
b9306711eb Jean*1141 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1142 errCount = errCount + 1
d4701cb6da Alis*1143 ENDIF
b9306711eb Jean*1144 IF ( nTimeSteps .NE. NINT((endTime-startTime)/deltaTClock)
1145 & ) THEN
1146 WRITE(msgBuf,'(A)')
1147 & 'S/R INI_PARMS: both endTime and nTimeSteps have been set'
1148 CALL PRINT_ERROR( msgBuf, myThid )
1149 WRITE(msgBuf,'(A)')
1150 & 'S/R INI_PARMS: but are inconsistent'
1151 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1152 errCount = errCount + 1
7a7a4899b4 Chri*1153 ENDIF
ea212a0263 Alis*1154
dc684458c1 Alis*1155
1156 IF (monitorFreq.LT.0.) THEN
1157 monitorFreq=0.
7abc7f9caf Alis*1158 IF (dumpFreq.NE.0.) monitorFreq=dumpFreq
3a279374db Alis*1159 IF (diagFreq.NE.0..AND.diagFreq.LT.monitorFreq)
1160 & monitorFreq=diagFreq
7abc7f9caf Alis*1161 IF (taveFreq.NE.0..AND.taveFreq.LT.monitorFreq)
1162 & monitorFreq=taveFreq
1163 IF (chkPtFreq.NE.0..AND.chkPtFreq.LT.monitorFreq)
1164 & monitorFreq=chkPtFreq
1165 IF (pChkPtFreq.NE.0..AND.pChkPtFreq.LT.monitorFreq)
1166 & monitorFreq=pChkPtFreq
db8d49beca Jean*1167 IF (monitorFreq.EQ.0.) monitorFreq=deltaTClock
dc684458c1 Alis*1168 ENDIF
f804abbd25 Jean*1169 IF ( monitorSelect.EQ.UNSET_I ) THEN
1170 monitorSelect = 2
1171 IF ( fluidIsWater ) monitorSelect = 3
1172 ENDIF
dc684458c1 Alis*1173
924557e60a Chri*1174
1175
7cd7754949 Jean*1176 DO k =1,Nr
1177 delZ(k) = UNSET_RL
1178 delP(k) = UNSET_RL
1179 delR(k) = UNSET_RL
924557e60a Chri*1180 ENDDO
1181
0127add478 Alis*1182 dxSpacing = UNSET_RL
1183 dySpacing = UNSET_RL
bb2fd3f1ad Jean*1184
1185 interViscAr_pCell = .FALSE.
1186 interDiffKr_pCell = .FALSE.
1187 pCellMix_select = 0
1188 pCellMix_maxFac = 1. _d 4
1189 pCellMix_delR = 0.
1190 DO k=1,Nr
1191 pCellMix_viscAr(k) = viscArNr(k)
1192 pCellMix_diffKr(k) = diffKrNrT(k)
1193 ENDDO
1194
aa8f58c33d Jean*1195 WRITE(msgBuf,'(A)') ' INI_PARMS ; starts to read PARM04'
6b6f69d6b0 Jean*1196 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*1197 & SQUEEZE_RIGHT, myThid )
b6e98c1b91 Jean*1198 READ(UNIT=iUnit,NML=PARM04)
c206b07c1d Chri*1199 IF ( errIO .LT. 0 ) THEN
924557e60a Chri*1200 WRITE(msgBuf,'(A)')
aa8f58c33d Jean*1201 & 'S/R INI_PARMS: Error reading model parameter file "data"'
b9306711eb Jean*1202 CALL PRINT_ERROR( msgBuf, myThid )
aa8f58c33d Jean*1203 WRITE(msgBuf,'(A)') 'S/R INI_PARMS: Problem in namelist PARM04'
b9306711eb Jean*1204 CALL PRINT_ERROR( msgBuf, myThid )
924557e60a Chri*1205 STOP 'ABNORMAL END: S/R INI_PARMS'
41f1ac6de1 Jean*1206 ELSE
aa8f58c33d Jean*1207 WRITE(msgBuf,'(A)') ' INI_PARMS ; read PARM04 : OK'
41f1ac6de1 Jean*1208 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*1209 & SQUEEZE_RIGHT, myThid )
809bdccbfc Jean*1210 ENDIF
9bec0e0933 Alis*1211
da40681a8e Jean*1212
2ad79bdf32 Jean*1213 IF ( Ro_SeaLevel .NE. UNSET_RL ) THEN
1214
1215 IF ( usingPCoords ) THEN
1216 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: ',
1217 & '"Ro_SeaLevel" (P @ bottom) depreciated (backward compat'
1218 ELSEIF ( usingZCoords ) THEN
1219 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: ',
1220 & '"Ro_SeaLevel" (Z @ top) depreciated (backward compat'
1221 ENDIF
1222 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
1223 & SQUEEZE_RIGHT, myThid )
1224 IF ( usingPCoords ) THEN
1225 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: ',
1226 & ' only). To set vert. axis, use instead "top_Pres".'
1227 ELSEIF ( usingZCoords ) THEN
1228 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: ',
1229 & ' only). To set vert. axis, use instead "seaLev_Z".'
1230 ENDIF
1231 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
1232 & SQUEEZE_RIGHT, myThid )
1233 ENDIF
da40681a8e Jean*1234 IF ( rkFac .NE. UNSET_RL ) THEN
1235 nRetired = nRetired+1
1236 WRITE(msgBuf,'(A,A)')
1237 & 'S/R INI_PARMS: "rkFac" has been replaced by -rkSign ',
1238 & ' and is no longer allowed in file "data".'
b9306711eb Jean*1239 CALL PRINT_ERROR( msgBuf, myThid )
da40681a8e Jean*1240 ENDIF
1241 IF ( groundAtK1 ) THEN
1242
1243 WRITE(msgBuf,'(A,A)')
1244 & 'S/R INI_PARMS: "groundAtK1" is set according to vertical ',
1245 & ' coordinate and is no longer allowed in file "data".'
b9306711eb Jean*1246 CALL PRINT_ERROR( msgBuf, myThid )
da40681a8e Jean*1247 ENDIF
9780090eaa Jean*1248 IF ( thetaMin .NE. UNSET_RL ) THEN
1249 nRetired = nRetired+1
1250 WRITE(msgBuf,'(A,A)')
1251 & 'S/R INI_PARMS: "thetaMin" no longer allowed,',
1252 & ' has been replaced by "xgOrigin"'
1253 CALL PRINT_ERROR( msgBuf, myThid )
1254 ENDIF
1255 IF ( phiMin .NE. UNSET_RL ) THEN
1256 nRetired = nRetired+1
1257 WRITE(msgBuf,'(A,A)')
1258 & 'S/R INI_PARMS: "phiMin" no longer allowed,',
1259 & ' has been replaced by "ygOrigin"'
1260 CALL PRINT_ERROR( msgBuf, myThid )
1261 ENDIF
da40681a8e Jean*1262
f76021331f Jean*1263
1264 goptCount = 0
1265 IF ( delX(1) .NE. UNSET_RL ) goptCount = goptCount + 1
1266 IF ( dxSpacing .NE. UNSET_RL ) goptCount = goptCount + 1
1267 IF ( delXFile .NE. ' ' ) goptCount = goptCount + 1
1268 IF ( goptCount.GT.1 ) THEN
b9306711eb Jean*1269 WRITE(msgBuf,'(A,A)') 'Too many specifications for delX:',
9bec0e0933 Alis*1270 & 'Specify only one of delX, dxSpacing or delXfile'
b9306711eb Jean*1271 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1272 errCount = errCount + 1
9bec0e0933 Alis*1273 ENDIF
0127add478 Alis*1274 IF ( dxSpacing .NE. UNSET_RL ) THEN
b9b591469d Jean*1275 DO i=1,gridNx
0127add478 Alis*1276 delX(i) = dxSpacing
1277 ENDDO
1278 ENDIF
f76021331f Jean*1279
1280 goptCount = 0
1281 IF ( delY(1) .NE. UNSET_RL ) goptCount = goptCount + 1
1282 IF ( dySpacing .NE. UNSET_RL ) goptCount = goptCount + 1
1283 IF ( delYFile .NE. ' ' ) goptCount = goptCount + 1
1284 IF ( goptCount.GT.1 ) THEN
b9306711eb Jean*1285 WRITE(msgBuf,'(A,A)') 'Too many specifications for delY:',
9bec0e0933 Alis*1286 & 'Specify only one of delY, dySpacing or delYfile'
b9306711eb Jean*1287 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1288 errCount = errCount + 1
9bec0e0933 Alis*1289 ENDIF
0127add478 Alis*1290 IF ( dySpacing .NE. UNSET_RL ) THEN
b9b591469d Jean*1291 DO j=1,gridNy
f76021331f Jean*1292 delY(j) = dySpacing
0127add478 Alis*1293 ENDDO
1294 ENDIF
2c7e43108f Jean*1295
910f05e765 Chri*1296
924557e60a Chri*1297 goptCount = 0
1298 IF ( usingCartesianGrid ) goptCount = goptCount+1
1299 IF ( usingSphericalPolarGrid ) goptCount = goptCount+1
aea29c8517 Alis*1300 IF ( usingCurvilinearGrid ) goptCount = goptCount+1
0ac260a803 Andr*1301 IF ( usingCylindricalGrid ) goptCount = goptCount+1
aea29c8517 Alis*1302 IF ( goptCount .GT. 1 ) THEN
924557e60a Chri*1303 WRITE(msgBuf,'(A)')
1304 & 'S/R INI_PARMS: More than one coordinate system requested'
b9306711eb Jean*1305 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1306 errCount = errCount + 1
924557e60a Chri*1307 ENDIF
aea29c8517 Alis*1308 IF ( goptCount .LT. 1 ) THEN
f2459ca318 Jean*1309
aea29c8517 Alis*1310 WRITE(msgBuf,'(A)')
f2459ca318 Jean*1311 & 'S/R INI_PARMS: No horizontal grid requested'
1312 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
b9306711eb Jean*1313 & SQUEEZE_RIGHT, myThid )
f2459ca318 Jean*1314 WRITE(msgBuf,'(A)')
1315 & 'S/R INI_PARMS: => Use Cartesian Grid as default'
1316 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
b9306711eb Jean*1317 & SQUEEZE_RIGHT, myThid )
f2459ca318 Jean*1318 usingCartesianGrid = .TRUE.
aea29c8517 Alis*1319 ENDIF
9744e36521 Jean*1320
1321 IF ( rSphere.EQ.UNSET_RL ) THEN
1322 IF ( usingCurvilinearGrid .AND.
1323 & radius_fromHorizGrid.NE.UNSET_RL ) THEN
1324 rSphere = radius_fromHorizGrid
1325 ELSE
1326 rSphere = 6370. _d 3
1327 ENDIF
1328 ENDIF
1329 IF ( radius_fromHorizGrid.EQ.UNSET_RL ) THEN
1330 radius_fromHorizGrid = rSphere
1331 ENDIF
1332 IF ( rSphere .NE. 0. ) THEN
1333 recip_rSphere = 1. _d 0/rSphere
1334 ELSE
1335 recip_rSphere = 0.
1336 ENDIF
427e24e121 Jean*1337
1338 IF ( phiEuler .NE. 0. _d 0 .OR. thetaEuler .NE. 0. _d 0
1339 & .OR. psiEuler .NE. 0. _d 0 ) rotateGrid = .TRUE.
1340
2ad79bdf32 Jean*1341
1342 IF ( Ro_SeaLevel .NE. UNSET_RL ) THEN
1343 IF ( usingPCoords .AND. top_Pres.NE.UNSET_RL ) THEN
ef7fa20135 Jean*1344 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS: ',
1345 & 'Cannot set both "Ro_SeaLevel" and "top_Pres"'
1346 CALL PRINT_ERROR( msgBuf, myThid )
1347 errCount = errCount + 1
2ad79bdf32 Jean*1348 ENDIF
1349 IF ( usingZCoords .AND. seaLev_Z.NE.UNSET_RL ) THEN
ef7fa20135 Jean*1350 WRITE(msgBuf,'(2A)') 'S/R INI_PARMS: ',
1351 & 'Cannot set both "Ro_SeaLevel" and "seaLev_Z"'
1352 CALL PRINT_ERROR( msgBuf, myThid )
1353 errCount = errCount + 1
2ad79bdf32 Jean*1354 ENDIF
1355 rF(1) = Ro_SeaLevel
1356 ELSE
1357 rF(1) = UNSET_RS
1358 ENDIF
1359 IF ( top_Pres.EQ.UNSET_RL ) top_Pres = 0.
1360 IF ( seaLev_Z.EQ.UNSET_RL ) seaLev_Z = 0.
9780090eaa Jean*1361
1362 IF ( xgOrigin .EQ. UNSET_RL ) xgOrigin = 0.
1363 IF ( ygOrigin .EQ. UNSET_RL ) THEN
1364 IF ( usingSphericalPolarGrid ) THEN
1365 ygOrigin = 0.
1366 ELSEIF ( usingCartesianGrid ) THEN
1367 ygOrigin = 0.
1368 ELSEIF ( usingCylindricalGrid ) THEN
1369 ygOrigin = 0.
1370 ELSEIF ( usingCurvilinearGrid ) THEN
1371 ygOrigin = 0.
1372 ELSE
1373 WRITE(msgBuf,'(A)')
1374 & 'S/R INI_PARMS: found no coordinate system to set ygOrigin'
1375 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1376 errCount = errCount + 1
9780090eaa Jean*1377 ENDIF
1378 ENDIF
7fbd8d8c1c Jean*1379
c28ce1627a Jean*1380
1381 setCenterDr = .FALSE.
7cd7754949 Jean*1382 DO k=1,Nr+1
1383 IF ( delRc(k).EQ.UNSET_RL ) THEN
b9306711eb Jean*1384 IF ( setCenterDr ) THEN
1385 WRITE(msgBuf,'(A,I4)')
7cd7754949 Jean*1386 & 'S/R INI_PARMS: No value for delRc at k =', k
b9306711eb Jean*1387 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1388 errCount = errCount + 1
b9306711eb Jean*1389 ENDIF
1390 ELSE
1391 IF ( k.EQ.1 ) setCenterDr = .TRUE.
1392 IF ( .NOT.setCenterDr ) THEN
1393 WRITE(msgBuf,'(A,I4)')
7cd7754949 Jean*1394 & 'S/R INI_PARMS: No value for delRc at k <', k
b9306711eb Jean*1395 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1396 errCount = errCount + 1
b9306711eb Jean*1397 ENDIF
1398 ENDIF
c28ce1627a Jean*1399 ENDDO
b9306711eb Jean*1400 IF ( setCenterDr ) rCoordInputData = .TRUE.
910f05e765 Chri*1401
b9306711eb Jean*1402 setInterFDr = .FALSE.
7cd7754949 Jean*1403 DO k = 1, Nr
1404 IF ( delZ(k) .NE. UNSET_RL ) zCoordInputData = .TRUE.
1405 IF ( delP(k) .NE. UNSET_RL ) pCoordInputData = .TRUE.
1406 IF ( delR(k) .NE. UNSET_RL ) rCoordInputData = .TRUE.
1407 IF ( delR(k) .EQ. UNSET_RL ) delR(k) = delZ(k)
1408 IF ( delR(k) .EQ. UNSET_RL ) delR(k) = delP(k)
1409 IF ( delR(k) .EQ. UNSET_RL ) THEN
b9306711eb Jean*1410 IF ( setInterFDr ) THEN
88830be691 Alis*1411 WRITE(msgBuf,'(A,I4)')
7cd7754949 Jean*1412 & 'S/R INI_PARMS: No value for delZ/delP/delR at k =', k
b9306711eb Jean*1413 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1414 errCount = errCount + 1
b9306711eb Jean*1415 ENDIF
1416 ELSE
1417 IF ( k.EQ.1 ) setInterFDr = .TRUE.
1418 IF ( .NOT.setInterFDr ) THEN
1419 WRITE(msgBuf,'(A,I4)')
7cd7754949 Jean*1420 & 'S/R INI_PARMS: No value for delZ/delP/delR at k <', k
b9306711eb Jean*1421 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1422 errCount = errCount + 1
b9306711eb Jean*1423 ENDIF
88830be691 Alis*1424 ENDIF
910f05e765 Chri*1425 ENDDO
1426
b9306711eb Jean*1427 coordsSet = 0
910f05e765 Chri*1428 IF ( zCoordInputData ) coordsSet = coordsSet + 1
1429 IF ( pCoordInputData ) coordsSet = coordsSet + 1
1430 IF ( rCoordInputData ) coordsSet = coordsSet + 1
1431 IF ( coordsSet .GT. 1 ) THEN
1432 WRITE(msgBuf,'(A)')
1433 & 'S/R INI_PARMS: Cannot mix z, p and r in the input data.'
b9306711eb Jean*1434 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1435 errCount = errCount + 1
910f05e765 Chri*1436 ENDIF
b9306711eb Jean*1437
1438 IF ( delRcFile.NE.' ' ) THEN
1439 IF ( setCenterDr ) THEN
1440 WRITE(msgBuf,'(A)')
1441 & 'S/R INI_PARMS: Cannot set both delRc and delRcFile'
1442 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1443 errCount = errCount + 1
b9306711eb Jean*1444 ENDIF
1445 setCenterDr = .TRUE.
1446 ENDIF
1447 IF ( delRFile.NE.' ' ) THEN
1448 IF ( setInterFDr ) THEN
1449 WRITE(msgBuf,'(A)')
1450 & 'S/R INI_PARMS: Cannot set both delR and delRFile'
1451 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1452 errCount = errCount + 1
b9306711eb Jean*1453 ENDIF
1454 setInterFDr = .TRUE.
1455 ENDIF
b05b067368 Chri*1456
81bc00c2f0 Chri*1457
aa8f58c33d Jean*1458 WRITE(msgBuf,'(A)') ' INI_PARMS ; starts to read PARM05'
6b6f69d6b0 Jean*1459 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*1460 & SQUEEZE_RIGHT, myThid )
88830be691 Alis*1461 READ(UNIT=iUnit,NML=PARM05)
809bdccbfc Jean*1462 IF ( errIO .LT. 0 ) THEN
81bc00c2f0 Chri*1463 WRITE(msgBuf,'(A)')
aa8f58c33d Jean*1464 & 'S/R INI_PARMS: Error reading model parameter file "data"'
b9306711eb Jean*1465 CALL PRINT_ERROR( msgBuf, myThid )
aa8f58c33d Jean*1466 WRITE(msgBuf,'(A)') 'S/R INI_PARMS: Problem in namelist PARM05'
b9306711eb Jean*1467 CALL PRINT_ERROR( msgBuf, myThid )
81bc00c2f0 Chri*1468 STOP 'ABNORMAL END: S/R INI_PARMS'
41f1ac6de1 Jean*1469 ELSE
aa8f58c33d Jean*1470 WRITE(msgBuf,'(A)') ' INI_PARMS ; read PARM05 : OK'
41f1ac6de1 Jean*1471 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
b9306711eb Jean*1472 & SQUEEZE_RIGHT, myThid )
809bdccbfc Jean*1473 ENDIF
8194bc4a99 Mart*1474
1475 IF ( shelfIceFile .NE. ' ' ) THEN
1476 nRetired = nRetired+1
1477 WRITE(msgBuf,'(A,A)')
1478 & 'S/R INI_PARMS: "shelfIceFile" is not allowed in "data", ',
1479 & 'substitute "SHELFICEtopoFile" in data.shelfice'
1480 CALL PRINT_ERROR( msgBuf, myThid )
1481 ENDIF
275ad78b78 Jean*1482 IF ( dQdTFile .NE. ' ' ) THEN
1483 nRetired = nRetired+1
1484 WRITE(msgBuf,'(A,A)')
1485 & 'S/R INI_PARMS: "dQdTFile" has been retired from file "data"'
1486 CALL PRINT_ERROR( msgBuf, myThid )
1487 ENDIF
cbc417d429 Mart*1488
1489 IF (mdsioLocalDir .NE. ' ' .AND. adTapeDir .NE. ' ') THEN
1490 WRITE(msgBuf,'(A)')
1491 & 'S/R INI_PARMS: mdsioLocalDir and adTapeDir cannot be'
1492 CALL PRINT_ERROR( msgBuf, myThid )
1493 WRITE(msgBuf,'(A)')
1494 & 'S/R INI_PARMS: specified at the same time'
1495 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1496 errCount = errCount + 1
cbc417d429 Mart*1497 ENDIF
81bc00c2f0 Chri*1498
7b50f71112 Gael*1499
1500 IF ( (tauThetaClimRelax.GT.0.).AND.
1501 & (thetaClimFile.EQ.' ') ) THEN
1502 WRITE(msgBuf,'(A)')
1503 & 'S/R INI_PARMS: tauThetaClimRelax > 0 but'
1504 CALL PRINT_ERROR( msgBuf, myThid )
1505 WRITE(msgBuf,'(A)')
1506 & 'S/R INI_PARMS: thetaClimFile is undefined'
1507 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1508 errCount = errCount + 1
7b50f71112 Gael*1509 ENDIF
1510 IF ( (tauSaltClimRelax.GT.0.).AND.
1511 & (saltClimFile.EQ.' ') ) THEN
1512 WRITE(msgBuf,'(A)')
1513 & 'S/R INI_PARMS: tauSaltClimRelax > 0 but'
1514 CALL PRINT_ERROR( msgBuf, myThid )
1515 WRITE(msgBuf,'(A)')
1516 & 'S/R INI_PARMS: saltClimFile is undefined'
1517 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1518 errCount = errCount + 1
7b50f71112 Gael*1519 ENDIF
13b4d13c98 Jean*1520
1521 #ifdef ALLOW_3D_DIFFKR
1522 IF ( specifiedDiffKrT ) THEN
3fcdbd6e66 Jean*1523 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: Ignores diffKr',
1524 & 'T (or Kp,Kz) setting in file "data" with ALLOW_3D_DIFFKR'
1525 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
1526 & SQUEEZE_RIGHT, myThid )
13b4d13c98 Jean*1527 ENDIF
1528 IF ( specifiedDiffKrS .AND. diffKrFile.NE.' ' ) THEN
3fcdbd6e66 Jean*1529 WRITE(msgBuf,'(2A)') '** WARNING ** INI_PARMS: Ignores diffKr',
1530 & 'S (or Kp,Kz) setting in file "data" and uses diffKrFile'
1531 CALL PRINT_MESSAGE( msgBuf, errorMessageUnit,
1532 & SQUEEZE_RIGHT, myThid )
13b4d13c98 Jean*1533 ENDIF
1534 #endif
7b50f71112 Gael*1535
599be48d9a Jean*1536
1537
427e24e121 Jean*1538
1539
e99e7fe4e6 Jean*1540 IF ( usingPCoords ) THEN
599be48d9a Jean*1541 mass2rUnit = gravity
1542 rUnit2mass = recip_gravity
1543 ELSE
1544 mass2rUnit = recip_rhoConst
1545 rUnit2mass = rhoConst
9800a0be87 Jean*1546 ENDIF
688d11fba8 Alis*1547
80d98e0151 Dimi*1548
b46f9da188 Jean*1549
80d98e0151 Dimi*1550 IF (temp_addMass .EQ. UNSET_RL) temp_addMass = temp_EvPrRn
1551 IF (salt_addMass .EQ. UNSET_RL) salt_addMass = salt_EvPrRn
1552
b46f9da188 Jean*1553
1554
1555
1556 IF ( cg2dUseMinResSol.EQ.UNSET_I ) THEN
1557 cg2dUseMinResSol = 0
1558 IF ( topoFile.EQ.' ' .AND. bathyFile.EQ.' '
1559 & .AND. usingCartesianGrid ) cg2dUseMinResSol = 1
1560 ENDIF
1561
aa8f58c33d Jean*1562
a6485b3244 Mart*1563 #ifdef SINGLE_DISK_IO
9ba4d7e4bd Chri*1564 CLOSE(iUnit)
a6485b3244 Mart*1565 #else
1566 CLOSE(iUnit,STATUS='DELETE')
1567 #endif /* SINGLE_DISK_IO */
9ba4d7e4bd Chri*1568
aa8f58c33d Jean*1569 WRITE(msgBuf,'(A)') ' INI_PARMS: finished reading file "data"'
1570 CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
1571 & SQUEEZE_RIGHT , myThid )
1572
ede2cfe569 Chri*1573
46979d29da Jean*1574 IF ( nRetired .GT. 0 ) THEN
ede2cfe569 Chri*1575 WRITE(msgBuf,'(A)')
90aed06079 Jean*1576 & 'S/R INI_PARMS: Error reading parameter file "data":'
b9306711eb Jean*1577 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1578 WRITE(msgBuf,'(I4,A)') nRetired,
1579 & ' out of date parameters were found in the namelist'
b9306711eb Jean*1580 CALL PRINT_ERROR( msgBuf, myThid )
90aed06079 Jean*1581 errCount = errCount + 1
1582 ENDIF
1583
1584 IF ( errCount .GE. 1 ) THEN
1585 WRITE(msgBuf,'(A,I3,A)')
1586 & 'S/R INI_PARMS: detected', errCount,' fatal error(s)'
1587 CALL PRINT_ERROR( msgBuf, myThid )
1588 CALL ALL_PROC_DIE( 0 )
1589 STOP 'ABNORMAL END: S/R INI_PARMS'
ede2cfe569 Chri*1590 ENDIF
1591
924557e60a Chri*1592 _END_MASTER(myThid)
1593
1594
1595 _BARRIER
2c7e43108f Jean*1596
924557e60a Chri*1597 RETURN
1598 END