MITgcm |
|
|
|
|||
File indexing completed on 2025-11-16 06:10:45 UTC
view on githubraw file Latest commit 204068be on 2025-11-16 04:31:18 UTCb58e51ce4e Jean*0001 C *==========================================================* 0002 C | SEAICE_PARAMS.h 0003 C | o Basic parameter header for sea ice model. 0004 C *==========================================================* 809c36b928 Patr*0005 0006 C-- COMMON /SEAICE_PARM_L/ Logical parameters of sea ice model. f82b316546 Jean*0007 C - dynamics: 425e8efc36 Jean*0008 C SEAICEuseDYNAMICS :: If false, do not use dynamics; 0009 C default is to use dynamics. 1c278edd09 Jean*0010 C SEAICEuseFREEDRIFT :: If True use free drift velocity instead of EVP fe80ae39de Mart*0011 C or LSR 1459e8fe5a Mart*0012 C SEAICEuseStrImpCpl:: If true use strongly implicit coupling formulation 210ee8461e jm-c 0013 C for LSR solver (Hutchings et al 2004, 1459e8fe5a Mart*0014 C Ocean Modelling, eq.44) 2bf4ed5c39 Mart*0015 C SEAICEuseEVP :: If true use elastic viscous plastic solver 210ee8461e jm-c 0016 C SEAICEuseEVPstar :: If true use modified elastic viscous plastic 2bf4ed5c39 Mart*0017 C solver (EVP*) by Lemieux et al (2012) 210ee8461e jm-c 0018 C SEAICEuseEVPrev :: If true use "revisited" elastic viscous plastic 0a27dce68e Mart*0019 C solver following Bouillon et al. (2013), very similar 0020 C to EVP*, but uses fewer implicit terms and drops 0021 C one 1/e^2 in equations for sigma2 and sigma12 425e8efc36 Jean*0022 C SEAICEuseEVPpickup :: Set to false in order to start EVP solver with d558d05cf5 Dimi*0023 C non-EVP pickup files. Default is true. 0024 C Applied only if SEAICEuseEVP=.TRUE. f82b316546 Jean*0025 C SEAICEuseMultiTileSolver :: in LSR, use full domain tri-diagonal solver c8739d4898 Mart*0026 C SEAICEuseLSR :: If true, use default Picard solver with Line- 0027 C Successive(-over)-Relaxation, can also be true 210ee8461e jm-c 0028 C if LSR is used as a preconditioner for the c8739d4898 Mart*0029 C non-linear JFNK solver 143d9ce879 Dami*0030 C SEAICEuseLSRflex :: If true, use default Picard solver with Line- 0031 C Successive(-over)-Relaxation, but determine the 0032 C number of non-linear iterations depends on the 0033 C residual resduction, similar to the Krylov and 0034 C JFNK solvers 210ee8461e jm-c 0035 C SEAICEusePicardAsPrecon :: If true, allow SEAICEuseLSR = .TRUE. as a c8739d4898 Mart*0036 C preconditioner for non-linear JFNK problem (def. = F) 1ac05d2425 Mart*0037 C SEAICEuseKrylov :: If true, use matrix-free Krylov solver with Picard 0038 C solver instead of LSR (default: false) 1aa3c3f516 Mart*0039 C SEAICEuseJFNK :: If true, use Jacobi-free Newton-Krylov solver 0040 C instead of LSR (default: false) 0041 C SEAICEuseIMEX :: use IMplicit/EXplicit scheme with JFNK c512e371cc drin*0042 C SEAICEuseTEM :: to use the truncated ellipse method (see Geiger et al. f82b316546 Jean*0043 C 1998) set this parameter to true, default is false c512e371cc drin*0044 C SEAICEuseMCS :: to use the Mohr-Coulomb yield curve with a shear 0045 C only flow rule (Ip et al 1991), set this parameter to 0046 C true, default is false 0047 C SEAICEuseMCE :: to use the Mohr-Coulomb yield curve with elliptical 0048 C plastic potential (similarly to Hibler and Schulson 0049 C 2000 without the elliptical cap) set this parameter 0050 C to true, default is false 0051 C SEAICEuseTD :: to use the teardrop yield curve (Zhang and Rothrock, 0052 C 2005) set this parameter to true, default is false 0053 C SEAICEusePL :: to use the parabolic lens yield curve (Zhang and 0054 C Rothrock, 2005) set this parameter to true, 0055 C default is false 1aa3c3f516 Mart*0056 C SEAICEuseTilt :: If true then include surface tilt term in dynamics 4935169c82 Mart*0057 C SEAICEuseMetricTerms :: use metric terms for dynamics solver 0058 C (default = .true. ) f82b316546 Jean*0059 C SEAICE_no_slip :: apply no slip boundary conditions to seaice velocity 5dac41bc68 Mart*0060 C SEAICE_2ndOrderBC :: apply 2nd order no slip boundary conditions (works 0061 C only with EVP, JFNK or KRYLOV solver, default=F) f82b316546 Jean*0062 C SEAICE_maskRHS :: mask the RHS of the solver where there is no ice 0063 C SEAICE_clipVelocities :: clip velocities to +/- 40cm/s 79022779f5 Mart*0064 C SEAICEaddSnowMass :: in computing seaiceMass, add snow contribution b8665dacca Mart*0065 C default is .TRUE. f82b316546 Jean*0066 C useHB87stressCoupling :: use an intergral over ice and ocean surface 0067 C layer to define surface stresses on ocean 0068 C following Hibler and Bryan (1987, JPO) 210ee8461e jm-c 0069 C SEAICEupdateOceanStress :: If TRUE, update ocean surface stress 0070 C accounting for seaice cover (default= T) e501eee760 Mart*0071 C SEAICEuseBDF2 :: use 2nd-order backward difference approach 0072 C for momentum equations as described in 0073 C Lemieux et al. 2014, JCP 1aa3c3f516 Mart*0074 C so far only implemented for JFNK-solver 210ee8461e jm-c 0075 C useHibler79IceStrength :: if true original ice strength parameterization 5fd900657f Mart*0076 C other use Rothrock (1975) parameterization based 0077 C on energetics and an ice thickness distribution 0078 C (default = .true.) 36d0480ad6 Mart*0079 C SEAICEscaleSurfStress :: if TRUE, scale ice-ocean and ice-atmosphere 0080 C stress on ice by concenration (AREA) following b8665dacca Mart*0081 C Connolley et al. (2004), JPO. (default = .TRUE.) 5fd900657f Mart*0082 C SEAICEsimpleRidging :: use Hibler(1979) ridging (default=.true.) 3bd2cd9e40 Mart*0083 C SEAICEuseLinRemapITD :: use linear remapping (Lipscomb et al. 2001) 0084 C .TRUE. by default f82b316546 Jean*0085 C - advection: 0086 C SEAICEuseFluxForm :: use flux form for advection and diffusion 0087 C of seaice 425e8efc36 Jean*0088 C SEAICEadvHeff :: turn on advection of effective thickness 09e3b53265 Mart*0089 C (default = .true.) 0090 C SEAICEadvArea :: turn on advection of fraction area 0091 C (default = .true.) eefc792916 Mart*0092 C SEAICEadvSnow :: turn on advection of snow (does not work with 0093 C non-default Leap-frog scheme for advection) fdfa8e151f Dimi*0094 C SEAICEadvSalt :: turn on advection of salt (does not work with 0095 C non-default Leap-frog scheme for advection) e0fa1cecbf Mart*0096 C SEAICEmultiDimAdvection:: internal flag, set to false if any sea ice 0097 C variable uses a non-multi-dimensional advection 0098 C scheme 2e05c0b826 Mart*0099 C SEAICEmomAdvection:: turn on advection of momentum (default = .false.) 0100 C SEAICEhighOrderVorticity :: momentum advection parameters analogous to 0101 C SEAICEupwindVorticity :: highOrderVorticity, upwindVorticity, 0102 C SEAICEuseAbsVorticity :: useAbsVorticity, useJamartMomAdv for vector 0103 C SEAICEuseJamartMomAdv :: invariant momentum in the ocean f82b316546 Jean*0104 C - thermodynamics: 09e3b53265 Mart*0105 C usePW79thermodynamics :: use "0-layer" thermodynamics as described in f414077a87 Mart*0106 C Parkinson and Washington (1979) and Hibler (1979) 210ee8461e jm-c 0107 C SEAICE_useMultDimSnow :: use same fixed pdf for snow as for b8665dacca Mart*0108 C multi-thickness-category ice (default=.TRUE.) f82b316546 Jean*0109 C SEAICEuseFlooding :: turn on scheme to convert submerged snow into ice 0110 C SEAICEheatConsFix :: If true then fix ocn<->seaice advective heat flux. b58e51ce4e Jean*0111 C useMaykutSatVapPoly :: use Maykut Polynomial for saturation vapor pressure 0112 C instead of extended temp-range exponential law; def=F. f82b316546 Jean*0113 C SEAICE_mcPheeStepFunc :: use step function (not linear tapering) in fe80ae39de Mart*0114 C ocean-ice turbulent flux 6ec4646d60 Gael*0115 C SEAICE_doOpenWaterGrowth :: use open water heat flux directly to grow ice fe80ae39de Mart*0116 C (when false cool ocean, and grow later if needed) 6ec4646d60 Gael*0117 C SEAICE_doOpenWaterMelt :: use open water heat flux directly to melt ice fe80ae39de Mart*0118 C (when false warm ocean, and melt later if needed) 210ee8461e jm-c 0119 C SEAICE_growMeltByConv :: grow/melt according to convergence of turbulence 62cc8945c8 Gael*0120 C and conduction, rather than in two steps (default) 1c278edd09 Jean*0121 C SEAICE_salinityTracer :: use SItracer to exchange and trace ocean fe80ae39de Mart*0122 C salt in ice f82b316546 Jean*0123 C SEAICE_ageTracer :: use SItracer to trace the age of ice f414077a87 Mart*0124 C SEAICErestoreUnderIce :: restore surface T/S also underneath ice 0125 C ( default is false ) f82b316546 Jean*0126 C - other (I/O, ...): 0127 C SEAICEwriteState :: If true, write sea ice state to file; 0128 C default is false. 07abca99f1 Jean*0129 C SEAICE_dump_mdsio :: write snap-shot output using MDSIO 0130 C SEAICE_mon_stdio :: write monitor to std-outp 0131 C SEAICE_dump_mnc :: write snap-shot output using MNC 0132 C SEAICE_mon_mnc :: write monitor to netcdf file a2d714f5b7 Jean*0133 LOGICAL 1459e8fe5a Mart*0134 & SEAICEuseDYNAMICS, SEAICEuseFREEDRIFT, SEAICEuseStrImpCpl, 0a27dce68e Mart*0135 & SEAICEuseEVP, SEAICEuseEVPstar, SEAICEuseEVPrev, 0136 & SEAICEuseEVPpickup, f82b316546 Jean*0137 & SEAICEuseMultiTileSolver, 143d9ce879 Dami*0138 & SEAICEuseLSR, SEAICEuseLSRflex, SEAICEuseKrylov, 1ac05d2425 Mart*0139 & SEAICEuseJFNK, SEAICEuseIMEX, SEAICEuseBDF2, c8739d4898 Mart*0140 & SEAICEusePicardAsPrecon, 5fd900657f Mart*0141 & useHibler79IceStrength, SEAICEsimpleRidging, c512e371cc drin*0142 & SEAICEuseLinRemapITD, SEAICEuseTD, SEAICEusePL, f82b316546 Jean*0143 & SEAICEuseTEM, SEAICEuseTilt, SEAICEuseMetricTerms, c512e371cc drin*0144 & SEAICEuseMCS, SEAICEuseMCE, 210ee8461e jm-c 0145 & SEAICE_no_slip, SEAICE_2ndOrderBC, 5dac41bc68 Mart*0146 & SEAICE_maskRHS, SEAICEscaleSurfStress, 210ee8461e jm-c 0147 & SEAICE_clipVelocities, SEAICEaddSnowMass, 0148 & useHB87stressCoupling, SEAICEupdateOceanStress, f82b316546 Jean*0149 & SEAICEuseFluxForm, SEAICEadvHeff, SEAICEadvArea, e0fa1cecbf Mart*0150 & SEAICEmultiDimAdvection, 2e05c0b826 Mart*0151 & SEAICEadvSnow, SEAICEadvSalt, SEAICEmomAdvection, 0152 & SEAICEhighOrderVorticity, SEAICEupwindVorticity, 0153 & SEAICEuseAbsVorticity, SEAICEuseJamartMomAdv, f82b316546 Jean*0154 & usePW79thermodynamics, 0155 & SEAICE_useMultDimSnow, SEAICEuseFlooding, SEAICEheatConsFix, 0156 & useMaykutSatVapPoly, SEAICE_mcPheeStepFunc, 6ec4646d60 Gael*0157 & SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt, be02c52974 Gael*0158 & SEAICE_salinityTracer, SEAICE_ageTracer, 62cc8945c8 Gael*0159 & SEAICErestoreUnderIce, SEAICE_growMeltByConv, f82b316546 Jean*0160 & SEAICEwriteState, b7411f1a84 Jean*0161 & SEAICE_dump_mdsio, SEAICE_mon_stdio, 0162 & SEAICE_dump_mnc, SEAICE_mon_mnc 809c36b928 Patr*0163 COMMON /SEAICE_PARM_L/ 1459e8fe5a Mart*0164 & SEAICEuseDYNAMICS, SEAICEuseFREEDRIFT, SEAICEuseStrImpCpl, 0a27dce68e Mart*0165 & SEAICEuseEVP, SEAICEuseEVPstar, SEAICEuseEVPrev, 0166 & SEAICEuseEVPpickup, f82b316546 Jean*0167 & SEAICEuseMultiTileSolver, 143d9ce879 Dami*0168 & SEAICEuseLSR, SEAICEuseLSRflex, SEAICEuseKrylov, 210ee8461e jm-c 0169 & SEAICEuseJFNK, SEAICEuseIMEX, SEAICEuseBDF2, c8739d4898 Mart*0170 & SEAICEusePicardAsPrecon, 5fd900657f Mart*0171 & useHibler79IceStrength, SEAICEsimpleRidging, c512e371cc drin*0172 & SEAICEuseLinRemapITD, SEAICEuseTD, SEAICEusePL, f82b316546 Jean*0173 & SEAICEuseTEM, SEAICEuseTilt, SEAICEuseMetricTerms, c512e371cc drin*0174 & SEAICEuseMCS, SEAICEuseMCE, 5dac41bc68 Mart*0175 & SEAICE_no_slip, SEAICE_2ndOrderBC, 0176 & SEAICE_maskRHS, SEAICEscaleSurfStress, 210ee8461e jm-c 0177 & SEAICE_clipVelocities, SEAICEaddSnowMass, 0178 & useHB87stressCoupling, SEAICEupdateOceanStress, f82b316546 Jean*0179 & SEAICEuseFluxForm, SEAICEadvHeff, SEAICEadvArea, 2e05c0b826 Mart*0180 & SEAICEadvSnow, SEAICEadvSalt, SEAICEmomAdvection, e0fa1cecbf Mart*0181 & SEAICEmultiDimAdvection, 2e05c0b826 Mart*0182 & SEAICEhighOrderVorticity, SEAICEupwindVorticity, 0183 & SEAICEuseAbsVorticity, SEAICEuseJamartMomAdv, f82b316546 Jean*0184 & usePW79thermodynamics, 0185 & SEAICE_useMultDimSnow, SEAICEuseFlooding, SEAICEheatConsFix, 0186 & useMaykutSatVapPoly, SEAICE_mcPheeStepFunc, 6ec4646d60 Gael*0187 & SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt, be02c52974 Gael*0188 & SEAICE_salinityTracer, SEAICE_ageTracer, 62cc8945c8 Gael*0189 & SEAICErestoreUnderIce, SEAICE_growMeltByConv, f82b316546 Jean*0190 & SEAICEwriteState, b7411f1a84 Jean*0191 & SEAICE_dump_mdsio, SEAICE_mon_stdio, 0192 & SEAICE_dump_mnc, SEAICE_mon_mnc 809c36b928 Patr*0193 0194 C-- COMMON /SEAICE_PARM_I/ Integer valued parameters of sea ice model. 210ee8461e jm-c 0195 C IMAX_TICE :: number of iterations for ice surface temp 4413db9420 Mart*0196 C (default=10) 210ee8461e jm-c 0197 C postSolvTempIter :: select flux calculation after surf. temp solver 4413db9420 Mart*0198 C iteration 0199 C 0 = none, i.e., from last iter 210ee8461e jm-c 0200 C 1 = use linearized approx (consistent with tsurf 4413db9420 Mart*0201 C finding) 0202 C 2 = full non-lin form 79df32c3f1 Mart*0203 C SOLV_NCHECK :: iteration interval for LSR-solver convergence test 0204 C SEAICEnonLinIterMax :: number of allowed non-linear solver iterations 0205 C for implicit solvers (JFNK and Picard) (>= 2) 0206 C SEAICElinearIterMax :: number of allowed linear solver iterations for 0207 C for implicit solvers (JFNK and Picard) C 0208 C SEAICEpreconNL_Iter :: number non-linear iterations in preconditioner 0209 C SEAICEpreconLinIter :: number linear iterations in preconditioner 2bf4ed5c39 Mart*0210 C SEAICEnEVPstarSteps :: number of evp*-steps e501eee760 Mart*0211 C SEAICEmomStartBDF :: number of previous u/vIce time levels available 0212 C to start (or restart) BDF2 scheme. a88fc20c6c Mart*0213 C SEAICE_JFNK_lsIter :: number of Newton iterations after which the 0214 C line search is started c704c5a1ef Mart*0215 C SEAICE_JFNK_lsLmax :: max. number line search iterations (default = 4) f6f4a9e227 Mart*0216 C SEAICE_JFNK_tolIter :: number of Newton iterations after which the 0217 C the tolerance is relaxed again (default = 100) b8665dacca Mart*0218 C SEAICE_OLx/y :: overlaps for LSR-solver and for the 0219 C LSR-preconditioner in JFNK and KRYLOV solver; 0220 C for 0 < SEAICE_OLx/y 0 <= OLx/y-2 the LSR solver 0221 C and preconditioner use a restricted additive 0222 C Schwarz method (default = OLx/y-2). 210ee8461e jm-c 0223 C LSR_mixIniGuess :: control mixing of free-drift sol. into LSR initial 4413db9420 Mart*0224 C guess 210ee8461e jm-c 0225 C :: =0 : nothing; =1 : no mix, but print free-drift 4413db9420 Mart*0226 C resid.; 604907fea6 Jean*0227 C :: =2,4 : mix with (1/local-residual)^2,4 factor 76c8933abb Gael*0228 C SEAICEpresPow0 :: HEFF exponent for ice strength below SEAICEpresH0 0229 C SEAICEpresPow1 :: HEFF exponent for ice strength above SEAICEpresH0 5fd900657f Mart*0230 C rigding parameters (only active when SEAICE_ITD is defined) 0231 C SEAICEpartFunc :: =0 use Thorndyke et al (1975) participation function 0232 C :: =1 use Lipscomb et al (2007) participation function 0233 C SEAICEredistFunc :: =0 assume ridged ice is uniformly distributed 0234 C (Hibler, 1980) 210ee8461e jm-c 0235 C =1 Following Lipscomb et al. (2007), ridged ice is 5fd900657f Mart*0236 C distributed following an exponentially 0237 C decaying function 0238 C SEAICEridgingIterMax :: maximum number of ridging iterations 0239 C end ridging parameters 2e05c0b826 Mart*0240 C SEAICEselectKEscheme :: momentum advection parameters analogous 0241 C SEAICEselectVortScheme :: to selectKEscheme and selectVortScheme 425e8efc36 Jean*0242 C SEAICEadvScheme :: sets the advection scheme for thickness and area b8665dacca Mart*0243 C (default = 77) 425e8efc36 Jean*0244 C SEAICEadvSchArea :: sets the advection scheme for area 0245 C SEAICEadvSchHeff :: sets the advection scheme for effective thickness 0246 C (=volume), snow thickness, and salt if available 0247 C SEAICEadvSchSnow :: sets the advection scheme for snow on sea-ice 0248 C SEAICEadvSchSalt :: sets the advection scheme for sea ice salinity 0249 C SEAICEadvSchSnow :: sets the advection scheme for snow on sea-ice 6ec4646d60 Gael*0250 C SEAICE_areaLossFormula :: selects formula for ice cover loss from melt 0251 C :: 1=from all but only melt conributions by ATM and OCN 0252 C :: 2=from net melt-growth>0 by ATM and OCN 0253 C :: 3=from predicted melt by ATM 210ee8461e jm-c 0254 C SEAICE_areaGainFormula :: selects formula for ice cover gain from open 4413db9420 Mart*0255 C water growth 6ec4646d60 Gael*0256 C :: 1=from growth by ATM 0257 C :: 2=from predicted growth by ATM 2349c46237 Mart*0258 C SEAICEetaZmethod :: determines how shear-viscosity eta is computed at 0259 C Z-points 0260 C 0=simple averaging from C-points (default and old) 210ee8461e jm-c 0261 C 3=weighted averaging of squares of strain rates c43a67fada Mart*0262 C (recommended for energy conservation) f5282c5b03 Gael*0263 C SEAICE_multDim :: number of ice categories 425e8efc36 Jean*0264 C SEAICE_debugPointI :: I,J index for seaice-specific debuggin 136908bfac Ian *0265 C SEAICE_debugPointJ 809c36b928 Patr*0266 C b58e51ce4e Jean*0267 INTEGER IMAX_TICE, postSolvTempIter 79df32c3f1 Mart*0268 INTEGER SOLV_NCHECK 0269 INTEGER SEAICEnonLinIterMax, SEAICElinearIterMax 0270 INTEGER SEAICEpreconLinIter, SEAICEpreconNL_Iter 3de1dea93d Jean*0271 INTEGER LSR_mixIniGuess 2bf4ed5c39 Mart*0272 INTEGER SEAICEnEVPstarSteps e501eee760 Mart*0273 INTEGER SEAICEmomStartBDF f6f4a9e227 Mart*0274 INTEGER SEAICE_JFNK_lsIter, SEAICE_JFNK_tolIter c704c5a1ef Mart*0275 INTEGER SEAICE_JFNK_lsLmax b3c9aef4a6 Mart*0276 INTEGER SEAICE_OLx, SEAICE_OLy 2e05c0b826 Mart*0277 INTEGER SEAICEselectKEscheme, SEAICEselectVortScheme 0f6ad6f523 Mart*0278 INTEGER SEAICEadvScheme bd4f8028de Mart*0279 INTEGER SEAICEadvSchArea 0280 INTEGER SEAICEadvSchHeff 0281 INTEGER SEAICEadvSchSnow fdfa8e151f Dimi*0282 INTEGER SEAICEadvSchSalt 3a3bf6419a Gael*0283 INTEGER SEAICEadjMODE 6ec4646d60 Gael*0284 INTEGER SEAICE_areaLossFormula 1c278edd09 Jean*0285 INTEGER SEAICE_areaGainFormula 2349c46237 Mart*0286 INTEGER SEAICEetaZmethod f5282c5b03 Gael*0287 INTEGER SEAICE_multDim 136908bfac Ian *0288 INTEGER SEAICE_debugPointI 0289 INTEGER SEAICE_debugPointJ 76c8933abb Gael*0290 INTEGER SEAICEpresPow0, SEAICEpresPow1 5fd900657f Mart*0291 INTEGER SEAICEpartFunc, SEAICEredistFunc 0292 INTEGER SEAICEridgingIterMax 425e8efc36 Jean*0293 COMMON /SEAICE_PARM_I/ 79df32c3f1 Mart*0294 & IMAX_TICE, postSolvTempIter, SOLV_NCHECK, 0295 & SEAICEnonLinIterMax, SEAICElinearIterMax, 0296 & SEAICEpreconLinIter, SEAICEpreconNL_Iter, 3de1dea93d Jean*0297 & LSR_mixIniGuess, 2bf4ed5c39 Mart*0298 & SEAICEnEVPstarSteps, e501eee760 Mart*0299 & SEAICEmomStartBDF, b3c9aef4a6 Mart*0300 & SEAICE_JFNK_lsIter, SEAICE_OLx, SEAICE_OLy, c704c5a1ef Mart*0301 & SEAICE_JFNK_lsLmax, SEAICE_JFNK_tolIter, 76c8933abb Gael*0302 & SEAICEpresPow0, SEAICEpresPow1, 210ee8461e jm-c 0303 & SEAICEpartFunc, SEAICEredistFunc, SEAICEridgingIterMax, 2e05c0b826 Mart*0304 & SEAICEselectKEscheme, SEAICEselectVortScheme, bd4f8028de Mart*0305 & SEAICEadvScheme, 0306 & SEAICEadvSchArea, 0307 & SEAICEadvSchHeff, fdfa8e151f Dimi*0308 & SEAICEadvSchSnow, f681b7f5d4 Dimi*0309 & SEAICEadvSchSalt, 65f34462d4 Gael*0310 & SEAICEadjMODE, 6ec4646d60 Gael*0311 & SEAICE_areaLossFormula, 0312 & SEAICE_areaGainFormula, f5282c5b03 Gael*0313 & SEAICE_multDim, 2349c46237 Mart*0314 & SEAICEetaZmethod, 136908bfac Ian *0315 & SEAICE_debugPointI, 0316 & SEAICE_debugPointJ 809c36b928 Patr*0317 0318 C-- COMMON /SEAICE_PARM_C/ Character valued sea ice model parameters. 5bb179ddc2 Mart*0319 C AreaFile :: File containing initial sea-ice concentration 0320 C HsnowFile :: File containing initial snow thickness 0321 C HsaltFile :: File containing initial sea ice salt content 0322 C HeffFile :: File containing initial sea-ice thickness 0323 C uIceFile :: File containing initial sea-ice U comp. velocity 0324 C vIceFile :: File containing initial sea-ice V comp. velocity 0325 C uCoastLineFile :: File containing the some measure of coastline 0326 C roughness length (in m) at the U-points in the 0327 C X-direction (i.e. for the U-equation). 0328 C vCoastLineFile :: Files containing the some measure of coastline 0329 C roughness length (in m) at the V-points in the 0330 C Y-direction (i.e. for the V-equation). 0331 C 0332 C This roughness length can be the subgrid 0333 C scale length of the coastline in a grid cell 0334 C projected in the direction normal to the u/v- 0335 C direction as in Liu et al. (2022), but it can 0336 C also be anything that is a good proxy of coast 0337 C line roughness. 0338 C 190e1777cb Dimi*0339 C !!! NOTE !!! Initial sea-ice thickness can also be set using 0340 C SEAICE_initialHEFF below. But a constant initial condition 0341 C can mean large artificial fluxes of heat and freshwater in 0342 C the surface layer during the first model time step. 809c36b928 Patr*0343 C 09066b09cb Mart*0344 CHARACTER*(MAX_LEN_FNAM) AreaFile de31ea8481 Dimi*0345 CHARACTER*(MAX_LEN_FNAM) HsnowFile fdfa8e151f Dimi*0346 CHARACTER*(MAX_LEN_FNAM) HsaltFile 6060ec2938 Dimi*0347 CHARACTER*(MAX_LEN_FNAM) HeffFile 425e8efc36 Jean*0348 CHARACTER*(MAX_LEN_FNAM) uIceFile 0349 CHARACTER*(MAX_LEN_FNAM) vIceFile 5bb179ddc2 Mart*0350 CHARACTER*(MAX_LEN_FNAM) uCoastLineFile 0351 CHARACTER*(MAX_LEN_FNAM) vCoastLineFile f681b7f5d4 Dimi*0352 COMMON /SEAICE_PARM_C/ 425e8efc36 Jean*0353 & AreaFile, HsnowFile, HsaltFile, HeffFile, 5bb179ddc2 Mart*0354 & uIceFile, vIceFile, uCoastLineFile, vCoastLineFile 809c36b928 Patr*0355 0356 C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model. 425e8efc36 Jean*0357 C SEAICE_deltaTtherm :: Seaice timestep for thermodynamic equations (s) 0358 C SEAICE_deltaTdyn :: Seaice timestep for dynamic solver (s) e45202e340 Mart*0359 C SEAICE_LSRrelaxU/V :: relaxation parameter for LSR-solver: U/V-component 425e8efc36 Jean*0360 C SEAICE_deltaTevp :: Seaice timestep for EVP solver (s) 0361 C SEAICE_elasticParm :: parameter that sets relaxation timescale 0362 C tau = SEAICE_elasticParm * SEAICE_deltaTdyn 0363 C SEAICE_evpTauRelax :: relaxation timescale tau (s) 3de1dea93d Jean*0364 C SEAICE_evpDampC :: evp damping constant (Hunke,JCP,2001) (kg/m^2) 2bf4ed5c39 Mart*0365 C SEAICE_evpAlpha :: dimensionless parameter 2*evpTauRelax/deltaTevp 0366 C SEAICE_evpBeta :: dimensionless parameter deltaTdyn/deltaTevp 1a0d5a451b Mart*0367 C SEAICEaEVPcoeff :: main coefficent for adaptive EVP (largest 0368 C stabilized frequency) 0369 C SEAICEaEVPcStar :: multiple of stabilty factor: alpha*beta=cstar*gamma 0370 C SEAICEaEVPalphaMin :: lower limit of alpha and beta, regularisation 210ee8461e jm-c 0371 C to prevent singularities of system matrix, 1a0d5a451b Mart*0372 C e.g. when ice concentration is too low. 204068be6a R. S*0373 C SEAICE_evpAreaReg :: Specifies a minimun ice fraction for the purposes 0374 C of regularizations in the calculation of denomU/V, 0375 C to enhance the stability of EVP; off by default, 0376 C turn on with a sensible value, e.g. 1e-5 79df32c3f1 Mart*0377 C SEAICEnonLinTol :: non-linear tolerance parameter for implicit solvers 812a8fb151 Mart*0378 C JFNKgamma_lin_min/max :: tolerance parameters for linear JFNK solver 0379 C JFNKres_t :: tolerance parameter for FGMRES residual 2e75855dde Mart*0380 C JFNKres_tFac :: if set, JFNKres_t=JFNKres_tFac*(initial residual) c8f5b4180e Mart*0381 C SEAICE_JFNKepsilon :: step size for the FD-gradient in s/r seaice_jacvec c704c5a1ef Mart*0382 C SEAICE_JFNK_lsGamma:: reduction factor for line search (default 0.5) 1f3ad2d627 Mart*0383 C SEAICE_JFNKphi :: [0,1] parameter for inexact Newton Method (def = 1) 0384 C SEAICE_JFNKalpha :: (1,2] parameter for inexact Newton Method (def = 1) 3de1dea93d Jean*0385 C SEAICE_zetaMaxFac :: factor determining the maximum viscosity (s) bb06934ea0 Mart*0386 C (default = 5.e+12/2.e4 = 2.5e8) fe80ae39de Mart*0387 C SEAICE_zetaMin :: lower bound for viscosity (default = 0) (N s/m^2) 76c8933abb Gael*0388 C SEAICEpresH0 :: HEFF threshold for ice strength (m) 425e8efc36 Jean*0389 C SEAICE_monFreq :: SEAICE monitor frequency. (s) 0390 C SEAICE_dumpFreq :: SEAICE dump frequency. (s) 0391 C SEAICE_initialHEFF :: initial sea-ice thickness (m) 0392 C SEAICE_rhoAir :: density of air (kg/m^3) 0393 C SEAICE_rhoIce :: density of sea ice (kg/m^3) 0394 C SEAICE_rhoSnow :: density of snow (kg/m^3) 0395 C ICE2WATR :: ratio of sea ice density to water density 0396 C SEAICE_cpAir :: specific heat of air (J/kg/K) f834b21bef Dimi*0397 C 5867b94c2f Mart*0398 C OCEAN_drag :: unitless air-ocean drag coefficient (default 0.001) 0399 C SEAICE_drag :: unitless air-ice drag coefficient (default 0.001) 0400 C SEAICE_waterDrag :: unitless water-ice drag coefficient (default 0.0055) 0401 C SEAICEdWatMin :: minimum linear water-ice drag applied to DWATN 0402 C (default 0.25 m/s) 0403 C 425e8efc36 Jean*0404 C SEAICE_dryIceAlb :: winter albedo 0405 C SEAICE_wetIceAlb :: summer albedo 0406 C SEAICE_drySnowAlb :: dry snow albedo 0407 C SEAICE_wetSnowAlb :: wet snow albedo 1c278edd09 Jean*0408 C HO :: AKA "lead closing parameter", demarcation thickness fe80ae39de Mart*0409 C between thin and thick ice. Alternatively, HO (in 0410 C meters) can be interpreted as the thickness of ice 0411 C formed in open water. 0412 C HO is a key ice-growth parameter that determines 0413 C the partition between vertical and lateral growth. 0414 C The default is 0.5m, increasing this value leads 0415 C slower formation of a closed ice cover and thus to 0416 C more ice (and thicker) ice, decreasing to faster 0417 C formation of a closed ice cover (leads are closing 0418 C faster) and thus less (thinner) ice. f834b21bef Dimi*0419 C 425e8efc36 Jean*0420 C SEAICE_drag_south :: Southern Ocean SEAICE_drag 0421 C SEAICE_waterDrag_south :: Southern Ocean SEAICE_waterDrag 0422 C SEAICE_dryIceAlb_south :: Southern Ocean SEAICE_dryIceAlb 0423 C SEAICE_wetIceAlb_south :: Southern Ocean SEAICE_wetIceAlb 0424 C SEAICE_drySnowAlb_south :: Southern Ocean SEAICE_drySnowAlb 0425 C SEAICE_wetSnowAlb_south :: Southern Ocean SEAICE_wetSnowAlb 0426 C HO_south :: Southern Ocean HO f834b21bef Dimi*0427 C 210ee8461e jm-c 0428 C Parameters for basal drag of grounded ice following abb637800a Mart*0429 C Lemieux et al. (2015), doi:10.1002/2014JC010678 0430 C SEAICE_cBasalStar (default = SEAICE_cStar) 0431 C SEAICEbasalDragU0 (default = 5e-5) 0432 C SEAICEbasalDragK1 (default = 8) 210ee8461e jm-c 0433 C SEAICEbasalDragK2 :: if > 0, turns on basal drag abb637800a Mart*0434 C (default = 0, Lemieux suggests 15) 5bb179ddc2 Mart*0435 C SEAICEsideDrag :: if > 0, turns on lateral static drag 0436 C if < 0, turns on lateral quadratic drag 0437 C both are different landfast ice parameterisations 0438 C (Liu et al 2022 use 2e-4, 0439 C the default = 0 turns off the parameterisations) abb637800a Mart*0440 C b58e51ce4e Jean*0441 C SEAICE_wetAlbTemp :: Temp (deg.C) above which wet-albedo values are used 425e8efc36 Jean*0442 C SEAICE_waterAlbedo :: water albedo 0443 C SEAICE_strength :: sea-ice strength Pstar ba6cfc5714 Mart*0444 C SEAICE_cStar :: sea-ice strength paramter C* (def: 20) 0445 C SEAICE_tensilFac :: sea-ice tensile strength factor, values in [0,1] d104051171 Mart*0446 C SEAICE_tensilDepth :: crtical depth for sea-ice tensile strength (def 0.) 0adbdb4edd Mart*0447 C SEAICEpressReplFac :: interpolator between PRESS0 and regularized PRESS 0448 C 1. (default): pure pressure replace method (PRESS) 0449 C 0. : pure Hibler (1979) method (PRESS0) 425e8efc36 Jean*0450 C SEAICE_eccen :: sea-ice eccentricity of the elliptical yield curve c512e371cc drin*0451 C SEAICE_eccfr :: sea-ice eccentricity of the elliptical flow rule 425e8efc36 Jean*0452 C SEAICE_lhFusion :: latent heat of fusion for ice and snow (J/kg) 0453 C SEAICE_lhEvap :: latent heat of evaporation for water (J/kg) 0454 C SEAICE_dalton :: Dalton number (= sensible heat transfer coefficient) 0455 C SEAICE_iceConduct :: sea-ice conductivity 0456 C SEAICE_snowConduct :: snow conductivity 0d2761b0d1 Mart*0457 C SEAICE_emissivity :: longwave ocean-surface emissivity (-) 0458 C SEAICE_ice_emiss :: longwave ice-surface emissivity (-) 0459 C SEAICE_snow_emiss :: longwave snow-surface emissivity (-) d778130a13 Mart*0460 C SEAICE_boltzmann :: Stefan-Boltzman constant (not a run time parameter) b58e51ce4e Jean*0461 C SEAICE_snowThick :: cutoff snow thickness (for snow-albedo) 0462 C SEAICE_shortwave :: ice penetration shortwave radiation factor 1c278edd09 Jean*0463 C SEAICE_saltFrac :: salinity of newly formed seaice defined as a fe80ae39de Mart*0464 C fraction of the ocean surface salinity at the time 0465 C of freezing 1c278edd09 Jean*0466 C SEAICE_salt0 :: prescribed salinity of seaice (in g/kg). fe80ae39de Mart*0467 C facOpenGrow :: 0./1. version of logical SEAICE_doOpenWaterGrowth 0468 C facOpenMelt :: 0./1. version of logical SEAICE_doOpenWaterMelt 1c278edd09 Jean*0469 C SEAICE_mcPheePiston:: ocean-ice turbulent flux "piston velocity" (m/s) fe80ae39de Mart*0470 C that sets melt efficiency. 1c278edd09 Jean*0471 C SEAICE_mcPheeTaper :: tapering down of turbulent flux term with ice 0472 C concentration. The 100% cover turb. flux is fe80ae39de Mart*0473 C multiplied by 1.-SEAICE_mcPheeTaper 1c278edd09 Jean*0474 C SEAICE_frazilFrac :: Fraction of surface level negative heat content 0475 C anomalies (relative to the local freezing point) fe80ae39de Mart*0476 C may contribute as frazil over one time step. 840c7fba30 Gael*0477 C SEAICE_tempFrz0 :: sea water freezing point is fe80ae39de Mart*0478 C SEAICE_dTempFrz_dS :: tempFrz = SEAICE_tempFrz0 + salt*SEAICE_dTempFrz_dS 210ee8461e jm-c 0479 C SEAICE_PDF :: prescribed sea-ice distribution within grid box 425e8efc36 Jean*0480 C SEAICEstressFactor :: factor by which ice affects wind stress (default=1) 0481 C LSR_ERROR :: sets accuracy of LSR solver 0482 C DIFF1 :: parameter used in advect.F c512e371cc drin*0483 C SEAICEtdMU :: slope parameter for the teardrop and parabolic lens 0484 C yield curves 3daf25222c Mart*0485 C SEAICE_deltaMin :: small number used to reduce singularities of Delta fe80ae39de Mart*0486 C SEAICE_area_max :: usually set to 1. Seeting areaMax below 1 specifies 1c278edd09 Jean*0487 C the minimun amount of leads (1-areaMax) in the fe80ae39de Mart*0488 C ice pack. 425e8efc36 Jean*0489 C SEAICE_area_floor :: usually set to 1x10^-5. Specifies a minimun 0490 C ice fraction in the ice pack. 0491 C SEAICE_area_reg :: usually set to 1x10^-5. Specifies a minimun 0492 C ice fraction for the purposes of regularization 0493 C SEAICE_hice_reg :: usually set to 5 cm. Specifies a minimun 0494 C ice thickness for the purposes of regularization 0495 C SEAICEdiffKhArea :: sets the diffusivity for area (m^2/s) 0496 C SEAICEdiffKhHeff :: sets the diffusivity for effective thickness (m^2/s) 0497 C SEAICEdiffKhSnow :: sets the diffusivity for snow on sea-ice (m^2/s) 0498 C SEAICEdiffKhSalt :: sets the diffusivity for sea ice salinity (m^2/s) 0499 C SEAICE_airTurnAngle :: turning angles of air-ice interfacial stress 0500 C SEAICE_waterTurnAngle :: and ice-water interfacial stress (in degrees) f82b316546 Jean*0501 C SEAICE_tauAreaObsRelax :: Timescale of relaxation to observed d32fe07ad8 Patr*0502 C sea ice concentration (s), default=unset 5fd900657f Mart*0503 C ridging parameters (Lipscomb et al, 2007, Bitz et al. 2001): 0504 C SEAICE_cf :: ratio of total energy sinks to gravitational sink 0505 C (scales ice strength, suggested values: 2 to 17) 0506 C SEAICEgStar :: maximum ice concentration that participates in ridging 0507 C SEAICEhStar :: empirical thickness (ridging parameter) 0508 C SEAICEaStar :: ice concentration parameter similar to gStar for 0509 C exponential distribution (Lipscomb et al 2007) 0510 C SEAICEshearParm :: <=1 reduces amount of energy lost to ridge building 0511 C SEAICEmuRidging :: tuning parameter similar to hStar for Lipcomb et al 0512 C (2007)-scheme 0513 C SEAICEmaxRaft :: regularization parameter (default=1) 3ac108d8cf Mart*0514 C SEAICEsnowFracRidge :: fraction of snow that remains on ridged dc54d31829 Ian *0515 C SINegFac :: SIADV over/undershoot factor in FW/Adjoint c512e371cc drin*0516 C SEAICEmcMu :: parameter for MC yield curve for useMCE, useMCS and 0517 C useTEM options, default is one 809c36b928 Patr*0518 C 6e2f4e58fa Mart*0519 _RL SEAICE_deltaTtherm, SEAICE_deltaTdyn, SEAICE_deltaTevp e45202e340 Mart*0520 _RL SEAICE_LSRrelaxU, SEAICE_LSRrelaxV b7411f1a84 Jean*0521 _RL SEAICE_monFreq, SEAICE_dumpFreq a2d714f5b7 Jean*0522 _RL SEAICE_initialHEFF 07b4a12853 Mart*0523 _RL SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR fff6be1885 Mart*0524 _RL SEAICE_cpAir 5867b94c2f Mart*0525 _RL SEAICE_drag, SEAICE_waterDrag, SEAICEdWatMin 0526 _RL SEAICE_dryIceAlb, SEAICE_wetIceAlb 0527 _RL SEAICE_drySnowAlb, SEAICE_wetSnowAlb, HO f834b21bef Dimi*0528 _RL SEAICE_drag_south, SEAICE_waterDrag_south 0529 _RL SEAICE_dryIceAlb_south, SEAICE_wetIceAlb_south f81e0cd482 Dimi*0530 _RL SEAICE_drySnowAlb_south, SEAICE_wetSnowAlb_south, HO_south abb637800a Mart*0531 _RL SEAICE_cBasalStar, SEAICEbasalDragU0 0532 _RL SEAICEbasalDragK1, SEAICEbasalDragK2 5bb179ddc2 Mart*0533 _RL SEAICEsideDrag b58e51ce4e Jean*0534 _RL SEAICE_wetAlbTemp, SEAICE_waterAlbedo 0adbdb4edd Mart*0535 _RL SEAICE_strength, SEAICE_cStar, SEAICEpressReplFac c512e371cc drin*0536 _RL SEAICE_tensilFac, SEAICE_tensilDepth 0537 _RL SEAICE_eccen, SEAICE_eccfr 0538 _RL SEAICEmcMu, SEAICEtdMU 136908bfac Ian *0539 _RL SEAICE_lhFusion, SEAICE_lhEvap fff6be1885 Mart*0540 _RL SEAICE_dalton 0d2761b0d1 Mart*0541 _RL SEAICE_iceConduct, SEAICE_snowConduct 0542 _RL SEAICE_emissivity, SEAICE_ice_emiss, SEAICE_snow_emiss d778130a13 Mart*0543 _RL SEAICE_boltzmann 1c278edd09 Jean*0544 _RL SEAICE_snowThick, SEAICE_shortwave 0545 _RL SEAICE_saltFrac, SEAICE_salt0, SEAICEstressFactor 0546 _RL SEAICE_mcPheeTaper, SEAICE_mcPheePiston 0547 _RL SEAICE_frazilFrac, SEAICE_availHeatFrac 840c7fba30 Gael*0548 _RL facOpenGrow, facOpenMelt 0549 _RL SEAICE_tempFrz0, SEAICE_dTempFrz_dS 4b6d456764 Mart*0550 _RL SEAICE_PDF(nITD) 136908bfac Ian *0551 _RL OCEAN_drag, LSR_ERROR, DIFF1 79df32c3f1 Mart*0552 _RL SEAICEnonLinTol, JFNKres_t, JFNKres_tFac c8f5b4180e Mart*0553 _RL JFNKgamma_lin_min, JFNKgamma_lin_max, SEAICE_JFNKepsilon c704c5a1ef Mart*0554 _RL SEAICE_JFNK_lsGamma 1f3ad2d627 Mart*0555 _RL SEAICE_JFNKphi, SEAICE_JFNKalpha 3daf25222c Mart*0556 _RL SEAICE_deltaMin 425e8efc36 Jean*0557 _RL SEAICE_area_reg, SEAICE_hice_reg 136908bfac Ian *0558 _RL SEAICE_area_floor, SEAICE_area_max 16f85413ea Mart*0559 _RL SEAICE_airTurnAngle, SEAICE_waterTurnAngle 4ece3a93a6 Mart*0560 _RL SEAICE_elasticParm, SEAICE_evpTauRelax 2bf4ed5c39 Mart*0561 _RL SEAICE_evpAlpha, SEAICE_evpBeta bb06934ea0 Mart*0562 _RL SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac 1a0d5a451b Mart*0563 _RL SEAICEaEVPcoeff, SEAICEaEVPcStar, SEAICEaEVPalphaMin 204068be6a R. S*0564 _RL SEAICE_evpAreaReg, SEAICEpresH0 6d78fc5463 Gael*0565 _RL SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow e54fe3e1f9 Gael*0566 _RL SEAICEdiffKhSalt d32fe07ad8 Patr*0567 _RL SEAICE_tauAreaObsRelax 5fd900657f Mart*0568 _RL SEAICEgStar, SEAICEhStar, SEAICEaStar, SEAICEshearParm 0569 _RL SEAICEmuRidging, SEAICEmaxRaft, SEAICE_cf 3ac108d8cf Mart*0570 _RL SEAICEsnowFracRidge dc54d31829 Ian *0571 _RL SINegFac bb06934ea0 Mart*0572 809c36b928 Patr*0573 COMMON /SEAICE_PARM_RL/ e501eee760 Mart*0574 & SEAICE_deltaTtherm, SEAICE_deltaTdyn, e45202e340 Mart*0575 & SEAICE_LSRrelaxU, SEAICE_LSRrelaxV, 4ece3a93a6 Mart*0576 & SEAICE_deltaTevp, SEAICE_elasticParm, SEAICE_evpTauRelax, 2bf4ed5c39 Mart*0577 & SEAICE_evpAlpha, SEAICE_evpBeta, 1a0d5a451b Mart*0578 & SEAICEaEVPcoeff, SEAICEaEVPcStar, SEAICEaEVPalphaMin, bb06934ea0 Mart*0579 & SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac, 204068be6a R. S*0580 & SEAICE_evpAreaReg, SEAICEpresH0, b7411f1a84 Jean*0581 & SEAICE_monFreq, SEAICE_dumpFreq, a2d714f5b7 Jean*0582 & SEAICE_initialHEFF, 07b4a12853 Mart*0583 & SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR, 5867b94c2f Mart*0584 & SEAICE_drag, SEAICE_waterDrag, SEAICEdWatMin, 0585 & SEAICE_dryIceAlb, SEAICE_wetIceAlb, 0586 & SEAICE_drySnowAlb, SEAICE_wetSnowAlb, HO, f834b21bef Dimi*0587 & SEAICE_drag_south, SEAICE_waterDrag_south, 0588 & SEAICE_dryIceAlb_south, SEAICE_wetIceAlb_south, f81e0cd482 Dimi*0589 & SEAICE_drySnowAlb_south, SEAICE_wetSnowAlb_south, HO_south, abb637800a Mart*0590 & SEAICE_cBasalStar, SEAICEbasalDragU0, 0591 & SEAICEbasalDragK1, SEAICEbasalDragK2, 5bb179ddc2 Mart*0592 & SEAICEsideDrag, b58e51ce4e Jean*0593 & SEAICE_wetAlbTemp, SEAICE_waterAlbedo, c512e371cc drin*0594 & SEAICE_strength, SEAICE_cStar, SEAICE_eccen, SEAICE_eccfr, 0595 & SEAICEtdMU, SEAICEmcMu, 210ee8461e jm-c 0596 & SEAICEpressReplFac, SEAICE_tensilFac, SEAICE_tensilDepth, 425e8efc36 Jean*0597 & SEAICE_lhFusion, SEAICE_lhEvap, fff6be1885 Mart*0598 & SEAICE_dalton, SEAICE_cpAir, 0d2761b0d1 Mart*0599 & SEAICE_iceConduct, SEAICE_snowConduct, 0600 & SEAICE_emissivity, SEAICE_ice_emiss, SEAICE_snow_emiss, d778130a13 Mart*0601 & SEAICE_boltzmann, 1c278edd09 Jean*0602 & SEAICE_snowThick, SEAICE_shortwave, 0603 & SEAICE_saltFrac, SEAICE_salt0, SEAICEstressFactor, 0604 & SEAICE_mcPheeTaper, SEAICE_mcPheePiston, 0605 & SEAICE_frazilFrac, SEAICE_availHeatFrac, 6ec4646d60 Gael*0606 & facOpenGrow, facOpenMelt, 4b6d456764 Mart*0607 & SEAICE_tempFrz0, SEAICE_dTempFrz_dS, SEAICE_PDF, 425e8efc36 Jean*0608 & OCEAN_drag, LSR_ERROR, DIFF1, 79df32c3f1 Mart*0609 & SEAICEnonLinTol, JFNKres_t, JFNKres_tFac, c8f5b4180e Mart*0610 & JFNKgamma_lin_min, JFNKgamma_lin_max, SEAICE_JFNKepsilon, c704c5a1ef Mart*0611 & SEAICE_JFNK_lsGamma, SEAICE_JFNKphi, SEAICE_JFNKalpha, 3daf25222c Mart*0612 & SEAICE_deltaMin, SEAICE_area_reg, SEAICE_hice_reg, 136908bfac Ian *0613 & SEAICE_area_floor, SEAICE_area_max, 6d78fc5463 Gael*0614 & SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow, d32fe07ad8 Patr*0615 & SEAICEdiffKhSalt, SEAICE_tauAreaObsRelax, 5fd900657f Mart*0616 & SEAICE_airTurnAngle, SEAICE_waterTurnAngle, 0617 & SEAICEgStar, SEAICEhStar, SEAICEaStar, SEAICEshearParm, 0618 & SEAICEmuRidging, SEAICEmaxRaft, SEAICE_cf, c512e371cc drin*0619 & SINegFac, 0620 & SEAICEsnowFracRidge 809c36b928 Patr*0621 0622 C-- COMMON /SEAICE_BOUND_RL/ Various bounding values 425e8efc36 Jean*0623 C MIN_ATEMP :: minimum air temperature (deg C) 0624 C MIN_LWDOWN :: minimum downward longwave (W/m^2) 0625 C MIN_TICE :: minimum ice temperature (deg C) 3daf25222c Mart*0626 C SEAICE_EPS :: small number 425e8efc36 Jean*0627 C SEAICE_EPS_SQ :: small number square 809c36b928 Patr*0628 C 050eb90cc6 Gael*0629 _RL MIN_ATEMP, MIN_LWDOWN, MIN_TICE 460cb5f999 Dimi*0630 _RL SEAICE_EPS, SEAICE_EPS_SQ 809c36b928 Patr*0631 COMMON /SEAICE_BOUND_RL/ 050eb90cc6 Gael*0632 & MIN_ATEMP, MIN_LWDOWN, MIN_TICE, 45a6dc0d5a Dimi*0633 & SEAICE_EPS, SEAICE_EPS_SQ 809c36b928 Patr*0634 86b84a92fc Patr*0635 #ifdef SEAICE_ITD 210ee8461e jm-c 0636 C Hlimit :: ice thickness category limits (m), array of 4413db9420 Mart*0637 C size nITD+1 210ee8461e jm-c 0638 C Hlimit_c1,_c2,_c3 :: coefficients set in seaice_readparams.F to 4413db9420 Mart*0639 C calculate Hlimit in seaice_init_fixed.F 86b84a92fc Patr*0640 _RL Hlimit(0:nITD) 0641 _RL Hlimit_c1, Hlimit_c2, Hlimit_c3 0642 COMMON /SEAICE_BOUND_ITD_RL/ 0643 & Hlimit, 0644 & Hlimit_c1,Hlimit_c2,Hlimit_c3 0645 #endif /* SEAICE_ITD */ 0646 baa476eeba Dimi*0647 C-- Constants used by sea-ice model 0648 _RL ZERO , ONE , TWO 3de1dea93d Jean*0649 PARAMETER ( ZERO = 0.0 _d 0, ONE = 1.0 _d 0, TWO = 2.0 _d 0 ) baa476eeba Dimi*0650 _RL QUART , HALF 3de1dea93d Jean*0651 PARAMETER ( QUART = 0.25 _d 0, HALF = 0.5 _d 0 ) f50f58ec54 Gael*0652 _RL siEps 3de1dea93d Jean*0653 PARAMETER ( siEps = 1. _d -5 ) 9a87ea5e43 Ed H*0654 1ed503f8a3 Gael*0655 C-- Constants needed by McPhee formulas for turbulent ocean fluxes : 1c278edd09 Jean*0656 C Stanton number (dimensionless), typical friction velocity 1ed503f8a3 Gael*0657 C beneath sea ice (m/s), and tapering factor (dimensionless) 0658 _RL STANTON_NUMBER, USTAR_BASE, MCPHEE_TAPER_FAC 1c278edd09 Jean*0659 PARAMETER ( MCPHEE_TAPER_FAC = 12.5 _d 0 , STANTON_NUMBER = 1ed503f8a3 Gael*0660 & 0.0056 _d 0, USTAR_BASE = 0.0125 _d 0 ) 0661 47852c9c0c Mart*0662 C-- identifiers for advected properties f681b7f5d4 Dimi*0663 INTEGER GAD_HEFF,GAD_AREA,GAD_QICE1,GAD_QICE2,GAD_SNOW f50f58ec54 Gael*0664 INTEGER GAD_SALT,GAD_SITR 425e8efc36 Jean*0665 PARAMETER ( GAD_HEFF = 1, bb3037b865 Gael*0666 & GAD_AREA = 2, 0667 & GAD_SNOW = 3, 0668 & GAD_SALT = 4, ccaa3c61f4 Patr*0669 & GAD_QICE1 = 5, 0670 & GAD_QICE2 = 6, f50f58ec54 Gael*0671 & GAD_SITR = 7) 0f6ad6f523 Mart*0672 9a87ea5e43 Ed H*0673 CEH3 ;;; Local Variables: *** 0674 CEH3 ;;; mode:fortran *** 0675 CEH3 ;;; End: ***
| [ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
|
This page was automatically generated from https://github.com/MITgcm/MITgcm by the 2.2.1-MITgcm-0.1 LXR engine. The LXR team |
|