Setup transient simulations to use "limit" rather than "cycle" or "extend" for historical stream files...

[ekluzek]

This was a problem for the CMIP5 simulations with CESM2 for the streams files: urbantv, pop-dens, and ndep for the future scenarios. With CESM2 the start date was 2015, but to be compatible with CMIP5 it needed to be 2006. So urbantv and pop-dens were were to extend so used 2015 values, while ndep cycled back to the future year of 2092. cycle and extend make sense for fixed 1850 or 2000 conditions, but for transient we should use "limit" so that it would abort if you got the start or end years wrong. "cycle" is especially bad in these cases.

This would also apply for lai-streams and prescribed soil moisture. Lightning right now is always a climatology so should always be cycle.

[ekluzek]

Only ndep_taxmode is currently in the namelist (the others are hardhired in the code), so it'll only be easy to change it for ndep. But, I think doing that should help to find these types of problems. But, longer term we should add taxmode variables to the namelist.

[ekluzek]

OK, at least with NUOPC this ends up being a problem.

Running SMS_D_Vnuopc_Mmpi-serial.CLM_USRDAT.I1PtClm51Bgc.izumi_nag.clm-NEON_NIWO

I get the following error which says that it can't use limit on the ndep file (valid between 1849-2100) when you tell it to start in 2018 because Jan/1/2018 is before the first valid 2018 in the time series.

(shr_stream_findBounds) ERROR: limit on and rDateIn lt rDatelvd 2.0180101000000000E+07 2.0180116500000000E+07
ERROR: (shr_stream_findBounds) ERROR: rDateIn lt rDatelvd limit true

@mvertens what I'd like NUOPC to do in this case is to use the Dec/2017 data in order to properly interpolate in time to get valid Jan/1st/2018 data. Is that something that can be done?

[ekluzek]

I had tested a case that started in the middle of the series which worked fine. But, now I've tested a MCT case and see the same problem.

SMS_D_Ld1_Mmpi-serial.1x1_smallvilleIA.IHistClm45BgcCropQianRs.cheyenne_intel.clm-cropMonthOutput

There error I see for the MCT case is...

/glade/p/cesmdata/cseg/inputdata/lnd/clm2/ndepdata/fndep_clm_hist_b.e21.BWHIST.
f09_g17.CMIP6-historical-WACCM.ensmean_1849-2015_monthly_0.9x1.25_c180926.nc
(shr_strdata_init_streams) timeName = 1 time
(shr_strdata_init_streams) lonName = 1 lon
(shr_strdata_init_streams) latName = 1 lat
(shr_strdata_init_streams) hgtName = 1 undefined
(shr_strdata_init_streams) maskName = 1 mask
(shr_strdata_init_streams) areaName = 1 area
(shr_strdata_init_mapping) calling shr_dmodel_mapSet for remap
(shr_strdata_print) ----------------------------------------------------------
(shr_strdata_print) name = CLMNDEP data
.
.
.
(shr_strdata_print) streams ( 1) = null
(shr_strdata_print) taxMode ( 1) = limit
(shr_strdata_print) dtlimit ( 1) = 1.500000E+00
.
.
.
(shr_strdata_print) readmode( 1) = single
(shr_strdata_print)
(shr_strdata_print) nvectors = 0
(shr_strdata_print) ----------------------------------------------------------
(shr_stream_findBounds) ERROR: limit on and rDateIn lt rDatelvd
18500101.0000000 18500116.5000000
ERROR: (shr_stream_findBounds) ERROR: rDateIn lt rDatelvd limit true

[ekluzek]

This is blocked by the CDEPS issue.