ctsm5.3 f19 spinup for PPE

[linniahawkins]

Description:
For the CLM6-PPE, we need a 2-deg f19 full grid cold start spinup.

We will start from these initial conditions and run our own spin-up protocol for each ensemble member.

For reference, this is the restart we used in the CLM5-PPE, but it's been moved and I don't know where it came from:
'/glade/p/cgd/tss/people/oleson/CLM5_restarts/clm51_PPEn02ctsm51d021_2deg_GSWP3V1_leafbiomassesai_PPE3_1850pAD.clm2.r.2041-01-01-00000.nc'

NOTE: CLM6-PPE is actually a ctsm5.3.0 PPE. It's okay if this spinup is done with a 5.2.02* tag, as long as it is close enough.

Timeframe: We would like to be able to start testing with the new restart by ~Oct.1 so we can kick off a full PPE by Oct. 15

---

Case directory:
Locally:
/glade/work/oleson/ctsm5.3.0/cime/scripts/ctsm530_f09_PPE_AD
/glade/work/oleson/ctsm5.3.0/cime/scripts/ctsm530_f09_PPE_SASU
/glade/work/oleson/ctsm5.3.0/cime/scripts/ctsm530_f09_PPE_pSASU
/glade/work/oleson/ctsm5.3.0/cime/scripts/ctsm530_f09_PPE_hist

---

Sandbox:
Locally:
/glade/work/oleson/ctsm5.3.0

---

user_nl_ changes:

user_nl_clm:
AD

! The 78pft version comes out of the box with this setup, so point explicitly to a 16pft version
fsurdat = '/glade/derecho/scratch/slevis/temp_work/new_rawdata/tools/mksurfdata_esmf/surfdata_1.9x2.5_hist_1850_16pfts_c240926.nc'

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

! turn off shifting cultivation
do_grossunrep = .false.

SASU:

finidat = '/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm530_f19_PPE_AD.clm2.r.0281-01-01-00000.nc'
reseed_dead_plants = .true.

! The 78pft version comes out of the box with this setup, so point explicitly to a 16pft version
fsurdat = '/glade/derecho/scratch/slevis/temp_work/new_rawdata/tools/mksurfdata_esmf/surfdata_1.9x2.5_hist_1850_16pfts_c240926.nc'

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

! turn off shifting cultivation
do_grossunrep = .false.

pSASU:

finidat = '/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm530_f19_PPE_SASU.clm2.r.0281-01-01-00000.nc'
reseed_dead_plants = .true.

! The 78pft version comes out of the box with this setup, so point explicitly to a 16pft version
fsurdat = '/glade/derecho/scratch/slevis/temp_work/new_rawdata/tools/mksurfdata_esmf/surfdata_1.9x2.5_hist_1850_16pfts_c240926.nc'

! turn off shifting cultivation
do_grossunrep = .false.

user_nl_datm_streams:

CLMCRUNCEPv7.Solar:datafiles=/glade/derecho/scratch/wwieder/TRENDY2024/inputs/three_stream/clmforc.CRUJRAv2.5_0.5x0.5.Solr.1901.nc, ...through 2023
CLMCRUNCEPv7.Solar:meshfile=/glade/derecho/scratch/slevis/temp_work/TRENDY2024/inputs/three_stream/mesh_cdf5.nc
CLMCRUNCEPv7.Solar:datavars= FSDS Faxa_swdn

CLMCRUNCEPv7.Precip:datafiles=/glade/derecho/scratch/wwieder/TRENDY2024/inputs/three_stream/clmforc.CRUJRAv2.5_0.5x0.5.Prec.1901.nc, ...through 2023
CLMCRUNCEPv7.Precip:meshfile=/glade/derecho/scratch/slevis/temp_work/TRENDY2024/inputs/three_stream/mesh_cdf5.nc

CLMCRUNCEPv7.TPQW:datafiles=/glade/derecho/scratch/wwieder/TRENDY2024/inputs/three_stream/clmforc.CRUJRAv2.5_0.5x0.5.TPQWL.1901.nc, ...through 2023
CLMCRUNCEPv7.TPQW:meshfile=/glade/derecho/scratch/slevis/temp_work/TRENDY2024/inputs/three_stream/mesh_cdf5.nc
CLMCRUNCEPv7.TPQW:datavars = QBOT Sa_shum,PSRF Sa_pbot,TBOT Sa_tbot,WIND Sa_wind,FLDS Faxa_lwdn

