/[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.209 by dimitri, Wed Nov 28 09:26:16 2007 UTC revision 1.217 by gforget, Fri May 30 02:46:19 2008 UTC
# Line 84  C     surfQnetFile    :: File containing Line 84  C     surfQnetFile    :: File containing
84  C     surfQswFile     :: File containing surface shortwave radiation  C     surfQswFile     :: File containing surface shortwave radiation
85  C     dQdTfile        :: File containing thermal relaxation coefficient  C     dQdTfile        :: File containing thermal relaxation coefficient
86  C     EmPmRfile       :: File containing surface fresh water flux  C     EmPmRfile       :: File containing surface fresh water flux
87    C           NOTE: for backward compatibility EmPmRfile is specified in
88    C                 m/s when using external_fields_load.F.  It is converted
89    C                 to kg/m2/s by multiplying by rhoConstFresh.
90  C     saltFluxFile    :: File containing surface salt flux  C     saltFluxFile    :: File containing surface salt flux
91  C     pLoadFile       :: File containing pressure loading  C     pLoadFile       :: File containing pressure loading
92  C     eddyTauxFile    :: File containing zonal Eddy stress data  C     eddyPsiXFile    :: File containing zonal Eddy streamfunction data
93  C     eddyTauyFile    :: File containing meridional Eddy stress data  C     eddyPsiYFile    :: File containing meridional Eddy streamfunction data
94  C     the_run_name    :: string identifying the name of the model "run"  C     the_run_name    :: string identifying the name of the model "run"
95        COMMON /PARM_C/        COMMON /PARM_C/
96       &                buoyancyRelation, eosType,       &                buoyancyRelation, eosType,
# Line 104  C     the_run_name    :: string identify Line 107  C     the_run_name    :: string identify
107       &                lambdaThetaFile, lambdaSaltFile,       &                lambdaThetaFile, lambdaSaltFile,
108       &                uVelInitFile, vVelInitFile, pSurfInitFile,       &                uVelInitFile, vVelInitFile, pSurfInitFile,
109       &                dQdTfile, ploadFile,       &                dQdTfile, ploadFile,
110       &                eddyTauxFile, eddyTauyFile,       &                eddyPsiXFile, eddyPsiYFile,
111       &                the_run_name       &                the_run_name
112        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
113        CHARACTER*(6)  eosType        CHARACTER*(6)  eosType
# Line 135  C     the_run_name    :: string identify Line 138  C     the_run_name    :: string identify
138        CHARACTER*(MAX_LEN_FNAM) pSurfInitFile        CHARACTER*(MAX_LEN_FNAM) pSurfInitFile
139        CHARACTER*(MAX_LEN_FNAM) dQdTfile        CHARACTER*(MAX_LEN_FNAM) dQdTfile
140        CHARACTER*(MAX_LEN_FNAM) ploadFile        CHARACTER*(MAX_LEN_FNAM) ploadFile
141        CHARACTER*(MAX_LEN_FNAM) eddyTauxFile        CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile
142        CHARACTER*(MAX_LEN_FNAM) eddyTauyFile        CHARACTER*(MAX_LEN_FNAM) eddyPsiYFile
143        CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile        CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile
144        CHARACTER*(MAX_LEN_FNAM) lambdaSaltFile        CHARACTER*(MAX_LEN_FNAM) lambdaSaltFile
145        CHARACTER*(MAX_LEN_PREC/2) the_run_name        CHARACTER*(MAX_LEN_PREC/2) the_run_name
# Line 170  C     tempVertAdvScheme   :: Temp. Vert. Line 173  C     tempVertAdvScheme   :: Temp. Vert.
173  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector
174  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector
175  C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)  C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)
176    C     selectVortScheme    :: Scheme selector for Vorticity term (Vector Inv.)
177    C     monitorSelect       :: select group of variables to monitor
178    C                            =1 : dynvars ; =2 : + vort ; =3 : + surface
179  C     debugLevel          :: debug level selector: higher -> more writing  C     debugLevel          :: debug level selector: higher -> more writing
180    
181        COMMON /PARM_I/        COMMON /PARM_I/
# Line 184  C     debugLevel          :: debug level Line 190  C     debugLevel          :: debug level
190       &        momForcingOutAB, tracForcingOutAB,       &        momForcingOutAB, tracForcingOutAB,
191       &        tempAdvScheme, tempVertAdvScheme,       &        tempAdvScheme, tempVertAdvScheme,
192       &        saltAdvScheme, saltVertAdvScheme,       &        saltAdvScheme, saltVertAdvScheme,
193       &        selectKEscheme,       &        selectKEscheme, selectVortScheme,
194       &        debugLevel       &        monitorSelect, debugLevel
195        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
196        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
197        INTEGER cg2dPreCondFreq        INTEGER cg2dPreCondFreq
# Line 203  C     debugLevel          :: debug level Line 209  C     debugLevel          :: debug level
209        INTEGER tempAdvScheme, tempVertAdvScheme        INTEGER tempAdvScheme, tempVertAdvScheme
210        INTEGER saltAdvScheme, saltVertAdvScheme        INTEGER saltAdvScheme, saltVertAdvScheme
211        INTEGER selectKEscheme        INTEGER selectKEscheme
212          INTEGER selectVortScheme
213          INTEGER monitorSelect
214        INTEGER debugLevel        INTEGER debugLevel
215    
216  C  C
# Line 214  C Line 222  C
222        PARAMETER(debLevB=2)        PARAMETER(debLevB=2)
223    
224  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
225    C- Coordinate + Grid params:
226    C     fluidIsAir       :: Set to indicate that the fluid major constituent
227    C                        is air
228    C     fluidIsWater     :: Set to indicate that the fluid major constituent
229    C                        is water
230    C     usingPCoords     :: Set to indicate that we are working in a pressure
231    C                        type coordinate (p or p*).
232    C     usingZCoords     :: Set to indicate that we are working in a height
233    C                        type coordinate (z or z*)
234    C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
235    C                         this requires specific code for restart & exchange
236  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
237  C                          coordinate frame.  C                          coordinate frame.
238  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a
239  C                               spherical polar frame.  C                               spherical polar frame.
240  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical  C     rotateGrid      :: rotate grid coordinates to geographical coordinates
241    C                        according to Euler angles phiEuler, thetaEuler, psiEuler
242  C     usingCurvilinearGrid :: If TRUE, use a curvilinear grid (to be provided)  C     usingCurvilinearGrid :: If TRUE, use a curvilinear grid (to be provided)
243    C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
244  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)
245  C     no_slip_sides :: Impose "no-slip" at lateral boundaries.  C     setInterFDr    :: set Interface depth (put cell-Center at the middle)
246    C     setCenterDr    :: set cell-Center depth (put Interface at the middle)
247    C- Momentum params:
248    C     no_slip_sides  :: Impose "no-slip" at lateral boundaries.
249  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.
250    C     useFullLeith   :: Set to true to use full Leith viscosity(may be unstable
251    C                       on irregular grids)
252    C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
253    C     useAreaViscLength :: Set to true to use old scaling for viscous lengths,
254    C                          e.g., L2=Raz.  May be preferable for cube sphere.
255  C     momViscosity  :: Flag which turns momentum friction terms on and off.  C     momViscosity  :: Flag which turns momentum friction terms on and off.
256  C     momAdvection  :: Flag which turns advection of momentum on and off.  C     momAdvection  :: Flag which turns advection of momentum on and off.
257  C     momForcing    :: Flag which turns external forcing of momentum on  C     momForcing    :: Flag which turns external forcing of momentum on
258  C                     and off.  C                      and off.
259  C     momPressureForcing :: Flag which turns pressure term in momentum equation  C     momPressureForcing :: Flag which turns pressure term in momentum equation
260  C                          on and off.  C                          on and off.
261  C     metricTerms   :: Flag which turns metric terms on or off.  C     metricTerms   :: Flag which turns metric terms on or off.
262  C     useNHMTerms   :: If TRUE use non-hydrostatic metric terms.  C     useNHMTerms   :: If TRUE use non-hydrostatic metric terms.
263  C     useCoriolis   :: Flag which turns the coriolis terms on and off.  C     useCoriolis   :: Flag which turns the coriolis terms on and off.
264  C     tempAdvection :: Flag which turns advection of temperature on  C     use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - off
265  C                     and off.  C     useConstantF  :: Coriolis parameter set to f0
266  C     tempForcing   :: Flag which turns external forcing of temperature on  C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y
267  C                     and off.  C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)
268  C     saltAdvection :: Flag which turns advection of salinity on  C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
269  C                     and off.  C     vectorInvariantMomentum :: use Vector-Invariant form (mom_vecinv package)
270  C     saltForcing   :: Flag which turns external forcing of salinity on  C                                (default = F = use mom_fluxform package)
271  C                     and off.  C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart & Ozer 1986)
272    C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term
273    C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
274    C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
275    C     useAbsVorticity :: work with f+zeta in Coriolis terms
276    C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
277    C     momStepping    :: Turns momentum equation time-stepping off
278    C- Temp. & Salt params:
279    C     tempStepping   :: Turns temperature equation time-stepping off
280    C     saltStepping   :: Turns salinity equation time-stepping off
281    C     tempAdvection  :: Flag which turns advection of temperature on and off.
282    C     tempIsActiveTr :: Pot.Temp. is a dynamically active tracer
283    C     tempForcing    :: Flag which turns external forcing of temperature on
284    C                       and off.
285    C     saltAdvection  :: Flag which turns advection of salinity on and off.
286    C     saltIsActiveTr :: Salinity  is a dynamically active tracer
287    C     saltForcing    :: Flag which turns external forcing of salinity on
288    C                       and off.
289  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
290  C                         as a real Fresh Water (=> changes the Sea Level)  C                         as a real Fresh Water (=> changes the Sea Level)
291  C                         if F, converts P+R-E to salt flux (no SL effect)  C                         if F, converts P+R-E to salt flux (no SL effect)
292  C     useFullLeith   :: Set to true to use full Leith viscosity(may be unstable  C- Time-stepping params:
 C                       on irregular grids)  
 C     useAreaViscLength :: Set to true to use old scaling for viscous  
 C              lengths, e.g., L2=Raz.  May be preferable for cube sphere.  
 C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms  
