ISSP regional run error: FATAL--length of vector y different from row count of sMat.Length of y

[Maggie Xia]

Hi all,

I am new to CESM and I am trying to run composet ISSP126Clm50BgcCrop, res f09_g16_gl4. Global runs seem fine with all the default data, however, problems occured when I am trying to perform regional runs, refering the CLM_USRDAT_NAME method (1.6.3. Running Single Point Configurations — ctsm release-clm5.0 documentation).

Basic steps I did were:
1. Generated regional domain and surfdata files using getregional_datasets.pl
2. create_newcase
3. ./case.setup
4.

./xmlchange DIN_LOC_ROOT=/home/user/cesmlab/inputdata
./xmlchange CLM_USRDAT_NAME=QTP
./xmlchange ATM_DOMAIN_PATH=/home/user/cesmlab/inputdata/share/domains, LND_DOMAIN_PATH=/home/user/cesmlab/inputdata/share/domainsa
./xmlchange ATM_DOMAIN_FILE=domain.lnd.fv0.9x1.25_gx1v6.QTP.090309.nc, LND_DOMAIN_FILE=domain.lnd.fv0.9x1.25_gx1v6.QTP.090309.nc

5. add lines in user_nl_clm:

fsurdat='/home/user/cesmlab/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.9x1.25_hist_78pfts_CMIP6_simyr1850_QTP_c190214.nc'
flanduse_timeseries='landuse.timeseries_0.9x1.25_SSP1-2.6_78pfts_CMIP6_simyr1850-2100_QTP_c190214.nc'
hist_fincl1='GPP','NPP','TLAI'
hist_mfilt=12
hist_nhtfrq=0

6. several maps need to be directed:

./xmlchange LND2ROF_FMAPNAME=/home/user/cesmlab/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.9x1.25_nomask_to_0.5x0.5_nomask_aave_da_c120522.nc
./xmlchange ROF2LND_FMAPNAME=/home/user/cesmlab/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_nomask_to_0.9x1.25_nomask_aave_da_c121019.nc
./xmlchange LND2GLC_FMAPNAME=/home/user/cesmlab/inputdata/cpl/gridmaps/fv0.9x1.25/map_fv0.9x1.25_TO_gland4km_aave.170429.nc
./xmlchange LND2GLC_SMAPNAME=/home/user/cesmlab/inputdata/cpl/gridmaps/fv0.9x1.25/map_fv0.9x1.25_TO_gland4km_blin.170429.nc
./xmlchange GLC2LND_FMAPNAME=/home/user/cesmlab/inputdata/cpl/gridmaps/gland4km/map_gland4km_TO_fv0.9x1.25_aave.170429.nc
./xmlchange GLC2LND_SMAPNAME=/home/user/cesmlab/inputdata/cpl/gridmaps/gland4km/map_gland4km_TO_fv0.9x1.25_aave.170429.nc

7. After I did ./preview_namelists, the 'datm.streams.txt.presaero.SSP1-2.6' and the 'datm.input_data_list' seemed off. The <fileNames> in the streams file and domain1, file1 in datm.input_data_list appeared to be 'aerosoldep_SSP1-2.6_monthly_1849-2101_QTP.nc', seemed like something created automatically based on my CLM_USRDAT_NAME. However, I want this to be 'aerodep_clm_SSP126_b.e21.BSSP126cmip6.f09_g17.CMIP6-SSP1-2.6.001_2014-2101_monthly_0.9x1.25_c190523.nc' and changed the <fileNames> accordingly. However, even though the streams file has been changed, everytime I re-ran ./preview_namelists, the streams file and the datm.input_data_list authomatically updated to the original. Therefore, I renamed the file to 'aerosoldep_SSP1-2.6_monthly_1849-2101_QTP.nc' (I know this must be improper, but just wanna give it a try).

8. ./case build, ./case.submit
Error occured saying 'FATAL--length of vector y different from row count of sMat.Length of y'. Here I attached log and streams.txt.presaero files. Hope anyone know what happened here? Thanks in advance!

P.S., I am running in the container on Linux OS (Linux version 3.10.0-957.10.1.el7.x86_64 and Docker version 1.13.1, build 7d71120/1.13.1).

 

[oleson]

Your error could be related to this:

I1850Clm50Bgc error

Hello, when I run I1850Clm50Bgc with CESM2.0, an error, Writing to file yang.cism.initial_hist.0001-01-01-00000.nc at time 0.0000000000000000 MCT::m_SparseMatrixPlus:: FATAL--length of vector y different from row count of sMat.Length of y = 25 Number of rows in sMat = 13824...

bb.cgd.ucar.edu

You can use the longname version of ISSP126Clm50BgcCrop to specify the stub glacier model, i.e.,

SSP126_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV

You might also need to specify the stub runoff model as well:

SSP126_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_SROF_SGLC_SWAV

Regarding the aerosol deposition stream, you should copy the datm.streams.txt.presaero.SSP1-2.6 file from your CaseDocs directory into your case directory and prepend user_ to it, i.e.,

user_datm.streams.txt.presaero.SSP1-2.6

 

[Maggie Xia]

Your error could be related to this:

I1850Clm50Bgc error

Hello, when I run I1850Clm50Bgc with CESM2.0, an error, Writing to file yang.cism.initial_hist.0001-01-01-00000.nc at time 0.0000000000000000 MCT::m_SparseMatrixPlus:: FATAL--length of vector y different from row count of sMat.Length of y = 25 Number of rows in sMat = 13824...

bb.cgd.ucar.edu

You can use the longname version of ISSP126Clm50BgcCrop to specify the stub glacier model, i.e.,

SSP126_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_MOSART_SGLC_SWAV

You might also need to specify the stub runoff model as well:

SSP126_DATM%GSWP3v1_CLM50%BGC-CROP_SICE_SOCN_SROF_SGLC_SWAV

Regarding the aerosol deposition stream, you should copy the datm.streams.txt.presaero.SSP1-2.6 file from your CaseDocs directory into your case directory and prepend user_ to it, i.e.,

user_datm.streams.txt.presaero.SSP1-2.6

Thanks, Keith, the model runs fine now. However, I encountered another problem. The TLAI, GPP, NPP in my regional monthly outputs are basically all 0s. I used the same datasets (default data) for global runs, and the variables seems reasonable over my region of interests.