/[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.47 by adcroft, Mon Aug 30 18:29:26 1999 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        IMPLICIT NONE  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"
 #ifdef ALLOW_KPP  
 #include "KPPMIX.h"  
 #endif  
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 45  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 62  C     dependency. Line 70  C     dependency.
70        _RL hFacMinDrDefault        _RL hFacMinDrDefault
71        _RL delRDefault(Nr)        _RL delRDefault(Nr)
72        _RS rkFacDefault        _RS rkFacDefault
73  C     zCoordInputData - These are used to select between different coordinate systems.  C     zCoordInputData :: Variables used to select between different coordinate systems.
74  C     pCoordInputData   The vertical coordinate system in the rest of the model is  C     pCoordInputData :: The vertical coordinate system in the rest of the model is
75  C     rCoordInputData   written in terms of r. In the model "data" file input data can  C     rCoordInputData :: written in terms of r. In the model "data" file input data can
76  C     coordsSet         be interms of z, p or r.  C     coordsSet       :: be interms of z, p or r.
77  C                       e.g. delZ or delP or delR for the vertical grid spacing.  C                     :: e.g. delZ or delP or delR for the vertical grid spacing.
78  C                       The following rules apply:  C                     :: The following rules apply:
79  C                       All parameters must use the same vertical coordinate system.  C                     :: All parameters must use the same vertical coordinate system.
80  C                       e.g. delZ and viscAz is legal but  C                     ::  e.g. delZ and viscAz is legal but
81  C                            delZ and viscAr is an error.  C                     ::       delZ and viscAr is an error.
82  C                       Similarly specifyinh delZ and delP is an error.  C                     :: Similarly specifyinh delZ and delP is an error.
83  C                       zCoord..., pCoord..., rCoord... are used to flag when z, p or r are  C                     :: zCoord..., pCoord..., rCoord... are used to flag when z, p or r are
84  C                       used. coordsSet counts how many vertical coordinate systems have been  C                     :: used. coordsSet counts how many vertical coordinate systems have been
85  C                       used to specify variables. coordsSet > 1 is an error.  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, cosPower,       & 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,       & 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,       & exactConserv,uniformLin_PhiSurf,nonlinFreeSurf,hFacInf,hFacSup,
117         & staggerTimeStep,
118         & tempStepping, saltStepping, momStepping, tr1Stepping,
119       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
120       & viscAr, diffKrT, diffKrS, hFacMinDr,       & viscAr, diffKrT, diffKrS, hFacMinDr,
121       & viscAp, diffKpT, diffKpS, hFacMinDp,       & viscAp, diffKpT, diffKpS, hFacMinDp,
122       & rhoConst, buoyancyRelation, HeatCapacity_Cp,       & rhoConst, buoyancyRelation, HeatCapacity_Cp,
123       & writeBinaryPrec, readBinaryPrec, writeStatePrec,       & writeBinaryPrec, readBinaryPrec, writeStatePrec,
124       & openBoundaries, nonHydrostatic, globalFiles,       & nonHydrostatic, globalFiles,
125       & allowFreezing, ivdc_kappa       & 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       & cg3dMaxIters, cg3dChkResFreq, cg3dTargetResidual
141    
142  C--   Time stepping parammeters  C--   Time stepping parammeters
143        NAMELIST /PARM03/        NAMELIST /PARM03/
144       & nIter0, nTimeSteps, nEndIter, deltaT, deltaTmom, deltaTtracer,       & nIter0, nTimeSteps, nEndIter,
145       & abEps, tauCD, rCD,       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf,
146       & startTime, endTime, chkPtFreq, dumpFreq, taveFreq, deltaTClock,       & forcing_In_AB, abEps, tauCD, rCD,
147       & pChkPtFreq, cAdjFreq, tauThetaClimRelax, tauSaltClimRelax,       & 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,       & zonalWindFile, meridWindFile,
165       & thetaClimFile, saltClimFile,       & thetaClimFile, saltClimFile,
166       & surfQfile, EmPmRfile,       & surfQfile, EmPmRfile, surfQswfile,
167       & uVelInitFile, vVelInitFile       & uVelInitFile, vVelInitFile, pSurfInitFile,
168         & dQdTFile, ploadFile
 C--   Open Boundaries  
       NAMELIST /PARM06/  
      & OB_Jnorth, OB_Jsouth, OB_Ieast, OB_Iwest  
   
 #ifdef ALLOW_KPP  
 C--   KPP vertical mixing parameters  
       NAMELIST /PARM07/  
      & usingKPPmixing, KPPmixingMaps, KPPwriteState,  
      & minKPPviscAz, maxKPPviscAz, minKPPghat, maxKPPghat,  
      & minKPPdiffKzT, maxKPPdiffKzT, minKPPdiffKzS, maxKPPdiffKzS,  
      & epsln, epsilon, vonk, conc1, conam, concm, conc2, zetam,  
      & conas, concs, conc3, zetas,  
      & Ricr, cekman, cmonob, concv, hbf, Vtc,  
      & zmin, zmax, umin, umax,  
      & num_v_smooth_Ri, num_v_smooth_BV,  
      & num_z_smooth_sh, num_m_smooth_sh,  
      & Riinfty, BVSQcon, difm0, difs0, dift0,  
      & difmiw, difsiw, diftiw, difmcon, difscon, diftcon,  
      & cstar, cg  
 #endif  
169    
170  C  C
171        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
# Line 169  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')
# Line 245  C--   Set default "physical" parameters Line 264  C--   Set default "physical" parameters
264        hFacMinDr           = UNSET_RL        hFacMinDr           = UNSET_RL
265        hFacMinDz           = UNSET_RL        hFacMinDz           = UNSET_RL
266        hFacMinDp           = UNSET_RL        hFacMinDp           = UNSET_RL
267          convertFW2Salt = UNSET_RL
268          tAlpha              = UNSET_RL
269          sBeta               = UNSET_RL
270        READ(UNIT=iUnit,NML=PARM01) !,IOSTAT=errIO)        READ(UNIT=iUnit,NML=PARM01) !,IOSTAT=errIO)
271        IF ( errIO .LT. 0 ) THEN        IF ( errIO .LT. 0 ) THEN
272         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
# Line 261  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          ELSE
287           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM01 : OK'
288           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
289         &                    SQUEEZE_RIGHT , 1)
290        ENDIF        ENDIF
291        IF ( implicitFreeSurface ) freeSurfFac = 1.D0        IF ( implicitFreeSurface ) freeSurfFac = 1.D0
292        IF ( rigidLid            ) freeSurfFac = 0.D0        IF ( rigidLid            ) freeSurfFac = 0.D0
# Line 300  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 327  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 (convertFW2Salt.EQ.UNSET_RL) THEN
366            convertFW2Salt = 35.
367            IF (useRealFreshWaterFlux) convertFW2Salt=-1
368          ENDIF
369    
370        IF ( ivdc_kappa .NE. 0. .AND. .NOT. implicitDiffusion ) THEN        IF ( ivdc_kappa .NE. 0. .AND. .NOT. implicitDiffusion ) THEN
371         WRITE(msgBuf,'(A,A)')          WRITE(msgBuf,'(A,A)')
372       &  'S/R INI_PARMS: To use ivdc_kappa you must enable implicit',       &  'S/R INI_PARMS: To use ivdc_kappa you must enable implicit',
373       &  ' vertical diffusion.'       &  ' 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    
       IF ( implicitFreeSurface .AND.  rigidLid ) THEN  
        WRITE(msgBuf,'(A,A)')  
      & 'S/R INI_PARMS: Cannot select both implicitFreeSurface',  
      & ' and rigidLid.'  
        CALL PRINT_ERROR( msgBuf , myThid)  
        STOP 'ABNORMAL END: S/R INI_PARMS'  
       ENDIF  
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 395  C     Set globalFiles flag for READ_WRIT Line 424  C     Set globalFiles flag for READ_WRIT
424  C     Set globalFiles flag for READ_WRITE_REC package  C     Set globalFiles flag for READ_WRITE_REC package
425        CALL SET_WRITE_GLOBAL_PICKUP( globalFiles )        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        READ(UNIT=iUnit,NML=PARM02) !,IOSTAT=errIO)        READ(UNIT=iUnit,NML=PARM02) !,IOSTAT=errIO)
443        IF ( errIO .LT. 0 ) THEN        IF ( errIO .LT. 0 ) THEN
# Line 412  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          ELSE
459           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM02 : OK'
460           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
461         &                    SQUEEZE_RIGHT , 1)
462        ENDIF            ENDIF    
463    
464  C--   Time stepping parameters  C--   Time stepping parameters
# Line 432  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          ELSE
483           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM03 : OK'
484           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
485         &                    SQUEEZE_RIGHT , 1)
486        ENDIF          ENDIF  
487  C     Process "timestepping" params  C     Process "timestepping" params
488  C     o Time step size  C     o Time step size
# Line 440  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 501  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 Tracer 1 climatology relaxation time scale
562          IF ( tauTr1ClimRelax .EQ. 0.D0 ) THEN
563           doTr1ClimRelax     = .FALSE.
564           lambdaTr1ClimRelax = 0.D0
565          ELSE
566           doTr1ClimRelax     = .TRUE.
567           lambdaTr1ClimRelax = 1./tauTr1ClimRelax
568          ENDIF
569    
570  C     o Start time  C     o Start time
571        IF ( nIter0 .NE. 0 .AND. startTime .EQ. 0. )        IF ( nIter0 .NE. 0 .AND. startTime .EQ. 0. )
# Line 549  C     o Consistent? Line 614  C     o Consistent?
614          STOP 'ABNORMAL END: S/R INI_PARMS'          STOP 'ABNORMAL END: S/R INI_PARMS'
615        ENDIF        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
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 INCLUDE_DIAGNOSTICS_INTERFACE_CODE  #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 560  C       code is being compiled Line 640  C       code is being compiled
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,A)')          WRITE(msgBuf,'(A,A)')
643       &  'Re-compile with:  #define INCLUDE_DIAGNOSTICS_INTERFACE_CODE',       &  'Re-compile with:  #define ALLOW_TIMEAVE',
644       &  '              or  -DINCLUDE_DIAGNOSTICS_INTERFACE_CODE'       &  '              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 579  C     In spherical polar distances are i Line 659  C     In spherical polar distances are i
659        recip_rSphere  = 1.D0/rSphere        recip_rSphere  = 1.D0/rSphere
660        dxSpacing = UNSET_RL        dxSpacing = UNSET_RL
661        dySpacing = UNSET_RL        dySpacing = UNSET_RL
662        READ(UNIT=iUnit,NML=PARM04) !,IOSTAT=errIO)        delXfile = ' '
663          delYfile = ' '
664          READ(UNIT=iUnit,NML=PARM04,IOSTAT=errIO)
665        IF ( errIO .LT. 0 ) THEN        IF ( errIO .LT. 0 ) THEN
666         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
667       &  'S/R INI_PARMS'       &  'S/R INI_PARMS'
# Line 595  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'
680          ELSE
681           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM04 : OK'
682           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
683         &                    SQUEEZE_RIGHT , 1)
684        ENDIF            ENDIF    
685  C  
686    C     X coordinate
687          IF ( delXfile .NE. ' ' ) THEN
688           IF ( delX(1) .NE. UNSET_RL .OR. dxSpacing .NE. UNSET_RL ) THEN
689             WRITE(msgBuf,'(A,A)') 'Too many specifications for delX:',
690         &   'Specify only one of delX, dxSpacing or delXfile'
691            CALL PRINT_ERROR( msgBuf , 1)
692            STOP 'ABNORMAL END: S/R INI_PARMS'
693           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        IF ( dxSpacing .NE. UNSET_RL ) THEN
716         DO i=1,Nx         DO i=1,Nx
717          delX(i) = dxSpacing          delX(i) = dxSpacing
718         ENDDO         ENDDO
719        ENDIF        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)
726            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
748          ENDIF
749        IF ( dySpacing .NE. UNSET_RL ) THEN        IF ( dySpacing .NE. UNSET_RL ) THEN
750         DO j=1,Ny         DO i=1,Ny
751          delY(j) = dySpacing          delY(i) = dySpacing
752         ENDDO         ENDDO
753        ENDIF        ENDIF
754    C
755        IF ( rSphere .NE. 0 ) THEN        IF ( rSphere .NE. 0 ) THEN
756         recip_rSphere = 1.D0/rSphere         recip_rSphere = 1.D0/rSphere
757        ELSE        ELSE
758         recip_rSphere = 0.         recip_rSphere = 0.
759        ENDIF        ENDIF
 C--   Initialize EOS coefficients (3rd order polynomial)  
       IF (eostype.eq.'POLY3') THEN  
        OPEN(37,FILE='POLY3.COEFFS',STATUS='OLD',FORM='FORMATTED')  
        READ(37,*) I  
        IF (I.NE.Nr) THEN  
         WRITE(msgBuf,'(A)')  
      &  'ini_parms: attempt to read POLY3.COEFFS failed'  
         CALL PRINT_ERROR( msgBuf , 1)  
         WRITE(msgBuf,'(A)')  
      &  '           because bad # of levels in data'  
         CALL PRINT_ERROR( msgBuf , 1)  
         STOP 'Bad data in POLY3.COEFFS'  
        ENDIF  
        READ(37,*) (eosRefT(K),eosRefS(K),eosSig0(K),K=1,Nr)  
        DO K=1,Nr  
         READ(37,*) (eosC(I,K),I=1,9)  
        ENDDO  
        CLOSE(37)  
       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 665  C--   p, z, r coord parameters Line 804  C--   p, z, r coord parameters
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(K)         IF ( delR(K) .EQ. UNSET_RL ) delR(K) = delRDefault(K)
807         IF ( delR(K) .EQ. 0. ) THEN         IF (.NOT.setCenterDr .AND. delR(K).EQ.delRDefault(K) ) THEN
808           WRITE(msgBuf,'(A,I4)')           WRITE(msgBuf,'(A,I4)')
809       &  'S/R INI_PARMS: No value for delZ/delP/delR at K = ',K       &  'S/R INI_PARMS: No value for delZ/delP/delR at K = ',K
810           CALL PRINT_ERROR( msgBuf , 1)           CALL PRINT_ERROR( msgBuf , 1)
811           STOP 'ABNORMAL END: S/R INI_PARMS'           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         ENDIF
818        ENDDO        ENDDO
819  C     Check for multiple coordinate systems  C     Check for multiple coordinate systems
# Line 688  C--   Input files Line 832  C--   Input files
832        READ(UNIT=iUnit,NML=PARM05) !,IOSTAT=errIO)        READ(UNIT=iUnit,NML=PARM05) !,IOSTAT=errIO)
833        IF ( errIO .LT. 0 ) THEN            IF ( errIO .LT. 0 ) THEN    
834         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
      &  'S/R INI_PARMS'  
        CALL PRINT_ERROR( msgBuf , 1)  
        WRITE(msgBuf,'(A)')  
