clm model run error: IHist compset at 50km (the error: mismatch of input dimension)

[MS_CLM]

I've compiled the CLM model at --res f05_g17 (50km) with --compset IHistClm50BgcCropGs. However, model run gives the following error:

#### error #####
Attempting to read crop dynamic landuse data .....
(GETFIL): attempting to find local file
landuse.timeseries_0.47x0.63_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c180508.nc
(GETFIL): using
/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/surfdata_map/landuse.times
eries_0.47x0.63_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c180508.nc
check_dim ERROR: mismatch of input dimension 2 with expected value
64 for variable cft
ERROR:
ERROR in /gpfs/fs1/home/c/cgf/msol/cesm2.1.3_clm/components/clm/src/main/ncdio_pio.F90.in at line 368
#################################

The above error is referred to the ncdio_pio.F90.in, which those lines attached below:

And, full details about the error are attached.

########
I have included the followings in my user_nl_clm:
irrigate = .false.
use_init_interp = .true.
fsurdat = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/surfdata_map/surfdata_0.47x0.63_78pfts_CMIP6_simyr1850_c180508.nc'
flanduse_timeseries = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/surfdata_map/landuse.timeseries_0.47x0.63_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c180508.nc'
######################

Also, I have the followings in my lnd_in
fatmlndfrc = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/share/domains/domain.lnd.fv0.47x0.63_gx1v7.180521.nc'
finidat = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/initdata_map/clmi.I2000Clm50BgcCrop.2011-01-01.1.9x2.5_gx1v7_gl4_simyr2000_c190312.nc'
fsnowaging = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/snicardata/snicar_drdt_bst_fit_60_c070416.nc'
fsnowoptics = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/snicardata/snicar_optics_5bnd_c090915.nc'
fsurdat = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/surfdata_map/surfdata_0.47x0.63_78pfts_CMIP6_simyr1850_c180508.nc'

flanduse_timeseries = '/scratch/c/cgf/msol/cesm2.1.5_clm/inputdata/lnd/clm2/surfdata_map/landuse.timeseries_0.47x0.63_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c180508.nc'

glacier_region_behavior = 'single_at_atm_topo','virtual','virtual','multiple'

#### Note: for finidat, I also did use ''clmi.I1850Clm50BgcCrop.1366-01-01.0.9x1.25_gx1v6_simyr1850_c171213.nc'', but the same error I got.


I have to say that, I have already successfully compiled and run the clm model with the same compset: --compset IHistClm50BgcCropGs at --res f19_g16 (200km) and --res f09_g16 (100km). But, it gives an error at 50km related to mismatch of input dimension.

Which datasets for fsurdat and flanduse_timeseries should be used to run the clm model with ''IHist'' compsets at ''50km''? is it possible at all? So far I have tested it with different IHist-related compsets at 50km resolutions, but all failed.


I would appreciate your suggestions/thoughts/recommendations to resolve this issue.
Thank you.

 

[oleson]

That error is occurring because you've invoked the crop model in your compset (BgcCrop) and you are using a 16pft dataset which does not have specific crop types on it. You either need a 78pft surface dataset or turn off the crop model.

 

[MS_CLM]

That error is occurring because you've invoked the crop model in your compset (BgcCrop) and you are using a 16pft dataset which does not have specific crop types on it. You either need a 78pft surface dataset or turn off the crop model.

Thank you very much for the reply. Using 78pft is not possible for the landuse.timeseries at 50km resolution (0.47x63): https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/lnd/clm2/surfdata_map/ so I'll try the option to turn off the crop model. To confirm, is it enough to 'false' the followings in lnd_in: (create_crop_landunit = .true. ; use_crop = .true. ; do_transient_crops = .true.). Thank you.

 

[oleson]

That might work. Another way is to create a new case and specify BGC in the compset longname, not BGC-CROP.

 

[MS_CLM]

That might work. Another way is to create a new case and specify BGC in the compset longname, not BGC-CROP.

Thanks. I'll try it.

 

[MS_CLM]

As an update: the clm model (50km) is not running if you just turn off the crop model in the namelist using the compset IHistClm50BgcCropGs; instead, I did compile the compset without the crop-model i.e. IHistClm50Bgc, then it worked and running.