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

[linniahawkins]

Description:
Testing a new parameter set for CESM3. Documentation, paramfile details, and evaluation on sparse grid here

Paramfile change log from ctsm60_c240208.nc (i.e., CLM5.2 params)

• Add Nfix params for BNF parameterization (nfix_tmin, nfix_topt, nfix_tmax)

• Fangs new fire parameterization (rswf_max, rswf_min, fsr_pft, others…)

• Undo dead veg tuning of FUN_fracfixers in arctic (FUN_fracfixers=0.25 for all pfts again)

• Reduce LAI in the tropics without inducing water stress (pft 4,6)

• Reduce LAI and GPP in needleleaf forests (pft 1,2,3)

• Increase mid-latitude evapotranspiration (pft 1,2)

• Reduce arctic LAI for C3 grasses and arctic shrubs (pft 11,12)

• Increase C3 grass LAI (survivability)

• Relax nitrogen limitation in C4 grass to improve survivability (KCN parameters)

• Use t232 ocean mask

Thanks!

---

Case directory:

./create_newcase --case ~slevis/cases_LMWG_dev/ctsm53019_f09_BNF_AD --compset I1850Clm60BgcCropCru --res f09_t232 --run-unsupported
./create_clone --case ~slevis/cases_LMWG_dev/ctsm53019_f09_BNF_SASU --clone ~slevis/cases_LMWG_dev/ctsm53019_f09_BNF_AD
./create_clone --case ~slevis/cases_LMWG_dev/ctsm53019_f09_BNF_pSASU --clone ~slevis/cases_LMWG_dev/ctsm53019_f09_BNF_SASU
./create_newcase --case ~slevis/cases_LMWG_dev/ctsm53019_f09_BNF_hist --compset HIST_DATM%CRUv7_CLM60%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV --res f09_t232 --run-unsupported

---

Sandbox:
/glade/work/slevis/git_ctsm_tags/ctsm5.3.019 (currently the latest ctsm tag)
In #88 and #89 (here) I had to update cime6.1.49 to cime6.1.56 to avoid the rpointer.cpl bug at restart.
In #89 I had to include a change to ccs_config/modelgrid_aliases_nuopc.xml for f09_t232 to work.

---

user_nl_ changes:

1. AD

nfix_method = 'Bytnerowicz'
paramfile = /glade/u/home/linnia/ctsm6_ppe/gen_ensembles/jobscripts/test/gen_paramfiles/ctsm5.3.012.Nfix_params.v9.c250112.nc' (feel free to rename this something 'sensible)

hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS','H2OSNO','NFIX'

do_grossunrep=.false.

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

During the spin-up, we also set MOSART_MODE = NULL in env_build.xml

2. SASU and pSASU same as AD plus

reseed_dead_plants = .true.

3. hist see CTSM 5.3 Spinup and HIST on f09 grid (2nd Attempt - with li2024crujra fire ON this time and using filled Antarctica datm data)  #69 and Sam R's post below and include pft survivability output

user_nl_clm ../ctsm53019_f09_BNF_pSASU/.
19a20,21
> reseed_dead_plants = .true.
> 
22,23d23
< use_init_interp = .true.
< finidat = '/glade/derecho/scratch/slevis/archive/ctsm53019_f09_BNF_pSASU/rest/0161-01-01-00000/ctsm53019_f09_BNF_pSASU.clm2.r.0161-01-01-00000.nc'
25c25
< ndep_taxmode = 'extend'
---
> hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS','H2OSNO','NFIX'
35,78d34
< 
< ! "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'
< 
< ! Requested by Sam Rabin
< hist_dov2xy(6) = .false.
< hist_mfilt(6) = 1
< hist_nhtfrq(6) = -8760
< hist_type1d_pertape(6) = ' '
< hist_fincl6 = 'SDATES:I', 'SDATES_PERHARV:I', 'SYEARS_PERHARV:I', 'HDATES:I', 'GRAINC_TO_FOOD_PERHARV:I', 'GRAINC_TO_FOOD_ANN:I', 'GRAINN_TO_FOOD_PERHARV:I', 'GRAINN_TO_FOOD_ANN:I', 'GRAINC_TO_SEED_PERHARV:I', 'GRAINC_TO_SEED_ANN:I', 'GRAINN_TO_SEED_PERHARV:I', 'GRAINN_TO_SEED_ANN:I', 'HDATES:I', 'GDDHARV_PERHARV:I', 'GDDACCUM_PERHARV:I', 'HUI_PERHARV:I', 'SOWING_REASON_PERHARV:I', 'HARVEST_REASON_PERHARV:I', 'SWINDOW_STARTS:I', 'SWINDOW_ENDS:I', 'GDD20_BASELINE:I', 'GDD20_SEASON_START:I', 'GDD20_SEASON_END:I'

