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