Land/ocean mask

Raising sea level

Re-assigning land to ocean is complicated process because it requires setting land points to ocean, which then must be initialized.  The ocean requires that all points be initialized at start up.  If you configure your experiment as a startup run, the new ocean cells will automatically be filled with Levitus (Levitus et al., 1998, Steel et al., 2001) values.  However, it is generally desirable to initialize the ocean from a branch or hybrid control experiment, because of the long time scales required to spin up the deep ocean.  In order to use a previous ocean restart file, you will need to modify the ocean source code to initialize the temperature and salinity fields for the new ocean cells.  This process can be complex and is beyond the scope of this document.

Lowering sea level

Changing ocean points to land is more straightforward than removing land points, but it still requires multiple steps.   

  1. Change the ocean bathymetry (KMT) file.   
  2. Remake the coupler mapping files that drive communication within CESM1.   
  3. Create a new initialization (CAMI) file for CAM.
  4. Create a new topography file (topo_bnd) file for CAM. 
  5. Create a new surface dataset:
    1. Create 0.5o resolution raw ‘mksrf’ datafiles that are consistent with your land/sea configuration. 
    2. Create mapping files for the raw mksrf datafiles that reflect the new land/sea mask for mapping the modified raw datasets to the new surface dataset.
  6.  Additional steps may be required if you are restarting from a previous case, and you have changed land cover (e.g., added land ice, changed vegetation distributions, etc,).
  7. In contrast to the case of simulated sea level rise, sea level lowering simulations can generally restart the ocean and sea ice models using restart files from present-day simulations, as long as no new ocean points have been added.


change_kmt.ncl (for Quaternary applications)

The NCL script change_kmt.ncl can be used to reassign ocean points to land.  For example, it is a straightforward way to remove the ocean cells representing Hudson Bay for Pre-Pleistocene simulations.  Or you may use change_kmt.ncl to modify the land/sea mask globally (e.g., to extend the continental margins for sea level lowstands).  The NCL code requires the default binary CESM1 topography (KMT) and region mask files as input, and produces two new binary files with the user-defined changes to the ocean map. The new binary files are used by the ocean model (user_nl_pop2) and are input to the paleo tool, mk_remap.csh, to create new coupler mapping files that map the new land/ocean map to the atmosphere. 

Hint:  Check the region mask to be sure you have not eliminated any ocean regions.  Eliminating ocean regions requires renumbering the ocean region in the region mask file and changing the region_ids file (e.g., gx1v6_region_ids) correspondingly.  By retaining at least one active ocean cell in the default regions, you can avoid having to change the region_ids.

Hint:  The NCL code we provide is a template for making minor changes to the default binary KMT map.  For best ocean model results, carefully examine your new KMT ocean/land mask for newly emergent islands in the Pacific, open closed basins or change your marginal seas, and widen or eliminate narrow channels.   These changes will require hand-editing of the change_kmt.ncl code to fix these problems.  We recommend using ncview on the netCDF output file from change_kmt.ncl to identify cells you would like to modify.