293  C     rigidLid            :: Set to true to use rigid lid  C     rigidLid            :: Set to true to use rigid lid
294  C     implicitFreeSurface :: Set to true to use implicit free surface  C     implicitFreeSurface :: Set to true to use implicit free surface
295  C     exactConserv        :: Set to true to conserve exactly the total Volume  C     exactConserv        :: Set to true to conserve exactly the total Volume
296  C     linFSConserveTr     :: Set to true to correct source/sink of tracer  C     linFSConserveTr     :: Set to true to correct source/sink of tracer
297  C                            at the surface due to Linear Free Surface  C                            at the surface due to Linear Free Surface
298  C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the  C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the
299  C                           linear relation Phi_surf = Bo_surf*eta  C                            linear relation Phi_surf = Bo_surf*eta
300  C     use3Dsolver   :: set to true to use 3-D pressure solver  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
301    C     nonHydrostatic   :: Using non-hydrostatic algorithm
302    C     use3Dsolver      :: set to true to use 3-D pressure solver
303  C     implicitIntGravWave :: treat Internal Gravity Wave implicitly  C     implicitIntGravWave :: treat Internal Gravity Wave implicitly
304  C     staggerTimeStep :: enable a Stagger time stepping T,S Rho then U,V  C     staggerTimeStep   :: enable a Stagger time stepping U,V (& W) then T,S
 C     momStepping   :: Turns momentum equation time-stepping off  
 C     tempStepping  :: Turns temperature equation time-stepping off  
 C     saltStepping  :: Turns salinity equation time-stepping off  
 C     useConstantF  :: Coriolis parameter set to f0  
 C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y  
 C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)  
 C     use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - off  
 C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.  
 C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)  
 C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term  
 C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny  
 C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term  
 C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)  
 C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)  
 C     useAbsVorticity :: work with f+zeta in Coriolis terms  
