7 Changing the Number of Ice Thickness Categories

The number of ice thickness categories affects ice model input files in three places:

• $NCAT in the run script
• The input template files ice_model_size.F.$RESx$NCAT in $CSIMDIR/input_templates
• The initial condition (restart) file in the input file directory

The number of ice thickness categories is set in ccsm3/scripts/$CASE/Buildexe/csim.buildexe.csh (coupled) or csim_run (uncoupled) using the variable called $NCAT. The default value is 5 categories. $NCAT is used to determine which module ice_model_size.F.$RESx$NCAT is copied to /obj before compilation where it is renamed ice_model_size.F. $RES is the resolution of the grid, 100x116 (gx3v5) and 320x384 (gx1v3) for low and high resolution grids, respectively. The input templates are located in $CSIMDIR/input_templates for 1, 3, 5 and 10 thickness categories on the default grid sizes of 100x116 and 320x384.

NOTE

To use one ice thickness category, the following changes will need to be made in the namelist:

  , kitd          = 0
  , kstrength     = 0

With these settings, the model will use the delta scheme instead of linear remapping and a strength parameterization based on open water area and mean ice thickness.

To create an input template with a number of categories other than the above values, copy an existing template that has the desired grid. Change ncat to the appropriate number of categories.

The information in the initial restart file is dependent on the number of ice thickness categories and the total number of layers in the ice distribution. An initial condition file exists only for the default case of 5 ice thickness categories, with four layers in each category. To create an initial condition file for a different number of categories or layers, these steps should be followed:

• Set $NCAT to the desired number of categories in csim_run (uncoupled) or
  ccsm3/scripts/$CASE/Buildexe/csim.buildexe.csh (coupled).
• Set the namelist variable dumpfreq = 'm' in csim.setup.csh (uncoupled) or
  ccsm3/scripts/$CASE/Buildnml_Prestage/csim.buildnml_prestage.csh
  (coupled) to print out restart files monthly.
• Set the namelist variable restart = .false. in csim.setup.csh (uncoupled) or
  ccsm3/scripts/$CASE/Buildnml_Prestage/csim.buildnml_prestage.csh
  (coupled) to use the initial conditions within the ice model.
• Create an input template (or use an existing one) with the desired number of categories.
• Run the model to equilibrium.
• The last restart file can be used as an initial condition file.
• Change the name of the last restart file to iced.0001-01-01.$GRID.
• Copy the file into the input data directory or directly into the /rest directory under the executable directory.

Note that the date printed inside the binary restart file will not be the same as 0001-01-01. For coupled runs, $BASEDATE will be the starting date and the date inside the file will not be used. For uncoupled runs, the starting date may have to be set in the code or changed in the restart file.