Create new land only initial conditions

[wwieder]

As we move to CTSM5.2 datasets and introduce CLM6_0 physics options, it seems like it would also be helpful to provide new initial conditions files. @dlawrenncar noted, this may be most important for coupled model simulations, so we can provide more realistic land states in F and B cases as we move forward with CESM3 development runs.

Given very high productivity biases in CRU-JRA forced runs, I'd suggest we provide 1850 initial conditions forced with GSWP3 using CTSM5.2 surface data. It may be most straight forward to do this after #2348 comes to main?

Looking through the LMWG-Dev issues, I don't think we have a GSWP3 forced run that's similar to this case, #54, do we @olyson? If not, we can create a new issue on LWMG-Dev for this.

We should make sure that isotopes, tillage, and residue removal are all on for these spin ups.

Definition of done:

• [ ] Decide on which LND_TUNING_MODES will have initial conditions (clm5_0_GSWP3, clm6_0_TRENDY, clm5_0_cam6, clm6_0_cam7?)
• [ ] Which resolutions will be run? ne30pg3, f09, f19, f45? What about NEON sites?
• [ ] Run all the configurations needed
• [x] Preliminary finidat file in place for clm5_1/clm6_0 all forcing options: lnd/clm2/initdata_esmf/ctsm5.2/clmi.I1850Clm60BgcCrop-ciso.1361-01-01.ne30pg3_mg17_c240317.nc
• [ ] Put files in place and in XML and into a tag

[ekluzek]

Since, this requires CTSM5.2 it actually should wait for #2372. If we need it sooner it could be added to CESM3_dev. But, I'd like to avoid that if we can. This was also on the CTSM5.2 board as a post-5.2 activity, so I'm linking this there.

Also this relates to #554, so linking it here.

[olyson]

@wwieder and I discussed this and thought that we could just clone the most recent deadveg branch simulation and switch from CRUJRA to GSWP3. In that recent simulation, tillage and residue removal was on. And I have been using CTSM5.2 surface datasets, i.e., for that run I used

/glade/campaign/cesm/cesmdata/cseg/inputdata/lnd/clm2/surfdata_esmf/ctsm5.2.0/surfdata_0.9x1.25_hist_78pfts_CMIP6_1850_c230517.nc

I know that @slevis-lmwg has recently re-generated all of the datasets, so I was planning to use that version which appears to be:

/glade/campaign/cesm/cesmdata/cseg/inputdata/lnd/clm2/surfdata_esmf/ctsm5.2.0/surfdata_0.9x1.25_hist_1850_78pfts_c240216.nc

Is there another reason to wait for #2372 ?

[ekluzek]

Nope that makes perfect sense @olyson. Better to not have to wait when you don't have to.

[olyson]

See https://github.com/NCAR/LMWG_dev/issues/57

[olyson]

I initially set this up as an f09 (1deg) simulation, but I think it should actually be an ne30 since we'll use it as initial conditions for the ne30 F- and B-cases coming up. I don't think the finidat interpolation will handle going from 1deg to ne30. Is that a correct assumption @ekluzek and @billsacks ?

[ekluzek]

@olyson it should be fine to interpolate from f09 to ne30np4.pg3 grid. The interpolation is just nearest neighbor so it will work between any type of grid.

But, I agree that we should spinup the land with the workhorse ne30 grid rather than f09. We'll interpolate to f09 from it for the land only simulations.

[billsacks]

I agree with @ekluzek .

[olyson]

Here is an ne30 spunup initial file for 1850 with configuration as noted above:

/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm51_ctsm51d166deadveg_ne30pg3ne30pg3mg17_GSWP3V1_ABsnoCDE_blk_A5BCD_1850pAD.clm2.r.1361-01-01-00000.nc

[wwieder]

Thanks @olyson should this be the default finidat file we point to for simulations with 'modern' tags? Here I'm especially thinking about the CESM alpha (or beta) 17 tag that's upcoming?

[olyson]

Yes, I think it will be suitable for the upcoming CESM tags.

[wwieder]

@ekluzek and @slevis-lmwg, can we make it so this initial conditions file is used out of the box with the next round of coupled model experiments?