835       &  'Error reading numerical model '       &  'Error reading numerical model '
836         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
837         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
# Line 701  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          ELSE
846           WRITE(msgBuf,'(A)') 'S/R INI_PARMS ; read PARM05 : OK'
847           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
848         &                    SQUEEZE_RIGHT , 1)
849        ENDIF              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 724  C     being combined and a single frame Line 886  C     being combined and a single frame
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        IF ( rkFac .EQ. UNSET_RS ) rkFac=rkFacDefault
894        recip_rkFac = 1.D0 / rkFac        recip_rkFac = 1.D0 / rkFac
895        recip_horiVertRatio = 1./horiVertRatio        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-- OBCS  C
900        IF (openBoundaries) THEN        CLOSE(iUnit)
        READ(UNIT=iUnit,NML=PARM06)  
        DO J=1,Ny  
         if (OB_Ieast(J).lt.0) OB_Ieast(J)=OB_Ieast(J)+Nx+1  
        ENDDO  
        DO I=1,Nx  
         if (OB_Jnorth(I).lt.0) OB_Jnorth(I)=OB_Jnorth(I)+Ny+1  
        ENDDO  
        write(0,*) 'OB Jn =',OB_Jnorth  
        write(0,*) 'OB Js =',OB_Jsouth  
        write(0,*) 'OB Ie =',OB_Ieast  
        write(0,*) 'OB Iw =',OB_Iwest  
       ENDIF  
