Dear scientists,
I am currently generating a surface dataset for a CMFD custom grid using CLM5.0, and I encountered a persistent error related to the domain and mapping file masks.
When I run mksurfdata_map, the following error occurs:
Attempting to make %fmax .....
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDMASK
Open soil fmax file:
/data/clm_test/inputdata/lnd/clm2/rawdata/mksrf_fmax_3x
3min_USGS_c120911.nc
(gridmap_map_read) reading mapping matrix data...
(gridmap_map_read) * file name : ../mkmapdata/map_3x3min_USGS_to_CMFD_nomask_aave_da_c260322.nc
* matrix dimensions rows x cols : 4147200 x 280000
* number of non-zero elements: 685461
domain_checksame ERROR: input domain mask and gridmap mask are not the same at
ni = 2421230
domain mask= 0
gridmap mask= 1
abort:
ERROR in mksurfdata_map: 34304
It appears that there is a mismatch between the domain mask in domain.lnd.CMFD.nc and the mapping file map_3x3min_USGS_to_CMFD_nomask_aave_da_c260322.nc. Is there a recommended procedure for aligning the masks when using custom grids to avoid this error?
Any guidance on resolving this issue would be greatly appreciated.
Thank you very much for your help!
Quote Reply
Report
Here is my code：
module load intel/ncl/6.6.2/gcc8.5.0_openmpi5.0.3-wqh7quu
export GRIDNAME=CMFD
export CDATE=`date +%y%m%d`
#export CDATE=260318
export CENLON=105
export CENLAT=35
export DX=70
export DY=40
export NX=700
export NY=400
./mknoocnmap.pl -p $CENLAT,$CENLON -n $GRIDNAME -dy $DY -ny $NY -dx $DX -nx $NX
export GRIDFILE=${HOME}/data/clm_test/clm5new/tools/mkmapgrids/SCRIPgrid_${GRIDNAME}_nomask_c${CDATE}.nc
export MAPFILE=${HOME}/data/clm_test/clm5new/tools/mkmapgrids/map_${GRIDNAME}_noocean_to_${GRIDNAME}_nomask_aave_da_${CDATE}.nc
export ESMFBIN_PATH=/apps/support/intel_spr_rocky8.9/spack/0.22.0/opt/spack/linux-rocky8-icelake/gcc-8.5.0/esmf-8.6.0-bzofgqehyzd6va4cbmujyaexnwvji5kh/bin
export CSMDATA=${HOME}/data/clm_test/inputdata
./mkmapdata.sh -r $GRIDNAME -f $GRIDFILE -t regional
export OCNDOM=noocean
export ATMDOM=$GRIDNAME
cd /tools/mapping/gen_domain_files/
./gen_domain -m $MAPFILE -o $OCNDOM -l $ATMDOM
cd clm5new/tools/mksurfdata_map
#第一次使用，要进入他的src文件里面，make一下，在mksurfdata_map/文件夹下产生mksurfdata_map可执行文件
#cd src
#export INC_NETCDF=/apps/support/intel_spr_rocky8.9/netcdf/4.8.1/intel2024.2_impi2024.2/include
#export LIB_NETCDF=/apps/support/intel_spr_rocky8.9/netcdf/4.8.1/intel2024.2_impi2024.2/lib
#export PATH=$PATH:$LIB_NETCDF
#cd ../
./mksurfdata.pl -r usrspec -usr_gname $GRIDNAME -usr_gdate $CDATE -l $CSMDATA -no-crop
I am currently generating a surface dataset for a CMFD custom grid using CLM5.0, and I encountered a persistent error related to the domain and mapping file masks.
When I run mksurfdata_map, the following error occurs:
Attempting to make %fmax .....
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDMASK
Open soil fmax file:
/data/clm_test/inputdata/lnd/clm2/rawdata/mksrf_fmax_3x
3min_USGS_c120911.nc
(gridmap_map_read) reading mapping matrix data...
(gridmap_map_read) * file name : ../mkmapdata/map_3x3min_USGS_to_CMFD_nomask_aave_da_c260322.nc
* matrix dimensions rows x cols : 4147200 x 280000
* number of non-zero elements: 685461
domain_checksame ERROR: input domain mask and gridmap mask are not the same at
ni = 2421230
domain mask= 0
gridmap mask= 1
abort:
ERROR in mksurfdata_map: 34304
It appears that there is a mismatch between the domain mask in domain.lnd.CMFD.nc and the mapping file map_3x3min_USGS_to_CMFD_nomask_aave_da_c260322.nc. Is there a recommended procedure for aligning the masks when using custom grids to avoid this error?
Any guidance on resolving this issue would be greatly appreciated.
Thank you very much for your help!
Quote Reply
Report
Here is my code：
module load intel/ncl/6.6.2/gcc8.5.0_openmpi5.0.3-wqh7quu
export GRIDNAME=CMFD
export CDATE=`date +%y%m%d`
#export CDATE=260318
export CENLON=105
export CENLAT=35
export DX=70
export DY=40
export NX=700
export NY=400
./mknoocnmap.pl -p $CENLAT,$CENLON -n $GRIDNAME -dy $DY -ny $NY -dx $DX -nx $NX
export GRIDFILE=${HOME}/data/clm_test/clm5new/tools/mkmapgrids/SCRIPgrid_${GRIDNAME}_nomask_c${CDATE}.nc
export MAPFILE=${HOME}/data/clm_test/clm5new/tools/mkmapgrids/map_${GRIDNAME}_noocean_to_${GRIDNAME}_nomask_aave_da_${CDATE}.nc
export ESMFBIN_PATH=/apps/support/intel_spr_rocky8.9/spack/0.22.0/opt/spack/linux-rocky8-icelake/gcc-8.5.0/esmf-8.6.0-bzofgqehyzd6va4cbmujyaexnwvji5kh/bin
export CSMDATA=${HOME}/data/clm_test/inputdata
./mkmapdata.sh -r $GRIDNAME -f $GRIDFILE -t regional
export OCNDOM=noocean
export ATMDOM=$GRIDNAME
cd /tools/mapping/gen_domain_files/
./gen_domain -m $MAPFILE -o $OCNDOM -l $ATMDOM
cd clm5new/tools/mksurfdata_map
#第一次使用，要进入他的src文件里面，make一下，在mksurfdata_map/文件夹下产生mksurfdata_map可执行文件
#cd src
#export INC_NETCDF=/apps/support/intel_spr_rocky8.9/netcdf/4.8.1/intel2024.2_impi2024.2/include
#export LIB_NETCDF=/apps/support/intel_spr_rocky8.9/netcdf/4.8.1/intel2024.2_impi2024.2/lib
#export PATH=$PATH:$LIB_NETCDF
#cd ../
./mksurfdata.pl -r usrspec -usr_gname $GRIDNAME -usr_gdate $CDATE -l $CSMDATA -no-crop