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