1850 with ctsm5.3 fsurdat, excess ice ON, A/C ON, and modified Sturm (ctsm52018_f09_xsiceON_acON_SturmExp1_[AD, SASU, pSASU, hist])

[slevis-lmwg]

Description:
Same as #62 but with modified Sturm.

• 1850 AD DONE
• 1850 SASU DONE
• 1850 pSASU DONE
• 1850-2023 DONE

---

Case directory:
Locally:
/glade/work/oleson/ctsm5.2.018/cime/scripts/ctsm52018_f09_xsiceON_acON_SturmExp1_AD
/glade/work/oleson/ctsm5.2.018/cime/scripts/ctsm52018_f09_xsiceON_acON_SturmExp1_SASU
/glade/work/oleson/ctsm5.2.018/cime/scripts/ctsm52018_f09_xsiceON_acON_SturmExp1_pSASU
/glade/work/oleson/ctsm5.2.018/cime/scripts/ctsm52018_f09_xsiceON_acON_SturmExp1_hist

---

Sandbox:
Locally:
/glade/work/oleson/ctsm5.2.018
git describe returns ctsm5.2.018

---

user_nl_ changes:
AD

urban_explicit_ac = .true.
use_excess_ice = .true.

fsurdat='/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_0.9x1.25_hist_1850_78pfts_c240807.nc'

! turn off shifting cultivation
do_grossunrep = .false.

! do this for spinup
reseed_dead_plants = .true.

! do this for isotopes
use_c13 = .true.
use_c14 = .true.
use_c13_timeseries = .true.
use_c14_bombspike = .true.

SASU also includes missing excess ice stream stuff and updated nnegcrit

stream_fldfilename_exice = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/paramdata/exice_init_0.125x0.125_c20220516.nc'
stream_mapalgo_exice = 'bilinear'
stream_meshfile_exice = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/paramdata/exice_init_0.125x0.125_ESMFmesh_cdf5_c20220802.nc'
use_excess_ice_streams = .true.

finidat = '/glade/derecho/scratch/oleson/archive/ctsm52018_f09_xsiceON_acON_SturmExp1_AD/rest/0301-01-01-00000/ctsm52018_f09_xsiceON_acON_SturmExp1_AD.clm2.r.0301-01-01-00000.nc'

nnegcrit = -7.d0

pSASU includes

finidat = '/glade/derecho/scratch/oleson/archive/ctsm52018_f09_xsiceON_acON_SturmExp1_SASU/rest/0301-01-01-00000/ctsm52018_f09_xsiceON_acON_SturmExp1_SASU.clm2.r.0301-01-01-00000.nc'

hist vs pSASU