305  C     implicitDiffusion :: Turns implicit vertical diffusion on  C     implicitDiffusion :: Turns implicit vertical diffusion on
306  C     implicitViscosity :: Turns implicit vertical viscosity on  C     implicitViscosity :: Turns implicit vertical viscosity on
307  C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature  C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature
# Line 284  C     momDissip_In_AB   :: if False, put Line 313  C     momDissip_In_AB   :: if False, put
313  C                          out off Adams-Bashforth time stepping.  C                          out off Adams-Bashforth time stepping.
314  C     doAB_onGtGs       :: if the Adams-Bashforth time stepping is used, always  C     doAB_onGtGs       :: if the Adams-Bashforth time stepping is used, always
315  C                          apply AB on tracer tendencies (rather than on Tracer)  C                          apply AB on tracer tendencies (rather than on Tracer)
316  C     pickupStrictlyMatch :: check and stop if pickup-file do not stricly match  C- Other forcing params -
317  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup  C     balanceEmPmR    :: substract global mean of EmPmR at every time step
318  C     usePickupBeforeC54 :: start from old-pickup files, generated with code from  C     balanceQnet     :: substract global mean of Qnet at every time step
319  C                           before checkpoint-54a, Jul 06, 2004.  C     balancePrintMean:: print substracted global means to STDOUT
320  C     doThetaClimRelax :: Set true if relaxation to temperature  C     doThetaClimRelax :: Set true if relaxation to temperature
321  C                        climatology is required.  C                        climatology is required.
322  C     doSaltClimRelax  :: Set true if relaxation to salinity  C     doSaltClimRelax  :: Set true if relaxation to salinity
323  C                        climatology is required.  C                        climatology is required.
324    C     allowFreezing  :: Allows surface water to freeze and form ice
325    C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
326  C     periodicExternalForcing :: Set true if forcing is time-dependant  C     periodicExternalForcing :: Set true if forcing is time-dependant
327  C     usingPCoords     :: Set to indicate that we are working in a pressure  C- I/O parameters -
 C                        type coordinate (p or p*).  
 C     usingZCoords     :: Set to indicate that we are working in a height  
 C                        type coordinate (z or z*)  
 C     fluidIsAir       :: Set to indicate that the fluid major constituent  
 C                        is air  
 C     fluidIsWater     :: Set to indicate that the fluid major constituent  
 C                        is water  
 C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)  
 C                         this requires specific code for restart & exchange  
 C     setInterFDr    :: set Interface depth (put cell-Center at the middle)  
 C     setCenterDr    :: set cell-Center depth (put Interface at the middle)  
 C     nonHydrostatic :: Using non-hydrostatic terms  
 C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm  
