TRENDY2024 yr-1700 + historical cases with ctsm5.2 fsurdat (TRENDY2024_f09_ctsm52_AD_20240726 followed by SASU, pSASU, HIST)

[slevis-lmwg]

Description:
"CONTROL" year-1700 spin-up and
"CONTROL" 1700-2023 historical simulations
using ctsm5.2 fsurdat/landuse files for an evaluation of the equivalent ctsm5.3 TRENDY2024 intercomparison simulations:

1700 AD DONE (cold start)
1700 SASU DONE (hybrid)
1700 pSASU DONE (hybrid)
1701-1920-2023 HIST DONE (hybrid)

---

Case directory:
Locally:
/glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_AD_20240726
/glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_SASU_20240801
/glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_pSASU_20240805
/glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_1700-2023_20240813

---

Sandbox:
Locally:
/glade/work/afoster/TRENDY_2024/CTSM_TRENDY
MODEL_VERSION = ctsm5.2.015 in env_case.xml

---

user_nl_ changes:
Diffs in AD from #60:

23c23
< fsurdat='/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_0.9x1.25_hist_1700_78pfts_c240726.nc'
---
> fsurdat='/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_0.9x1.25_hist_1700_78pfts_c240731.nc'

Diffs between SASU and AD similar to #58 but not identical

13c13
<  finidat = '/glade/derecho/scratch/slevis/archive/TRENDY2024_f09_ctsm52_AD_20240726/rest/0321-01-01-00000/TRENDY2024_f09_ctsm52_AD_20240726.clm2.r.0321-01-01-00000.nc'
---
>  finidat = ' '
31d30
<  iloop_avg = -999
38,39d36
<  nyr_forcing = 20
<  nyr_sasu = 20
53,54c50
<  spinup_matrixcn = .true.
<  spinup_state = 0
---
>  spinup_state = 2
77c73
<  use_matrixcn = .true.
---
>  use_matrixcn = .false.
81c77
<  use_soil_matrixcn = .true.
---
>  use_soil_matrixcn = .false.
256c252
<  nnegcrit = -7.d0
---
>  nnegcrit = -6.d+0

Diffs between pSASU and SASU:

<  finidat = '/glade/derecho/scratch/slevis/archive/TRENDY2024_f09_ctsm52_SASU_20240801/rest/0241-01-01-00000/TRENDY2024_f09_ctsm52_SASU_20240801.clm2.r.0241-01-01-00000.nc'
---
>  finidat = '/glade/derecho/scratch/slevis/archive/TRENDY2024_f09_ctsm52_AD_20240726/rest/0321-01-01-00000/TRENDY2024_f09_ctsm52_AD_20240726.clm2.r.0321-01-01-00000.nc'
24a25
>  hist_empty_htapes = .true.
25a27,29
>  hist_fincl1 = 'TOTECOSYSC', 'TOTECOSYSN', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TLAI', 'GPP', 'CPOOL', 'NPP', 'TWS'
>  hist_mfilt = 20
>  hist_nhtfrq = -8760
26a31
>  iloop_avg = -999
32a38,39
>  nyr_forcing = 20
>  nyr_sasu = 20
45a53
>  spinup_matrixcn = .true.
69c77
<  use_matrixcn = .false.
---
>  use_matrixcn = .true.
73c81
<  use_soil_matrixcn = .false.
---
>  use_soil_matrixcn = .true.
248c256
<  nnegcrit = -6.d+0
---
>  nnegcrit = -7.d0

Diffs between Hist and pSASU:

