wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_[AD, SASU, pSASU, Hist]

[slevis-lmwg]

Description:
New spinup using latest wetland branch tag (see git describe below).

---

Case directory:
Locally:
~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_AD
~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_SASU
~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_pSASU
~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_Hist

./create_newcase --case ~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_AD --compset I1850Clm60BgcCropCrujra --res f09_t232 --run-unsupported
./create_clone --clone ~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_AD --case ~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_SASU
./create_clone --clone ~/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_SASU --case ~/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_pSASU
./create_newcase --case ~slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_Hist --compset IHistClm60BgcCropCrujra --res f09_t232 --run-unsupported

---

Sandbox:
Locally:
/glade/work/slevis/git/wetlands

git describe: branch_tags/wetlands.n04.ctsm5.4.044

---

user_nl_ changes:
All spin-up cases (AD, SASU, pSASU), not the historical (Hist):
From meeting with @ekluzek

./xmlchange DATM_PRESNDEP=clim_1850_cmip7
./xmlchange CLM_NDEP_FROM_CPL=TRUE

user_nl_datm_streams

presndep.clim_1850_cmip7:datafiles = $DIN_LOC_ROOT/atm/cam/chem/trop_strat_mam5_ts4_aero_spinup/ndep/fndep_clm_B1850C_MT4s.1850monthly.ne30_c251220.nc
presndep.clim_1850_cmip7:meshfile = $DIN_LOC_ROOT/share/meshes/ne30pg3_ESMFmesh_cdf5_c20211018.nc
presndep.clim_1850_cmip7:datavars = dry_deposition_NHx_as_N Faxa_ndep_nhx_dry, wet_deposition_NHx_as_N Faxa_ndep_nhx_wet, dry_deposition_NOy_as_N Faxa_ndep_noy_dry, wet_deposition_NOy_as_N Faxa_ndep_noy_wet

to use the correct ndep (not cmip7, yet). The simulation will use constant co2 and cmip6 aerosols.

SASU
user_nl_clm
finidat = '/glade/derecho/scratch/slevis/archive/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_AD/rest/0401-01-01-00000/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_AD.clm2.r.0401-01-01-00000.nc'

reseed_dead_plants has been a discussion issue in the past, so I checked:
It went from AD .true. to SASU and pSASU .false. as intended by the code in namelist_defaults. In some past cases we have kept it .true. in SASU, and a test showed "tiny" diffs between SASU .true. and SASU .false. (#60).

pSASU
user_nl_clm
finidat = '/glade/derecho/scratch/slevis/archive/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_SASU/rest/0241-01-01-00000/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_SASU.clm2.r.0241-01-01-00000.nc'

Hist
user_nl_clm
finidat = '/glade/derecho/scratch/slevis/archive/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_pSASU/rest/0161-01-01-00000/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_pSASU.clm2.r.0161-01-01-00000.nc'
use_init_interp = .true.

...and after 1920 want ./xmlchange DATM_SKIP_RESTART_READ=TRUE for CONTINUE to work.

---

Diagnostics:
LDF 174 vs. 124 (last 20 yrs pSASU)
ILAMB 174 vs. 124 (last 20 yrs pSASU)

LDF 174 vs. 124 (timeseries 1850-2023 and climo last 20 yrs IHist)

---

Output:
/glade/derecho/scratch/slevis/archive/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_AD
/glade/derecho/scratch/slevis/archive/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_SASU
/glade/derecho/scratch/slevis/archive/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_pSASU
/glade/derecho/scratch/slevis/archive/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_Hist

---

Contacts:
@slevis-lmwg

---

Extra details:
AD 400 yrs: 2.90% of the land sfc not in eqm
SASU 240 yrs based on my past records, which originated from Daniel's guidance: 21.25% of land not in eqm
pSASU 160 yrs: 2.77% of the land sfc not in eqm

[slevis-lmwg]

@wwieder based on most of your wetlands cases, I went with noCrop, but should I do Crop?

[wwieder]

Yea, let's have this be a full BGC crop run

[slevis-lmwg]

@ekluzek the case
/glade/u/home/slevis/cases_dev/wetlands.n04.ctsm5.4.044_f09BgcCrop_flat_174_AD
failed to submit due to not finding the ndep file. I will look closer on Monday...

UPDATE We solved by adding the correct ndep filename in user_nl_datm_streams. I added this info above.

[slevis-lmwg]

Spin-up complete.

[wwieder]

Seems like we can compare this case to #131. I'll look and see if I have any pSASU results left and run LDF. We could also look at more recent FLAT cases from wetland branch and WIEMIP runs, but these have a different DATM file and NDEP.

[wwieder]

Looks like I didn't archive results from #131 😒
@slevis-lmwg can you run diagnostics on #124 instead? results are here /glade/campaign/cgd/tss/common/Land_Only_Simulations/CTSM54_DEV/ctsm5.4.CMIP7_ciso_ctsm5.3.075_f09_124_pSASU

[slevis-lmwg]

Posted diagnostics (above in the diags section) for the last 20 yrs of 174 vs. 124.

[wwieder]

Thanks Sam.
Do you think you can try running ILAMB on these two cases too? If you get it to work we can try adding more cases to the ilamb analysis later

[slevis-lmwg]

Thanks Sam. Do you think you can try running ILAMB on these two cases too? If you get it to work we can try adding more cases to the ilamb analysis later

I have ILAMB running. It seems to have completed the "initialized" phase successfully. Next it's "Running model-confrontation pairs..." and reporting "VarNotInModel" a lot. I hope that will not be a showstopper...

Also I hope that it completes in the next 19' because my interactive session ends in 20': Didn't. Starting over:

• Waiting for advice from Keith.
• Also I had --filter .clm2.h0 and I wonder, should I change this to .clm2.h0a? Seems likely to me.

[wwieder]

you may have to do something with time in the txt file that ILAMB reads in to get this to work as expected, since ILAMB is assuming a historical run. It looks like there are column in the model_setup.txt file that allow you to define Shift From and Shift To.

@olyson can you advise here on how to run pSASU cases through ILAMB?