328  C     globalFiles    :: Selects between "global" and "tiled" files  C     globalFiles    :: Selects between "global" and "tiled" files
329  C     useSingleCpuIO :: On SGI platforms, option globalFiles is either  C     useSingleCpuIO :: On SGI platforms, option globalFiles is either
330  C                       slow (f77) or does not work (f90).  When  C                       slow (f77) or does not work (f90).  When
331  C                       useSingleCpuIO is set, mdsio_writefield.F  C                       useSingleCpuIO is set, mdsio_writefield.F
332  C                       outputs from master mpi process only.  C                       outputs from master mpi process only.
333  C     allowFreezing  :: Allows surface water to freeze and form ice  C     pickupStrictlyMatch :: check and stop if pickup-file do not stricly match
334  C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
335    C     usePickupBeforeC54 :: start from old-pickup files, generated with code from
336    C                           before checkpoint-54a, Jul 06, 2004.
337  C     pickup_write_mdsio :: use mdsio to write pickups  C     pickup_write_mdsio :: use mdsio to write pickups
338  C     pickup_read_mdsio  :: use mdsio to read  pickups  C     pickup_read_mdsio  :: use mdsio to read  pickups
339  C     pickup_write_immed :: echo the pickup immediately (for conversion)  C     pickup_write_immed :: echo the pickup immediately (for conversion)
340    C     writePickupAtEnd   :: write pickup at the last timestep
341  C     timeave_mdsio      :: use mdsio for timeave output  C     timeave_mdsio      :: use mdsio for timeave output
342  C     snapshot_mdsio     :: use mdsio for "snapshot" (dumpfreq/diagfreq) output  C     snapshot_mdsio     :: use mdsio for "snapshot" (dumpfreq/diagfreq) output
343  C     monitor_stdio      :: use stdio for monitor output  C     monitor_stdio      :: use stdio for monitor output
344  C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)  C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)
345  C                        & Last iteration, in addition multiple of dumpFreq iter.  C                        & Last iteration, in addition multiple of dumpFreq iter.
 C     balanceEmPmR    :: substract global mean of EmPmR at every time step  
 C     balanceQnet     :: substract global mean of Qnet at every time step  
 C     balancePrintMean:: print substracted global means to STDOUT  
346    
347        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/
348         & fluidIsAir, fluidIsWater,
349         & usingPCoords, usingZCoords, useDynP_inEos_Zc,
350         & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,
351       & usingCurvilinearGrid, usingCylindricalGrid,       & usingCurvilinearGrid, usingCylindricalGrid,
352       & deepAtmosphere, setInterFDr, setCenterDr,       & deepAtmosphere, setInterFDr, setCenterDr,
353       & no_slip_sides,no_slip_bottom,       & no_slip_sides, no_slip_bottom,
354       & momViscosity, momAdvection, momForcing, useCoriolis,       & useFullLeith, useStrainTensionVisc, useAreaViscLength,
355       & momPressureForcing, vectorInvariantMomentum,       & momViscosity, momAdvection, momForcing,
356       & tempAdvection, tempForcing,       & momPressureForcing, metricTerms, useNHMTerms,
357       & saltAdvection, saltForcing,       & useCoriolis, use3dCoriolis,
358         & useConstantF, useBetaPlaneF, useSphereF,
359         & useCDscheme, vectorInvariantMomentum,
360         & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
361         & upwindVorticity, highOrderVorticity,
362         & useAbsVorticity, upwindShear,
363         & momStepping, tempStepping, saltStepping,
364         & tempAdvection, tempIsActiveTr, tempForcing,
365         & saltAdvection, saltIsActiveTr, saltForcing,
366       & useRealFreshWaterFlux,       & useRealFreshWaterFlux,
      & useFullLeith, useStrainTensionVisc,  
      & useAreaViscLength,  
367       & rigidLid, implicitFreeSurface, exactConserv, linFSConserveTr,       & rigidLid, implicitFreeSurface, exactConserv, linFSConserveTr,
368       & uniformLin_PhiSurf,       & uniformLin_PhiSurf,
369         & quasiHydrostatic, nonHydrostatic,
370       & use3Dsolver, implicitIntGravWave, staggerTimeStep,       & use3Dsolver, implicitIntGravWave, staggerTimeStep,
      & momStepping, tempStepping, saltStepping,  
      & metricTerms, useNHMTerms,  
      & useConstantF, useBetaPlaneF, useSphereF, use3dCoriolis,  
      & useCDscheme,  
      & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,  
      & SadournyCoriolis, upwindVorticity, highOrderVorticity,  
      & useAbsVorticity, upwindShear,  
371       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
372       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
373       & multiDimAdvection, useMultiDimAdvec,       & multiDimAdvection, useMultiDimAdvec,
374       & momDissip_In_AB, doAB_onGtGs,       & momDissip_In_AB, doAB_onGtGs,
375       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,       & balanceEmPmR, balanceQnet, balancePrintMean,
376       & periodicExternalForcing,       & doThetaClimRelax, doSaltClimRelax,
      & fluidIsAir, fluidIsWater,  
      & usingPCoords, usingZCoords, useDynP_inEos_Zc,  
      & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,  
