wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_ [AD, SASU, pSASU]

[wwieder]

Description:
Flat 1 degree case for WIEMIP spinup

---

Case directory:
Locally:
/glade/work/wwieder/wiemip/wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_AD

---

Sandbox:
Locally:
/glade/u/home/wwieder/CTSM (wetlands)
branch_tags/wetlands.n01.ctsm5.4.029

./create_newcase --case /glade/work/wwieder/wiemip/wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_AD --compset I1850Clm60BgcCrujra --res f09_t232 --run-unsupported --user-mods-dirs /glade/u/home/wwieder/CTSM/cime_config/usermods_dirs/clm/WIEMIP

---

user_nl_ changes:

fsurdat = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/surfdata_esmf/ctsm5.4.0/surfdata_0.9x1.25_hist_1850_78pfts_WIEMIP_c260408.nc'
paramfile = '/glade/work/wwieder/modified_inputs/ctsm60_params.genericCropFix.ch4cal_c.c260410.nc'
use_c13 = .false.
use_c14 = .false.

---

SourceMods:

1. AD

./xmlchange MOSART_MODE=NULL
./xmlchange CLM_ACCELERATED_SPINUP=on,CLM_FORCE_COLDSTART=on
./xmlchange CCSM_CO2_PPMV=280.
./xmlchange DATM_YR_START=1850,DATM_YR_END=1869
./xmlchange STOP_N=80,STOP_OPTION=nyears
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=3
./xmlchange JOB_PRIORITY=regular
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=12:00:00

2. SASU

export REFCASE=wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_AD

./xmlchange MOSART_MODE=NULL
./xmlchange RUN_TYPE=hybrid
./xmlchange STOP_N=80,REST_N=40,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=sasu,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=3
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFDATE=0321-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0321-01-01-00000
./xmlchange RUN_REFCASE=$REFCASE
./xmlchange RUN_REFDIR=/glade/derecho/scratch/wwieder/archive/$REFCASE/rest/0321-01-01-00000
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=12:00:00

3. pSASU

export REFCASE=wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_SASU

/xmlchange MOSART_MODE=NULL
./xmlchange STOP_N=80,STOP_OPTION=nyears,REST_N=40
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=2
./xmlchange CLM_ACCELERATED_SPINUP=off,CLM_FORCE_COLDSTART=off
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFDATE=0321-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0321-01-01-00000
./xmlchange RUN_REFCASE=$REFCASE
./xmlchange RUN_REFDIR=/glade/derecho/scratch/$USER/archive/$REFCASE/rest/0321-01-01-00000
./xmlchange JOB_PRIORITY=premium
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=10:00:00

---

Diagnostics:
Diags (if available)

https://webext.cgd.ucar.edu/I1850/$CASE/lnd/

---

Output:
Output (if still available):
/glade/scratch/$USER/archive/$CASE/

---

Contacts:
@wwieder

---

Extra details:

[wwieder]

Also started a 158_test run to look at DATM data we're getting from WIEMIP.
/glade/work/wwieder/wiemip/wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_climoTest

./create_clone --clone /glade/work/wwieder/wiemip/wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_AD \
  --case /glade/work/wwieder/wiemip/wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_climoTest \
  --keepexe

finidat='/glade/derecho/scratch/wwieder/archive/wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_AD/rest/0081-01-01-00000/wetlands.n01.ctsm5.4.029_f09noCrop_flat_158_AD.clm2.r.0081-01-01-00000.nc'

./xmlchange STOP_N=20,RESUBMIT=0,CONTINUE_RUN=FALSE
./xmlchange CLM_ACCELERATED_SPINUP=off,CLM_FORCE_COLDSTART=off
./xmlchange JOB_PRIORITY=premium
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=03:30:00

Diags that compare DATM with another 1 deg CRUJRA pSASU case (#147).

[wwieder]

It's kind of tricky to compare DATM from:

• WIEMIP (1850-1869) from TRENDY2024 vs.
• CRUJRA (1901-1920), from TRENDY2025* which also data replaced over Antarctica and Greenland.

Some variables seem to be just:

• Out of phase (e.g. TBOT, PREC), which might be expected but others seem to have
• Different baselines (PBOT, QBOT, WIND),

[wwieder]

Also using this flat case to test out BGC and RAD usermods

• /glade/u/home/wwieder/wetlands/cime_config/usermods_dirs/clm/WIEMIP/BGC

• /glade/u/home/wwieder/wetlands/cime_config/usermods_dirs/clm/WIEMIP/RAD

• For both cases what do we do for aerosols and ozone?

• Transient N dep with BGC, but cycles in RAD case?

• Transient PopDens? Urban?

• Fixed LULCC

• !use_init_interp = .true., incompatible with hillslope hydrology?

./create_newcase --case /glade/work/wwieder/wiemip/wetlands.n03.ctsm5.4.029_f09noCrop_flat_158_testBGC \
  --compset HIST_DATM%CRUJRA2024b_CLM60%BGC_SICE_SOCN_MOSART_SGLC_SWAV \
  --res f09_t232 --run-unsupported \
  --user-mods-dirs /glade/u/home/wwieder/wetlands/cime_config/usermods_dirs/clm/WIEMIP/BGC

./xmlchange STOP_N=5,RESUBMIT=0,STOP_OPTION=ndays
./xmlchange RUN_TYPE=startup
./xmlchange JOB_PRIORITY=premium
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=00:10:00

[slevis-lmwg]

@wwieder I noticed the following (may or may not be of concern):

1. In BGC and RAD, should co2 align be 1850 and last be 1869 or 2000?

co2tseries.cmip6_20tr:year_first=1850
co2tseries.cmip6_20tr:year_last=1850
co2tseries.cmip6_20tr:year_align=1

2. I have one or two more comments but glade has frozen on me...