domain_checksame ERROR: output domain lon and gridmap lon not the same

[zqzq]

Dear all,
When I execute the command “./mksurfdata.pl -r usrspec -usr_gname 1x1_PAILONG -usr_gdate $CDATE”, it prompts the following questions:

Attempting to make %lake and %wetland .....
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDMASK
Open lake file:
/es01/paratera/sce0489/cesm2.2/my_cesm_sandbox/inputdata/lnd/clm2/rawdata/mksrf
_LakePnDepth_3x3min_simyr2004_csplk_c151015.nc
(gridmap_map_read) reading mapping matrix data...
(gridmap_map_read) * file name : ../mkmapdata/map_3x3min_MODIS-wCsp_to_1x1_PAILONG_nomask_aave_da_c211014.nc
* matrix dimensions rows x cols : 25920000 x 1
* number of non-zero elements: 9
domain_checksame ERROR: output domain lon and gridmap lon not the same at ni =
1
domain lon= 95.3200000000000
gridmap lon= 95.6100000000000
abort:
ERROR in mksurfdata_map: 34304


however，I checked the corresponding domain file, the latitude and longitude is no problem.


Hope to get your help, thank you！

 

[oleson]

Check to see whether "xc_b" in the file pointed to by mksrf_fgrid in your namelist is the same as "xc_b" in this file:

../mkmapdata/map_3x3min_MODIS-wCsp_to_1x1_PAILONG_nomask_aave_da_c211014.nc

 

[zqzq]

Check to see whether "xc_b" in the file pointed to by mksrf_fgrid in your namelist is the same as "xc_b" in this file:

../mkmapdata/map_3x3min_MODIS-wCsp_to_1x1_PAILONG_nomask_aave_da_c211014.nc

Thank you for your reply，Can you tell me which domain file refers to in the model error? What is the mksrf_fgrid you pointed out above? Looking forward to your reply.

 

[oleson]

There is a namelist file created that is used by mksurfdata.pl in the mksurfdata_map directory. It should be something like:

surdata_............namelist

In that file there should be a definition for mksrf_fgrid

 

[zqzq]

There is a namelist file created that is used by mksurfdata.pl in the mksurfdata_map directory. It should be something like:

surdata_............namelist

In that file there should be a definition for mksrf_fgrid

Thank you for your patience，According to your prompt, I have found the corresponding list: mksrf_fgrid ='$griddata', but I still cannot find the corresponding file and its corresponding xc_b accurately. Thank you again for your reply.

 

[oleson]

When you run mksurfdata.pl there should be some output like this telling you what file it's using:

Attempting to create surface boundary data .....
------------------------------------------------------------
mksrf_fgrid =
../mkmapdata/map_0.25x0.25_MODIS_to_1x1_boulderCO_nomask_aave_da_c200522.nc

 

[zqzq]

When you run mksurfdata.pl there should be some output like this telling you what file it's using:

Attempting to create surface boundary data .....
------------------------------------------------------------
mksrf_fgrid =
../mkmapdata/map_0.25x0.25_MODIS_to_1x1_boulderCO_nomask_aave_da_c200522.nc

Indeed, the xc_b in the two files are different. What should I do next? thank you for your reply

When you run mksurfdata.pl there should be some output like this telling you what file it's using:

Attempting to create surface boundary data .....
------------------------------------------------------------
mksrf_fgrid =
../mkmapdata/map_0.25x0.25_MODIS_to_1x1_boulderCO_nomask_aave_da_c200522.nc

 

[oleson]

Presumably the map_*.nc files were created with mkmapdata.pl. I would check the process you used to create those.