<  co2_type = 'diagnostic'
---
>  co2_ppmv = 284.7
>  co2_type = 'constant'
13c14
<  finidat = '/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm52018_f09_xsiceON_acON_SturmExp1_pSASU.clm2.r.0181-01-01-00000.nc'
---
>  finidat = '/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm52018_f09_xsiceON_acON_SturmExp1_SASU.clm2.r.0301-01-01-00000.nc'
26,34d26
<  hist_dov2xy(1) = .true.
<  hist_dov2xy(2) = .false.
<  hist_dov2xy(3) = .false.
<  hist_dov2xy(4) = .true.
<  hist_dov2xy(5) = .false.
<  hist_fexcl1 = 'PCT_GLC_MEC', 'C14_SOILC_vr', 'SOM_ACT_C_vr', 'SOM_ACT_N_vr', 'SOM_SLO_C_vr', 'SOM_SLO_N_vr', 'SOM_PAS_C_vr',
<          'SOM_PAS_N_vr', 'SOILC_vr', 'SOILN_vr', 'CWD_C_vr', 'LIT_MET_C_vr', 'LIT_CEL_C_vr', 'LIT_LIG_C_vr',
<          'LIT_MET_N_vr', 'LIT_CEL_N_vr', 'LIT_LIG_C_vr', 'CWD_N_vr', 'SMIN_NO3_vr', 'CONC_O2_UNSAT', 'CONC_O2_SAT',
<          'SMIN_NH4_vr', 'SMINN_vr', 'PCT_CFT'
36,73d27
<  hist_fincl1 = 'FSDSND', 'FSDSVD', 'FSDSNI', 'FSDSVI', 'RAIN_FROM_ATM', 'SNOW_FROM_ATM', 'FLDS_NOT_DOWNSCALED',
<          'ZBOT', 'Tair_from_atm', 'Thair_from_atm', 'QBOT_NOT_DOWNSCALED', 'PBOT_NOT_DOWNSCALED', 'Rho_from_atm', 'PCO2',
<          'UWIND', 'VWIND', 'BCPHIDRY', 'BCPHODRY', 'BCPHIWET', 'OCPHIDRY', 'OCPHODRY',
<          'OCPHIWET', 'DSTWET1', 'DSTDRY1', 'DSTWET2', 'DSTDRY2', 'DSTWET3', 'DSTDRY3',
<          'DSTWET4', 'DSTDRY4', 'ATM_TOPO', 'EFLX_LH_TOT_ICE', 'FIRE_ICE', 'FLDS_ICE', 'FSH_ICE',
<          'FSNO_ICE', 'FSR_ICE', 'QFLX_SOLIDEVAP_FROM_TOP_LAYER_ICE', 'QRUNOFF_ICE', 'QSNOFRZ_ICE', 'QSNOMELT_ICE', 'RAIN_ICE',
<          'SNOW_ICE', 'SNOWICE_ICE', 'SNOWLIQ_ICE', 'SNOTXMASS_ICE', 'TG_ICE', 'TOPO_COL_ICE', 'TSA_ICE',
<          'TSOI_ICE', 'LEAFC_TO_LITTER', 'FROOTC_TO_LITTER', 'LIT_MET_C_TO_SOM_ACT_C', 'LIT_MET_N_TO_SOM_ACT_N', 'LIT_CEL_C_TO_SOM_ACT_C', 'LIT_CEL_N_TO_SOM_ACT_N',
<          'LIT_LIG_C_TO_SOM_SLO_C', 'LIT_LIG_N_TO_SOM_SLO_N', 'DWT_WOOD_PRODUCTC_GAIN_PATCH'
<  hist_fincl2 = 'TLAI', 'TSA', 'TREFMNAV', 'TREFMXAV', 'BTRANMN', 'Vcmx25Z', 'FSH', 'VEGWP', 'FCTR', 'FCEV', 'FGEV',
<          'FIRE', 'FSR', 'FIRA', 'FSA', 'GSSUNLN', 'GSSHALN', 'TSKIN', 'GPP', 'NPP', 'AGNPP', 'TOTVEGC',
<          'NPP_NUPTAKE', 'AR', 'HR', 'HTOP', 'GRAINC_TO_FOOD', 'NFERTILIZATION'
<  hist_fincl3 = 'FSR', 'H2OSNO', 'Q2M', 'SNOWDP', 'TSA', 'TREFMNAV', 'TREFMXAV', 'TG', 'QRUNOFF',
<          'FSH', 'FIRE', 'FIRA', 'FGR', 'EFLX_LH_TOT', 'RH2M', 'TLAI', 'SOILWATER_10CM', 'TOTSOILLIQ',
<          'TOTSOILICE', 'U10', 'TSOI_10CM', 'QIRRIG_DRIP', 'URBAN_HEAT', 'WASTEHEAT', 'TSKIN', 'GPP', 'NPP',
<          'AR', 'HR', 'DWT_CONV_CFLUX_PATCH', 'WOOD_HARVESTC', 'DWT_WOOD_PRODUCTC_GAIN_PATCH', 'SLASH_HARVESTC', 'COL_FIRE_CLOSS', 'DWT_SLASH_CFLUX_PATCH', 'FROOTC:I',
<          'HTOP', 'GRAINC_TO_FOOD'
<  hist_fincl4 = 'PCT_GLC_MEC', 'QICE_FORC', 'TSRF_FORC', 'TOPO_FORC', 'PCT_NAT_PFT', 'PCT_LANDUNIT', 'FSNO_ICE',
<          'C14_SOILC_vr', 'SOILC_vr', 'SOILN_vr', 'CWD_C_vr', 'LIT_MET_C_vr', 'LIT_CEL_C_vr', 'LIT_LIG_C_vr',
<          'LIT_MET_N_vr', 'LIT_CEL_N_vr', 'LIT_LIG_N_vr', 'CWD_N_vr', 'TOTLITC:I', 'TOT_WOODPRODC:I', 'TOTSOMC:I',
<          'TOTVEGC:I', 'PCT_CFT', 'CROPPROD1C:I'
<  hist_fincl5 = 'TOTSOMC:I', 'TOTSOMC_1m:I', 'TOTECOSYSC:I', 'TOTVEGC:I', 'WOODC:I', 'TOTLITC:I', 'LIVECROOTC:I',
<          'DEADCROOTC:I', 'FROOTC:I'
<  hist_mfilt(1) = 1
<  hist_mfilt(2) = 1
<  hist_mfilt(3) = 1
<  hist_mfilt(4) = 1
<  hist_mfilt(5) = 1
<  hist_nhtfrq(1) = 0
<  hist_nhtfrq(2) = 0
<  hist_nhtfrq(3) = 0
<  hist_nhtfrq(4) = -8760
<  hist_nhtfrq(5) = -8760
<  hist_type1d_pertape(1) = ' '
<  hist_type1d_pertape(2) = ' '
<  hist_type1d_pertape(3) = 'LAND'
<  hist_type1d_pertape(4) = ' '
<  hist_type1d_pertape(5) = 'LAND'
75c29
<  irrigate = .true.
---
>  irrigate = .false.
118d71
<  use_init_interp = .true.
133,134c86
<  model_year_align_ndep = 1850
<  ndep_taxmode = 'extend'
---
>  ndep_taxmode = 'cycle'
137c89
<  stream_fldfilename_ndep = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/ndepdata/fndep_clm_hist_b.e21.BWHIST.f09_g17.CMIP6-historical-WACCM.ensmean_1849-2015_monthly_0.9x1.25_c180926.nc'
---
>  stream_fldfilename_ndep = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/ndepdata/fndep_clm_WACCM6_CMIP6piControl001_y21-50avg_1850monthly_0.95x1.25_c180802.nc'
140c92
<  stream_year_last_ndep = 2015
---
>  stream_year_last_ndep = 1850
143d94
<  model_year_align_popdens = 1850
148c99
<  stream_year_last_popdens = 2016
---
>  stream_year_last_popdens = 1850
151d101
<  model_year_align_urbantv = 1850
155c105
<  stream_year_last_urbantv = 2106
---
>  stream_year_last_urbantv = 1850
215,220d164
<  do_harvest = .true.
<  do_transient_crops = .true.
<  do_transient_lakes = .true.
<  do_transient_pfts = .true.
<  do_transient_urban = .true.
<  flanduse_timeseries = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/landuse.timeseries_0.9x1.25_hist_1700-2023_78pfts_c240731.nc'
229d172
<  check_dynpft_consistency = .false.
232d174
<  init_interp_method = 'general'
378a321
>  reseed_dead_plants = .true.
402a346,349
>  stream_fldfilename_exice = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/paramdata/exice_init_0.125x0.125_c20220516.nc'
>  stream_mapalgo_exice = 'bilinear'
>  stream_meshfile_exice = '/glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/paramdata/exice_init_0.125x0.125_ESMFmesh_cdf5_c20220802.nc'
>  use_excess_ice_streams = .true.

