/[MITgcm]/MITgcm/model/inc/PARAMS.h
ViewVC logotype

Diff of /MITgcm/model/inc/PARAMS.h

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

revision 1.39 by adcroft, Tue Jun 29 18:33:21 1999 UTC revision 1.52 by adcroft, Tue May 29 14:01:36 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3  C  C
4  C     /==========================================================\  C     /==========================================================\
5  C     | PARAMS.h                                                 |  C     | PARAMS.h                                                 |
# Line 56  C                       in relaxation te Line 57  C                       in relaxation te
57  C     saltClimFile    - File containing salt climataology used  C     saltClimFile    - File containing salt climataology used
58  C                       in relaxation term -lambda(salt-salt*)  C                       in relaxation term -lambda(salt-salt*)
59  C     surfQfile       - File containing surface heat flux  C     surfQfile       - File containing surface heat flux
60    C     surfQswfile     - File containing surface shortwave radiation
61    C     dQdTfile        - File containing thermal relaxation coefficient
62  C     EmPmRfile       - File containing surface fresh water flux  C     EmPmRfile       - File containing surface fresh water flux
63  C     buoyancyRelation - Flag used to indicate which relation to use to  C     buoyancyRelation - Flag used to indicate which relation to use to
64  C                        get buoyancy.  C                        get buoyancy.
# Line 63  C                        get buoyancy. Line 66  C                        get buoyancy.
66       &                bathyFile, hydrogThetaFile, hydrogSaltFile,       &                bathyFile, hydrogThetaFile, hydrogSaltFile,
67       &                zonalWindFile, meridWindFile, thetaClimFile,       &                zonalWindFile, meridWindFile, thetaClimFile,
68       &                saltClimFile, buoyancyRelation,       &                saltClimFile, buoyancyRelation,
69       &                EmPmRfile,surfQfile       &                EmPmRfile, surfQfile, surfQswfile,
70         &                uVelInitFile, vVelInitFile, pSurfInitFile,
71         &                dQdTfile
72        CHARACTER*(5) checkPtSuff(maxNoChkptLev)        CHARACTER*(5) checkPtSuff(maxNoChkptLev)
73        CHARACTER*(MAX_LEN_FNAM) bathyFile        CHARACTER*(MAX_LEN_FNAM) bathyFile
74        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
# Line 73  C                        get buoyancy. Line 78  C                        get buoyancy.
78        CHARACTER*(MAX_LEN_FNAM) thetaClimFile        CHARACTER*(MAX_LEN_FNAM) thetaClimFile
79        CHARACTER*(MAX_LEN_FNAM) saltClimFile        CHARACTER*(MAX_LEN_FNAM) saltClimFile
80        CHARACTER*(MAX_LEN_FNAM) surfQfile        CHARACTER*(MAX_LEN_FNAM) surfQfile
81          CHARACTER*(MAX_LEN_FNAM) surfQswfile
82        CHARACTER*(MAX_LEN_FNAM) EmPmRfile        CHARACTER*(MAX_LEN_FNAM) EmPmRfile
83        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
84          CHARACTER*(MAX_LEN_FNAM) uVelInitFile
85          CHARACTER*(MAX_LEN_FNAM) vVelInitFile
86          CHARACTER*(MAX_LEN_FNAM) pSurfInitFile
87          CHARACTER*(MAX_LEN_FNAM) dQdTfile
88    
89  C--   COMMON /PARM_I/ Integer valued parameters used by the model.  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
90  C     cg2dMaxIters        - Maximum number of iterations in the  C     cg2dMaxIters        - Maximum number of iterations in the
# Line 93  C     writeStatePrec      - Precision us Line 103  C     writeStatePrec      - Precision us
103  C     writeBinaryPrec     - Precision used for writing binary files  C     writeBinaryPrec     - Precision used for writing binary files
104  C     readBinaryPrec      - Precision used for reading binary files  C     readBinaryPrec      - Precision used for reading binary files
105  C     nCheckLev           - Holds current checkpoint level  C     nCheckLev           - Holds current checkpoint level
 C     nShap               - "Order" of shapiro filter to apply to  
 C                           model prognositic fields.  
 C                           nShap == 1 => del2  
 C                           nShap == 2 => del4  
 C                           nShap == 3 => del6  
 C                           etc...  
