Spinning up NorESM

[mvertens]

The spinup process for NorESM involves running a fully coupled run and setting mediator namelist to create auxiliary history output that can then be ingested by data components to force the target prognostic components.

---

1. Creating mediator history files

• create_newcase commnad

 ./create_newcase --case <case> --compset 1850_CAM%DEV%LT%NORESM%CAMoslo_CLM60%FATES-SP_CICE_BLOM%ECO_MOSART_DGLC%NOEVOLVE_SWAV_SESP  \\
 --res ne30pg3_tn14 \\
 --project <project> --run-unsupported --mach betzy 

• user_nl_cpl

 histaux_atm2med_file1_enabled = .true.
 histaux_atm2med_file2_enabled = .true.
 histaux_atm2med_file3_enabled = .true.
 histaux_atm2med_file4_enabled = .true.
 histaux_atm2med_file5_enabled = .true.
 histaux_rof2med_file1_enabled = .true.
 histaux_atm2med_file5_ntperfile = 1
 histaux_atm2med_file5_history_n = 1
 histaux_atm2med_file5_history_option = 'ndays'
 histaux_rof2med_file1_ntperfile = 1
 histaux_rof2med_file1_history_n = 1
 histaux_rof2med_file1_history_option = 'ndays'
 histaux_l2x1yrg = .true.

Using the generated mediatory history files as input to DATM for forcing CLM, input to DATM and DROF for forcing BLOM and input to DLND for forcing CISM.

The output from the mediator history files is as follows:

• DATM output to spinup CLM and BLOM - that contains the following variables in the input file nuopc.runconfig

  • instantaneous fields

     histaux_atm2med_file1_auxname = atm.1h.inst
     histaux_atm2med_file1_doavg = .true.
     histaux_atm2med_file1_enabled = .false.
     histaux_atm2med_file1_flds = Faxa_swndr:Faxa_swvdr:Faxa_swndf:Faxa_swvdf
     histaux_atm2med_file1_history_n = 1
     histaux_atm2med_file1_history_option = nhours
     histaux_atm2med_file1_ntperfile = 24
          result in files
     <case>cpl.hx.atm.1h.inst.0001-01-01-03600.nc
     <case>cpl.hx.atm.1h.inst.0001-01-02-03600.nc
     ..... 
    

  • hourly averaged fields

     histaux_atm2med_file2_auxname = atm.1h.avrg
     histaux_atm2med_file2_doavg = .true.
     histaux_atm2med_file2_enabled = .false.
     histaux_atm2med_file2_flds = Sa_u:Sa_v
     histaux_atm2med_file2_history_n = 1
     histaux_atm2med_file2_history_option = nhours
     histaux_atm2med_file2_ntperfile = 24
          results in the files
     <case>.cpl.hx.atm.1h.avrg.0001-01-01-03600.nc
     <case>.cpl.hx.atm.1h.avrg.0001-01-02-03600.nc
      ....
    

  • 3-hourly averaged fields precipitation fields

    histaux_atm2med_file3_auxname = atm.3hprec.avrg
    histaux_atm2med_file3_doavg = .true.
    histaux_atm2med_file3_enabled = .true.
    histaux_atm2med_file3_flds = Faxa_rainc:Faxa_rainl:Faxa_snowc:Faxa_snowl
    histaux_atm2med_file3_history_n = 3
    histaux_atm2med_file3_history_option = nhours
    histaux_atm2med_file3_ntperfile = 8
        results in files
    <case>.cpl.hx.atm.3hprec.avrg.0001-01-01-10800.nc
    <case>.cpl.hx.atm.3hprec.avrg.0001-01-02-10800.nc
     ...
    

  • 3-hourly averaged fields other than precipitation

    histaux_atm2med_file4_auxname = atm.3h.avrg
    histaux_atm2med_file4_doavg = .true.
    histaux_atm2med_file4_enabled = .true.
    histaux_atm2med_file4_flds = Sa_z:Sa_topo:Sa_u:Sa_v:Sa_tbot:Sa_ptem:Sa_shum:Sa_dens:Sa_pbot:Sa_pslv:Faxa_lwdn:Faxa_rainc:Faxa_rainl:Faxa_snowc:Faxa_snowl:Faxa_swndr:Faxa_swvdr:Faxa_swndf:Faxa_swvdf:Sa_co2diag:Sa_co2prog
    histaux_atm2med_file4_history_n = 3
    histaux_atm2med_file4_history_option = nhours
    histaux_atm2med_file4_ntperfile = 8
         results in files 
    cpl.hx.atm.3h.avrg.0001-01-01-10800.nc
    cpl.hx.atm.3h.avrg.0001-01-02-10800.nc
    .....
    

  • daily averaged fields

    histaux_atm2med_file5_auxname = atm.24h.avrg
    histaux_atm2med_file5_doavg = .true.
    histaux_atm2med_file5_enabled = .false.
    histaux_atm2med_file5_flds = Faxa_bcph:Faxa_ocph:Faxa_dstwet:Faxa_dstdry:Sa_co2prog:Sa_co2diag
    histaux_atm2med_file5_history_n = 3
    histaux_atm2med_file5_history_option = nhours
    histaux_atm2med_file5_ntperfile = 2
           results in files
    <case>.cpl.hx.atm.24h.avrg.0001-01-02-00000.nc
    <case>.cpl.hx.atm.24h.avrg.0001-01-03-00000.nc
    <case>.cpl.hx.atm.24h.avrg.0003-01-01-00000.nc
    

