The namelist variables for the ice model physics are listed in Table 3. restart is almost always true since most run types begin by reading in a binary restart file. See section 5 for a description of the run types and about using restart files and internally generated model data as initial conditions. kcolumn is a flag that will run the model as a single column if is set to 1. This option has not been thoroughly tested and is not supported.
The calculation of the ice velocities is subcycled ndte times per timestep so that the elastic waves are damped before the next timestep. The subcycling timestep is calculated as dte = dt/ndte and must be sufficiently smaller than the damping timescale T, which needs to be sufficiently shorter than dt.
| dte < T < dt | (2) |
This relationship is discussed in (Hunke(2001)); also see (Hunke and Lipscomb(2008)), section 4.4. The best ratio for [dte : T : dt] is [1 : 40 : 120]. Typical combinations of dt and ndte are (3600., 120), (7200., 240) (10800., 120). The default ndte is 120 as set in ice_init.F90.
kitd determines the scheme used to redistribute sea ice within the ice thickness distribution (ITD) as the ice grows and melts. The linear remapping scheme is the default and approximates the thickness distribution in each category as a linear function ((Lipscomb(2001))). The delta function method represents g(h) in each category as a delta function ((Bitz et al.(2001))). This method can leave some categories mostly empty at any given time and cause jumps in the properties of g(h).
kdyn determines the ice dynamics used in the model. The default is the elastic-viscous-plastic (EVP) dynamics (Hunke and Dukowicz(1997)). If kdyn is set to o 0, the ice dynamics is inactive. In this case, ice velocities are not computed and ice is not transported. Since the initial ice velocities are read in from the restart file, the maximum and minimum velocities written to the log file will be non-zero in this case, but they are not used in any calculations.
The value of kstrength determines which formulation is used to calculate the strength of the pack ice. The (Hibler(1979)) calculation depends on mean ice thickness and open water fraction. The calculation of (Rothrock(1975)) is based on energetics and should not be used if the ice that participates in ridging is not well resolved.
evp_damping is used to control the damping of elastic waves in the ice dynamics. It is typically set to .true. for high-resolution simulations where the elastic waves are not sufficiently damped out in a small timestep without a significant amount of subcycling. This procedure works by reducing the effective ice strength that's used by the dynamics and is not a supported option.
advection determines the horizontal transport scheme used. The default scheme is the incremental remapping method ((Lipscomb and Hunke(2004))). This method is less diffusive and is computationally efficient for large numbers of categories or tracers. The upwind scheme is also available. The upwind scheme is only first order accurate.
The base values of the snow and ice albedos for the CCSM3 shortwave option are set in the namelist. The ice albedos are those for ice thicker than ahmax, which is currently set at 0.5 m. This thickness is a parameter that can be changed in ice_shortwave.F90. The snow albedos are for cold snow.
For the new delta-Eddington shortwave radiative transfer scheme (Briegleb and Light(2007)), the base albedos are computed based on the inherent optical properties of snow, sea ice, and melt ponds. These albedos are tunable through adjustments to the snow grain radius, R_snw, temperature to transition to melting snow, and maximum snow grain radius.