38 |
C rhoRefFile :: File containing reference density profile rhoRef (1.D) |
C rhoRefFile :: File containing reference density profile rhoRef (1.D) |
39 |
C delRFile :: File containing vertical grid spacing delR (1.D array) |
C delRFile :: File containing vertical grid spacing delR (1.D array) |
40 |
C delRcFile :: File containing vertical grid spacing delRc (1.D array) |
C delRcFile :: File containing vertical grid spacing delRc (1.D array) |
41 |
|
C hybSigmFile :: File containing hybrid-sigma vertical coord. coeff. (2x 1.D) |
42 |
C delXFile :: File containing X-spacing grid definition (1.D array) |
C delXFile :: File containing X-spacing grid definition (1.D array) |
43 |
C delYFile :: File containing Y-spacing grid definition (1.D array) |
C delYFile :: File containing Y-spacing grid definition (1.D array) |
44 |
C horizGridFile :: File containing horizontal-grid definition |
C horizGridFile :: File containing horizontal-grid definition |
80 |
& buoyancyRelation, eosType, |
& buoyancyRelation, eosType, |
81 |
& pickupSuff, mdsioLocalDir, |
& pickupSuff, mdsioLocalDir, |
82 |
& tRefFile, sRefFile, rhoRefFile, |
& tRefFile, sRefFile, rhoRefFile, |
83 |
& delRFile, delRcFile, |
& delRFile, delRcFile, hybSigmFile, |
84 |
& delXFile, delYFile, horizGridFile, |
& delXFile, delYFile, horizGridFile, |
85 |
& bathyFile, topoFile, |
& bathyFile, topoFile, |
86 |
& viscAhDfile, viscAhZfile, |
& viscAhDfile, viscAhZfile, |
87 |
& viscA4Dfile, viscA4Zfile, |
& viscA4Dfile, viscA4Zfile, |
88 |
& hydrogThetaFile, hydrogSaltFile, diffKrFile, |
& hydrogThetaFile, hydrogSaltFile, diffKrFile, |
89 |
& zonalWindFile, meridWindFile, thetaClimFile, |
& zonalWindFile, meridWindFile, thetaClimFile, |
90 |
& saltClimFile, |
& saltClimFile, |
104 |
CHARACTER*(MAX_LEN_FNAM) rhoRefFile |
CHARACTER*(MAX_LEN_FNAM) rhoRefFile |
105 |
CHARACTER*(MAX_LEN_FNAM) delRFile |
CHARACTER*(MAX_LEN_FNAM) delRFile |
106 |
CHARACTER*(MAX_LEN_FNAM) delRcFile |
CHARACTER*(MAX_LEN_FNAM) delRcFile |
107 |
|
CHARACTER*(MAX_LEN_FNAM) hybSigmFile |
108 |
CHARACTER*(MAX_LEN_FNAM) delXFile |
CHARACTER*(MAX_LEN_FNAM) delXFile |
109 |
CHARACTER*(MAX_LEN_FNAM) delYFile |
CHARACTER*(MAX_LEN_FNAM) delYFile |
110 |
CHARACTER*(MAX_LEN_FNAM) horizGridFile |
CHARACTER*(MAX_LEN_FNAM) horizGridFile |
151 |
C writeStatePrec :: Precision used for writing model state. |
C writeStatePrec :: Precision used for writing model state. |
152 |
C writeBinaryPrec :: Precision used for writing binary files |
C writeBinaryPrec :: Precision used for writing binary files |
153 |
C readBinaryPrec :: Precision used for reading binary files |
C readBinaryPrec :: Precision used for reading binary files |
154 |
|
C selectCoriMap :: select setting of Coriolis parameter map: |
155 |
|
C =0 f-Plane (Constant Coriolis, = f0) |
156 |
|
C =1 Beta-Plane Coriolis (= f0 + beta.y) |
157 |
|
C =2 Spherical Coriolis (= 2.omega.sin(phi)) |
158 |
|
C =3 Read Coriolis 2-d fields from files. |
159 |
|
C selectSigmaCoord :: option related to sigma vertical coordinate |
160 |
C nonlinFreeSurf :: option related to non-linear free surface |
C nonlinFreeSurf :: option related to non-linear free surface |
161 |
C =0 Linear free surface ; >0 Non-linear |
C =0 Linear free surface ; >0 Non-linear |
162 |
C select_rStar :: option related to r* vertical coordinate |
C select_rStar :: option related to r* vertical coordinate |
178 |
C selectVortScheme :: Scheme selector for Vorticity term (Vector Inv.) |
C selectVortScheme :: Scheme selector for Vorticity term (Vector Inv.) |
179 |
C monitorSelect :: select group of variables to monitor |
C monitorSelect :: select group of variables to monitor |
180 |
C =1 : dynvars ; =2 : + vort ; =3 : + surface |
C =1 : dynvars ; =2 : + vort ; =3 : + surface |
181 |
C debugLevel :: debug level selector: higher -> more writing |
C- debugLevel :: controls printing of algorithm intermediate results |
182 |
|
C and statistics ; higher -> more writing |
183 |
|
|
184 |
COMMON /PARM_I/ |
COMMON /PARM_I/ |
185 |
& cg2dMaxIters, |
& cg2dMaxIters, |
189 |
& nIter0, nTimeSteps, nEndIter, |
& nIter0, nTimeSteps, nEndIter, |
190 |
& writeStatePrec, |
& writeStatePrec, |
191 |
& writeBinaryPrec, readBinaryPrec, |
& writeBinaryPrec, readBinaryPrec, |
192 |
|
& selectCoriMap, |
193 |
|
& selectSigmaCoord, |
194 |
& nonlinFreeSurf, select_rStar, |
& nonlinFreeSurf, select_rStar, |
195 |
& selectNHfreeSurf, |
& selectNHfreeSurf, |
196 |
& selectAddFluid, |
& selectAddFluid, |
210 |
INTEGER writeStatePrec |
INTEGER writeStatePrec |
211 |
INTEGER writeBinaryPrec |
INTEGER writeBinaryPrec |
212 |
INTEGER readBinaryPrec |
INTEGER readBinaryPrec |
213 |
|
INTEGER selectCoriMap |
214 |
|
INTEGER selectSigmaCoord |
215 |
INTEGER nonlinFreeSurf |
INTEGER nonlinFreeSurf |
216 |
INTEGER select_rStar |
INTEGER select_rStar |
217 |
INTEGER selectNHfreeSurf |
INTEGER selectNHfreeSurf |
265 |
C useNHMTerms :: If TRUE use non-hydrostatic metric terms. |
C useNHMTerms :: If TRUE use non-hydrostatic metric terms. |
266 |
C useCoriolis :: Flag which turns the coriolis terms on and off. |
C useCoriolis :: Flag which turns the coriolis terms on and off. |
267 |
C use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - off |
C use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - 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) |
|
268 |
C useCDscheme :: use CD-scheme to calculate Coriolis terms. |
C useCDscheme :: use CD-scheme to calculate Coriolis terms. |
269 |
C vectorInvariantMomentum :: use Vector-Invariant form (mom_vecinv package) |
C vectorInvariantMomentum :: use Vector-Invariant form (mom_vecinv package) |
270 |
C (default = F = use mom_fluxform package) |
C (default = F = use mom_fluxform package) |
333 |
C useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12) |
C useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12) |
334 |
C periodicExternalForcing :: Set true if forcing is time-dependant |
C periodicExternalForcing :: Set true if forcing is time-dependant |
335 |
C- I/O parameters - |
C- I/O parameters - |
336 |
C globalFiles :: Selects between "global" and "tiled" files |
C globalFiles :: Selects between "global" and "tiled" files. |
337 |
C useSingleCpuIO :: On SGI platforms, option globalFiles is either |
C On some platforms with MPI, option globalFiles is either |
338 |
C slow (f77) or does not work (f90). When |
C slow or does not work. Use useSingleCpuIO instead. |
339 |
C useSingleCpuIO is set, mdsio_writefield.F |
C useSingleCpuIO :: moved to EEPARAMS.h |
|
C outputs from master mpi process only. |
|
340 |
C pickupStrictlyMatch :: check and stop if pickup-file do not stricly match |
C pickupStrictlyMatch :: check and stop if pickup-file do not stricly match |
341 |
C startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup |
C startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup |
342 |
C usePickupBeforeC54 :: start from old-pickup files, generated with code from |
C usePickupBeforeC54 :: start from old-pickup files, generated with code from |
350 |
C monitor_stdio :: use stdio for monitor output |
C monitor_stdio :: use stdio for monitor output |
351 |
C dumpInitAndLast :: dumps model state to files at Initial (nIter0) |
C dumpInitAndLast :: dumps model state to files at Initial (nIter0) |
352 |
C & Last iteration, in addition multiple of dumpFreq iter. |
C & Last iteration, in addition multiple of dumpFreq iter. |
353 |
|
C debugMode :: controls printing of debug msg (sequence of S/R calls). |
354 |
|
|
355 |
COMMON /PARM_L/ |
COMMON /PARM_L/ |
356 |
& fluidIsAir, fluidIsWater, |
& fluidIsAir, fluidIsWater, |
363 |
& momViscosity, momAdvection, momForcing, |
& momViscosity, momAdvection, momForcing, |
364 |
& momPressureForcing, metricTerms, useNHMTerms, |
& momPressureForcing, metricTerms, useNHMTerms, |
365 |
& useCoriolis, use3dCoriolis, |
& useCoriolis, use3dCoriolis, |
|
& useConstantF, useBetaPlaneF, useSphereF, |
|
366 |
& useCDscheme, vectorInvariantMomentum, |
& useCDscheme, vectorInvariantMomentum, |
367 |
& useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv, |
& useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv, |
368 |
& upwindVorticity, highOrderVorticity, |
& upwindVorticity, highOrderVorticity, |
384 |
& doThetaClimRelax, doSaltClimRelax, |
& doThetaClimRelax, doSaltClimRelax, |
385 |
& allowFreezing, useOldFreezing, |
& allowFreezing, useOldFreezing, |
386 |
& periodicExternalForcing, |
& periodicExternalForcing, |
387 |
& globalFiles, useSingleCpuIO, |
& globalFiles, |
388 |
& pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2, |
& pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2, |
389 |
& pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed, |
& pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed, |
390 |
& writePickupAtEnd, |
& writePickupAtEnd, |
442 |
LOGICAL tempStepping |
LOGICAL tempStepping |
443 |
LOGICAL saltStepping |
LOGICAL saltStepping |
444 |
LOGICAL metricTerms |
LOGICAL metricTerms |
|
LOGICAL useConstantF |
|
|
LOGICAL useBetaPlaneF |
|
|
LOGICAL useSphereF |
|
445 |
LOGICAL use3dCoriolis |
LOGICAL use3dCoriolis |
446 |
LOGICAL useCDscheme |
LOGICAL useCDscheme |
447 |
LOGICAL useEnergyConservingCoriolis |
LOGICAL useEnergyConservingCoriolis |
469 |
LOGICAL useOldFreezing |
LOGICAL useOldFreezing |
470 |
LOGICAL periodicExternalForcing |
LOGICAL periodicExternalForcing |
471 |
LOGICAL globalFiles |
LOGICAL globalFiles |
|
LOGICAL useSingleCpuIO |
|
472 |
LOGICAL pickupStrictlyMatch |
LOGICAL pickupStrictlyMatch |
473 |
LOGICAL usePickupBeforeC54 |
LOGICAL usePickupBeforeC54 |
474 |
LOGICAL startFromPickupAB2 |
LOGICAL startFromPickupAB2 |
545 |
C f0 :: Reference coriolis parameter ( 1/s ) |
C f0 :: Reference coriolis parameter ( 1/s ) |
546 |
C ( Southern edge f for beta plane ) |
C ( Southern edge f for beta plane ) |
547 |
C beta :: df/dy ( s^-1.m^-1 ) |
C beta :: df/dy ( s^-1.m^-1 ) |
548 |
|
C fPrime :: Second Coriolis parameter ( 1/s ), related to Y-component |
549 |
|
C of rotation (reference value = 2.Omega.Cos(Phi)) |
550 |
C omega :: Angular velocity ( rad/s ) |
C omega :: Angular velocity ( rad/s ) |
551 |
C rotationPeriod :: Rotation period (s) (= 2.pi/omega) |
C rotationPeriod :: Rotation period (s) (= 2.pi/omega) |
552 |
C viscArNr :: vertical profile of Eddy viscosity coeff. |
C viscArNr :: vertical profile of Eddy viscosity coeff. |
690 |
C ivdc_kappa :: implicit vertical diffusivity for convection [m^2/s] |
C ivdc_kappa :: implicit vertical diffusivity for convection [m^2/s] |
691 |
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 |
692 |
C starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel) |
C starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel) |
693 |
|
C rSigmaBnd :: vertical position (in r-unit) of r/sigma transition (Hybrid-Sigma) |
694 |
C sideDragFactor :: side-drag scaling factor (used only if no_slip_sides) |
C sideDragFactor :: side-drag scaling factor (used only if no_slip_sides) |
695 |
C (default=2: full drag ; =1: gives half-slip BC) |
C (default=2: full drag ; =1: gives half-slip BC) |
696 |
C bottomDragLinear :: Linear bottom-drag coefficient (units of [r]/s) |
C bottomDragLinear :: Linear bottom-drag coefficient (units of [r]/s) |
709 |
& delR, delRc, delX, delY, |
& delR, delRc, delX, delY, |
710 |
& deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock, |
& deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock, |
711 |
& abEps, alph_AB, beta_AB, |
& abEps, alph_AB, beta_AB, |
712 |
& xgOrigin, ygOrigin, rSphere, recip_RSphere, f0, beta, |
& xgOrigin, ygOrigin, rSphere, recip_rSphere, |
713 |
& viscAh, viscAhW, viscAhMax, |
& f0, beta, fPrime, omega, rotationPeriod, |
714 |
|
& viscFacAdj, viscAh, viscAhW, viscAhMax, |
715 |
& viscAhGrid, viscAhGridMax, viscAhGridMin, |
& viscAhGrid, viscAhGridMax, viscAhGridMin, |
716 |
& viscC2leith, viscC2leithD, |
& viscC2leith, viscC2leithD, |
717 |
& viscC2smag, viscC4smag, |
& viscC2smag, viscC4smag, |
737 |
& chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq, |
& chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq, |
738 |
& diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq, |
& diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq, |
739 |
& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom, |
& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom, |
740 |
& cosPower, cAdjFreq, omega, rotationPeriod, |
& cosPower, cAdjFreq, |
741 |
& tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax, |
& tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax, |
742 |
& externForcingCycle, externForcingPeriod, |
& externForcingCycle, externForcingPeriod, |
743 |
& convertFW2Salt, temp_EvPrRn, salt_EvPrRn, |
& convertFW2Salt, temp_EvPrRn, salt_EvPrRn, |
744 |
& temp_addMass, salt_addMass, hFacMinDr, hFacMinDp, |
& temp_addMass, salt_addMass, hFacMinDr, hFacMinDp, |
745 |
& ivdc_kappa, hMixCriteria, dRhoSmall, hMixSmooth, Ro_SeaLevel, |
& ivdc_kappa, hMixCriteria, dRhoSmall, hMixSmooth, |
746 |
|
& Ro_SeaLevel, rSigmaBnd, |
747 |
& sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2, |
& sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2, |
748 |
& smoothAbsFuncRange, |
& smoothAbsFuncRange, |
749 |
& tCylIn, tCylOut, |
& tCylIn, tCylOut, |
768 |
_RL rSphere |
_RL rSphere |
769 |
_RL recip_rSphere |
_RL recip_rSphere |
770 |
_RL f0 |
_RL f0 |
771 |
|
_RL beta |
772 |
|
_RL fPrime |
773 |
|
_RL omega |
774 |
|
_RL rotationPeriod |
775 |
_RL freeSurfFac |
_RL freeSurfFac |
776 |
_RL implicSurfPress |
_RL implicSurfPress |
777 |
_RL implicDiv2Dflow |
_RL implicDiv2Dflow |
782 |
_RL hFacMinDr |
_RL hFacMinDr |
783 |
_RL hFacInf |
_RL hFacInf |
784 |
_RL hFacSup |
_RL hFacSup |
|
_RL beta |
|
785 |
_RL viscArNr(Nr) |
_RL viscArNr(Nr) |
786 |
|
_RL viscFacAdj |
787 |
_RL viscAh |
_RL viscAh |
788 |
_RL viscAhW |
_RL viscAhW |
789 |
_RL viscAhD |
_RL viscAhD |
855 |
_RL mtFacMom |
_RL mtFacMom |
856 |
_RL cosPower |
_RL cosPower |
857 |
_RL cAdjFreq |
_RL cAdjFreq |
|
_RL omega |
|
|
_RL rotationPeriod |
|
858 |
_RL tauThetaClimRelax |
_RL tauThetaClimRelax |
859 |
_RL tauSaltClimRelax |
_RL tauSaltClimRelax |
860 |
_RL latBandClimRelax |
_RL latBandClimRelax |
870 |
_RL dRhoSmall |
_RL dRhoSmall |
871 |
_RL hMixSmooth |
_RL hMixSmooth |
872 |
_RL Ro_SeaLevel |
_RL Ro_SeaLevel |
873 |
|
_RL rSigmaBnd |
874 |
_RL sideDragFactor |
_RL sideDragFactor |
875 |
_RL bottomDragLinear |
_RL bottomDragLinear |
876 |
_RL bottomDragQuadratic |
_RL bottomDragQuadratic |