377       & allowFreezing, useOldFreezing,       & allowFreezing, useOldFreezing,
378         & periodicExternalForcing,
379         & globalFiles, useSingleCpuIO,
380       & pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2,       & pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2,
381       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
382         & writePickupAtEnd,
383       & timeave_mdsio, snapshot_mdsio, monitor_stdio,       & timeave_mdsio, snapshot_mdsio, monitor_stdio,
384       & outputTypesInclusive, dumpInitAndLast, debugMode,       & outputTypesInclusive, dumpInitAndLast, debugMode,
385       & inAdMode, inAdTrue, inAdFalse, inAdExact,       & inAdMode, inAdTrue, inAdFalse, inAdExact
      & balanceEmPmR, balanceQnet, balancePrintMean  
386    
387          LOGICAL fluidIsAir
388          LOGICAL fluidIsWater
389          LOGICAL usingPCoords
390          LOGICAL usingZCoords
391          LOGICAL useDynP_inEos_Zc
392        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
393        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid, rotateGrid
394        LOGICAL usingCylindricalGrid        LOGICAL usingCylindricalGrid
395        LOGICAL usingCurvilinearGrid        LOGICAL usingCurvilinearGrid
396        LOGICAL deepAtmosphere        LOGICAL deepAtmosphere
# Line 381  C     balancePrintMean:: print substract Line 406  C     balancePrintMean:: print substract
406        LOGICAL useCoriolis        LOGICAL useCoriolis
407        LOGICAL vectorInvariantMomentum        LOGICAL vectorInvariantMomentum
408        LOGICAL tempAdvection        LOGICAL tempAdvection
409          LOGICAL tempIsActiveTr
410        LOGICAL tempForcing        LOGICAL tempForcing
411        LOGICAL saltAdvection        LOGICAL saltAdvection
412          LOGICAL saltIsActiveTr
413        LOGICAL saltForcing        LOGICAL saltForcing
414        LOGICAL useRealFreshWaterFlux        LOGICAL useRealFreshWaterFlux
415        LOGICAL useFullLeith        LOGICAL useFullLeith
# Line 393  C     balancePrintMean:: print substract Line 420  C     balancePrintMean:: print substract
420        LOGICAL exactConserv        LOGICAL exactConserv
421        LOGICAL linFSConserveTr        LOGICAL linFSConserveTr
422        LOGICAL uniformLin_PhiSurf        LOGICAL uniformLin_PhiSurf
423          LOGICAL quasiHydrostatic
424          LOGICAL nonHydrostatic
425        LOGICAL use3Dsolver        LOGICAL use3Dsolver
426        LOGICAL implicitIntGravWave        LOGICAL implicitIntGravWave
427        LOGICAL staggerTimeStep        LOGICAL staggerTimeStep
# Line 408  C     balancePrintMean:: print substract Line 437  C     balancePrintMean:: print substract
437        LOGICAL useEnergyConservingCoriolis        LOGICAL useEnergyConservingCoriolis
438        LOGICAL useJamartWetPoints        LOGICAL useJamartWetPoints
439        LOGICAL useJamartMomAdv        LOGICAL useJamartMomAdv
       LOGICAL SadournyCoriolis  
440        LOGICAL upwindVorticity        LOGICAL upwindVorticity
441        LOGICAL highOrderVorticity        LOGICAL highOrderVorticity
442        LOGICAL useAbsVorticity        LOGICAL useAbsVorticity
# Line 422  C     balancePrintMean:: print substract Line 450  C     balancePrintMean:: print substract
450        LOGICAL useMultiDimAdvec        LOGICAL useMultiDimAdvec
451        LOGICAL momDissip_In_AB        LOGICAL momDissip_In_AB
452        LOGICAL doAB_onGtGs        LOGICAL doAB_onGtGs
453          LOGICAL balanceEmPmR
454          LOGICAL balanceQnet
455          LOGICAL balancePrintMean
456        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
457        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
458        LOGICAL doTr1ClimRelax        LOGICAL allowFreezing
459          LOGICAL useOldFreezing
460        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
       LOGICAL fluidIsAir  
       LOGICAL fluidIsWater  
       LOGICAL usingPCoords  
       LOGICAL usingZCoords  
       LOGICAL useDynP_inEos_Zc  
       LOGICAL nonHydrostatic  
       LOGICAL quasiHydrostatic  
461        LOGICAL globalFiles        LOGICAL globalFiles
462        LOGICAL useSingleCpuIO        LOGICAL useSingleCpuIO
       LOGICAL allowFreezing  
       LOGICAL useOldFreezing  
463        LOGICAL pickupStrictlyMatch        LOGICAL pickupStrictlyMatch
464        LOGICAL usePickupBeforeC54        LOGICAL usePickupBeforeC54
465        LOGICAL startFromPickupAB2        LOGICAL startFromPickupAB2
       LOGICAL dumpInitAndLast  
       LOGICAL debugMode  
