ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_ [AD, SASU, pSASU, HIST]

[wwieder]

Description:
Clone of #123, but on the f09 grid.
Also turning on some isotope history output for AD to visually check results
Mainly doing this because we're thinking this setup will provide default initial conditions for the CTSM5.4 tag

HIST diags from LDF posted here:
https://webext.cgd.ucar.edu/I20TR/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_HIST/lnd/website/

also ILAMB for CLM6_BGC_CRUJRA vs. CLM5_BGC_GSWP3

Not using the tethering script for this... will do this later.

---

Case directory:
Locally:
/glade/work/wwieder/ctsm5.4.0/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_AD

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_AD --compset I1850Clm60BgcCropCrujra --res f09_t232 --run-unsupported --project P08010000

./create_clone --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_AD --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_SASU --keepexe

./create_clone --clone /glade/work/wwieder/ctsm5.4.0/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_AD --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_pSASU --keepexe

./create_newcase --case /glade/work/wwieder/ctsm5.4.0/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_HIST --compset IHISTClm60BgcCropCrujra --res f09_t232 --run-unsupported

---

Sandbox:
Locally:
/glade/u/home/wwieder/ctsm5.4.CMIP7_ciso

---

1. AD

xml changes

./xmlchange MOSART_MODE=NULL
./xmlchange STOP_N=40,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=on,CLM_FORCE_COLDSTART=on
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=8

user_nl_ changes:

hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN',
              'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS', 'H2OSNO',
              'C13_GPP','C14_GPP','C13_TOTVEGC','C14_TOTVEGC'

2. SASU

./xmlchange MOSART_MODE=NULL
./xmlchange RUN_TYPE=hybrid
./xmlchange STOP_N=80,STOP_OPTION=nyears
./xmlchange CLM_ACCELERATED_SPINUP=sasu,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=3
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFCASE=ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_AD
./xmlchange RUN_REFDATE=0361-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0361-01-01-00000
./xmlchange RUN_REFDIR=/glade/derecho/scratch/wwieder/archive/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_AD/rest/0361-01-01-00000

3. pSASU

./xmlchange MOSART_MODE=NULL
./xmlchange STOP_N=80,STOP_OPTION=nyears,REST_N=40
./xmlchange CLM_ACCELERATED_SPINUP=off,CLM_FORCE_COLDSTART=off
./xmlchange CONTINUE_RUN=FALSE,RESUBMIT=1
./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFDATE=0321-01-01
./xmlchange DRV_RESTART_POINTER=rpointer.cpl.0321-01-01-00000
./xmlchange RUN_REFCASE=ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_SASU
./xmlchange RUN_REFDIR=/glade/derecho/scratch/wwieder/archive/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_SASU/rest/0321-01-01-00000

4. HIST
   user_nl_datm_streams

!! manually update datm_streams to get through 2023
co2tseries.20tr:year_last = 2024
co2tseries.20tr:datafiles = /glade/campaign/cesm/cesmdata/inputdata/atm/datm7/CO2/fco2_datm_globalSSP3-7.0_simyr_1750-2501_CMIP6_c201101.nc

presndep.hist:datafiles = /glade/campaign/cesm/cesmdata/inputdata/lnd/clm2/ndepdata/fndep_clm_SSP370_b.e21.BWSSP370cmip6.f09_g17.CMIP6-SSP3-7.0-WACCM.002_1849-2101_monthly_0.9x1.25_c211216.nc
presndep.hist:year_last = 2024

! No CMIP6 aero files? 
! We'll just use default CMIP6 HIST files and cycle after 2014 
!presaero.hist:datafiles = /glade/campaign/cesm/cesmdata/inputdata/atm/cam/chem/trop_mozart_aero/aero/aerosoldep_rcp4.5_monthly_1849-2104_0.9x1.25_c100407.nc
!presaero.hist:year_last = 2024

user_nl_clm

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

use_init_interp = .true.
init_interp_fill_missing_urban_with_HD = .true.
ndep_taxmode = 'extend'

! "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'

1850-1871

./xmlchange RUN_TYPE=startup
./xmlchange RUN_STARTDATE=1850-01-01
./xmlchange STOP_OPTION=nyears
./xmlchange STOP_N=21
./xmlchange CONTINUE_RUN=FALSE
./xmlchange RESUBMIT=0
./xmlchange DATM_YR_ALIGN=1901
./xmlchange DATM_YR_START=1901
./xmlchange DATM_YR_END=1920

1871-1901

./xmlchange STOP_N=30
./xmlchange CONTINUE_RUN=TRUE

1901-2023

./xmlchange STOP_N=31
./xmlchange DATM_YR_END=2023
./xmlchange RESUBMIT=3
./xmlchange DATM_SKIP_RESTART_READ=TRUE