• ROF output to spinup BLOM (daily averaged fields ) - that contains the following variables in the input file nuopc.runconfig

  histaux_rof2med_file1_auxname = rof.24h.avrg
  histaux_rof2med_file1_doavg = .true.
  histaux_rof2med_file1_enabled = .true.
  histaux_rof2med_file1_flds = all
  histaux_rof2med_file1_history_n = 1
  histaux_rof2med_file1_history_option = ndays
  histaux_rof2med_file1_ntperfile = 1
     results in files 
  <case>.cpl.hx.rof.24h.avrg.0002-01-01-00000.nc
  <case>.cpl.hx.rof.24h.avrg.0003-01-01-00000.nc
  .....
  

---

2. Spinning up CTSM using above auxiliary files

NOTE: this assumes that noresm2_5_alpha09 will be used - since changes to cdeps had to be made to have this working
Assumptions for this example:

1. spinup data is in directory /cluster/work/users/mvertens/noresm/spinup_tests/spinup_data
2. the casename for the spinup data is n1850.ne30_tn14.hybrid_fatessp_cplhist.202401031

> ./create_newcase --case <casepath> --compset 1850_DATM%CPLHIST_CLM60%BGC-CROP_SICE_SOCN_SROF_SGLC_SWAV --res ne30pg3_ne30pg3_mtn14 --project nn9039k --run-unsupported --mach betzy

> cd <casepath>

> ./xmlchange NTASKS=512
> ./xmlchange NTASKS_OCN=354
> ./xmlchange DATM_CPLHIST_CASE=n1850.ne30_tn14.hybrid_fatessp_cplhist.202401031
> ./xmlchange DATM_CPLHIST_DIR=/cluster/work/users/mvertens/noresm/spinup_tests/spinup_data
> ./xmlchange DATM_YR_START=1
> ./xmlchange DATM_YR_END=1
> ./xmlchange DATM_PRESNDEP=none

> ./case.setup
> ./caset.build
> ./case.submit

---

3. Spinning up BLOM using above auxiliary files

NOTE: this assumes that noresm2_5_alpha09 will be used - since changes to cdeps had to be made to have this working
Assumptions for this example:

1. spinup data is in dirctory /cluster/work/users/mvertens/noresm/spinup_tests/spinup_data
2. the casename for the spinup data is n1850.ne30_tn14.hybrid_fatessp_cplhist.202401031

>  ./create_newcase --case <casepath> --compset 1850_DATM%CPLHIST_SLND_CICE_BLOM%ECO_DROF%CPLHIST_SGLC_SWAV --res ne30pg3_tn14 --project <project> --run-unsupported --mach betzy

> cd <casepath>

> ./xmlchange NTASKS=512
> ./xmlchange NTASKS_OCN=354
> ./xmlchange DATM_CPLHIST_CASE=n1850.ne30_tn14.hybrid_fatessp_cplhist.202401031
> ./xmlchange DATM_CPLHIST_DIR=/cluster/work/users/mvertens/noresm/spinup_tests/spinup_data
> ./xmlchange DATM_YR_START=1
> ./xmlchange DATM_YR_END=1
> ./xmlchange DATM_PRESNDEP=clim_1850
> ./xmlchange DROF_CPLHIST_CASE=n1850.ne30_tn14.hybrid_fatessp_cplhist.202401031
> ./xmlchange DROF_CPLHIST_DIR=/cluster/work/users/mvertens/noresm/spinup_tests/spinup_data

> ./case.setup
> ./case.build
> ./case.submit

---

4. Spinning up CISM using above auxiliary files

NOTE: this assumes that noresm2_5_alpha09 will be used - since changes to cdeps had to be made to have this working
Assumptions for this example:

1. spinup data is in dirctory /cluster/work/users/mvertens/noresm/spinup_tests/spinup_data
2. the casename for the spinup data is n1850.ne30_tn14.hybrid_fatessp_cplhist.202401031

NOTE: changes must still be implemented in CDEPS/CMEPS in order to have docn multiple levels be written to auxiliary stream fiels

[adagj]

@mvertens I suggest to move this recipe to the wiki

[mvertens]

Right. I agree. My plan is to add the instructions for the second part here and then resolve any questions then move it to the wiki once everything was clear. Does that sound reasonable?