466        LOGICAL pickup_read_mdsio, pickup_write_mdsio        LOGICAL pickup_read_mdsio, pickup_write_mdsio
467        LOGICAL pickup_write_immed        LOGICAL pickup_write_immed, writePickupAtEnd
468        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio
469        LOGICAL outputTypesInclusive        LOGICAL outputTypesInclusive
470          LOGICAL dumpInitAndLast
471          LOGICAL debugMode
472        LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact        LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact
473    
       LOGICAL balanceEmPmR  
       LOGICAL balanceQnet  
       LOGICAL balancePrintMean  
   
474  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
475  C     cg2dTargetResidual  C     cg2dTargetResidual
476  C          :: Target residual for cg2d solver; no unit (RHS normalisation)  C          :: Target residual for cg2d solver; no unit (RHS normalisation)
# Line 516  C                 ( Southern edge f for Line 535  C                 ( Southern edge f for
535  C     beta      :: df/dy ( s^-1.m^-1 )  C     beta      :: df/dy ( s^-1.m^-1 )
536  C     omega     :: Angular velocity ( rad/s )  C     omega     :: Angular velocity ( rad/s )
537  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
538    C     viscAr    :: Eddy viscosity coeff. for mixing of
539    C                 momentum vertically ( units of r^2/s )
540  C     viscAh    :: Eddy viscosity coeff. for mixing of  C     viscAh    :: Eddy viscosity coeff. for mixing of
541  C                 momentum laterally ( m^2/s )  C                 momentum laterally ( m^2/s )
542  C     viscAhW   :: Eddy viscosity coeff. for mixing of vertical  C     viscAhW   :: Eddy viscosity coeff. for mixing of vertical
543  C                 momentum laterally, no effect for hydrostatic  C                 momentum laterally, no effect for hydrostatic
544  C                 model, defaults to viscAh if unset ( m^2/s )  C                 model, defaults to viscAh if unset ( m^2/s )
545  C                 Not used if variable horiz. viscosity is used.  C                 Not used if variable horiz. viscosity is used.
 C     viscAr    :: Eddy viscosity coeff. for mixing of  
 C                 momentum vertically ( units of r^2/s )  
546  C     viscA4    :: Biharmonic viscosity coeff. for mixing of  C     viscA4    :: Biharmonic viscosity coeff. for mixing of
547  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
548  C     viscA4W   :: Biharmonic viscosity coeff. for mixing of vertical  C     viscA4W   :: Biharmonic viscosity coeff. for mixing of vertical
# Line 538  C     viscA4D   :: Biharmonic viscosity Line 557  C     viscA4D   :: Biharmonic viscosity
557  C                  (act on Divergence part) ( m^4/s )  C                  (act on Divergence part) ( m^4/s )
558  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally
559  C                  (act on Vorticity  part) ( m^4/s )  C                  (act on Vorticity  part) ( m^4/s )
560  C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))  C     viscC2leith  :: Leith non-dimensional viscosity factor (grad(vort))
561  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))
562  C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor (harmonic)  C     viscC4leith  :: Leith non-dimensional viscosity factor (grad(vort))
563  C     viscC4smag  :: Smagorinsky non-dimensional viscosity factor (biharmonic)  C     viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div))
564  C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of  C     viscC2smag   :: Smagorinsky non-dimensional viscosity factor (harmonic)
565  C                   momentum laterally ( m^2/s )  C     viscC4smag   :: Smagorinsky non-dimensional viscosity factor (biharmonic)
566  C     viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity  C     viscAhMax    :: Maximum eddy viscosity coeff. for mixing of
567  C                   coeff. for mixing of momentum laterally (non-dim)  C                    momentum laterally ( m^2/s )
568    C     viscAhReMax  :: Maximum gridscale Reynolds number for eddy viscosity
569    C                     coeff. for mixing of momentum laterally (non-dim)
570    C     viscAhGrid   :: non-dimensional grid-size dependent viscosity
571  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...
572  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent visc.  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent visc.
573  C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of  C     viscA4Max    :: Maximum biharmonic viscosity coeff. for mixing of
574  C                 momentum laterally ( m^4/s )  C                     momentum laterally ( m^4/s )
575  C     viscA4ReMax :: Maximum Gridscale Reynolds number for  C     viscA4ReMax  :: Maximum Gridscale Reynolds number for
576  C                  biharmonic viscosity coeff. momentum laterally (non-dim)  C                     biharmonic viscosity coeff. momentum laterally (non-dim)
577  C     viscAhGrid:: non-dimensional grid-size dependent viscosity  C     viscA4Grid   :: non-dimensional grid-size dependent bi-harmonic viscosity
 C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity  
578  C     viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ...  C     viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ...
579  C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity  C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity
 C     viscC4leith :: Leith non-dimensional viscosity factor (grad(vort))  
 C     viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div))  