901    
902  #ifdef ALLOW_KPP  C--   Check whether any retired parameters were found.
903  C--   KPP vertical mixing parameters  C--   Stop if they were
904        usingKPPmixing      = .FALSE.        IF ( nRetired .GT. 0 ) THEN    
       KPPmixingMaps       = .FALSE.  
       KPPwriteState       = .FALSE.  
       minKPPghat    = 0.  
       maxKPPghat    = 1.e10  
       minKPPviscAz  = 0.  
       DO K =1,Nr  
        maxKPPviscAz(K) = 100.  
       ENDDO  
       minKPPdiffKzT = 0.  
       maxKPPdiffKzT = 100.  
       minKPPdiffKzS = 0.  
       maxKPPdiffKzS = 100.  
   
 c-----------------------------------------------------------------------  
 c define some non-dimensional constants and  
 c the vertical mixing coefficients in m-k-s units  
 c-----------------------------------------------------------------------  
   
       epsln   = 1.e-20      !  
       epsilon = 0.1         ! nondimensional extent of surface layer  
       vonk    = 0.40        ! von Karmans constant  
       conc1   = 5.0         ! scalar coefficients  
       conam   = 1.257       ! .  
       concm   = 8.380       ! .  
       conc2   = 16.0        ! .  
       zetam   = -0.2        ! .  
       conas   = -28.86      ! .  
       concs   = 98.96       ! .  
       conc3   = 16.0        ! .  
       zetas   = -1.0        ! .  
   
 c     parameters for subroutine "bldepth"  
   
       Ricr    = 0.30        ! critical bulk Richardson Number  
       cekman  = 0.7         ! coefficient for ekman depth  
       cmonob  = 1.0         ! coefficient for Monin-Obukhov depth  
       concv   = 1.8         ! ratio of interior buoyancy frequency to  
                             ! buoyancy frequency at entrainment depth  
       hbf     = 1.0         ! fraction of bounadry layer depth to  
                             ! which absorbed solar radiation  
                             ! contributes to surface buoyancy forcing  
       Vtc     = concv * sqrt(0.2/concs/epsilon) / vonk**2 / Ricr  
                             ! non-dimensional coefficient for velocity  
                             ! scale of turbulant velocity shear  
   
 c     parameters and common arrays for subroutines "kmixinit" and "wscale"  
   
       zmin    = -4.e-7      ! minimum limit for zehat in table (m3/s3)  
       zmax    = 0.0         ! maximum limit for zehat in table (m3/s3)  
       umin    = 0.0         ! minimum limit for ustar in table (m/s)  
       umax    = .04         ! maximum limit for ustar in table (m/s)  
   
 c     parameters for subroutine "Ri_iwmix"  
   
       num_v_smooth_Ri = 1   ! number of times Ri is vertically smoothed  
       num_v_smooth_BV = 1   ! number of times BV is vertically smoothed  
       num_z_smooth_sh = 0   ! number of times shear is zonally smoothed  
       num_m_smooth_sh = 0   ! number of times shear is meridionally smoothed  
       Riinfty = 0.7         ! local Richardson Number limit for shear instability  
       BVSQcon = -0.2e-4     ! Brunt-Vaisala squared                   (1/s^2)  
       difm0   = 0.005       ! max visc due to shear instability       (m^2/s)  
       difs0   = 0.005       ! max tracer diffusivity ..               (m^2/s)  
       dift0   = 0.005       ! max heat diffusivity ..                 (m^2/s)  
       difmiw  = 0.001       ! viscosity from background internal wave (m^2/s)  
       difsiw  = 0.00003     ! tracer diffusivity ..                   (m^2/s)  
       diftiw  = 0.00003     ! heat diffusivity ..                     (m^2/s)  
       difmcon = 0.1         ! viscosity due to convective instability (m^2/s)  
       difscon = 0.1         ! tracer diffusivity ..                   (m^2/s)  
       diftcon = 0.1         ! heat diffusivity ..                     (m^2/s)  
   
 c     parameters for subroutine "blmix"  
   
       cstar   = 10.       ! proportionality coefficient for nonlocal transport  
   
                           ! non-dimensional coefficient for counter-gradient term  
       cg      = cstar * vonk * (concs * vonk * epsilon)**(1./3.)  
   
       READ(UNIT=iUnit,NML=PARM07,IOSTAT=errIO,err=13)  
       IF ( errIO .GE. 0 ) GOTO 14  
    13 CONTINUE  
905         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
906       &  'S/R INI_PARMS'       &  'Error reading '
        CALL PRINT_ERROR( msgBuf , 1)  
        WRITE(msgBuf,'(A)')  
      &  'Error reading numerical model '  
907         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
908         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
909       &  'parameter file "data"'       &  'parameter file "data"'
910         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
911         WRITE(msgBuf,'(A)')         WRITE(msgBuf,'(A)')
912       &  'Problem in namelist PARM07'       &  'some out of date parameters were found in the namelist'
913         CALL PRINT_ERROR( msgBuf , 1)         CALL PRINT_ERROR( msgBuf , 1)
        CALL MODELDATA_EXAMPLE( myThid )  
914         STOP 'ABNORMAL END: S/R INI_PARMS'         STOP 'ABNORMAL END: S/R INI_PARMS'
915     14 CONTINUE        ENDIF
 #endif  
   
 C  
       CLOSE(iUnit)  
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.47  
changed lines
  Added in v.1.83
  ViewVC Help
Powered by ViewVC 1.1.22