14 |
|
|
15 |
CEOP |
CEOP |
16 |
|
|
|
C Macros for special grid options |
|
|
#include "PARAMS_MACROS.h" |
|
|
|
|
17 |
C-- Contants |
C-- Contants |
18 |
C Useful physical values |
C Useful physical values |
19 |
Real*8 PI |
Real*8 PI |
68 |
C (old version, kept for backward compatibility) |
C (old version, kept for backward compatibility) |
69 |
C surfQnetFile :: File containing surface net heat flux |
C surfQnetFile :: File containing surface net heat flux |
70 |
C surfQswFile :: File containing surface shortwave radiation |
C surfQswFile :: File containing surface shortwave radiation |
|
C dQdTfile :: File containing thermal relaxation coefficient |
|
71 |
C EmPmRfile :: File containing surface fresh water flux |
C EmPmRfile :: File containing surface fresh water flux |
72 |
C NOTE: for backward compatibility EmPmRfile is specified in |
C NOTE: for backward compatibility EmPmRfile is specified in |
73 |
C m/s when using external_fields_load.F. It is converted |
C m/s when using external_fields_load.F. It is converted |
94 |
& surfQfile, surfQnetFile, surfQswFile, |
& surfQfile, surfQnetFile, surfQswFile, |
95 |
& lambdaThetaFile, lambdaSaltFile, |
& lambdaThetaFile, lambdaSaltFile, |
96 |
& uVelInitFile, vVelInitFile, pSurfInitFile, |
& uVelInitFile, vVelInitFile, pSurfInitFile, |
97 |
& dQdTfile, pLoadFile, addMassFile, |
& pLoadFile, addMassFile, |
98 |
& eddyPsiXFile, eddyPsiYFile, |
& eddyPsiXFile, eddyPsiYFile, |
99 |
& the_run_name |
& the_run_name |
100 |
CHARACTER*(MAX_LEN_FNAM) buoyancyRelation |
CHARACTER*(MAX_LEN_FNAM) buoyancyRelation |
130 |
CHARACTER*(MAX_LEN_FNAM) uVelInitFile |
CHARACTER*(MAX_LEN_FNAM) uVelInitFile |
131 |
CHARACTER*(MAX_LEN_FNAM) vVelInitFile |
CHARACTER*(MAX_LEN_FNAM) vVelInitFile |
132 |
CHARACTER*(MAX_LEN_FNAM) pSurfInitFile |
CHARACTER*(MAX_LEN_FNAM) pSurfInitFile |
|
CHARACTER*(MAX_LEN_FNAM) dQdTfile |
|
133 |
CHARACTER*(MAX_LEN_FNAM) pLoadFile |
CHARACTER*(MAX_LEN_FNAM) pLoadFile |
134 |
CHARACTER*(MAX_LEN_FNAM) addMassFile |
CHARACTER*(MAX_LEN_FNAM) addMassFile |
135 |
CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile |
CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile |
145 |
C in con. grad solver. |
C in con. grad solver. |
146 |
C cg2dPreCondFreq :: Frequency for updating cg2d preconditioner |
C cg2dPreCondFreq :: Frequency for updating cg2d preconditioner |
147 |
C (non-linear free-surf.) |
C (non-linear free-surf.) |
148 |
|
C cg2dUseMinResSol :: =0 : use last-iteration/converged solution |
149 |
|
C =1 : use solver minimum-residual solution |
150 |
C cg3dMaxIters :: Maximum number of iterations in the |
C cg3dMaxIters :: Maximum number of iterations in the |
151 |
C three-dimensional con. grad solver. |
C three-dimensional con. grad solver. |
152 |
C cg3dChkResFreq :: Frequency with which to check residual |
C cg3dChkResFreq :: Frequency with which to check residual |
188 |
C and statistics ; higher -> more writing |
C and statistics ; higher -> more writing |
189 |
|
|
190 |
COMMON /PARM_I/ |
COMMON /PARM_I/ |
191 |
& cg2dMaxIters, cg2dChkResFreq, cg2dPreCondFreq, |
& cg2dMaxIters, cg2dChkResFreq, |
192 |
|
& cg2dPreCondFreq, cg2dUseMinResSol, |
193 |
& cg3dMaxIters, cg3dChkResFreq, |
& cg3dMaxIters, cg3dChkResFreq, |
194 |
& printResidualFreq, |
& printResidualFreq, |
195 |
& nIter0, nTimeSteps, nEndIter, |
& nIter0, nTimeSteps, nEndIter, |
208 |
INTEGER cg2dMaxIters |
INTEGER cg2dMaxIters |
209 |
INTEGER cg2dChkResFreq |
INTEGER cg2dChkResFreq |
210 |
INTEGER cg2dPreCondFreq |
INTEGER cg2dPreCondFreq |
211 |
|
INTEGER cg2dUseMinResSol |
212 |
INTEGER cg3dMaxIters |
INTEGER cg3dMaxIters |
213 |
INTEGER cg3dChkResFreq |
INTEGER cg3dChkResFreq |
214 |
INTEGER printResidualFreq |
INTEGER printResidualFreq |
259 |
C- Momentum params: |
C- Momentum params: |
260 |
C no_slip_sides :: Impose "no-slip" at lateral boundaries. |
C no_slip_sides :: Impose "no-slip" at lateral boundaries. |
261 |
C no_slip_bottom :: Impose "no-slip" at bottom boundary. |
C no_slip_bottom :: Impose "no-slip" at bottom boundary. |
262 |
|
C useSmag3D :: Use isotropic 3-D Smagorinsky |
263 |
C useFullLeith :: Set to true to use full Leith viscosity(may be unstable |
C useFullLeith :: Set to true to use full Leith viscosity(may be unstable |
264 |
C on irregular grids) |
C on irregular grids) |
265 |
C useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms |
C useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms |
374 |
& usingCartesianGrid, usingSphericalPolarGrid, rotateGrid, |
& usingCartesianGrid, usingSphericalPolarGrid, rotateGrid, |
375 |
& usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners, |
& usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners, |
376 |
& deepAtmosphere, setInterFDr, setCenterDr, |
& deepAtmosphere, setInterFDr, setCenterDr, |
377 |
& no_slip_sides, no_slip_bottom, |
& no_slip_sides, no_slip_bottom, useSmag3D, |
378 |
& useFullLeith, useStrainTensionVisc, useAreaViscLength, |
& useFullLeith, useStrainTensionVisc, useAreaViscLength, |
379 |
& momViscosity, momAdvection, momForcing, |
& momViscosity, momAdvection, momForcing, |
380 |
& momPressureForcing, metricTerms, useNHMTerms, |
& momPressureForcing, metricTerms, useNHMTerms, |
426 |
|
|
427 |
LOGICAL no_slip_sides |
LOGICAL no_slip_sides |
428 |
LOGICAL no_slip_bottom |
LOGICAL no_slip_bottom |
429 |
|
LOGICAL useSmag3D |
430 |
LOGICAL useFullLeith |
LOGICAL useFullLeith |
431 |
LOGICAL useStrainTensionVisc |
LOGICAL useStrainTensionVisc |
432 |
LOGICAL useAreaViscLength |
LOGICAL useAreaViscLength |
584 |
C momentum laterally ( m^2/s ) |
C momentum laterally ( m^2/s ) |
585 |
C viscAhW :: Eddy viscosity coeff. for mixing of vertical |
C viscAhW :: Eddy viscosity coeff. for mixing of vertical |
586 |
C momentum laterally, no effect for hydrostatic |
C momentum laterally, no effect for hydrostatic |
587 |
C model, defaults to viscAh if unset ( m^2/s ) |
C model, defaults to viscAhD if unset ( m^2/s ) |
588 |
C Not used if variable horiz. viscosity is used. |
C Not used if variable horiz. viscosity is used. |
589 |
C viscA4 :: Biharmonic viscosity coeff. for mixing of |
C viscA4 :: Biharmonic viscosity coeff. for mixing of |
590 |
C momentum laterally ( m^4/s ) |
C momentum laterally ( m^4/s ) |
591 |
C viscA4W :: Biharmonic viscosity coeff. for mixing of vertical |
C viscA4W :: Biharmonic viscosity coeff. for mixing of vertical |
592 |
C momentum laterally, no effect for hydrostatic |
C momentum laterally, no effect for hydrostatic |
593 |
C model, defaults to viscA4 if unset ( m^2/s ) |
C model, defaults to viscA4D if unset ( m^2/s ) |
594 |
C Not used if variable horiz. viscosity is used. |
C Not used if variable horiz. viscosity is used. |
595 |
C viscAhD :: Eddy viscosity coeff. for mixing of momentum laterally |
C viscAhD :: Eddy viscosity coeff. for mixing of momentum laterally |
596 |
C (act on Divergence part) ( m^2/s ) |
C (act on Divergence part) ( m^2/s ) |
600 |
C (act on Divergence part) ( m^4/s ) |
C (act on Divergence part) ( m^4/s ) |
601 |
C viscA4Z :: Biharmonic viscosity coeff. for mixing of momentum laterally |
C viscA4Z :: Biharmonic viscosity coeff. for mixing of momentum laterally |
602 |
C (act on Vorticity part) ( m^4/s ) |
C (act on Vorticity part) ( m^4/s ) |
603 |
|
C smag3D_coeff :: Isotropic 3-D Smagorinsky coefficient (-) |
604 |
C viscC2leith :: Leith non-dimensional viscosity factor (grad(vort)) |
C viscC2leith :: Leith non-dimensional viscosity factor (grad(vort)) |
605 |
C viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div)) |
C viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div)) |
606 |
C viscC4leith :: Leith non-dimensional viscosity factor (grad(vort)) |
C viscC4leith :: Leith non-dimensional viscosity factor (grad(vort)) |
747 |
& abEps, alph_AB, beta_AB, |
& abEps, alph_AB, beta_AB, |
748 |
& rSphere, recip_rSphere, radius_fromHorizGrid, |
& rSphere, recip_rSphere, radius_fromHorizGrid, |
749 |
& f0, beta, fPrime, omega, rotationPeriod, |
& f0, beta, fPrime, omega, rotationPeriod, |
750 |
& viscFacAdj, viscAh, viscAhW, viscAhMax, |
& viscFacAdj, viscAh, viscAhW, smag3D_coeff, |
751 |
& viscAhGrid, viscAhGridMax, viscAhGridMin, |
& viscAhMax, viscAhGrid, viscAhGridMax, viscAhGridMin, |
752 |
& viscC2leith, viscC2leithD, |
& viscC2leith, viscC2leithD, |
753 |
& viscC2smag, viscC4smag, |
& viscC2smag, viscC4smag, |
754 |
& viscAhD, viscAhZ, viscA4D, viscA4Z, |
& viscAhD, viscAhZ, viscA4D, viscA4Z, |
823 |
_RL viscAhW |
_RL viscAhW |
824 |
_RL viscAhD |
_RL viscAhD |
825 |
_RL viscAhZ |
_RL viscAhZ |
826 |
|
_RL smag3D_coeff |
827 |
_RL viscAhMax |
_RL viscAhMax |
828 |
_RL viscAhReMax |
_RL viscAhReMax |
829 |
_RL viscAhGrid, viscAhGridMax, viscAhGridMin |
_RL viscAhGrid, viscAhGridMax, viscAhGridMin |
978 |
LOGICAL useStreamIce |
LOGICAL useStreamIce |
979 |
LOGICAL useICEFRONT |
LOGICAL useICEFRONT |
980 |
LOGICAL useThSIce |
LOGICAL useThSIce |
981 |
|
LOGICAL useLand |
982 |
LOGICAL useATM2d |
LOGICAL useATM2d |
983 |
LOGICAL useAIM |
LOGICAL useAIM |
984 |
LOGICAL useLand |
LOGICAL useAtm_Phys |
985 |
LOGICAL useFizhi |
LOGICAL useFizhi |
986 |
LOGICAL useGridAlt |
LOGICAL useGridAlt |
987 |
LOGICAL useDiagnostics |
LOGICAL useDiagnostics |
999 |
& useGrdchk,useSMOOTH,usePROFILES,useECCO,useSBO, useFLT, |
& useGrdchk,useSMOOTH,usePROFILES,useECCO,useSBO, useFLT, |
1000 |
& usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX, |
& usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX, |
1001 |
& useFRAZIL, useSEAICE, useSALT_PLUME, useShelfIce, |
& useFRAZIL, useSEAICE, useSALT_PLUME, useShelfIce, |
1002 |
& useStreamIce, useICEFRONT, useThSIce, |
& useStreamIce, useICEFRONT, useThSIce, useLand, |
1003 |
& useATM2D, useAIM, useLand, useFizhi, useGridAlt, |
& useATM2D, useAIM, useAtm_Phys, useFizhi, useGridAlt, |
1004 |
& useDiagnostics, useREGRID, useLayers, useMNC, |
& useDiagnostics, useREGRID, useLayers, useMNC, |
1005 |
& useRunClock, useEMBED_FILES, |
& useRunClock, useEMBED_FILES, |
1006 |
& useMYPACKAGE |
& useMYPACKAGE |