106    
107        COMMON /PARM_I/        COMMON /PARM_I/
108       &        cg2dMaxIters,       &        cg2dMaxIters,
109       &        cg2dChkResFreq,       &        cg2dChkResFreq,
110       &        cg3dMaxIters,       &        cg3dMaxIters,
111       &        cg3dChkResFreq,       &        cg3dChkResFreq,
112       &        nIter0, nTimeSteps,       &        nIter0, nTimeSteps, nEndIter,
113       &        numStepsPerPickup,       &        numStepsPerPickup,
114       &        writeStatePrec, nCheckLev,       &        writeStatePrec, nCheckLev,
115       &        writeBinaryPrec, readBinaryPrec,       &        writeBinaryPrec, readBinaryPrec,
116       &        nShap       &        zonal_filt_sinpow, zonal_filt_cospow
117        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
118        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
119        INTEGER cg3dMaxIters        INTEGER cg3dMaxIters
120        INTEGER cg3dChkResFreq        INTEGER cg3dChkResFreq
121        INTEGER nIter0        INTEGER nIter0
122        INTEGER nTimeSteps        INTEGER nTimeSteps
123          INTEGER nEndIter
124        INTEGER numStepsPerPickup        INTEGER numStepsPerPickup
125        INTEGER writeStatePrec        INTEGER writeStatePrec
126        INTEGER writeBinaryPrec        INTEGER writeBinaryPrec
127        INTEGER readBinaryPrec        INTEGER readBinaryPrec
128        INTEGER nCheckLev        INTEGER nCheckLev
129        INTEGER nShap        INTEGER zonal_filt_sinpow
130          INTEGER zonal_filt_cospow
131    
132  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
133  C     usingCartesianGrid - If TRUE grid generation will be in a cartesian  C     usingCartesianGrid - If TRUE grid generation will be in a cartesian
# Line 130  C     usingSphericalPolarGrid - If TRUE Line 136  C     usingSphericalPolarGrid - If TRUE
136  C                               spherical polar frame.  C                               spherical polar frame.
137  C     no_slip_sides - Impose "no-slip" at lateral boundaries.  C     no_slip_sides - Impose "no-slip" at lateral boundaries.
138  C     no_slip_bottom- Impose "no-slip" at bottom boundary.  C     no_slip_bottom- Impose "no-slip" at bottom boundary.
139    C     staggerTimeStep - enable a Stagger time stepping T,S Rho then U,V
140  C     momViscosity  - Flag which turns momentum friction terms on and off.  C     momViscosity  - Flag which turns momentum friction terms on and off.
141  C     momAdvection  - Flag which turns advection of momentum on and off.  C     momAdvection  - Flag which turns advection of momentum on and off.
142  C     momForcing    - Flag which turns external forcing of momentum on  C     momForcing    - Flag which turns external forcing of momentum on
# Line 151  C     saltAdvection - Flag which turns a Line 158  C     saltAdvection - Flag which turns a
158  C                     and off.  C                     and off.
159  C     saltForcing   - Flag which turns external forcing of salinit on  C     saltForcing   - Flag which turns external forcing of salinit on
160  C                     and off.  C                     and off.
 C     implicitFreeSurface - Set to true to use implcit free surface  