---

SourceMods:

---

Diagnostics:
Diags (if available)

https://webext.cgd.ucar.edu/I20TR/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_HIST/lnd/website/

---

Output:
Output (if still available):
/glade/derecho/scratch/wwieder/archive/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_HIST

/glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM54_DEV/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_*

---

Contacts:
@wwieder

---

Extra details:
Spin-up workflow notes by @slevis-lmwg for use in ESCOMP/CTSM#3998.
Plots located in /glade/work/slevis/git/ans_chging_tags/tools/contrib/.

Last AD file is ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_AD.clm2.h0a.0342-01-01-00000.nc

(0)     EQUILIBRIUM SUMMARY
(0)     =======================================================================
(0)     TOTECOSYSC is in equilibrium. Eq. Yr. = 182
(0)     TOTSOMC is in equilibrium. Eq. Yr. = 122
(0)     TOTVEGC is in equilibrium. Eq. Yr. = 122
(0)     TLAI is in equilibrium. Eq. Yr. = 42
(0)     GPP is in equilibrium. Eq. Yr. = 122
(0)     TWS is in equilibrium. Eq. Yr. = 182
(0)     H2OSNO is in equilibrium. Eq. Yr. = 182
(0)     FATAL: Not enough of the land surface is in equilibrium (  3.26% >   3.00%)

Last SASU file is ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_SASU.clm2.h0a.0302-01-01-00000.nc

(0)     EQUILIBRIUM SUMMARY
(0)     =======================================================================
(0)     FATAL: TOTECOSYSC is NOT in equilibrium
(0)     FATAL: TOTSOMC is NOT in equilibrium
(0)     TOTVEGC is in equilibrium. Eq. Yr. = 162
(0)     TLAI is in equilibrium. Eq. Yr. = 2
(0)     GPP is in equilibrium. Eq. Yr. = 42
(0)     TWS is in equilibrium. Eq. Yr. = 2
(0)     H2OSNO is in equilibrium. Eq. Yr. = 2
(0)     FATAL: Not enough of the land surface is in equilibrium ( 16.12% >   3.00%)

Last pSASU file is ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_pSASU.clm2.h0a.0160-12.nc

(0)     EQUILIBRIUM SUMMARY
(0)     =======================================================================
(0)     TOTECOSYSC is in equilibrium. Eq. Yr. = 1
(0)     TOTSOMC is in equilibrium. Eq. Yr. = 1
(0)     TOTVEGC is in equilibrium. Eq. Yr. = 1
(0)     TLAI is in equilibrium. Eq. Yr. = 1
(0)     GPP is in equilibrium. Eq. Yr. = 1
(0)     TWS is in equilibrium. Eq. Yr. = 1
(0)     H2OSNO is in equilibrium. Eq. Yr. = 1
(0)     At least 97 percent of the land surface is in TOTECOSYSC equilibrium. Eq. Yr. = 61
(0)     Percent of the land surface not in equilibrium (  2.28% )
(0)     Congratulations! Your simulation is in equilibrium

[wwieder]

This historically case is finally done, and I wonder if I can get some help from @lifang0209 and @samsrabin to look at crop and fire results?

raw history data are here
/glade/derecho/scratch/wwieder/archive/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_HIST
with single variable time series and climo results for a few variables
/glade/derecho/scratch/wwieder/ADF_unstruct/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_HIST

Best comparison would be from #118, which uses the same surface dataset and (nearly) the same parameter file with raw data:
/glade/derecho/scratch/wwieder/archive/ctsm5.4_5.3.068_PPEcal115f09_118_HIST/lnd/hist

and post processed data:
/glade/derecho/scratch/wwieder/ADF_unstruct/ctsm5.4_5.3.068_PPEcal115f09_118_HIST

I'm running LDF and will post results above when they are finished.

[wwieder]

@slevis-lmwg once we get approval on the fire and crop results these restart files are now ready for the 5.4 tag.

[samsrabin]

Global yield looks fine to me.

Failed crop area is lower, which is good:

[lifang0209]

@wwieder I’ll check the fire simulations and get back to you tomorrow. I’ve been out of the office

[lifang0209]

Fire simulations look generally good. The global total is reasonable, and the global spatial pattern for run 124 is better than run 118 and PPE101–105.
Issues: in the Arctic–Boreal Zone (ABZ): overestimation west of 60°E, and underestimation over East Siberia and North America. The interannual variability in 124 is similar to 118, and both are still skillful but worse than the earlier version.
These issues should be not hard to address. Overall, I prefer 124. I’ll review in more detail and work on solutions tomorrow.

Burned area fraction (the first three are observations; global total and spatial correlation between model and observations are given in parentheses):

Fire Carbon emissions (the first two are satellite-based products)