/glade/campaign/cgd/tss/people/oleson/CLM5_restarts/ctsm51_ctsm51d166deadveg_ne30pg3ne30pg3mg17_GSWP3V1_ABsnoCDE_blk_A5BCD_1850pAD.clm2.r.1361-01-01-00000.nc

[wwieder]

looking at the project board, I wonder if this should come in with #2492, or actually #2501?

[ekluzek]

@wwieder it looks like that's GSWP3 forcing finidat file. Do you want that with I1850Clm60BgcCrop and also with B1850/F1850 or just the former because it's for GSWP3 forcing?

[wwieder]

all of the above B, I, & F1850 until we get new cpl.hist files to try a spinup for the B case. We want to provide living arctic vegetation as much as possible.

[ekluzek]

It looks like the fsurdat file used for that IC file was with a preliminary version of ctsm5.2.0 datasets, and is different by at least roundoff level. That means that it will likely require use_init_interp set to TRUE.

Also @olyson note that in CSEG Meeting Notes Mike Levy tells us that they will be going to the t232 mask for MOM. So we should start using it for our simulations as well.

[olyson]

Ok, sounds like it is a showstopper, so I'll start another spinup with the latest code and the latest surface dataset (which out of the box seems to be /glade/campaign/cesm/cesmdata/cseg/inputdata/lnd/clm2/surfdata_esmf/ctsm5.2.0/surfdata_ne30np4.pg3_hist_1850_78pfts_c231026.nc) and the t232 mask.

[ekluzek]

@olyson I don't mean to say it's a showstopper. I mean to say we need to transition over to it. I think we need this in place by end of June for the science capability freeze. But, actually we should talk as a group about this in CTSM as well as CESM.

[olyson]

Ah, ok, sorry, I didn't read your comments carefully enough.

[wwieder]

FWIW it also seems like @slevis-lmwg is creating some new initial condition files for matrix-CN work, although these also are unlikely to have the correct ocean mask.

[ekluzek]

New finidat file name will be:

lnd/clm2/initdata_esmf/ctsm5.2/clmi.I1850Clm60BgcCrop-ciso.1361-01-01.ne30pg3_mg17_c240317.nc

And used for standard 1850 control and hist compsets starting at 1850 with clm5_1 and clm6_0 physics.

[ekluzek]

The list of initial condition files is complex. I'm going to map out what we currently have, so we can start to modify the list:

Resolution	Tuning	BGC	Irrigate	Date
f09_g17	clm4_5_GSW	Bgc	T	18500101
f09_g17	clm4_5_CRU	Bgc	F	18500101
f09_g17	clm4_5_CAM6	BgcCrop	F	18500101
f09_g17	clm5_0_GSW	BgcCrop	F	18500101
f09_g17	clm5_0_GSW	Sp	T	18500101
f09_g17	clm5_0_CRU	BgcCrop	F	18500101
f09_g17	clm5_0_CAM6	BgcCrop	F	18500101
f19_g17	clm4_5_GSW	BgcCrop	T	20110101
f09_g17	clm5_0_CAM6	BgcCrop	T	20000101
f19_g17	clm5_0_CAM6	Sp	T	20030101
f09_g17	clm5_0_CAM6	Sp	T	19790101
f19_g17	clm5_0_CAM6	Sp	T	19790101
ne0np4.ARCTIC	clm5_0_CAM6	Sp	T	19790101
ne0np4.ARCTICGRIS	clm5_0_CAM6	Sp	T	19790101
ne0np4.CONUS	clm5_0_CAM6	Sp	T	20130101
ne120np4.pg3	clm5_0_CAM6	Sp	T	20000101
ne30np4.pg3	clm5_1_*	BgcCrop	F	18500101
ne30np4.pg3	clm6_0_*	BgcCrop	F	18500101

[ekluzek]

@adamrher we will likely want your help in getting initial conditions for the VR and ne120 grids for CESM3. We aren't ready to do this now, but wanted to give a heads up to even see if you are able to do that when the time comes. For these grids in CESM3.0 we could use the current finidat files, but that strikes at me as bad to do for the release.

[adamrher]