161  C     rigidLid            - Set to true to use rigid lid  C     rigidLid            - Set to true to use rigid lid
162    C     implicitFreeSurface - Set to true to use implcit free surface
163    C     exactConserv   -   Set to true to conserve exactly the total Volume
164    C     nonlinFreeSurf -   Set to true to use non-linear free surface
165  C     momStepping   - Turns momentum equation time-stepping off  C     momStepping   - Turns momentum equation time-stepping off
166  C     tempStepping  - Turns temperature equation time-stepping off  C     tempStepping  - Turns temperature equation time-stepping off
167  C     saltStepping  - Turns salinity equation time-stepping off  C     saltStepping  - Turns salinity equation time-stepping off
# Line 160  C     useConstantF  - Coriolis parameter Line 169  C     useConstantF  - Coriolis parameter
169  C     useBetaPlaneF - Coriolis parameter set to f0 + beta.y  C     useBetaPlaneF - Coriolis parameter set to f0 + beta.y
170  C     useSphereF    - Coriolis parameter set to 2.omega.sin(phi)  C     useSphereF    - Coriolis parameter set to 2.omega.sin(phi)
171  C     implicitDiffusion - Turns implicit vertical diffusion on  C     implicitDiffusion - Turns implicit vertical diffusion on
172    C     implicitViscosity - Turns implicit vertical viscosity on
173  C     doThetaClimRelax - Set true if relaxation to temperature  C     doThetaClimRelax - Set true if relaxation to temperature
174  C                        climatology is required.  C                        climatology is required.
175  C     doSaltClimRelax  - Set true if relaxation to salinity  C     doSaltClimRelax  - Set true if relaxation to salinity
# Line 169  C     usingPCoords     - Set to indicate Line 179  C     usingPCoords     - Set to indicate
179  C                        coords.  C                        coords.
180  C     usingZCoords     - Set to indicate that we are working in height  C     usingZCoords     - Set to indicate that we are working in height
181  C                        coords.  C                        coords.
 C     openBoundaries - Using open-boundaries  
182  C     nonHydrostatic - Using non-hydrostatic terms  C     nonHydrostatic - Using non-hydrostatic terms
183  C     globalFiles    - Selects between "global" and "tiled" files  C     globalFiles    - Selects between "global" and "tiled" files
184  C     allowFreezing  - Allows water to freeze and form ice  C     allowFreezing  - Allows water to freeze and form ice
185    C     groundAtK1  - put the surface(k=1) at the Lower Boundary (=ground)
186        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
187         & usingCurvilinearGrid,
188       & no_slip_sides,no_slip_bottom,       & no_slip_sides,no_slip_bottom,
189         & staggerTimeStep,
190       & momViscosity, momAdvection, momForcing, useCoriolis,       & momViscosity, momAdvection, momForcing, useCoriolis,
191       & momPressureForcing,tempDiffusion, tempAdvection, tempForcing,       & momPressureForcing,tempDiffusion, tempAdvection, tempForcing,
192       & saltDiffusion, saltAdvection, saltForcing,       & saltDiffusion, saltAdvection, saltForcing,
193       & implicitFreeSurface, rigidLid,       & rigidLid, implicitFreeSurface, exactConserv, nonlinFreeSurf,
194       & momStepping, tempStepping, saltStepping,       & momStepping, tempStepping, saltStepping,
195       & metricTerms, usingSphericalPolarMTerms,       & metricTerms, usingSphericalPolarMTerms,
196       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
197       & implicitDiffusion, doThetaClimRelax, doSaltClimRelax,       & implicitDiffusion, implicitViscosity,
198         & doThetaClimRelax, doSaltClimRelax,
199       & periodicExternalForcing, usingPCoords, usingZCoords,       & periodicExternalForcing, usingPCoords, usingZCoords,
200       & openBoundaries, nonHydrostatic, globalFiles,       & nonHydrostatic, globalFiles,
201       & allowFreezing       & allowFreezing, groundAtK1,
202         & usePickupBeforeC35
203        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
204        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid
205          LOGICAL usingCurvilinearGrid
206        LOGICAL usingSphericalPolarMTerms        LOGICAL usingSphericalPolarMTerms
207        LOGICAL no_slip_sides        LOGICAL no_slip_sides
208        LOGICAL no_slip_bottom        LOGICAL no_slip_bottom
209          LOGICAL staggerTimeStep
210        LOGICAL momViscosity        LOGICAL momViscosity
211        LOGICAL momAdvection        LOGICAL momAdvection
212        LOGICAL momForcing        LOGICAL momForcing
# Line 202  C     allowFreezing  - Allows water to f Line 218  C     allowFreezing  - Allows water to f
218        LOGICAL saltDiffusion        LOGICAL saltDiffusion
219        LOGICAL saltAdvection        LOGICAL saltAdvection
220        LOGICAL saltForcing        LOGICAL saltForcing
       LOGICAL implicitFreeSurface  
