POP: How do I initializing POP2 with a spun-up initial condition?

The recommended method for initializing the CESM active ocean model (POP2) in a CESM startup case is to use the default settings; these initialize the ocean model from a state of rest. Occasionally, however, researchers are interested in initializing the ocean model from a "spun up" ocean condition that was generated from an existing CESM run. To accommodate their request, a nonstandard method of initializing the model was developed. It is called the startup_spunup option. The startup_spunup initialization is a research option that is designed for use by expert users only.

Because of the complex interactions between the ocean-model parameterizations used to generate the spun-up case, $CASE_SP, and those used in the new startup case, it is impossible to provide a single recommended spun-up ocean initial condition for all cases. Instead, researchers must carefully select an existing solution whose case conditions closely match those in the new case. A mismatch of options between the spun-up case and the new case can result in scientifically invalid solutions.

When a startup_spunup case is necessary, use this procedure:

  1. Create a new case, $CASE. By default, $CASE will be a "startup run."

    
> create_newcase -case ~/cesm1/EXAMPLE_CASEocn \
                     -mach bluefire \
                     -compset B20TR \ 
                     -res 0.9x1.25_gx1v6 
    > cd ~/cesm1/EXAMPLE_CASEocn
    > ./cesm_setup
    
  2. The ocean restart filename is of the form ${CASE_SP}.pop.r.$date, where $date is the model date of your spun-up dataset. If the ocean restart files were written in binary format, a companion ascii-formatted restart "header" file will also exist. The companion header file will have the same name as the restart file, except that it will have the suffix ".hdr" appended at the end of the filename. You must copy both the binary restart file and the header file to your data directory.

  3. The spun-up ocean restart and restart header files must be available to your new case. Copy them directly into $RUNDIR. It is critically important to copy both the binary restart file and its companion header file to the $RUNDIR.

    
> cp ${CASE_SP}.pop.r.$date       $RUNDIR 
    > cp ${CASE_SP}.pop.r.${date}.hdr $RUNDIR
    
  4. Redefine the ocean-model initial-condition dataset by editing user_nl_pop2 and add the following lines at the end of the file (enter the resolved string for ${CASE_SP}).

    
set init_ts_suboption = 'spunup'
    init_ts_file = '${CASE_SP}.pop.r.$date
    
    Note that the model will automatically look for the ${CASE_SP}.pop.r.${date}.hdr file in $RUNDIR.

  5. Build and run as usual.