3c3,4
<  co2_type = 'diagnostic'
---
>  co2_ppmv = 277.75
>  co2_type = 'constant'
12c13
<  finidat = '/glade/derecho/scratch/slevis/archive/TRENDY2024_f09_ctsm52_pSASU_20240805/rest/0221-01-01-00000/TRENDY2024_f09_ctsm52_pSASU_20240805.clm2.r.0221-01-01-00000.nc'
---
>  finidat = '/glade/derecho/scratch/slevis/archive/TRENDY2024_f09_ctsm52_SASU_20240801/rest/0241-01-01-00000/TRENDY2024_f09_ctsm52_SASU_20240801.clm2.r.0241-01-01-00000.nc'
24d24
<  hist_dov2xy = .true.,.false.
26,51d25
<  hist_fincl1 = 'AR', 'COL_FIRE_CLOSS', 'CROPPROD1C_LOSS', 'CWD_C', 'DEADCROOTC', 'DENIT', 'DWT_CONV_CFLUX',
<          'F_N2O_DENIT', 'F_N2O_NIT', 'FAREA_BURNED', 'FFIX_TO_SMINN', 'FROOTC', 'FSDS', 'GPP',
<          'HR', 'LEAFC', 'LEAFN', 'LITFALL', 'LITFIRE', 'LIVECROOTC', 'NBP',
<          'NDEP_TO_SMINN', 'NET_NMIN', 'NFIX', 'NPP', 'QRUNOFF', 'QSOIL', 'QVEGE',
<          'QVEGT', 'RAIN', 'RR', 'SMIN_NO3_LEACHED', 'SMIN_NO3_RUNOFF', 'SMINN_TO_PLANT', 'SNOW',
<          'SOILICE', 'SOILLIQ', 'TLAI', 'TOT_WOODPRODC', 'TOT_WOODPRODC_LOSS', 'TOTLITC', 'TOTLITN',
<          'TOTSOILICE', 'TOTSOILLIQ', 'TOTSOMC', 'TOTSOMN', 'TOTVEGC', 'TOTVEGN', 'TSA',
<          'TSOI', 'WOODC', 'H2OSNO', 'FSR', 'FGEV', 'FCTR', 'FSRNI',
<          'FSRVI', 'FSDSNI', 'FSDSVI', 'FSH', 'FLDS', 'FIRE', 'SOMC_FIRE',
<          'SOM_ACT_C', 'SOM_ACT_C_1m', 'SOM_SLO_C', 'SOM_SLO_C_1m', 'SOM_PAS_C', 'SOM_PAS_C_1m', 'LEAFC_TO_LITTER',
<          'CWD_C_TO_LIT_CEL_C', 'CWD_C_TO_LIT_LIG_C', 'FROOTC_TO_LITTER', 'LIT_CEL_C_TO_SOM_ACT_C', 'LIT_LIG_C_TO_SOM_SLO_C', 'LIT_MET_C_TO_SOM_ACT_C', 'SOM_ACT_C_TO_SOM_PAS_C',
<          'SOM_ACT_C_TO_SOM_SLO_C', 'SOM_SLO_C_TO_SOM_ACT_C', 'SOM_SLO_C_TO_SOM_PAS_C', 'SOM_PAS_C_TO_SOM_ACT_C', 'SOMFIRE', 'QIRRIG_DRIP', 'DEADCROOTN',
<          'DEADSTEMN', 'LIVECROOTN', 'LIVESTEMN', 'FROOTN', 'SMINN', "TOT_WOODPRODN", 'M_DEADROOTC_STORAGE_TO_FIRE',
<          'M_DEADROOTC_TO_FIRE', 'M_DEADROOTC_XFER_TO_FIRE', 'M_DEADSTEMC_STORAGE_TO_FIRE', 'M_DEADSTEMC_TO_FIRE', 'M_DEADSTEMC_XFER_TO_FIRE', 'M_FROOTC_STORAGE_TO_FIRE', 'M_FROOTC_TO_FIRE',
<          'M_FROOTC_XFER_TO_FIRE', 'M_LEAFC_STORAGE_TO_FIRE', 'M_LEAFC_TO_FIRE', 'M_LEAFC_XFER_TO_FIRE', 'M_LIVEROOTC_STORAGE_TO_FIRE', 'M_LIVEROOTC_TO_FIRE', 'M_LIVEROOTC_XFER_TO_FIRE',
<          'M_LIVESTEMC_STORAGE_TO_FIRE', 'M_LIVESTEMC_TO_FIRE', 'M_LIVESTEMC_XFER_TO_FIRE', 'M_GRESP_STORAGE_TO_FIRE', 'M_GRESP_XFER_TO_FIRE', 'LEAFC_STORAGE', 'LEAFC_XFER',
<          'FROOTC_XFER', 'FROOTC_STORAGE', "LIVECROOTC_STORAGE", "LIVECROOTC_XFER", 'DEADCROOTC_STORAGE', 'DEADCROOTC_XFER', "LEAFN_STORAGE",
<          'LEAFN_XFER', 'DEADCROOTN_STORAGE', 'DEADCROOTN_XFER', 'DEADSTEMN_STORAGE', 'DEADSTEMN_XFER', 'LIVECROOTN_STORAGE', 'LIVECROOTN_XFER',
<          'LIVESTEMN_STORAGE', 'LIVESTEMN_XFER', 'FROOTN_STORAGE', 'FROOTN_XFER', "FSRND", "FSRVD", "FSDSND",
<          "FSDSVD"
<  hist_fincl2 = 'AGNPP', 'FCEV', 'FCTR', 'FGEV', 'FIRE', 'FLDS', 'FROOTC_ALLOC', 'FSDS', 'FSH', 'FSR', 'GPP',
<          'GRAINC_TO_FOOD', 'HR', 'HTOP', 'LEAFC_ALLOC', 'NBP', 'NPP', 'QIRRIG_DRIP', 'TLAI', 'TOTVEGC', 'TOTLITC', 'TOTSOMC',
<          'TV', 'QVEGE', 'QVEGT', 'QSOIL', 'FSRNI', 'FSRVI', 'FSDSNI', 'FSDSVI', 'WOODC_ALLOC', "FSRND", "FSRVD",
<          "FSDSND", "FSDSVD"
<  hist_mfilt = 1,1
<  hist_nhtfrq = 0,0
53c27
<  irrigate = .true.
---
>  irrigate = .false.
92d65
<  use_init_interp = .true.
108c81
<  ndep_taxmode = 'extend'
---
>  ndep_taxmode = 'cycle'
114c87
<  stream_year_last_ndep = 2099
---
>  stream_year_last_ndep = 1850
117d89
<  model_year_align_popdens = 1850
122c94
<  stream_year_last_popdens = 2016
---
>  stream_year_last_popdens = 1850
125d96
<  model_year_align_urbantv = 1850
129c100
<  stream_year_last_urbantv = 2106
---
>  stream_year_last_urbantv = 1850
178,183d148
<  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_c240801.nc'
194d158
<  init_interp_method = 'general'
340a305
>  reseed_dead_plants = .true.

