133 |
C SEAICE_areaGainFormula :: selects formula for ice cover gain from open water growth |
C SEAICE_areaGainFormula :: selects formula for ice cover gain from open water growth |
134 |
C :: 1=from growth by ATM |
C :: 1=from growth by ATM |
135 |
C :: 2=from predicted growth by ATM |
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 |
C SEAICE_multDim :: number of ice categories |
142 |
C SEAICE_debugPointI :: I,J index for seaice-specific debuggin |
C SEAICE_debugPointI :: I,J index for seaice-specific debuggin |
143 |
C SEAICE_debugPointJ |
C SEAICE_debugPointJ |
155 |
INTEGER SEAICEadjMODE |
INTEGER SEAICEadjMODE |
156 |
INTEGER SEAICE_areaLossFormula |
INTEGER SEAICE_areaLossFormula |
157 |
INTEGER SEAICE_areaGainFormula |
INTEGER SEAICE_areaGainFormula |
158 |
|
INTEGER SEAICEetaZmethod |
159 |
INTEGER SEAICE_multDim |
INTEGER SEAICE_multDim |
160 |
INTEGER SEAICE_debugPointI |
INTEGER SEAICE_debugPointI |
161 |
INTEGER SEAICE_debugPointJ |
INTEGER SEAICE_debugPointJ |
176 |
& SEAICE_areaLossFormula, |
& SEAICE_areaLossFormula, |
177 |
& SEAICE_areaGainFormula, |
& SEAICE_areaGainFormula, |
178 |
& SEAICE_multDim, |
& SEAICE_multDim, |
179 |
|
& SEAICEetaZmethod, |
180 |
& SEAICE_debugPointI, |
& SEAICE_debugPointI, |
181 |
& SEAICE_debugPointJ |
& SEAICE_debugPointJ |
182 |
|
|
205 |
C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model. |
C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model. |
206 |
C SEAICE_deltaTtherm :: Seaice timestep for thermodynamic equations (s) |
C SEAICE_deltaTtherm :: Seaice timestep for thermodynamic equations (s) |
207 |
C SEAICE_deltaTdyn :: Seaice timestep for dynamic solver (s) |
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) |
C SEAICE_deltaTevp :: Seaice timestep for EVP solver (s) |
210 |
C SEAICE_elasticParm :: parameter that sets relaxation timescale |
C SEAICE_elasticParm :: parameter that sets relaxation timescale |
211 |
C tau = SEAICE_elasticParm * SEAICE_deltaTdyn |
C tau = SEAICE_elasticParm * SEAICE_deltaTdyn |
308 |
C sea ice concentration (s), default=unset |
C sea ice concentration (s), default=unset |
309 |
C |
C |
310 |
_RL SEAICE_deltaTtherm, SEAICE_deltaTdyn, SEAICE_deltaTevp |
_RL SEAICE_deltaTtherm, SEAICE_deltaTdyn, SEAICE_deltaTevp |
311 |
|
_RL SEAICE_LSRrelaxU, SEAICE_LSRrelaxV |
312 |
_RL SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq |
_RL SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq |
313 |
_RL SEAICE_initialHEFF |
_RL SEAICE_initialHEFF |
314 |
_RL SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR |
_RL SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR |
346 |
|
|
347 |
COMMON /SEAICE_PARM_RL/ |
COMMON /SEAICE_PARM_RL/ |
348 |
& SEAICE_deltaTtherm, SEAICE_deltaTdyn, |
& SEAICE_deltaTtherm, SEAICE_deltaTdyn, |
349 |
|
& SEAICE_LSRrelaxU, SEAICE_LSRrelaxV, |
350 |
& SEAICE_deltaTevp, SEAICE_elasticParm, SEAICE_evpTauRelax, |
& SEAICE_deltaTevp, SEAICE_elasticParm, SEAICE_evpTauRelax, |
351 |
& SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac, |
& SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac, |
352 |
& SEAICEpresH0, |
& SEAICEpresH0, |
387 |
C MIN_TICE :: minimum ice temperature (deg C) |
C MIN_TICE :: minimum ice temperature (deg C) |
388 |
C SEAICE_EPS :: small number used to reduce derivative singularities |
C SEAICE_EPS :: small number used to reduce derivative singularities |
389 |
C SEAICE_EPS_SQ :: small number square |
C SEAICE_EPS_SQ :: small number square |
|
CToM<<< |
|
|
C Hlimit :: ice thickness category limits (m), array of size nITD+1 |
|
|
C Hlimit_c1,_c2,_c3 :: coefficients set in seaice_readparams.F to calculate Hlimit in seaice_init_fixed.F |
|
|
C>>>ToM |
|
390 |
C |
C |
391 |
_RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE |
_RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE |
392 |
_RL SEAICE_EPS, SEAICE_EPS_SQ |
_RL SEAICE_EPS, SEAICE_EPS_SQ |
|
CToM<<< |
|
|
_RL Hlimit(0:nITD) |
|
|
_RL Hlimit_c1, Hlimit_c2, Hlimit_c3 |
|
|
C>>>ToM |
|
393 |
COMMON /SEAICE_BOUND_RL/ |
COMMON /SEAICE_BOUND_RL/ |
394 |
& MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE, |
& MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE, |
395 |
& SEAICE_EPS, SEAICE_EPS_SQ |
& SEAICE_EPS, SEAICE_EPS_SQ |
396 |
|
|
397 |
|
#ifdef SEAICE_ITD |
398 |
CToM<<< |
CToM<<< |
399 |
C & SEAICE_EPS, SEAICE_EPS_SQ |
C Hlimit :: ice thickness category limits (m), array of size nITD+1 |
400 |
& SEAICE_EPS, SEAICE_EPS_SQ, |
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, |
& Hlimit, |
405 |
& Hlimit_c1,Hlimit_c2,Hlimit_c3 |
& Hlimit_c1,Hlimit_c2,Hlimit_c3 |
406 |
C>>>ToM |
C>>>ToM |
407 |
|
#endif /* SEAICE_ITD */ |
408 |
|
|
409 |
C-- Constants used by sea-ice model |
C-- Constants used by sea-ice model |
410 |
_RL ZERO , ONE , TWO |
_RL ZERO , ONE , TWO |