ERROR: (shr_dmodel_readstrm) ERROR in data sizes

Dear All,

I am currently running CLM5 for my region using GSWP3v1 forcing data. I have successfully created the domain and surface datasets at a 0.5x0.5 resolution and built the model without issues. However, when I submit the case, an error occurs. Could anyone provide guidance on how to resolve this problem?

I would greatly appreciate any help or suggestions you can offer.

Thank you in advance!
Best regards,
Belasoviet
-Campset: I2000Clm50BgcCropGs

./xmlchange STOP_N=1
./xmlchange DIN_LOC_ROOT=/home/phouta/cesm/inputdata
./xmlchange ATM_DOMAIN_FILE=domain.lnd.TRSR_0.5_noocean.250221.nc
./xmlchange RUN_REFDATE=2010-01-01
./xmlchange RUN_STARTDATE=2010-01-01
./xmlchange DATM_CLMNCEP_YR_START=2010
./xmlchange DATM_CLMNCEP_YR_END=2010
./xmlchange ATM_DOMAIN_PATH=/home/phouta/cesm/inputdata/share/domains
./xmlchange LND_DOMAIN_FILE=domain.lnd.TRSR_0.5_noocean.250221.nc
./xmlchange LND_DOMAIN_PATH=/home/phouta/cesm/inputdata/share/domains
echo "fsurdat = '/home/phouta/cesm/inputdata/lnd/clm2/surfdata_map/surfdata_TRSR_0.5_78pfts_CMIP6_simyr2000_c250221.nc'" >> user_nl_clm
echo "mapalgo = 'nn','nn','nn','nn','nn'" >> user_nl_datm
echo "use_init_interp = .true." >> user_nl_clm
echo "hist_fincl1 = 'TLAI', 'QRUNOFF'" >> user_nl_clm
echo "check_dynpft_consistency = .false." >> user_nl_clm
./xmlchange DIN_LOC_ROOT_CLMFORC=/home/phouta/cesm/inputdata/lmwg/atm_forcing.datm7.GSWP3.0.5d.v1.c170516
echo "hist_nhtfrq = -24" >> user_nl_clm

There seem to be some mismatches in your surface dataset and domain file. Surface dataset has lsmlon=19, lsmlat=13. While domain file has 18 and 13.
Also, you should not be using domain.lnd.TRSR_0.5_noocean.250221.nc as the domain file in your datm.streams files. You should be using the default GSWP3 domain file.

Dear Oleson,
thank for your replied!
Base on your suggestion earlier, I have change domain file to default GSWP3 and I masked the file by using ncl_scripts. But I am still facing a problem when I submit case (ERROR: (shr_ncread_domain) ERROR var does not exist frac). Do you have any ideas about this?
Thank you!

I think the error is saying that the variable "frac" does not exist on the domain file you are using:

fatmlndfrc = '/home/phouta/cesm/inputdata/share/domains/domain.lnd.360x720_TRSR.c170606.nc'

Thank you for your response!
I am so confusing about how to do regional running by using the default GSWP3 forcing data. Previously, I created the domain and surface data by using tools in clm5 (mkmapdata, gen_domain_files, mksurfdata_map) but it had a problem with data sizes. After that I used a domain file from the default GSWP3 and it's still error with "frac" does not exist. So my question are :

1. Do I need to create the domain and surface data by using (tools/site_and_regional) instead of clm tools?
2. If I use the default GSWP3 domain data, where can I get surface data from? Since The error above I used surfdata_360x720cru_78pfts_simyr2000_c170428.nc and I masked it together with domain file by using ncl_scripts.
If none of the above reasons apply, could you please provide me with the flow step for running the model in my region using GSWP3 forcing data?

thank you!