221        LOGICAL rigidLid        LOGICAL rigidLid
222          LOGICAL implicitFreeSurface
223          LOGICAL exactConserv
224          LOGICAL nonlinFreeSurf
225        LOGICAL momStepping        LOGICAL momStepping
226        LOGICAL tempStepping        LOGICAL tempStepping
227        LOGICAL saltStepping        LOGICAL saltStepping
# Line 212  C     allowFreezing  - Allows water to f Line 230  C     allowFreezing  - Allows water to f
230        LOGICAL useBetaPlaneF        LOGICAL useBetaPlaneF
231        LOGICAL useSphereF        LOGICAL useSphereF
232        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
233          LOGICAL implicitViscosity
234        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
235        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
236        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
237        LOGICAL usingPCoords        LOGICAL usingPCoords
238        LOGICAL usingZCoords        LOGICAL usingZCoords
       LOGICAL openBoundaries  
239        LOGICAL nonHydrostatic        LOGICAL nonHydrostatic
240        LOGICAL globalFiles        LOGICAL globalFiles
241        LOGICAL allowFreezing        LOGICAL allowFreezing
242          LOGICAL groundAtK1
243          LOGICAL usePickupBeforeC35
244    
245  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
246  C     cg2dTargetResidual  C     gg2dTargetResidual
247  C               - Target residual for cg2d solver.  C          - Target residual for cg2d solver; no unit (RHS normalisation)
248    C     cg2dTargetResWunit
249    C          - Target residual for cg2d solver; W unit (No RHS normalisation)
250  C     cg3dTargetResidual  C     cg3dTargetResidual
251  C               - Target residual for cg3d solver.  C               - Target residual for cg3d solver.
252  C     cg2dpcOffDFac - Averaging weight for preconditioner off-diagonal.  C     cg2dpcOffDFac - Averaging weight for preconditioner off-diagonal.
# Line 304  C     deltaTtracer - Timestep for tracer Line 326  C     deltaTtracer - Timestep for tracer
326  C     freesurfFac  - Parameter to turn implicit free surface term on or off  C     freesurfFac  - Parameter to turn implicit free surface term on or off
327  C                    freesurfac = 1. uses implicit free surface  C                    freesurfac = 1. uses implicit free surface
328  C                    freesurfac = 0. uses rigid lid  C                    freesurfac = 0. uses rigid lid
329    C     implicSurfPress - parameter of the Crank-Nickelson time stepping :
330    C                     Implicit part of Surface Pressure Gradient ( 0-1 )
331    C     implicDiv2Dflow - parameter of the Crank-Nickelson time stepping :
332    C                     Implicit part of barotropic flow Divergence ( 0-1 )
333  C     hFacMin      - Minimum fraction size of a cell (affects hFacC etc...)  C     hFacMin      - Minimum fraction size of a cell (affects hFacC etc...)
334  C     hFacMinDz    - Minimum dimesional size of a cell (affects hFacC etc..., m)  C     hFacMinDz    - Minimum dimesional size of a cell (affects hFacC etc..., m)
335  C     hFacMinDp    - Minimum dimesional size of a cell (affects hFacC etc..., Pa)  C     hFacMinDp    - Minimum dimesional size of a cell (affects hFacC etc..., Pa)
336  C     hFacMinDr    - Minimum dimesional size of a cell (affects hFacC etc..., units of r)  C     hFacMinDr    - Minimum dimesional size of a cell (affects hFacC etc..., units of r)
337    C     hFacInf      - Threshold (inf and sup) for fraction size of surface cell
338    C     hFacSup        that control vanishing and creating levels
339  C     tauCD        - CD scheme coupling timescale ( 1/s )  C     tauCD        - CD scheme coupling timescale ( 1/s )
340  C     rCD          - CD scheme normalised coupling parameter ( 0-1 )  C     rCD          - CD scheme normalised coupling parameter ( 0-1 )
 C     GMmaxslope    - max. slope allowed in GM/Redi tensor  
 C     GMlength      - Length to use in Visbeck et al. formula for K (m)  
 C     GMalpha       - alpha to use in Visbeck et al. formula for K  
 C     GMdepth       - Depth over which to integrate Richardson # (Visbeck et al.)  
 C     GMkbackground - background value of GM/Redi diffusion coefficient ( m^2/s )  
 C     GMmaxval      - max. value of KapGM allowed in GM/Redi scheme ( m^2/s )  
