Setting up WIEMIP Transient runs

[wwieder]

How do we initialize transient hillslope hydrology runs from 1850 initial conditions?

I'm trying to figure out how to set up our transient experiments for WIEMIP and need advice.

My plan was to use this compset HIST_DATM%CRUJRA2024b_CLM60%BGC_SICE_SOCN_MOSART_SGLC_SWAV with modifications for the following experiments (Note, there is no landuse change in any of these experiments, flanduse_timeseries = '').

• control (CTRL), which is basically a continuation of our spinup, cycling over 20 years of ndep and climate with constant CO2
• Biogeochemically coupled (BGC) with time evolving CO2 and Ndep, but cycling climate
• Radiatiatively coupled (RAD) with time evolving climate, but cycling NDEP and constant CO2
• Coupled (COU) with time evolving CO2, climate & Ndep

I have a sandbox of usermods where this is working for non-hillslope cases here /glade/u/home/wwieder/wetlands/cime_config/usermods_dirs/clm/WIEMIP

Smoke tests with with a flat case seem to work (e.g. /glade/work/wwieder/wiemip/wetlands.n03.ctsm5.4.029_f09noCrop_flat_158_testBGC), but require setting use_init_interp = .true..

Similar tests with a hillslope hydrology case fail because we can't use use_init_interp = .true. with hillslope hydrology, but it also seems like it's required for starting a historical run (even though flanduse_timeseries = '')

 Reading restart dataset:
 /glade/derecho/scratch/wwieder/archive/wetlands.n03.ctsm5.4.029_f09noCrop_HH_16
 5_AD/rest/0181-01-01-00000/wetlands.n03.ctsm5.4.029_f09noCrop_HH_165_AD.clm2.r.
 0181-01-01-00000.nc
 check_dim_size ERROR: mismatch of input dimension       176267
  with expected value       178441  for variable column
Did you mean to set use_init_interp = .true. in user_nl_clm?
(Setting use_init_interp = .true. is needed when doing a
transient run using an initial conditions file from a non-transient run,
or a non-transient run using an initial conditions file from a transient run,
or when running a resolution or configuration that differs from the initial conditions.)

Suggestions are welcome here @slevis-lmwg, @olyson, @ekluzek, @swensosc (thanks in advance).

[swensosc]

we could run for a month with use_init_interp on, get a restart file, then write a script to overwrite the values on the new restart file using the correct / spun-up restart file.

…

On Mon, Apr 27, 2026 at 9:17 AM will wieder ***@***.***> wrote: *wwieder* created an issue (NCAR/LMWG_dev#165) <#165> How do we initialize *transient* hillslope hydrology runs from 1850 initial conditions? I'm trying to figure out how to set up our transient experiments for WIEMIP and need advice. My plan was to use this compset HIST_DATM%CRUJRA2024b_CLM60%BGC_SICE_SOCN_MOSART_SGLC_SWAV with modifications for the following experiments (Note, there is *no landuse change in any of these experiments*, flanduse_timeseries = ''). - control (CON), which is basically a continuation of our spinup, cycling over 20 years of ndep and climate with constant CO2 - Biogeochemically coupled (BGC) with time evolving CO2 and Ndep, but cycling climate - Radiatiatively coupled (RAD) with time evolving climate, but cycling NDEP and constant CO2 - Coupled (COU) with time evolving CO2, climate & Ndep I have a sandbox of usermods where this is working for non-hillslope cases here /glade/u/home/wwieder/wetlands/cime_config/usermods_dirs/clm/WIEMIP Smoke tests with with a *flat* case seem to work (e.g. /glade/work/wwieder/wiemip/wetlands.n03.ctsm5.4.029_f09noCrop_flat_158_testBGC), but require setting use_init_interp = .true.. Similar tests with a hillslope hydrology case fail because we can't use use_init_interp = .true. with hillslope hydrology, but it also seems like it's required for starting a historical run (even though flanduse_timeseries = '') Reading restart dataset: /glade/derecho/scratch/wwieder/archive/wetlands.n03.ctsm5.4.029_f09noCrop_HH_16 5_AD/rest/0181-01-01-00000/wetlands.n03.ctsm5.4.029_f09noCrop_HH_165_AD.clm2.r. 0181-01-01-00000.nc check_dim_size ERROR: mismatch of input dimension 176267 with expected value 178441 for variable column Did you mean to set use_init_interp = .true. in user_nl_clm? (Setting use_init_interp = .true. is needed when doing a transient run using an initial conditions file from a non-transient run, or a non-transient run using an initial conditions file from a transient run, or when running a resolution or configuration that differs from the initial conditions.) Suggestions are welcome here @slevis-lmwg <https://github.com/slevis-lmwg>, @olyson <https://github.com/olyson>, @ekluzek <https://github.com/ekluzek>, @swensosc <https://github.com/swensosc> (thanks in advance). — Reply to this email directly, view it on GitHub <#165>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AGRN57EM2WQNWDSPUJBVR3T4X52XZAVCNFSM6AAAAACYICBPNSVHI2DSMVQWIX3LMV43ASLTON2WKOZUGMZTMOBSHEZDQMY> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you were mentioned.Message ID: ***@***.***>

[wwieder]

Here's a list of things I'm trying:

1. Out of the box (see lnd.log.5984915.desched1.260427-084731)

ERROR: mismatch of input dimension       176267
 with expected value       178441  for variable column

2. Turn off transient urban stream_fldfilename_urbantv = ''
   • doesn't seem to work, is there another way to do this? See lnd.log.5985052.desched1.260427-085614
3. reset_dynbal_baselines = .true. (same "mismatch of input dimension" error as in 1)
4. stream_year_last_urbantv=1850
5. irrigate = .false.. This solves the issue!

This is an irrigated crop issue!
Do we need to allow for irrigated columns in the spinup with hillslope hydrology to avoid changing column counts? Alternatively, we could turn off irrigated crops in the HIST case?

[wwieder]

Usermods

Working development of usermods here:

• /glade/u/home/wwieder/wetlands/cime_config/usermods_dirs/clm/WIEMIP/BGC

• /glade/u/home/wwieder/wetlands/cime_config/usermods_dirs/clm/WIEMIP/RAD

• Constant N dep with BGC & RAD case?

• Fixed LULCC (flanduse_timeseries = '')

• Urban? (how do I turn this off? automatic with flanduse_timeseries is null) (stream_year_last_urbantv=1850)?

• Aerosols (hold constant, 1850)

• Ozone (does it even matter, hold constant for consistency)

• PopDens (hold constant)

• Irrigation maybe also turn off (irrigate = .false. for HIST cases)

• megan soil moisture (megan_use_gamma_sm = .false., which is default for CLM6)