| 1 |
torge |
1.6 |
C $Header: /u/gcmpack/MITgcm/pkg/seaice/SEAICE_PARAMS.h,v 1.109 2013/02/28 16:21:40 mlosch Exp $ |
| 2 |
dimitri |
1.1 |
C $Name: $ |
| 3 |
|
|
|
| 4 |
|
|
C *==========================================================* |
| 5 |
|
|
C | SEAICE_PARAMS.h |
| 6 |
|
|
C | o Basic parameter header for sea ice model. |
| 7 |
|
|
C *==========================================================* |
| 8 |
|
|
|
| 9 |
|
|
C-- COMMON /SEAICE_PARM_L/ Logical parameters of sea ice model. |
| 10 |
torge |
1.6 |
C - dynamics: |
| 11 |
dimitri |
1.1 |
C SEAICEuseDYNAMICS :: If false, do not use dynamics; |
| 12 |
|
|
C default is to use dynamics. |
| 13 |
|
|
C SEAICEuseFREEDRIFT :: If True use free drift velocity instead of EVP |
| 14 |
|
|
C or LSR |
| 15 |
torge |
1.6 |
C SEAICEuseEVP :: If true use elastic viscous plastic solver |
| 16 |
|
|
C SEAICEuseEVPstar :: If true use modified elastic viscous plastic |
| 17 |
|
|
C solver (EVP*) by Lemieux et al (2012) |
| 18 |
dimitri |
1.1 |
C SEAICEuseEVPpickup :: Set to false in order to start EVP solver with |
| 19 |
|
|
C non-EVP pickup files. Default is true. |
| 20 |
|
|
C Applied only if SEAICEuseEVP=.TRUE. |
| 21 |
torge |
1.6 |
C SEAICEuseMultiTileSolver :: in LSR, use full domain tri-diagonal solver |
| 22 |
|
|
C SEAICEuseJFNK :: If true, use Jacobi-free Newton-Krylov solver |
| 23 |
|
|
C instead of LSR (default: false) |
| 24 |
|
|
C SEAICEuseTEM :: to use truncated ellipse method (see Geiger et al. |
| 25 |
|
|
C 1998) set this parameter to true, default is false |
| 26 |
|
|
C SEAICEuseTilt :: If true then include surface tilt term in dynamics |
| 27 |
|
|
C SEAICEuseMetricTerms :: use metric terms for dynamics solver |
| 28 |
|
|
C (default = .true. ) |
| 29 |
|
|
C SEAICE_no_slip :: apply no slip boundary conditions to seaice velocity |
| 30 |
|
|
C SEAICE_maskRHS :: mask the RHS of the solver where there is no ice |
| 31 |
|
|
C SEAICE_clipVelocities :: clip velocities to +/- 40cm/s |
| 32 |
|
|
C useHB87stressCoupling :: use an intergral over ice and ocean surface |
| 33 |
|
|
C layer to define surface stresses on ocean |
| 34 |
|
|
C following Hibler and Bryan (1987, JPO) |
| 35 |
|
|
C - advection: |
| 36 |
dimitri |
1.1 |
C SEAICEuseFluxForm :: use flux form for advection and diffusion |
| 37 |
|
|
C of seaice |
| 38 |
|
|
C SEAICEadvHeff :: turn on advection of effective thickness |
| 39 |
|
|
C (default = .true.) |
| 40 |
|
|
C SEAICEadvArea :: turn on advection of fraction area |
| 41 |
|
|
C (default = .true.) |
| 42 |
|
|
C SEAICEadvSnow :: turn on advection of snow (does not work with |
| 43 |
|
|
C non-default Leap-frog scheme for advection) |
| 44 |
|
|
C SEAICEadvSalt :: turn on advection of salt (does not work with |
| 45 |
|
|
C non-default Leap-frog scheme for advection) |
| 46 |
torge |
1.6 |
C - thermodynamics: |
| 47 |
dimitri |
1.1 |
C usePW79thermodynamics :: use "0-layer" thermodynamics as described in |
| 48 |
|
|
C Parkinson and Washington (1979) and Hibler (1979) |
| 49 |
torge |
1.6 |
C SEAICE_useMultDimSnow :: use same fixed pdf for snow as for |
| 50 |
|
|
C MULITCATEGORY ice |
| 51 |
|
|
C SEAICEuseFlooding :: turn on scheme to convert submerged snow into ice |
| 52 |
|
|
C SEAICEheatConsFix :: If true then fix ocn<->seaice advective heat flux. |
| 53 |
dimitri |
1.1 |
C useMaykutSatVapPoly :: use Maykut Polynomial for saturation vapor pressure |
| 54 |
|
|
C instead of extended temp-range exponential law; def=F. |
| 55 |
torge |
1.6 |
C SEAICE_mcPheeStepFunc :: use step function (not linear tapering) in |
| 56 |
dimitri |
1.1 |
C ocean-ice turbulent flux |
| 57 |
|
|
C SEAICE_doOpenWaterGrowth :: use open water heat flux directly to grow ice |
| 58 |
|
|
C (when false cool ocean, and grow later if needed) |
| 59 |
|
|
C SEAICE_doOpenWaterMelt :: use open water heat flux directly to melt ice |
| 60 |
|
|
C (when false warm ocean, and melt later if needed) |
| 61 |
|
|
C SEAICE_salinityTracer :: use SItracer to exchange and trace ocean |
| 62 |
|
|
C salt in ice |
| 63 |
torge |
1.6 |
C SEAICE_ageTracer :: use SItracer to trace the age of ice |
| 64 |
dimitri |
1.1 |
C SEAICErestoreUnderIce :: restore surface T/S also underneath ice |
| 65 |
|
|
C ( default is false ) |
| 66 |
torge |
1.6 |
C - other (I/O, ...): |
| 67 |
|
|
C SEAICEwriteState :: If true, write sea ice state to file; |
| 68 |
|
|
C default is false. |
| 69 |
dimitri |
1.1 |
C SEAICE_tave_mdsio :: write TimeAverage output using MDSIO |
| 70 |
|
|
C SEAICE_dump_mdsio :: write snap-shot output using MDSIO |
| 71 |
|
|
C SEAICE_mon_stdio :: write monitor to std-outp |
| 72 |
|
|
C SEAICE_tave_mnc :: write TimeAverage output using MNC |
| 73 |
|
|
C SEAICE_dump_mnc :: write snap-shot output using MNC |
| 74 |
|
|
C SEAICE_mon_mnc :: write monitor to netcdf file |
| 75 |
|
|
LOGICAL |
| 76 |
torge |
1.6 |
& SEAICEuseDYNAMICS, SEAICEuseFREEDRIFT, |
| 77 |
|
|
& SEAICEuseEVP, SEAICEuseEVPstar, SEAICEuseEVPpickup, |
| 78 |
|
|
& SEAICEuseMultiTileSolver, |
| 79 |
torge |
1.3 |
& SEAICEuseJFNK, |
| 80 |
torge |
1.6 |
& SEAICEuseTEM, SEAICEuseTilt, SEAICEuseMetricTerms, |
| 81 |
|
|
& SEAICE_no_slip, SEAICE_maskRHS, |
| 82 |
|
|
& SEAICE_clipVelocities, useHB87stressCoupling, |
| 83 |
|
|
& SEAICEuseFluxForm, SEAICEadvHeff, SEAICEadvArea, |
| 84 |
dimitri |
1.1 |
& SEAICEadvSnow, SEAICEadvSalt, |
| 85 |
torge |
1.6 |
& usePW79thermodynamics, |
| 86 |
|
|
& SEAICE_useMultDimSnow, SEAICEuseFlooding, SEAICEheatConsFix, |
| 87 |
|
|
& useMaykutSatVapPoly, SEAICE_mcPheeStepFunc, |
| 88 |
dimitri |
1.1 |
& SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt, |
| 89 |
|
|
& SEAICE_salinityTracer, SEAICE_ageTracer, |
| 90 |
|
|
& SEAICErestoreUnderIce, |
| 91 |
torge |
1.6 |
& SEAICEwriteState, |
| 92 |
dimitri |
1.1 |
& SEAICE_tave_mdsio, SEAICE_dump_mdsio, SEAICE_mon_stdio, |
| 93 |
|
|
& SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc |
| 94 |
|
|
COMMON /SEAICE_PARM_L/ |
| 95 |
torge |
1.6 |
& SEAICEuseDYNAMICS, SEAICEuseFREEDRIFT, |
| 96 |
|
|
& SEAICEuseEVP, SEAICEuseEVPstar, SEAICEuseEVPpickup, |
| 97 |
|
|
& SEAICEuseMultiTileSolver, |
| 98 |
torge |
1.3 |
& SEAICEuseJFNK, |
| 99 |
torge |
1.6 |
& SEAICEuseTEM, SEAICEuseTilt, SEAICEuseMetricTerms, |
| 100 |
|
|
& SEAICE_no_slip, SEAICE_maskRHS, |
| 101 |
|
|
& SEAICE_clipVelocities, useHB87stressCoupling, |
| 102 |
|
|
& SEAICEuseFluxForm, SEAICEadvHeff, SEAICEadvArea, |
| 103 |
dimitri |
1.1 |
& SEAICEadvSnow, SEAICEadvSalt, |
| 104 |
torge |
1.6 |
& usePW79thermodynamics, |
| 105 |
|
|
& SEAICE_useMultDimSnow, SEAICEuseFlooding, SEAICEheatConsFix, |
| 106 |
|
|
& useMaykutSatVapPoly, SEAICE_mcPheeStepFunc, |
| 107 |
dimitri |
1.1 |
& SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt, |
| 108 |
|
|
& SEAICE_salinityTracer, SEAICE_ageTracer, |
| 109 |
|
|
& SEAICErestoreUnderIce, |
| 110 |
torge |
1.6 |
& SEAICEwriteState, |
| 111 |
dimitri |
1.1 |
& SEAICE_tave_mdsio, SEAICE_dump_mdsio, SEAICE_mon_stdio, |
| 112 |
|
|
& SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc |
| 113 |
|
|
|
| 114 |
|
|
C-- COMMON /SEAICE_PARM_I/ Integer valued parameters of sea ice model. |
| 115 |
torge |
1.6 |
C IMAX_TICE :: number of iterations for ice surface temp |
| 116 |
|
|
C (default=10) |
| 117 |
|
|
C postSolvTempIter :: select flux calculation after surf. temp solver |
| 118 |
|
|
C iteration |
| 119 |
|
|
C 0 = none, i.e., from last iter |
| 120 |
|
|
C 1 = use linearized approx (consistent with tsurf |
| 121 |
|
|
C finding) |
| 122 |
|
|
C 2 = full non-lin form |
| 123 |
dimitri |
1.1 |
C SOLV_MAX_ITERS :: maximum number of allowed LSR-solver iterations |
| 124 |
|
|
C SOLV_NCHECK :: iteration interval for solver convergence test |
| 125 |
|
|
C NPSEUDOTIMESTEPS :: number of extra pseudo time steps (>= 2) |
| 126 |
torge |
1.6 |
C SEAICEnEVPstarSteps :: number of evp*-steps |
| 127 |
|
|
C SEAICEnewtonIterMax :: maximum number of allowed Newton iterations |
| 128 |
torge |
1.3 |
C in JFNK-solver |
| 129 |
|
|
C SEAICEkrylovIterMax :: maximum number of allowed Krylov iterations |
| 130 |
|
|
C in JFNK-solver |
| 131 |
torge |
1.6 |
C SEAICE_JFNK_lsIter :: number of Newton iterations after which the |
| 132 |
|
|
C line search is started |
| 133 |
|
|
C SEAICE_JFNK_tolIter :: number of Newton iterations after which the |
| 134 |
|
|
C the tolerance is relaxed again (default = 100) |
| 135 |
|
|
C SEAICE_OLx/y :: overlaps for LSR-preconditioner in JFNK solver; |
| 136 |
|
|
C for 0 < SEAICE_OLx/y 0 <= OLx/y-2 |
| 137 |
|
|
C the preconditioner is a restricted additive |
| 138 |
|
|
C Schwarz method (default = OLx/y-2). |
| 139 |
|
|
C LSR_mixIniGuess :: control mixing of free-drift sol. into LSR initial |
| 140 |
|
|
C guess |
| 141 |
|
|
C :: =0 : nothing; =1 : no mix, but print free-drift |
| 142 |
|
|
C resid.; |
| 143 |
|
|
C :: =2,4 : mix with (1/local-residual)^2,4 factor |
| 144 |
torge |
1.3 |
C SEAICEpresPow0 :: HEFF exponent for ice strength below SEAICEpresH0 |
| 145 |
|
|
C SEAICEpresPow1 :: HEFF exponent for ice strength above SEAICEpresH0 |
| 146 |
dimitri |
1.1 |
C SEAICEadvScheme :: sets the advection scheme for thickness and area |
| 147 |
|
|
C SEAICEadvSchArea :: sets the advection scheme for area |
| 148 |
|
|
C SEAICEadvSchHeff :: sets the advection scheme for effective thickness |
| 149 |
|
|
C (=volume), snow thickness, and salt if available |
| 150 |
|
|
C SEAICEadvSchSnow :: sets the advection scheme for snow on sea-ice |
| 151 |
|
|
C SEAICEadvSchSalt :: sets the advection scheme for sea ice salinity |
| 152 |
|
|
C SEAICEadvSchSnow :: sets the advection scheme for snow on sea-ice |
| 153 |
|
|
C SEAICE_areaLossFormula :: selects formula for ice cover loss from melt |
| 154 |
|
|
C :: 1=from all but only melt conributions by ATM and OCN |
| 155 |
|
|
C :: 2=from net melt-growth>0 by ATM and OCN |
| 156 |
|
|
C :: 3=from predicted melt by ATM |
| 157 |
torge |
1.6 |
C SEAICE_areaGainFormula :: selects formula for ice cover gain from open |
| 158 |
|
|
C water growth |
| 159 |
dimitri |
1.1 |
C :: 1=from growth by ATM |
| 160 |
|
|
C :: 2=from predicted growth by ATM |
| 161 |
torge |
1.5 |
C SEAICEetaZmethod :: determines how shear-viscosity eta is computed at |
| 162 |
|
|
C Z-points |
| 163 |
|
|
C 0=simple averaging from C-points (default and old) |
| 164 |
|
|
C 1=linear averaging of strain rates to Z-points |
| 165 |
|
|
C 2=averaging of squares of strain rates |
| 166 |
dimitri |
1.1 |
C SEAICE_multDim :: number of ice categories |
| 167 |
|
|
C SEAICE_debugPointI :: I,J index for seaice-specific debuggin |
| 168 |
|
|
C SEAICE_debugPointJ |
| 169 |
|
|
C |
| 170 |
|
|
INTEGER IMAX_TICE, postSolvTempIter |
| 171 |
|
|
INTEGER SOLV_MAX_ITERS, SOLV_NCHECK |
| 172 |
|
|
INTEGER NPSEUDOTIMESTEPS |
| 173 |
|
|
INTEGER LSR_mixIniGuess |
| 174 |
torge |
1.6 |
INTEGER SEAICEnEVPstarSteps |
| 175 |
torge |
1.3 |
INTEGER SEAICEnewtonIterMax, SEAICEkrylovIterMax |
| 176 |
torge |
1.6 |
INTEGER SEAICE_JFNK_lsIter, SEAICE_JFNK_tolIter |
| 177 |
|
|
INTEGER SEAICE_OLx, SEAICE_OLy |
| 178 |
dimitri |
1.1 |
INTEGER SEAICEadvScheme |
| 179 |
|
|
INTEGER SEAICEadvSchArea |
| 180 |
|
|
INTEGER SEAICEadvSchHeff |
| 181 |
|
|
INTEGER SEAICEadvSchSnow |
| 182 |
|
|
INTEGER SEAICEadvSchSalt |
| 183 |
|
|
INTEGER SEAICEadjMODE |
| 184 |
|
|
INTEGER SEAICE_areaLossFormula |
| 185 |
|
|
INTEGER SEAICE_areaGainFormula |
| 186 |
torge |
1.5 |
INTEGER SEAICEetaZmethod |
| 187 |
dimitri |
1.1 |
INTEGER SEAICE_multDim |
| 188 |
|
|
INTEGER SEAICE_debugPointI |
| 189 |
|
|
INTEGER SEAICE_debugPointJ |
| 190 |
torge |
1.3 |
INTEGER SEAICEpresPow0, SEAICEpresPow1 |
| 191 |
dimitri |
1.1 |
COMMON /SEAICE_PARM_I/ |
| 192 |
|
|
& IMAX_TICE, postSolvTempIter, |
| 193 |
|
|
& SOLV_MAX_ITERS, SOLV_NCHECK, |
| 194 |
|
|
& NPSEUDOTIMESTEPS, |
| 195 |
|
|
& LSR_mixIniGuess, |
| 196 |
torge |
1.6 |
& SEAICEnEVPstarSteps, |
| 197 |
torge |
1.3 |
& SEAICEnewtonIterMax, SEAICEkrylovIterMax, |
| 198 |
torge |
1.6 |
& SEAICE_JFNK_lsIter, SEAICE_OLx, SEAICE_OLy, |
| 199 |
|
|
& SEAICE_JFNK_tolIter, |
| 200 |
torge |
1.3 |
& SEAICEpresPow0, SEAICEpresPow1, |
| 201 |
dimitri |
1.1 |
& SEAICEadvScheme, |
| 202 |
|
|
& SEAICEadvSchArea, |
| 203 |
|
|
& SEAICEadvSchHeff, |
| 204 |
|
|
& SEAICEadvSchSnow, |
| 205 |
|
|
& SEAICEadvSchSalt, |
| 206 |
|
|
& SEAICEadjMODE, |
| 207 |
|
|
& SEAICE_areaLossFormula, |
| 208 |
|
|
& SEAICE_areaGainFormula, |
| 209 |
|
|
& SEAICE_multDim, |
| 210 |
torge |
1.5 |
& SEAICEetaZmethod, |
| 211 |
dimitri |
1.1 |
& SEAICE_debugPointI, |
| 212 |
|
|
& SEAICE_debugPointJ |
| 213 |
|
|
|
| 214 |
|
|
C-- COMMON /SEAICE_PARM_C/ Character valued sea ice model parameters. |
| 215 |
|
|
C AreaFile :: File containing initial sea-ice concentration |
| 216 |
|
|
C HsnowFile :: File containing initial snow thickness |
| 217 |
|
|
C HsaltFile :: File containing initial sea ice salt content |
| 218 |
|
|
C HeffFile :: File containing initial sea-ice thickness |
| 219 |
|
|
C uIceFile :: File containing initial sea-ice U comp. velocity |
| 220 |
|
|
C vIceFile :: File containing initial sea-ice V comp. velocity |
| 221 |
|
|
C !!! NOTE !!! Initial sea-ice thickness can also be set using |
| 222 |
|
|
C SEAICE_initialHEFF below. But a constant initial condition |
| 223 |
|
|
C can mean large artificial fluxes of heat and freshwater in |
| 224 |
|
|
C the surface layer during the first model time step. |
| 225 |
|
|
C |
| 226 |
|
|
CHARACTER*(MAX_LEN_FNAM) AreaFile |
| 227 |
|
|
CHARACTER*(MAX_LEN_FNAM) HsnowFile |
| 228 |
|
|
CHARACTER*(MAX_LEN_FNAM) HsaltFile |
| 229 |
|
|
CHARACTER*(MAX_LEN_FNAM) HeffFile |
| 230 |
|
|
CHARACTER*(MAX_LEN_FNAM) uIceFile |
| 231 |
|
|
CHARACTER*(MAX_LEN_FNAM) vIceFile |
| 232 |
|
|
COMMON /SEAICE_PARM_C/ |
| 233 |
|
|
& AreaFile, HsnowFile, HsaltFile, HeffFile, |
| 234 |
|
|
& uIceFile, vIceFile |
| 235 |
|
|
|
| 236 |
|
|
C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model. |
| 237 |
|
|
C SEAICE_deltaTtherm :: Seaice timestep for thermodynamic equations (s) |
| 238 |
|
|
C SEAICE_deltaTdyn :: Seaice timestep for dynamic solver (s) |
| 239 |
torge |
1.5 |
C SEAICE_LSRrelaxU/V :: relaxation parameter for LSR-solver: U/V-component |
| 240 |
dimitri |
1.1 |
C SEAICE_deltaTevp :: Seaice timestep for EVP solver (s) |
| 241 |
|
|
C SEAICE_elasticParm :: parameter that sets relaxation timescale |
| 242 |
|
|
C tau = SEAICE_elasticParm * SEAICE_deltaTdyn |
| 243 |
|
|
C SEAICE_evpTauRelax :: relaxation timescale tau (s) |
| 244 |
|
|
C SEAICE_evpDampC :: evp damping constant (Hunke,JCP,2001) (kg/m^2) |
| 245 |
torge |
1.6 |
C SEAICE_evpAlpha :: dimensionless parameter 2*evpTauRelax/deltaTevp |
| 246 |
|
|
C SEAICE_evpBeta :: dimensionless parameter deltaTdyn/deltaTevp |
| 247 |
torge |
1.3 |
C JFNKgamma_nonlin :: non-linear tolerance parameter for JFNK solver |
| 248 |
|
|
C JFNKgamma_lin_min/max :: tolerance parameters for linear JFNK solver |
| 249 |
|
|
C JFNKres_t :: tolerance parameter for FGMRES residual |
| 250 |
torge |
1.6 |
C JFNKres_tFac :: if set, JFNKres_t=JFNKres_tFac*(initial residual) |
| 251 |
|
|
C SEAICE_JFNKepsilon :: step size for the FD-gradient in s/r seaice_jacvec |
| 252 |
dimitri |
1.1 |
C SEAICE_zetaMaxFac :: factor determining the maximum viscosity (s) |
| 253 |
|
|
C (default = 5.e+12/2.e4 = 2.5e8) |
| 254 |
|
|
C SEAICE_zetaMin :: lower bound for viscosity (default = 0) (N s/m^2) |
| 255 |
torge |
1.3 |
C SEAICEpresH0 :: HEFF threshold for ice strength (m) |
| 256 |
dimitri |
1.1 |
C SEAICE_monFreq :: SEAICE monitor frequency. (s) |
| 257 |
|
|
C SEAICE_dumpFreq :: SEAICE dump frequency. (s) |
| 258 |
|
|
C SEAICE_taveFreq :: SEAICE time-averaging frequency. (s) |
| 259 |
|
|
C SEAICE_initialHEFF :: initial sea-ice thickness (m) |
| 260 |
|
|
C SEAICE_rhoAir :: density of air (kg/m^3) |
| 261 |
|
|
C SEAICE_rhoIce :: density of sea ice (kg/m^3) |
| 262 |
|
|
C SEAICE_rhoSnow :: density of snow (kg/m^3) |
| 263 |
|
|
C ICE2WATR :: ratio of sea ice density to water density |
| 264 |
|
|
C OCEAN_drag :: air-ocean drag coefficient |
| 265 |
|
|
C SEAICE_cpAir :: specific heat of air (J/kg/K) |
| 266 |
|
|
C |
| 267 |
|
|
C SEAICE_drag :: air-ice drag coefficient |
| 268 |
|
|
C SEAICE_waterDrag :: water-ice drag coefficient * water density |
| 269 |
|
|
C SEAICE_dryIceAlb :: winter albedo |
| 270 |
|
|
C SEAICE_wetIceAlb :: summer albedo |
| 271 |
|
|
C SEAICE_drySnowAlb :: dry snow albedo |
| 272 |
|
|
C SEAICE_wetSnowAlb :: wet snow albedo |
| 273 |
|
|
C HO :: AKA "lead closing parameter", demarcation thickness |
| 274 |
|
|
C between thin and thick ice. Alternatively, HO (in |
| 275 |
|
|
C meters) can be interpreted as the thickness of ice |
| 276 |
|
|
C formed in open water. |
| 277 |
|
|
C HO is a key ice-growth parameter that determines |
| 278 |
|
|
C the partition between vertical and lateral growth. |
| 279 |
|
|
C The default is 0.5m, increasing this value leads |
| 280 |
|
|
C slower formation of a closed ice cover and thus to |
| 281 |
|
|
C more ice (and thicker) ice, decreasing to faster |
| 282 |
|
|
C formation of a closed ice cover (leads are closing |
| 283 |
|
|
C faster) and thus less (thinner) ice. |
| 284 |
|
|
C |
| 285 |
|
|
C SEAICE_drag_south :: Southern Ocean SEAICE_drag |
| 286 |
|
|
C SEAICE_waterDrag_south :: Southern Ocean SEAICE_waterDrag |
| 287 |
|
|
C SEAICE_dryIceAlb_south :: Southern Ocean SEAICE_dryIceAlb |
| 288 |
|
|
C SEAICE_wetIceAlb_south :: Southern Ocean SEAICE_wetIceAlb |
| 289 |
|
|
C SEAICE_drySnowAlb_south :: Southern Ocean SEAICE_drySnowAlb |
| 290 |
|
|
C SEAICE_wetSnowAlb_south :: Southern Ocean SEAICE_wetSnowAlb |
| 291 |
|
|
C HO_south :: Southern Ocean HO |
| 292 |
|
|
C |
| 293 |
|
|
C SEAICE_wetAlbTemp :: Temp (deg.C) above which wet-albedo values are used |
| 294 |
|
|
C SEAICE_waterAlbedo :: water albedo |
| 295 |
|
|
C SEAICE_strength :: sea-ice strength Pstar |
| 296 |
|
|
C SEAICE_eccen :: sea-ice eccentricity of the elliptical yield curve |
| 297 |
|
|
C SEAICE_lhFusion :: latent heat of fusion for ice and snow (J/kg) |
| 298 |
|
|
C SEAICE_lhEvap :: latent heat of evaporation for water (J/kg) |
| 299 |
|
|
C SEAICE_dalton :: Dalton number (= sensible heat transfer coefficient) |
| 300 |
|
|
C SEAICE_iceConduct :: sea-ice conductivity |
| 301 |
|
|
C SEAICE_snowConduct :: snow conductivity |
| 302 |
|
|
C SEAICE_emissivity :: longwave ocean-surface emissivity (-) |
| 303 |
|
|
C SEAICE_ice_emiss :: longwave ice-surface emissivity (-) |
| 304 |
|
|
C SEAICE_snow_emiss :: longwave snow-surface emissivity (-) |
| 305 |
|
|
C SEAICE_boltzmann :: Stefan-Boltzman constant (not a run time parameter) |
| 306 |
|
|
C SEAICE_snowThick :: cutoff snow thickness (for snow-albedo) |
| 307 |
|
|
C SEAICE_shortwave :: ice penetration shortwave radiation factor |
| 308 |
|
|
C SEAICE_saltFrac :: salinity of newly formed seaice defined as a |
| 309 |
|
|
C fraction of the ocean surface salinity at the time |
| 310 |
|
|
C of freezing |
| 311 |
|
|
C SEAICE_salt0 :: prescribed salinity of seaice (in g/kg). |
| 312 |
|
|
C facOpenGrow :: 0./1. version of logical SEAICE_doOpenWaterGrowth |
| 313 |
|
|
C facOpenMelt :: 0./1. version of logical SEAICE_doOpenWaterMelt |
| 314 |
|
|
C SEAICE_mcPheePiston:: ocean-ice turbulent flux "piston velocity" (m/s) |
| 315 |
|
|
C that sets melt efficiency. |
| 316 |
|
|
C SEAICE_mcPheeTaper :: tapering down of turbulent flux term with ice |
| 317 |
|
|
C concentration. The 100% cover turb. flux is |
| 318 |
|
|
C multiplied by 1.-SEAICE_mcPheeTaper |
| 319 |
|
|
C SEAICE_frazilFrac :: Fraction of surface level negative heat content |
| 320 |
|
|
C anomalies (relative to the local freezing point) |
| 321 |
|
|
C may contribute as frazil over one time step. |
| 322 |
|
|
C SEAICE_tempFrz0 :: sea water freezing point is |
| 323 |
|
|
C SEAICE_dTempFrz_dS :: tempFrz = SEAICE_tempFrz0 + salt*SEAICE_dTempFrz_dS |
| 324 |
|
|
C SEAICEstressFactor :: factor by which ice affects wind stress (default=1) |
| 325 |
|
|
C LSR_ERROR :: sets accuracy of LSR solver |
| 326 |
|
|
C DIFF1 :: parameter used in advect.F |
| 327 |
|
|
C SEAICE_area_max :: usually set to 1. Seeting areaMax below 1 specifies |
| 328 |
|
|
C the minimun amount of leads (1-areaMax) in the |
| 329 |
|
|
C ice pack. |
| 330 |
|
|
C SEAICE_area_floor :: usually set to 1x10^-5. Specifies a minimun |
| 331 |
|
|
C ice fraction in the ice pack. |
| 332 |
|
|
C SEAICE_area_reg :: usually set to 1x10^-5. Specifies a minimun |
| 333 |
|
|
C ice fraction for the purposes of regularization |
| 334 |
|
|
C SEAICE_hice_reg :: usually set to 5 cm. Specifies a minimun |
| 335 |
|
|
C ice thickness for the purposes of regularization |
| 336 |
|
|
C SEAICEdiffKhArea :: sets the diffusivity for area (m^2/s) |
| 337 |
|
|
C SEAICEdiffKhHeff :: sets the diffusivity for effective thickness (m^2/s) |
| 338 |
|
|
C SEAICEdiffKhSnow :: sets the diffusivity for snow on sea-ice (m^2/s) |
| 339 |
|
|
C SEAICEdiffKhSalt :: sets the diffusivity for sea ice salinity (m^2/s) |
| 340 |
|
|
C SEAICE_airTurnAngle :: turning angles of air-ice interfacial stress |
| 341 |
|
|
C SEAICE_waterTurnAngle :: and ice-water interfacial stress (in degrees) |
| 342 |
torge |
1.6 |
C SEAICE_tauAreaObsRelax :: Timescale of relaxation to observed |
| 343 |
dimitri |
1.2 |
C sea ice concentration (s), default=unset |
| 344 |
dimitri |
1.1 |
C |
| 345 |
|
|
_RL SEAICE_deltaTtherm, SEAICE_deltaTdyn, SEAICE_deltaTevp |
| 346 |
torge |
1.5 |
_RL SEAICE_LSRrelaxU, SEAICE_LSRrelaxV |
| 347 |
dimitri |
1.1 |
_RL SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq |
| 348 |
|
|
_RL SEAICE_initialHEFF |
| 349 |
|
|
_RL SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR |
| 350 |
|
|
_RL SEAICE_cpAir |
| 351 |
|
|
_RL SEAICE_drag, SEAICE_waterDrag, SEAICE_dryIceAlb |
| 352 |
|
|
_RL SEAICE_wetIceAlb, SEAICE_drySnowAlb, SEAICE_wetSnowAlb, HO |
| 353 |
|
|
_RL SEAICE_drag_south, SEAICE_waterDrag_south |
| 354 |
|
|
_RL SEAICE_dryIceAlb_south, SEAICE_wetIceAlb_south |
| 355 |
|
|
_RL SEAICE_drySnowAlb_south, SEAICE_wetSnowAlb_south, HO_south |
| 356 |
|
|
_RL SEAICE_wetAlbTemp, SEAICE_waterAlbedo |
| 357 |
|
|
_RL SEAICE_strength, SEAICE_eccen |
| 358 |
|
|
_RL SEAICE_lhFusion, SEAICE_lhEvap |
| 359 |
|
|
_RL SEAICE_dalton |
| 360 |
|
|
_RL SEAICE_iceConduct, SEAICE_snowConduct |
| 361 |
|
|
_RL SEAICE_emissivity, SEAICE_ice_emiss, SEAICE_snow_emiss |
| 362 |
|
|
_RL SEAICE_boltzmann |
| 363 |
|
|
_RL SEAICE_snowThick, SEAICE_shortwave |
| 364 |
|
|
_RL SEAICE_saltFrac, SEAICE_salt0, SEAICEstressFactor |
| 365 |
|
|
_RL SEAICE_mcPheeTaper, SEAICE_mcPheePiston |
| 366 |
|
|
_RL SEAICE_frazilFrac, SEAICE_availHeatFrac |
| 367 |
|
|
_RL facOpenGrow, facOpenMelt |
| 368 |
|
|
_RL SEAICE_tempFrz0, SEAICE_dTempFrz_dS |
| 369 |
|
|
_RL OCEAN_drag, LSR_ERROR, DIFF1 |
| 370 |
torge |
1.6 |
_RL JFNKgamma_nonlin, JFNKres_t, JFNKres_tFac |
| 371 |
|
|
_RL JFNKgamma_lin_min, JFNKgamma_lin_max, SEAICE_JFNKepsilon |
| 372 |
dimitri |
1.1 |
_RL SEAICE_area_reg, SEAICE_hice_reg |
| 373 |
|
|
_RL SEAICE_area_floor, SEAICE_area_max |
| 374 |
|
|
_RL SEAICE_airTurnAngle, SEAICE_waterTurnAngle |
| 375 |
|
|
_RL SEAICE_elasticParm, SEAICE_evpTauRelax |
| 376 |
torge |
1.6 |
_RL SEAICE_evpAlpha, SEAICE_evpBeta |
| 377 |
dimitri |
1.1 |
_RL SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac |
| 378 |
torge |
1.3 |
_RL SEAICEpresH0 |
| 379 |
dimitri |
1.1 |
_RL SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow |
| 380 |
|
|
_RL SEAICEdiffKhSalt |
| 381 |
dimitri |
1.2 |
_RL SEAICE_tauAreaObsRelax |
| 382 |
dimitri |
1.1 |
|
| 383 |
|
|
COMMON /SEAICE_PARM_RL/ |
| 384 |
|
|
& SEAICE_deltaTtherm, SEAICE_deltaTdyn, |
| 385 |
torge |
1.5 |
& SEAICE_LSRrelaxU, SEAICE_LSRrelaxV, |
| 386 |
dimitri |
1.1 |
& SEAICE_deltaTevp, SEAICE_elasticParm, SEAICE_evpTauRelax, |
| 387 |
torge |
1.6 |
& SEAICE_evpAlpha, SEAICE_evpBeta, |
| 388 |
dimitri |
1.1 |
& SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac, |
| 389 |
torge |
1.3 |
& SEAICEpresH0, |
| 390 |
dimitri |
1.1 |
& SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq, |
| 391 |
|
|
& SEAICE_initialHEFF, |
| 392 |
|
|
& SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR, |
| 393 |
|
|
& SEAICE_drag, SEAICE_waterDrag, SEAICE_dryIceAlb, |
| 394 |
|
|
& SEAICE_wetIceAlb, SEAICE_drySnowAlb, SEAICE_wetSnowAlb, HO, |
| 395 |
|
|
& SEAICE_drag_south, SEAICE_waterDrag_south, |
| 396 |
|
|
& SEAICE_dryIceAlb_south, SEAICE_wetIceAlb_south, |
| 397 |
|
|
& SEAICE_drySnowAlb_south, SEAICE_wetSnowAlb_south, HO_south, |
| 398 |
|
|
& SEAICE_wetAlbTemp, SEAICE_waterAlbedo, |
| 399 |
|
|
& SEAICE_strength, SEAICE_eccen, |
| 400 |
|
|
& SEAICE_lhFusion, SEAICE_lhEvap, |
| 401 |
|
|
& SEAICE_dalton, SEAICE_cpAir, |
| 402 |
|
|
& SEAICE_iceConduct, SEAICE_snowConduct, |
| 403 |
|
|
& SEAICE_emissivity, SEAICE_ice_emiss, SEAICE_snow_emiss, |
| 404 |
|
|
& SEAICE_boltzmann, |
| 405 |
|
|
& SEAICE_snowThick, SEAICE_shortwave, |
| 406 |
|
|
& SEAICE_saltFrac, SEAICE_salt0, SEAICEstressFactor, |
| 407 |
|
|
& SEAICE_mcPheeTaper, SEAICE_mcPheePiston, |
| 408 |
|
|
& SEAICE_frazilFrac, SEAICE_availHeatFrac, |
| 409 |
|
|
& facOpenGrow, facOpenMelt, |
| 410 |
|
|
& SEAICE_tempFrz0, SEAICE_dTempFrz_dS, |
| 411 |
|
|
& OCEAN_drag, LSR_ERROR, DIFF1, |
| 412 |
torge |
1.6 |
& JFNKgamma_nonlin, JFNKres_t, JFNKres_tFac, |
| 413 |
|
|
& JFNKgamma_lin_min, JFNKgamma_lin_max, SEAICE_JFNKepsilon, |
| 414 |
dimitri |
1.1 |
& SEAICE_area_reg, SEAICE_hice_reg, |
| 415 |
|
|
& SEAICE_area_floor, SEAICE_area_max, |
| 416 |
|
|
& SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow, |
| 417 |
dimitri |
1.2 |
& SEAICEdiffKhSalt, SEAICE_tauAreaObsRelax, |
| 418 |
dimitri |
1.1 |
& SEAICE_airTurnAngle, SEAICE_waterTurnAngle |
| 419 |
|
|
|
| 420 |
|
|
C-- COMMON /SEAICE_BOUND_RL/ Various bounding values |
| 421 |
|
|
C MIN_ATEMP :: minimum air temperature (deg C) |
| 422 |
|
|
C MIN_LWDOWN :: minimum downward longwave (W/m^2) |
| 423 |
|
|
C MIN_TICE :: minimum ice temperature (deg C) |
| 424 |
|
|
C SEAICE_EPS :: small number used to reduce derivative singularities |
| 425 |
|
|
C SEAICE_EPS_SQ :: small number square |
| 426 |
|
|
C |
| 427 |
torge |
1.6 |
_RL MIN_ATEMP, MIN_LWDOWN, MIN_TICE |
| 428 |
dimitri |
1.1 |
_RL SEAICE_EPS, SEAICE_EPS_SQ |
| 429 |
|
|
COMMON /SEAICE_BOUND_RL/ |
| 430 |
torge |
1.6 |
& MIN_ATEMP, MIN_LWDOWN, MIN_TICE, |
| 431 |
torge |
1.3 |
& SEAICE_EPS, SEAICE_EPS_SQ |
| 432 |
heimbach |
1.4 |
|
| 433 |
|
|
#ifdef SEAICE_ITD |
| 434 |
dimitri |
1.2 |
CToM<<< |
| 435 |
torge |
1.6 |
C Hlimit :: ice thickness category limits (m), array of |
| 436 |
|
|
C size nITD+1 |
| 437 |
|
|
C Hlimit_c1,_c2,_c3 :: coefficients set in seaice_readparams.F to |
| 438 |
|
|
C calculate Hlimit in seaice_init_fixed.F |
| 439 |
heimbach |
1.4 |
_RL Hlimit(0:nITD) |
| 440 |
|
|
_RL Hlimit_c1, Hlimit_c2, Hlimit_c3 |
| 441 |
|
|
COMMON /SEAICE_BOUND_ITD_RL/ |
| 442 |
dimitri |
1.2 |
& Hlimit, |
| 443 |
|
|
& Hlimit_c1,Hlimit_c2,Hlimit_c3 |
| 444 |
|
|
C>>>ToM |
| 445 |
heimbach |
1.4 |
#endif /* SEAICE_ITD */ |
| 446 |
dimitri |
1.1 |
|
| 447 |
|
|
C-- Constants used by sea-ice model |
| 448 |
|
|
_RL ZERO , ONE , TWO |
| 449 |
|
|
PARAMETER ( ZERO = 0.0 _d 0, ONE = 1.0 _d 0, TWO = 2.0 _d 0 ) |
| 450 |
|
|
_RL QUART , HALF |
| 451 |
|
|
PARAMETER ( QUART = 0.25 _d 0, HALF = 0.5 _d 0 ) |
| 452 |
|
|
_RL siEps |
| 453 |
|
|
PARAMETER ( siEps = 1. _d -5 ) |
| 454 |
|
|
INTEGER MPSEUDOTIMESTEPS |
| 455 |
|
|
PARAMETER (MPSEUDOTIMESTEPS=2) |
| 456 |
|
|
|
| 457 |
|
|
C-- Constants needed by McPhee formulas for turbulent ocean fluxes : |
| 458 |
|
|
C Stanton number (dimensionless), typical friction velocity |
| 459 |
|
|
C beneath sea ice (m/s), and tapering factor (dimensionless) |
| 460 |
|
|
_RL STANTON_NUMBER, USTAR_BASE, MCPHEE_TAPER_FAC |
| 461 |
|
|
PARAMETER ( MCPHEE_TAPER_FAC = 12.5 _d 0 , STANTON_NUMBER = |
| 462 |
|
|
& 0.0056 _d 0, USTAR_BASE = 0.0125 _d 0 ) |
| 463 |
|
|
|
| 464 |
|
|
C-- identifiers for advected properties |
| 465 |
|
|
INTEGER GAD_HEFF,GAD_AREA,GAD_QICE1,GAD_QICE2,GAD_SNOW |
| 466 |
|
|
INTEGER GAD_SALT,GAD_SITR |
| 467 |
|
|
PARAMETER ( GAD_HEFF = 1, |
| 468 |
|
|
& GAD_AREA = 2, |
| 469 |
|
|
& GAD_SNOW = 3, |
| 470 |
|
|
& GAD_SALT = 4, |
| 471 |
|
|
& GAD_QICE1 = 5, |
| 472 |
|
|
& GAD_QICE2 = 6, |
| 473 |
|
|
& GAD_SITR = 7) |
| 474 |
|
|
|
| 475 |
|
|
CEH3 ;;; Local Variables: *** |
| 476 |
|
|
CEH3 ;;; mode:fortran *** |
| 477 |
|
|
CEH3 ;;; End: *** |
Parent Directory
| 
Revision Log
| 
Revision Graph