341  C     startTime     - Starting time for this integration ( s ).  C     startTime     - Starting time for this integration ( s ).
342  C     endTime       - Ending time for this integration ( s ).  C     endTime       - Ending time for this integration ( s ).
343  C     chkPtFreq     - Frequency of rolling check pointing ( s ).  C     chkPtFreq     - Frequency of rolling check pointing ( s ).
# Line 346  C     recip_horiVertRatio  ( 1 if horiz Line 368  C     recip_horiVertRatio  ( 1 if horiz
368  C                          ( g*rho if horiz in m and vertical in Pa ).  C                          ( g*rho if horiz in m and vertical in Pa ).
369  C     latFFTFiltLo       - Low latitude for FFT filtering of latitude  C     latFFTFiltLo       - Low latitude for FFT filtering of latitude
370  C                          circles ( see filter*.F )  C                          circles ( see filter*.F )
371        COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac,  C     Ro_SeaLevel        - standard position of Sea-Level in "R" coordinate, used as
372       & cg3dTargetResidual,  C                          starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
373    C     bottomDragLinear   - Drag coefficient built in to core dynamics
374    C         "     Quadratic  ( linear: 1/s, quadratic: 1/m )
375          COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
376         & cg2dpcOffDFac, cg3dTargetResidual,
377       & delP, delZ, delR, delX, delY,       & delP, delZ, delR, delX, delY,
378       & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime,       & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime,
379       & phiMin, thetaMin, rSphere, recip_RSphere, f0, fCori, beta,       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
380         & fCori, fCoriG,
381       & viscAh,  viscAz,  viscA4,  viscAr,       & viscAh,  viscAz,  viscA4,  viscAr,
382       & diffKhT, diffKzT, diffK4T, diffKrT,       & diffKhT, diffKzT, diffK4T, diffKrT,
383       & diffKhS, diffKzS, diffK4S, diffKrS,       & diffKhS, diffKzS, diffK4S, diffKrS,
384       & delT, tauCD, rCD, freeSurfFac, hFacMin, hFacMinDz,       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
385       & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,GMmaxval,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
386       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
387       & recip_rhoConst, rhoConst, tRef, sRef,       & recip_rhoConst, rhoConst, tRef, sRef,
388       & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,       & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,
# Line 366  C                          circles ( see Line 393  C                          circles ( see
393       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
394       & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,       & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,
395       & theta_S, specVol_S, horiVertRatio, recip_horiVertRatio,       & theta_S, specVol_S, horiVertRatio, recip_horiVertRatio,
396       & latFFTFiltLo       & latFFTFiltLo, ivdc_kappa, Ro_SeaLevel, zonal_filt_lat,
397         & bottomDragLinear,bottomDragQuadratic
398    
399        _RL cg2dTargetResidual        _RL cg2dTargetResidual
400          _RL cg2dTargetResWunit
401        _RL cg3dTargetResidual        _RL cg3dTargetResidual
402        _RL cg2dpcOffDFac        _RL cg2dpcOffDFac
403        _RL delZ(Nr)        _RL delZ(Nr)
# Line 387  C                          circles ( see Line 416  C                          circles ( see
416        _RL recip_RSphere        _RL recip_RSphere
417        _RL f0        _RL f0
418        _RL freeSurfFac        _RL freeSurfFac
419          _RL implicSurfPress
420          _RL implicDiv2Dflow
421        _RL hFacMin        _RL hFacMin
422        _RL hFacMinDz        _RL hFacMinDz
423        _RL hFacMinDp        _RL hFacMinDp
424        _RL hFacMinDr        _RL hFacMinDr
425          _RL hFacInf
426          _RL hFacSup
427        _RL beta        _RL beta
428        _RL viscAh        _RL viscAh
429        _RL viscAz        _RL viscAz
# Line 410  C                          circles ( see Line 443  C                          circles ( see
443        _RL delt        _RL delt
444        _RL tauCD        _RL tauCD
445        _RL rCD        _RL rCD
       _RL GMmaxslope  
       _RL GMlength  
       _RL GMalpha  
       _RL GMdepth  
       _RL GMkbackground  
       _RL GMmaxval  