user_nl_datm_streams same as in #60 and #61.

---

SourceMods:

/glade/work/oleson/ctsm5.2.018/cime/scripts/ctsm52018_f09_xsiceON_acON_SturmExp1_(AD)(SASU)(pSASU)(hist)/SourceMods/src.clm/SoilTemperatureMod.F90: Modifications to split the difference between the snow thermal conductivity curves of Sturm1997 and Jordan1991. See the discussion in #10 beginning July 25, 2024.

---

Diagnostics:
Diags (if available)
ILAMB:
https://webext.cgd.ucar.edu/I20TR/ctsm52018_f09_xsiceON_acON_SturmExp1_hist/lnd/_build_ctsm52018_f09_xsiceON_SturmExp1/
Standard diagnostics:
https://webext.cgd.ucar.edu/I20TR/ctsm52018_f09_xsiceON_acON_SturmExp1_hist/lnd/ctsm52018_f09_xsiceON_acON_SturmExp1_hist_2004_2023-ctsm52018_f09_xsiceON_acON_hist_2004_2023/setsIndex.html

---

Output:
Output (if still available):
/glade/derecho/scratch/oleson/archive/ctsm52018_f09_xsiceON_acON_SturmExp1_AD/
/glade/derecho/scratch/oleson/archive/ctsm52018_f09_xsiceON_acON_SturmExp1_SASU
/glade/derecho/scratch/oleson/archive/ctsm52018_f09_xsiceON_acON_SturmExp1_pSASU
/glade/derecho/scratch/oleson/archive/ctsm52018_f09_xsiceON_acON_SturmExp1_hist

---

Contacts:
@olyson

---

Extra details:
Sequence of simulations:
1850AD 300 yrs --> 1850SASU 300 yrs --> 1850pSASU 180 yrs --> hist 1850-2023

AD spunup to 3.80% disequilibrium, 300 years
SASU spunup to 13.05% disequilibrium, 300 years
pSASU spunup to 2.91% disequilibrium, 180 years