580  C     diffKhT   :: Laplacian diffusion coeff. for mixing of  C     diffKhT   :: Laplacian diffusion coeff. for mixing of
581  C                 heat laterally ( m^2/s )  C                 heat laterally ( m^2/s )
582  C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.  C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.
583  C                 for mixing of heat vertically ( units of r^2/s )  C                 for mixing of heat vertically ( units of r^2/s )
584  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
585  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
586  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
587  C                 salt laterally ( m^2/s )  C                 salt laterally ( m^2/s )
588  C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.  C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.
589  C                 for mixing of salt vertically ( units of r^2/s ),  C                 for mixing of salt vertically ( units of r^2/s ),
590  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
591  C                 salt laterally ( m^4/s )  C                 salt laterally ( m^4/s )
592  C     diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979  C     diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979
# Line 626  C     mtFacMom      :: Metric terms trac Line 645  C     mtFacMom      :: Metric terms trac
645  C     cosPower      :: Power of cosine of latitude to multiply viscosity  C     cosPower      :: Power of cosine of latitude to multiply viscosity
646  C     cAdjFreq      :: Frequency of convective adjustment  C     cAdjFreq      :: Frequency of convective adjustment
647  C  C
648  C     taveFreq      :: Frequency with which time-averaged model state  C     taveFreq      :: Frequency with which time-averaged model state
649  C                      is written to post-processing files ( s ).  C                      is written to post-processing files ( s ).
650  C     tave_lastIter :: (for state variable only) fraction of the last time  C     tave_lastIter :: (for state variable only) fraction of the last time
651  C                      step (of each taveFreq period) put in the time average.  C                      step (of each taveFreq period) put in the time average.
652  C                      (fraction for 1rst iter = 1 - tave_lastIter)  C                      (fraction for 1rst iter = 1 - tave_lastIter)
653  C     tauThetaClimRelax :: Relaxation to climatology time scale ( s ).  C     tauThetaClimRelax :: Relaxation to climatology time scale ( s ).
654  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
# Line 641  C                          (note: extern Line 660  C                          (note: extern
660  C                           number times externForcingPeriod)  C                           number times externForcingPeriod)
661  C     convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux  C     convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux
662  C                       (use model surface (local) value if set to -1)  C                       (use model surface (local) value if set to -1)
663  C     temp_EvPrRn :: temperature of Rain & Evap.  C     temp_EvPrRn :: temperature of Rain & Evap.
664  C     salt_EvPrRn :: salinity of Rain & Evap.  C     salt_EvPrRn :: salinity of Rain & Evap.
665  C        (notes: a) tracer content of Rain/Evap only used if both  C        (notes: a) tracer content of Rain/Evap only used if both
666  C                     NonLin_FrSurf & useRealFreshWater are set.  C                     NonLin_FrSurf & useRealFreshWater are set.
667  C                b) use model surface (local) value if set to UNSET_RL)  C                b) use model surface (local) value if set to UNSET_RL)
668  C     hMixCrit    :: criteria for mixed-layer diagnostic  C     hMixCrit    :: criteria for mixed-layer diagnostic
669  C     ivdc_kappa  :: implicit vertical diffusivity for convection [m^2/s]  C     ivdc_kappa  :: implicit vertical diffusivity for convection [m^2/s]
670  C     Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as  C     Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as
671  C                    starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)  C                    starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
672  C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)  C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)
673  C                           (default=2: full drag ; =1: gives half-slip BC)  C                           (default=2: full drag ; =1: gives half-slip BC)
674  C     bottomDragLinear    :: Linear    bottom-drag coefficient (units of [r]/s)  C     bottomDragLinear    :: Linear    bottom-drag coefficient (units of [r]/s)
675  C     bottomDragQuadratic :: Quadratic bottom-drag coefficient (units of [r]/m)  C     bottomDragQuadratic :: Quadratic bottom-drag coefficient (units of [r]/m)
# Line 660  C                           is to be rep Line 679  C                           is to be rep
679  C                           (affects myabs, mymin, mymax)  C                           (affects myabs, mymin, mymax)
680  C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales  C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales
681  C                      (i.e. allows convection at different Rayleigh numbers)  C                      (i.e. allows convection at different Rayleigh numbers)
682    C     phiEuler      :: Euler angle, rotation about original z-axis
683    C     thetaEuler    :: Euler angle, rotation about new x-axis
684    C     psiEuler      :: Euler angle, rotation about new z-axis
685        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
686       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
687       & delR, delRc, delX, delY,       & delR, delRc, delX, delY,
# Line 671  C                      (i.e. allows conv Line 693  C                      (i.e. allows conv
693       & viscC2leith, viscC2leithD,       & viscC2leith, viscC2leithD,
694       & viscC2smag, viscC4smag,       & viscC2smag, viscC4smag,
695       & viscAhD, viscAhZ, viscA4D, viscA4Z,       & viscAhD, viscAhZ, viscA4D, viscA4Z,
696       & viscA4, viscA4W,       & viscA4, viscA4W, viscA4Max,
697       & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,       & viscA4Grid, viscA4GridMax, viscA4GridMin,
698       & viscAhRemax, viscA4Remax,       & viscAhReMax, viscA4ReMax,
699       & viscC4leith, viscC4leithD, viscAr,       & viscC4leith, viscC4leithD, viscAr,
700       & diffKhT, diffK4T, diffKrNrT,       & diffKhT, diffK4T, diffKrNrT,
701       & diffKhS, diffK4S, diffKrNrS,       & diffKhS, diffK4S, diffKrNrS,
# Line 692  C                      (i.e. allows conv Line 714  C                      (i.e. allows conv
714       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,
715       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
716       & cosPower, cAdjFreq, omega, rotationPeriod,       & cosPower, cAdjFreq, omega, rotationPeriod,
717       & tauThetaClimRelax, tauSaltClimRelax,       & tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax,
      & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,  