446        _RL gravity        _RL gravity
447        _RL recip_gravity        _RL recip_gravity
448        _RL gBaro        _RL gBaro
# Line 427  C                          circles ( see Line 454  C                          circles ( see
454        _RL tRef(Nr)        _RL tRef(Nr)
455        _RL theta_S(Nr)        _RL theta_S(Nr)
456        _RL sRef(Nr)        _RL sRef(Nr)
457        _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
458          _RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
459        _RL startTime        _RL startTime
460        _RL endTime        _RL endTime
461        _RL chkPtFreq        _RL chkPtFreq
# Line 452  C                          circles ( see Line 480  C                          circles ( see
480        _RL horiVertRatio        _RL horiVertRatio
481        _RL recip_horiVertRatio        _RL recip_horiVertRatio
482        _RL latFFTFiltLo        _RL latFFTFiltLo
483          _RL ivdc_kappa
484          _RL Ro_SeaLevel
485          _RL zonal_filt_lat
486          _RL bottomDragLinear
487          _RL bottomDragQuadratic
488    
489        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,
490       &                Lamba_theta       &                Lamba_theta
# Line 470  C     sBeta     - Linear EOS haline cont Line 503  C     sBeta     - Linear EOS haline cont
503        _RL sBeta        _RL sBeta
504        character*(6) eosType        character*(6) eosType
505    
506  C These are input arrays (of integers) that contain the *absolute*  C Logical flags for selecting packages
507  C computational index of an open-boundary (OB) point.        LOGICAL useKPP
508  C A zero (0) element means there is no corresponding OB in that column/row.        LOGICAL useGMRedi
509  C The computational coordinate refers to "tracer" cells.        LOGICAL useOBCS
510  C For a northern/southern OB, the OB V point is to the south/north.        LOGICAL useAIM
511  C For an eastern/western OB, the OB U point is to the west/east.        LOGICAL useECCO
512  C eg.        LOGICAL useSHAP_FILT
513  C     OB_Jnorth(3)=34  means that:        COMMON /PARM_PACKAGES/
514  C          T( 3 ,34) is a an OB point       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO, useSHAP_FILT
515  C          U(3:4,34) is a an OB point  
 C          V( 4 ,34) is a an OB point  
 C while  
 C     OB_Jsouth(3)=1  means that:  
 C          T( 3 ,1) is a an OB point  
 C          U(3:4,1) is a an OB point  
 C          V( 4 ,2) is a an OB point  
 C  
 C For convenience, negative values for Jnorth/Ieast refer to  
 C points relative to the Northern/Eastern edges of the model  
 C eg. OB_Jnorth(3)=-1  means that the point (3,Ny-1) is a northern O-B.  
 C  
       COMMON /PARM_OB/  
      & OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest  
       INTEGER OB_Jnorth(Nx)  
       INTEGER OB_Jsouth(Nx)  
       INTEGER OB_Ieast(Ny)  
       INTEGER OB_Iwest(Ny)  
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.52
  ViewVC Help
Powered by ViewVC 1.1.22