--- MITgcm/model/inc/PARAMS.h 2010/10/01 02:00:40 1.241 +++ MITgcm/model/inc/PARAMS.h 2010/12/06 10:45:51 1.245 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/inc/PARAMS.h,v 1.241 2010/10/01 02:00:40 heimbach Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/inc/PARAMS.h,v 1.245 2010/12/06 10:45:51 mlosch Exp $ C $Name: $ C @@ -33,6 +33,11 @@ C and read pickup files with this suffix (max 10 Char.) C mdsioLocalDir :: read-write tiled file from/to this directory name C (+ 4 digits Processor-Rank) instead of current dir. +C adTapeDir :: read-write checkpointing tape files from/to this +C directory name instead of current dir. Conflicts +C mdsioLocalDir, so only one of the two can be set. +C In contrast to mdsioLocalDir, if specified adTapeDir +C must exist before the model starts. C tRefFile :: File containing reference Potential Temperat. tRef (1.D) C sRefFile :: File containing reference salinity/spec.humid. sRef (1.D) C rhoRefFile :: File containing reference density profile rhoRef (1.D) @@ -78,7 +83,7 @@ C the_run_name :: string identifying the name of the model "run" COMMON /PARM_C/ & buoyancyRelation, eosType, - & pickupSuff, mdsioLocalDir, + & pickupSuff, mdsioLocalDir, adTapeDir, & tRefFile, sRefFile, rhoRefFile, & delRFile, delRcFile, hybSigmFile, & delXFile, delYFile, horizGridFile, @@ -99,6 +104,7 @@ CHARACTER*(6) eosType CHARACTER*(10) pickupSuff CHARACTER*(MAX_LEN_FNAM) mdsioLocalDir + CHARACTER*(MAX_LEN_FNAM) adTapeDir CHARACTER*(MAX_LEN_FNAM) tRefFile CHARACTER*(MAX_LEN_FNAM) sRefFile CHARACTER*(MAX_LEN_FNAM) rhoRefFile @@ -151,13 +157,18 @@ C writeStatePrec :: Precision used for writing model state. C writeBinaryPrec :: Precision used for writing binary files C readBinaryPrec :: Precision used for reading binary files +C selectCoriMap :: select setting of Coriolis parameter map: +C =0 f-Plane (Constant Coriolis, = f0) +C =1 Beta-Plane Coriolis (= f0 + beta.y) +C =2 Spherical Coriolis (= 2.omega.sin(phi)) +C =3 Read Coriolis 2-d fields from files. +C selectSigmaCoord :: option related to sigma vertical coordinate C nonlinFreeSurf :: option related to non-linear free surface C =0 Linear free surface ; >0 Non-linear C select_rStar :: option related to r* vertical coordinate C =0 (default) use r coord. ; > 0 use r* C selectNHfreeSurf :: option for Non-Hydrostatic (free-)Surface formulation: C =0 (default) hydrostatic surf. ; > 0 add NH effects. -C selectSigmaCoord :: option related to sigma vertical coordinate C selectAddFluid :: option to add mass source/sink of fluid in the interior C (3-D generalisation of oceanic real-fresh water flux) C =0 off ; =1 add fluid ; =-1 virtual flux (no mass added) @@ -173,7 +184,8 @@ C selectVortScheme :: Scheme selector for Vorticity term (Vector Inv.) C monitorSelect :: select group of variables to monitor C =1 : dynvars ; =2 : + vort ; =3 : + surface -C debugLevel :: debug level selector: higher -> more writing +C- debugLevel :: controls printing of algorithm intermediate results +C and statistics ; higher -> more writing COMMON /PARM_I/ & cg2dMaxIters, @@ -183,9 +195,10 @@ & nIter0, nTimeSteps, nEndIter, & writeStatePrec, & writeBinaryPrec, readBinaryPrec, + & selectCoriMap, + & selectSigmaCoord, & nonlinFreeSurf, select_rStar, & selectNHfreeSurf, - & selectSigmaCoord, & selectAddFluid, & momForcingOutAB, tracForcingOutAB, & tempAdvScheme, tempVertAdvScheme, @@ -203,10 +216,11 @@ INTEGER writeStatePrec INTEGER writeBinaryPrec INTEGER readBinaryPrec + INTEGER selectCoriMap + INTEGER selectSigmaCoord INTEGER nonlinFreeSurf INTEGER select_rStar INTEGER selectNHfreeSurf - INTEGER selectSigmaCoord INTEGER selectAddFluid INTEGER momForcingOutAB, tracForcingOutAB INTEGER tempAdvScheme, tempVertAdvScheme @@ -257,9 +271,6 @@ C useNHMTerms :: If TRUE use non-hydrostatic metric terms. C useCoriolis :: Flag which turns the coriolis terms on and 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) C useCDscheme :: use CD-scheme to calculate Coriolis terms. C vectorInvariantMomentum :: use Vector-Invariant form (mom_vecinv package) C (default = F = use mom_fluxform package) @@ -345,6 +356,7 @@ C monitor_stdio :: use stdio for monitor output C dumpInitAndLast :: dumps model state to files at Initial (nIter0) C & Last iteration, in addition multiple of dumpFreq iter. +C debugMode :: controls printing of debug msg (sequence of S/R calls). COMMON /PARM_L/ & fluidIsAir, fluidIsWater, @@ -357,7 +369,6 @@ & momViscosity, momAdvection, momForcing, & momPressureForcing, metricTerms, useNHMTerms, & useCoriolis, use3dCoriolis, - & useConstantF, useBetaPlaneF, useSphereF, & useCDscheme, vectorInvariantMomentum, & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv, & upwindVorticity, highOrderVorticity, @@ -437,9 +448,6 @@ LOGICAL tempStepping LOGICAL saltStepping LOGICAL metricTerms - LOGICAL useConstantF - LOGICAL useBetaPlaneF - LOGICAL useSphereF LOGICAL use3dCoriolis LOGICAL useCDscheme LOGICAL useEnergyConservingCoriolis @@ -543,6 +551,8 @@ C f0 :: Reference coriolis parameter ( 1/s ) C ( Southern edge f for beta plane ) C beta :: df/dy ( s^-1.m^-1 ) +C fPrime :: Second Coriolis parameter ( 1/s ), related to Y-component +C of rotation (reference value = 2.Omega.Cos(Phi)) C omega :: Angular velocity ( rad/s ) C rotationPeriod :: Rotation period (s) (= 2.pi/omega) C viscArNr :: vertical profile of Eddy viscosity coeff. @@ -705,7 +715,8 @@ & delR, delRc, delX, delY, & deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock, & abEps, alph_AB, beta_AB, - & xgOrigin, ygOrigin, rSphere, recip_RSphere, f0, beta, + & xgOrigin, ygOrigin, rSphere, recip_rSphere, + & f0, beta, fPrime, omega, rotationPeriod, & viscFacAdj, viscAh, viscAhW, viscAhMax, & viscAhGrid, viscAhGridMax, viscAhGridMin, & viscC2leith, viscC2leithD, @@ -732,7 +743,7 @@ & chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq, & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq, & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom, - & cosPower, cAdjFreq, omega, rotationPeriod, + & cosPower, cAdjFreq, & tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax, & externForcingCycle, externForcingPeriod, & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, @@ -763,6 +774,10 @@ _RL rSphere _RL recip_rSphere _RL f0 + _RL beta + _RL fPrime + _RL omega + _RL rotationPeriod _RL freeSurfFac _RL implicSurfPress _RL implicDiv2Dflow @@ -773,7 +788,6 @@ _RL hFacMinDr _RL hFacInf _RL hFacSup - _RL beta _RL viscArNr(Nr) _RL viscFacAdj _RL viscAh @@ -847,8 +861,6 @@ _RL mtFacMom _RL cosPower _RL cAdjFreq - _RL omega - _RL rotationPeriod _RL tauThetaClimRelax _RL tauSaltClimRelax _RL latBandClimRelax @@ -942,7 +954,6 @@ LOGICAL useMNC LOGICAL useRunClock LOGICAL useEMBED_FILES - LOGICAL useOASIS LOGICAL useMYPACKAGE COMMON /PARM_PACKAGES/ & useGAD, useOBCS, useSHAP_FILT, useZONAL_FILT, @@ -956,8 +967,8 @@ & useATM2D, useAIM, useLand, useFizhi, useGridAlt, & useDiagnostics, useREGRID, useLayers, useMNC, & useRunClock, useEMBED_FILES, - & useOASIS, & useMYPACKAGE + CEH3 ;;; Local Variables: *** CEH3 ;;; mode:fortran *** CEH3 ;;; End: ***