---

SourceMods:
None.

---

Diagnostics:
https://webext.cgd.ucar.edu/I20TR/TRENDY2024_f09_clm60_S3.2004_2023-TRENDY2024_f09_ctsm52_1700-2023_20240813.2004_2023/setsIndex.html

---

Output:
Output (if still available):
Initially: /glade/derecho/scratch/slevis/archive/TRENDY2024_f09_ctsm52_AD_20240726
later: /glade/campaign/cgd/tss/people/slevis/...
Equivalent paths for other cases.

---

Contacts:
@slevis-lmwg @wwieder @olyson @ekluzek

---

Extra details:

./create_clone --clone /glade/work/afoster/TRENDY_2024/simulations/TRENDY2024_f09_clm60_ADspinup_20240726-144515 --case /glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_AD_20240726 --cime-output-root /glade/derecho/scratch/slevis
./create_clone --clone /glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_AD_20240726 --case /glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_SASU_20240801
./create_clone --clone /glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_SASU_20240801 --case /glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_pSASU_20240805
./create_newcase --case /glade/u/home/slevis/cases_ctsm5.2/TRENDY2024_f09_ctsm52_1700-2023_20240813 --compset IHistClm60BgcCrop --res f09_g17 --run-unsupported

Planned sequence of simulations:
1700AD 320 yrs --> 1700SASU 240 yrs --> 1700pSASU 220 yrs --> hist 1701-2023

[slevis-lmwg]

@olyson
I don't want to jinx it, but the SASU simulation got past the error and is in progress now. The fix:
I upped the tolerance nnegcrit = -7.d0 from -6, which is what @ekluzek had recommended from the beginning, but I resisted...

[slevis-lmwg]

@adrifoster
I have created the 1700-2023 case.

1. I recommend that you clone it with ./create_clone from here /glade/work/afoster/TRENDY_2024/CTSM_TRENDY/cime/scripts
2. Now change the following:

• In user_nl_clm replace three lines
  fsurdat = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/surfdata_0.9x1.25_hist_1700_78pfts_c240731.nc'
  flanduse_timeseries = '/glade/work/slevis/git/mksurfdata_toolchain/tools/mksurfdata_esmf/landuse.timeseries_0.9x1.25_hist_1700-2023_78pfts_c240731.nc'
  finidat = '<your 221-01-01 file from pSASU>'
• Run ./case.build (or ./preview_namelists?) to look at the co2 stream in datm.streams.xml, and then modify user_nl_datm_streams to include TRENDY's co2 stream.
• Any other mods that you think you need...

3. At that point I will compare my case to yours and bring back over any changes I need, except your 5.3-related changes.

[adrifoster]

Thanks! Are you going to use TRENDY's co2 stream?

[slevis-lmwg]

If you do, then yes. Let me know.

[adrifoster]

Yes, I have to do that for TRENDY.