1-deg I-case spin-up and historical to test new CLM6 paramfile (ctsm53026_BNF_[AD, SASU, pSASU, hist])

[linniahawkins]

Description:

Testing a new parameter set for CESM3. This is nearly identical to #89 with two modifications (1) We reverted needleleaf evergreen boreal trees (PFT2) to the CLM5 default tuning in order to improve survivability and (2) revert C3 non-arctic grass (PFT13) to the CLM5 default to reduce overly productive grasses.

Anticipated changes from #89 are (1) increased LAI in Boreal forests, (2) higher global mean biomass (~650 PgC) and GPP (150 PgC/yr).

Parameter file:
/glade/u/home/linnia/ctsm6_ppe/gen_ensembles/jobscripts/test/gen_paramfiles/ctsm5.3.012.Nfix_params.v13.c250221.nc

I'll let @wwieder comment what simulation is needed and how it should be set up.

Thanks!

---

Case directory:
./create_clone --clone ~slevis/cases_LMWG_dev/ctsm53019_f09_BNF_AD/ --case ~slevis/cases_LMWG_dev/ctsm53026_BNF_AD (cloning seems to have worked @wwieder)

./create_clone --clone ~slevis/cases_LMWG_dev/ctsm53026_BNF_AD/ --case ~slevis/cases_LMWG_dev/ctsm53026_BNF_SASU
./create_clone --clone ~slevis/cases_LMWG_dev/ctsm53026_BNF_SASU/ --case ~slevis/cases_LMWG_dev/ctsm53026_BNF_pSASU
./create_newcase --case ~slevis/cases_LMWG_dev/ctsm53026_BNF_hist --compset HIST_DATM%CRUv7_CLM60%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV --res f09_t232 --run-unsupported

---

Sandbox:
Locally:
/glade/work/slevis/git_ctsm_tags/ctsm5.3.026 (currently the latest ctsm tag)

---

user_nl_ changes:

1. AD as in 1-deg I-case spin-up and historical to test new paramfile for CESM3 (ctsm53019_f09_BNF_[AD, SASU, pSASU, hist]) #89 plus
   paramfile = '/glade/u/home/linnia/ctsm6_ppe/gen_ensembles/jobscripts/test/gen_paramfiles/ctsm5.3.012.Nfix_params.v13.c250221.nc'

2. SASU and pSASU also

reseed_dead_plants = .true.

3. hist

diff ctsm53026_BNF_hist/user_nl_clm ctsm53026_BNF_pSASU/user_nl_clm
22c22,23
< finidat = '/glade/derecho/scratch/slevis/archive/ctsm53026_BNF_pSASU/rest/0161-01-01-00000/ctsm53026_BNF_pSASU.clm2.r.0161-01-01-00000.nc'
---
> 
> hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS','H2OSNO','NFIX'
26a28,29
> reseed_dead_plants = .true.
> 
32,71d34
<
< use_init_interp = .true.
< ndep_taxmode = 'extend'
<
< ! "Survivability" output
< hist_dov2xy(1) = .true.
< hist_dov2xy(2) = .false.
< hist_dov2xy(3) = .false.
< hist_dov2xy(4) = .true.
< hist_dov2xy(5) = .false.
< hist_fields_list_file = .false.
< 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'
< 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'
< hist_fincl4 = 'PCT_GLC_MEC', 'QICE_FORC', 'TSRF_FORC', 'TOPO_FORC', 'PCT_NAT_PFT', 'PCT_LANDUNIT', 'FSNO_ICE',
<         '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'
< 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'

---

SourceMods:
None

---

Diagnostics:
land diagnostics
https://webext.cgd.ucar.edu/I20TR/ctsm53026_BNF_hist/lnd/ctsm53026_BNF_hist.2004_2023-ctsm53019_f09_BNF_hist.2004_2023/setsIndex.html

ILAMB
https://webext.cgd.ucar.edu/I20TR/ctsm53026_BNF_hist/lnd/_build_newscore_CLM50CTSM51BGC

---

Output:
/glade/derecho/scratch/slevis/archive/ctsm53026_BNF_AD
/glade/derecho/scratch/slevis/archive/ctsm53026_BNF_SASU
/glade/derecho/scratch/slevis/archive/ctsm53026_BNF_pSASU
/glade/derecho/scratch/slevis/archive/ctsm53026_BNF_hist
Archived at:
/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM53_DEV/ctsm53026_BNF_pSASU
/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM53_DEV/ctsm53026_BNF_hist

---

Contacts:
@linniahawkins

---

Extra details:
AD startup 400 yrs, 1.80% disequilibrium
SASU hybrid 320 yrs, 11.77% disequilibrium
pSASU hybrid 160 yrs, 2.24% disequilibrium
hist 1850-2023

[slevis-lmwg]

AD DONE
SASU DONE
pSASU DONE
hist DONE
land diags DONE
ILAMB DONE

[slevis-lmwg]

@wwieder I assume that I should add "survivability" output to the first (or second?) half of the hist simulation.

@samsrabin do you also want the extra output that you requested for #89?

[samsrabin]

Is the only difference from #89 the new parameter file? If so, then no, don't worry about it—crops shouldn't be affected.

[wwieder]

Yes, at some point having h1 files with monthly LAI by PFT would be excellent.

doing this for the same time period as #89 makes sense and will make the evaluation cleaner.

[slevis-lmwg]

Yes, at some point having h1 files with monthly LAI by PFT would be excellent.

doing this for the same time period as #89 makes sense and will make the evaluation cleaner.

Right, I remember you saying that survivability output at the beginning of the run was most helpful, because you could look at it early. The 1850-1920 segment of the hist simulation is in progress.

[wwieder]

Survivability in this case (not weighted by pft weights on grid cells)

For comparison, this is where we were with #89

[wwieder]

No diagnostics yet, but it really only looks like PFT2 and C3 grasses have any significant changes (as intended).

[linniahawkins]

PFT2 looks great! C3 grasses are still weird (as expected), but I think they will require a physiologically informed tuning.