I can help when the time comes. I'm thinking that I could initialize a ne120 clm60 AMIP run with the clm50 ne120 finidat file we are currently using to initialize the ARCTIC and ARCTICGRIS cases, and run it out for a bit.

[slevis-lmwg]

I am happy to help, as well.

[adamrher]

HI @ekluzek, I'm setting up ne120 to create some land inic for another project, so I can start this work now. Based on your table there are 4 grids that need to be addressed (ARCTIC, ARCTICGRIS, CONUS, and ne120pg3). Three of them use the same finidat file -- the ne120pg3 finidat file I create a few years back. However, CONUS has a special 2013 findat file on the CONUS grid that I don't really know the backstory on. Ill have to do some digging to figure out how or whether to change this CONUS finidat file.

In the meantime, I'm going to generate a new ne120 finidat file to replace the current one used by those three grids. I am using cesm2_3_alpha17f and currently generating landuse_timeseries for the ne120 run. However I just noticed that our workhorse compsets FLTHIST and BLTHIST_v0c have CLM51%BGC-CROP in this latest tag. Should these be CLM52 or CLM60?

[ekluzek]

Thanks for that update @adamrher this sounds great. Using cesm2_3_beta17 will be great for this purpose. The one caution is that there will be tuning changes that come in that might mean we should wait until AFTER the tuning to do finidat files? Or should we do it now as well as later?

I think the CAM testing starts in 2013 for the CONUS grid, so it's probably best to have a finidat file for that date. We might need to brainstorm how to get it though...

Using CLM51%BGC-CROP is perfectly fine (it'll be identical to CLM60%BGC-CROP. So I'd just leave it like that. If it were CLM50 we'd need to change it, but not CLM51.

Also note I haven't filled out the full table, but think I should to get the best handle on the big picture. But, in terms of your contribution it's only the VS grids and ne120, so your part should be covered.

[adamrher]

OK thanks for clarifying my confusion regarding CLM51 and CLM60 (they're the same). I do agree that we may want to hold off on generating new inic until we're closer to a final tuned up version. Let's revisit sometime after the workshop.

To re-do the 2013 inic, I think we can just run the CONUS grid with nudging for a couple months before the 2013 inic date. Someone at ACOM must have done this so I will reach out to them.

[wwieder]

I'll add a comment here, but wanted to note that the initial conditions being created in LMWG_DEV_#65 that has:

• 5.3 surface datasets,
• all the CLM6 physics options turned on excess ice, new fire model, etc, and
• hand tuned parameters to address LAI and TOTVEGC biases in CRUJRA forced runs

Is also being forced by CRUJRAv2, which has no data over Antarctica (consistent with other TRENDY products). I worry that these these initial conditions on the ne30 grid will not be suitable for coupled model simulations. see /glade/campaign/cgd/tss/projects/TRENDY2024/inputs/three_stream

Here's my suggestion:

• 1a) short term fix, combine restart files from CLM5.3_CRUJRA (N of 60S) and CLM5.1_GSWP3 (S of 60S).
• 1b) Alternative short term fix: Do another spinup of CLM5.3 run on the ne30 grid with GSWP3 (so we have the same code base) and cut and paste initial conditions as above.
• 2 ) Longer term bandaid. Blend the forcing datasets so CRUJRA has GSWP3 climate over Antarctica with a climatology over the southern continent for the last decade of the CRUJRA record. Hopefully this isn't too difficult with xarray?

[samsrabin]

Discussed at CTSM SE meeting today.

Needed for coupled model runs that Cecile will start hopefully in the next few weeks. Longest-term fix is to have a version of CRUJRA we can support that includes Antarctica—just put GSWP3 data there.

Decision: Go ahead and do that; skip intermediate steps. @wwieder, @slevis-lmwg, and @adrifoster to work on this. Can be done on Izumi next week during Derecho/Casper downtime if the datasets needed get copied—Sam L. will do this. Even just 1901-1920 to have something.

[slevis-lmwg]

Keith shared his script from 2014 where we ran the model at 0.5 (resolution of CRUNCEP data) forced by Qian data (at T42) and saved 3-hourly coupler history files and then blended that into CRUNCEP forcing data. /glade/u/home/oleson/cruncep/blend_Qian_CRUNCEP.ncl