718       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
719       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
720       & hFacMinDr, hFacMinDp,       & hFacMinDr, hFacMinDp,
721       & ivdc_kappa, hMixCriteria, Ro_SeaLevel,       & ivdc_kappa, hMixCriteria, Ro_SeaLevel,
722       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
723       & smoothAbsFuncRange,       & smoothAbsFuncRange,
724       & tCylIn, tCylOut       & tCylIn, tCylOut,
725         & phiEuler, thetaEuler, psiEuler
726    
727        _RL cg2dTargetResidual        _RL cg2dTargetResidual
728        _RL cg2dTargetResWunit        _RL cg2dTargetResWunit
# Line 731  C                      (i.e. allows conv Line 753  C                      (i.e. allows conv
753        _RL hFacInf        _RL hFacInf
754        _RL hFacSup        _RL hFacSup
755        _RL beta        _RL beta
756          _RL viscAr
757        _RL viscAh        _RL viscAh
758        _RL viscAhW        _RL viscAhW
759        _RL viscAhD        _RL viscAhD
760        _RL viscAhZ        _RL viscAhZ
761        _RL viscAhMax        _RL viscAhMax
762        _RL viscAhReMax        _RL viscAhReMax
763        _RL viscAhGrid        _RL viscAhGrid, viscAhGridMax, viscAhGridMin
       _RL viscAhGridMax  
       _RL viscAhGridMin  
764        _RL viscC2leith        _RL viscC2leith
765        _RL viscC2leithD        _RL viscC2leithD
766        _RL viscC2smag        _RL viscC2smag
       _RL viscC4smag  
       _RL viscAr  
767        _RL viscA4        _RL viscA4
768        _RL viscA4W        _RL viscA4W
769        _RL viscA4D        _RL viscA4D
# Line 754  C                      (i.e. allows conv Line 773  C                      (i.e. allows conv
773        _RL viscA4Grid, viscA4GridMax, viscA4GridMin        _RL viscA4Grid, viscA4GridMax, viscA4GridMin
774        _RL viscC4leith        _RL viscC4leith
775        _RL viscC4leithD        _RL viscC4leithD
776          _RL viscC4smag
777        _RL diffKhT        _RL diffKhT
778        _RL diffKrNrT(Nr)        _RL diffKrNrT(Nr)
779        _RL diffK4T        _RL diffK4T
# Line 810  C                      (i.e. allows conv Line 830  C                      (i.e. allows conv
830        _RL rotationPeriod        _RL rotationPeriod
831        _RL tauThetaClimRelax        _RL tauThetaClimRelax
832        _RL tauSaltClimRelax        _RL tauSaltClimRelax
       _RL tauTr1ClimRelax  
       _RL lambdaTr1ClimRelax  
833        _RL latBandClimRelax        _RL latBandClimRelax
834        _RL externForcingCycle        _RL externForcingCycle
835        _RL externForcingPeriod        _RL externForcingPeriod
# Line 828  C                      (i.e. allows conv Line 846  C                      (i.e. allows conv
846        _RL nh_Am2        _RL nh_Am2
847        _RL tCylIn        _RL tCylIn
848        _RL tCylOut        _RL tCylOut
849          _RL phiEuler, thetaEuler, psiEuler
850    
851  C--   COMMON /PARM_A/ Thermodynamics constants ?  C--   COMMON /PARM_A/ Thermodynamics constants ?
852        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp
# Line 843  C     atm_kappa :: kappa = R/Cp (R: cons Line 862  C     atm_kappa :: kappa = R/Cp (R: cons
862  C     atm_Rq    :: water vapour specific volume anomaly relative to dry air  C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
863  C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])  C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
864  C     integr_GeoPot :: option to select the way we integrate the geopotential  C     integr_GeoPot :: option to select the way we integrate the geopotential
865  C                     (still a subject of discussions ...)  C                     (still a subject of discussions ...)
866  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
867  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
868        COMMON /PARM_ATM/        COMMON /PARM_ATM/
Legend:
Removed from v.1.209  
changed lines
  Added in v.1.217
  ViewVC Help
Powered by ViewVC 1.1.22