post-5.4.n02.ctsm5.4.000_HIST_test0

[linniahawkins]

Description:
Evaluate post5.4 changes #3626, which inlcudes:

• paramfile bug fixes from #3585
• fire calibration #3608
• crop calendar and parameter updates [#3618](updates post-5.4.n00: New crop calendar inputs and updated params (omni02) ESCOMP/CTSM#3618)

This case is a 34 year run starting from 1990 initial conditions from #123 designed for a sanity check on the results from the post5.4 branch. I will do a longer spinup after 5.4.004 tag comes in with additional crop model changes

Note: finidat came from a case that used the wrong datm DATM_YR_END, see #123

---

• new paramfile: /glade/work/linnia/CLM6-PPE/ctsm6_cal/paramfiles/clm60_params.ctsm6_li2024_xlbugfix.c251121.nc
  • This file is now pointing to lnd/clm2/paramdata/ctsm60_params.c251124.nc

paramfile changes:
Starting from clm60_params.ctsm6_li2024.c250822.nc
(1) Set xl parameters back to CLM5 default settings for PFTs 2,8,9,10,11.
Keep the calibrated xl value for PFT4.
(2) Raise bgc_rf_s2s1 from 0.38 to 0.625 (as per calibration).

Expected changes:
(1) Slight increase in albedo for PFTs: 2,8,9,10,11 (<0.01)
(2) Reduction in soil carbon stocks from ~1650PgC to ~1450 PgC

I will leave the choice of model configuration and simulation protocol to you @wwieder

Thanks to @samsrabin 's param_utils for making this easy!

---

Case directory:
Locally:
/glade/work/wwieder/ctsm5.4.0/post-5.4.n02.ctsm5.4.000_HIST_test0

---

Sandbox:
Locally:
/glade/u/home/wwieder/CTSM (post-5.4.n02.ctsm5.4.000 branch)

---

user_nl_ changes:

finidat='/glade/derecho/scratch/wwieder/archive/ctsm5.4.CMIP7_ciso_ctsm5.3.075_ne30_123_HIST/rest/1990-01-01-00000/ctsm5.4.CMIP7_ciso_ctsm5.3.075_ne30_123_HIST.clm2.r.1990-01-01-00000.nc'

ndep_taxmode = 'extend'
!Not sure if I need these?
!use_init_interp = .true.
!init_interp_fill_missing_urban_with_HD = .true.

! turning off isotopes for now...
use_c13 = .false.
use_c14 = .false.

! "Survivability" output
hist_dov2xy(1) = .true.
hist_dov2xy(2) = .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_mfilt(1) = 1
hist_mfilt(2) = 1
hist_nhtfrq(1) = 0
hist_nhtfrq(2) = 0
hist_type1d_pertape(1) = ' '
hist_type1d_pertape(2) = ' '

! Crop output
hist_dov2xy(3) = .false.
hist_mfilt(3) = 1
hist_nhtfrq(3) = -8760
hist_type1d_pertape(3) = ' '
hist_fincl3 = '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:

---

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]

@linniahawkins can you create your requested param file on top of this file from Sam?
...lnd/clm2/paramdata/clm60_params.crop_omni02.c251120.nc

It's OK for you to make it in your /glade/work/linnia/CLM6-PPE/ctsm6_cal/paramfiles/ directory.

[linniahawkins]

Updated paramfile here:
/glade/work/linnia/CLM6-PPE/ctsm6_cal/paramfiles/clm60_params.crop_omni02_xlbugfix.c251124.nc

[wwieder]

@ekluzek and @slevis-lmwg do you want to rimport this param file and I'll also update this in the namelist PR (see linked issue)? If so, do you want to rename the file too?

[slevis-lmwg]

I will ./rimport and relink.py in the next few minutes.

[slevis-lmwg]

Oh, sorry, first rename to ctsm60_params.c251124.nc?

[wwieder]

That name works for me, @slevis-lmwg

[slevis-lmwg]

Done in /glade/campaign/cesm/cesmdata/inputdata:

./rimport -file lnd/clm2/paramdata/ctsm60_params.c251124.nc
./relink.py

[wwieder]

Biggest changes I see are

• LAI declines in crop regions (especially wheat), which seems expected based on conversastion with @samsrabin
• LAI increases in higher latitudes, likely reflecting a temporary N fertilization from faster SOM turnover.
• Increases in burned area, mostly in boreal forests (not shown)
• slight increases in high latitude albedo (not shown)
• declines in global grain harvest (mostly from wheat?)

Generally, these results seem as expected based on the changes in the post 5.4 branch