---

SourceMods:

/glade/work/oleson/ctsm5.3.0/cime/scripts/ctsm530_f19_PPE_AD/SourceMods/src.clm/CanopyFluxesMod.F90 (sa_leaf = sa_leaf + esai; see PR#2788)

---

Diagnostics:
Diags (if available)

https://webext.cgd.ucar.edu/I1850/$CASE/lnd/

---

Output:
Output (if still available):
/glade/scratch/oleson/archive/ctsm530_f09_PPE_AD/
/glade/scratch/oleson/archive/ctsm530_f09_PPE_SASU/
/glade/scratch/oleson/archive/ctsm530_f09_PPE_pSASU/
/glade/scratch/oleson/archive/ctsm530_f09_PPE_hist/

---

Contacts:
@linniahawkins
@olyson

---

Extra details:
AD: 280 yrs, 2.08% disequilibrium
SASU: 280 yrs, 13.73% disequilibrium
pSASU: 160 yrs, 2.15% disequilibrium

[wwieder]

Thanks for creating this issue, Linnia. I'd suggest we use the same set up as in #69, but with f19. Maybe it makes sense for us to look at 69 first (at least in 1850) to make sure we don't need any quick fixes or additional adjustments?

Alternatively we can at least do this in the ctsm5.3 code base when it's tagged next week? What do you think @slevis-lmwg and @olyson?

[slevis-lmwg]

Alternatively we can at least do this in the ctsm5.3 code base when it's tagged next week? What do you think @slevis-lmwg and @olyson?

At the moment either seems fine, since #69 is using the latest ctsm5.3 tag (branch_tags/ctsm5.3.n04_ctsm5.2.028), though we should first have resolution regarding Keith's latest posts in issue 2567.

[olyson]

The f09 and ne30 spinups should be in their post-SASU phase by this afternoon, so we could probably start looking at results tomorrow afternoon. As I posted in ESCOMP/CTSM#2567, the differences in CNVegStructMod.F90 code were bfb, and the changes to sa_leaf in CanopyFluxesMod.F90 resulted in small differences, at least in I2000 cases as I posted in ESCOMP/CTSM#2777.

[wwieder]

Is it worth touching base on this at some point on Thursday? 930 or 100 times are better for me

[olyson]

Looks like we have a CTSM5.3 tag now, so it seems like we could start this. Could include the sa_leaf change as a SourceMod. One question I have is that our approach last time was to run without the crop model on since we weren't changing crop parameters. But we weren't tuning either. So I assume we'll run with the crop model on...? I can meet Thurs, but we have CTSM software meeting at 9:30?

[wwieder]

Good idea on the sa_leaf issue, @olyson. I'll let @linniahawkins or @djk2120 weigh in here if they want the crop model on for the 2 degree spinup?

[olyson]

A couple of other questions:

1. No need to have isotopes on right?
2. We have been setting do_grossunrep = .false. (turn off shifting cultivation). Do the same here?

[wwieder]

@slevis-lmwg let's make a 16PFT surface dataset at the f19 resolution in 1850 + a historical simulation.

[slevis-lmwg]

I submitted the job to generate new fsurdat/landuse files earlier this afternoon.
(Issue ESCOMP/CTSM#2791)

[olyson]

Crops off (use 16pft dataset).

[slevis-lmwg]

The files are ready (again, details in ESCOMP/CTSM#2791)

[olyson]

I think this is ready to start (beginning with the AD). But it might be good for anyone who is interested to check the setup:
/glade/work/oleson/ctsm5.3.0/cime/scripts/ctsm530_f19_PPE_AD
@wwieder @slevis-lmwg @linniahawkins @djk2120

[djk2120]

everything matches my expectations

[djk2120]

is it still appropriate to include: reseed_dead_plants = .true. edit: I can see that this is indeed the setting

[slevis-lmwg]

My understanding is: yes.

[slevis-lmwg]

Looks good to me.

[olyson]

reseed_dead_plants is set to .true. by default, but in a cold start (as this AD is) I don't think it actually does anything.

[olyson]

This spinup is complete. The restart file is:

/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm530_f19_PPE_pSASU.clm2.r.0161-01-01-00000.nc

I'm running a branch now to look at survivability. Pending that looking ok, I plan to proceed with the historical.

[olyson]