---

SourceMods:
None.

---

Diagnostics:
#89 versus #69
Standard land diagnostics
https://webext.cgd.ucar.edu/I20TR/ctsm53019_f09_BNF_hist/lnd/ctsm53019_f09_BNF_hist.2004_2023-ctsm53n04ctsm52028_f09_hist_2004_2023/setsIndex.html
ILAMB
https://webext.cgd.ucar.edu/I20TR/ctsm53019_f09_BNF_hist/lnd/_build_newscore_CLM50CTSM51BGC

---

Output:
Shown for AD and similar for all the cases.
Initially:
/glade/derecho/scratch/slevis/ctsm53019_f09_BNF_AD/run
Later:
/glade/derecho/scratch/slevis/archive/ctsm53019_f09_BNF_AD
Later:
/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM53_DEV/ctsm53019_f09_BNF_pSASU
/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM53_DEV/ctsm53019_f09_BNF_hist

---

Contacts:
@slevis-lmwg

---

Extra details:
AD startup 360 yrs, 2.18% disequilibrium
SASU hybrid 320 yrs, 11.78% disequilibrium
pSASU hybrid 160 yrs, 2.83% disequilibrium
hist 1850-2023

For reference, the spin-up steps AD + SASU + pSASU took about a week of real time to complete.

[olyson]

Regarding the ocean mask, I think the choice here for --res would be f09_t232.
I'm on PTO beginning in about a half-hour :), so I can either get this going next Tues or is @slevis-lmwg available/willing?

[wwieder]

For the ocean mask to work on the f09_t232 grid this has to be manually added to ccs_config
https://github.com/ESMCI/ccs_config_cesm/pull/205/files

@slevis-lmwg do you think you can start an AD/SASU/pSASU spinup before next week, or should I do this one?

• Namelist changes should basically be the same as CTSM 5.3 Spinup and HIST on f09 grid (2nd Attempt - with li2024crujra fire ON this time and using filled Antarctica datm data)  #69 for the setup + changes Linnia noted above (parameter file + nfix method)
• Create case with t232 grid

[slevis-lmwg]

@slevis-lmwg do you think you can start an AD/SASU/pSASU spinup before next week, or should I do this one? We can talk tomorrow at our 10:00?

Yes, absolutely!

[samsrabin]

@slevis-lmwg For my upcoming crop calibration work, could you save the following instantaneous variables at annual frequency during the historical period?

'SDATES', 'SDATES_PERHARV', 'SYEARS_PERHARV', 'HDATES', 'GRAINC_TO_FOOD_PERHARV', 'GRAINC_TO_FOOD_ANN', 'GRAINN_TO_FOOD_PERHARV', 'GRAINN_TO_FOOD_ANN', 'GRAINC_TO_SEED_PERHARV', 'GRAINC_TO_SEED_ANN', 'GRAINN_TO_SEED_PERHARV', 'GRAINN_TO_SEED_ANN', 'HDATES', 'GDDHARV_PERHARV', 'GDDACCUM_PERHARV', 'HUI_PERHARV', 'SOWING_REASON_PERHARV', 'HARVEST_REASON_PERHARV', 'SWINDOW_STARTS', 'SWINDOW_ENDS', 'GDD20_BASELINE', 'GDD20_SEASON_START', 'GDD20_SEASON_END'

It would also be great to have a restart available around the mid-20th century. Thanks!

[slevis-lmwg]

UPDATE
I've been monitoring the spin-up, and I just started the pSASU phase.

LATEST
I just started the hist phase.

[wwieder]

@samsrabin to evaluate the crop model:

1. I'm assuming you also want monthly CFT output to look at yields (vector format)? Is there other output you'd need (e.g. higher frequency crop model data)?
2. For the list you added above for "instantaneous variables at annual frequency during the historical period" do you want these as gridcell averages, or are CFT level output more helpful?

---

to further clarify will this be enough:

1. Standard crop model output
2. PLUS the list you added above added the h3 file?

[samsrabin]

1. I'm assuming you also want monthly CFT output to look at yields (vector format)? Is there other output you'd need (e.g. higher frequency crop model data)?

I don't need any monthly CFT output because of the work I did a while back to make annual yield outputs possible. The relevant variables in my above list are GRAINC_TO_FOOD_PERHARV and GRAINC_TO_FOOD_ANN.

2. For the list you added above for "instantaneous variables at annual frequency during the historical period" do you want these as gridcell averages, or are CFT level output more helpful?

Good point, sorry—these need to be CFT level to be of any use, so they'd need to be on a new history file (h4?) according to the setup you linked. Vector format is preferred, since that's what my scripts are set up to deal with.

I don't need (and will not use!) GRAINC_TO_FOOD, NFERTILIZATION, or CROPPROD1.