d_max, land only (ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_SturmSnowtk_dmax_1850)

[wwieder]

Description:
The coupled and land only PPE projects identified d_max (the dry surface layer (DSL) parameter) as a parameter that could help snow melt, although this maybe more important for coupled simulations...

It looks the the default value for d_max is 0.15m, and for this I think we want to use the high side change (0.6m) from the PPE? @linniahawkins or @aswann can you confirm this?

I'd build on the Sturm case, #10

• starting an AD spinup with initial conditions from the postAD Sturm run
• setting reseed=.true.

---

Case directory:
Locally:
/glade/work/oleson/cesm2_3_alpha16b/cime/scripts/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_SturmSnowtk_dmax_GSWP3V1_(1850AD, 1850pAD, 1850)

---

Sandbox:
Locally:
/glade/work/oleson/cesm2_3_alpha16b

---

user_nl_ changes:
hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS', 'H2OSNO','NDEP_TO_SMINN'
use_c13 = .true.
use_c14 = .true.
use_c13_timeseries = .true.
use_c14_bombspike = .true.
fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_ne30np4.pg3_SSP5-8.5_78pfts_CMIP6_1850-2100_c230227.nc'
finidat = '/glade/scratch/oleson/archive/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_SturmSnowtk_1850pAD/rest/0401-01-01-00000/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_GSWP3V1_SturmSnowtk_1850pAD.clm2.r.0401-01-01-00000.nc'
paramfile = '/glade/p/cgd/tss/people/oleson/modify_param/ctsm51_params.c211112.nc.dmax.nc' (d_max=60 cm)

---

SourceMods:

src.clm/atm2lndType.F90: Add some additional history fields
src.clm/clm_varpar.F90: nlevurb=5 -> nlevurb=10
src.clm/surfrdMod.F90: Remove read of PFTDATA_MASK
src.clm/Wateratm2lndType.F90: Add some additional history fields
src.clm/SoilTemperatureMod.F90 Sturm Snow tk

---

Diagnostics:
Diags (if available)

---

Output:
Output:
/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM52_DEV/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_SturmSnowtk_dmax_GSWP3V1_1850AD/
/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM52_DEV/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_SturmSnowtk_dmax_GSWP3V1_1850pAD/
/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM52_DEV/ctsm51_cesm23a16bctsm51d130_ne30pg3ne30pg3mg17_SturmSnowtk_dmax_GSWP3V1_1850/

---

Contacts:
@olyson
@wwieder

Extra details:
See discussion in #3
AD spunup to 1.88% disequilibrium, 160 years
pAD spunup to 4.78% disequilibrium, 440 years
Complete

[olyson]

I'll work on this when I'm back on Tues.

[aswann]

I confirmed with Claire that we used the values from Daniel's one at a time PPE, so whatever was used there is what we implemented in our coupled runs. We didn't save a separate version of that spreadsheet.

[olyson]