Create surface dataset with mksurfdata_map (CLM4.5)

[yfan1@uni-goettingen_de]

Dear all, I followed the whole process to create user surface dataset for a point simulation named "1x1pt_Bariri" on a user machine.There was no problem in the first 3 steps: 1. Create SCRIP grid datasets with mknoocnmap.pl , 2. Create domain dataset with gen_domain, 3. Create mapping files for mksurfdata_map with mkmapdata.shThe mapping files were successfully created and raw data files "mksrf_..." were downloaded but mksurfdata_map gave an error message when making %fmax with the mapping file "map_3x3min_USGS_to_1x1pt_Bariri_nomask_aave_da_c131213.nc" "* matrix dimensions rows x cols : 4147200 x 1
* number of non-zero elements: 4
domain_checksame ERROR: output domain size and gridmap destination size are not the same size"(see below near the end)This is strange because the other mapping file "map_0.5x0.5_MODIS_to_1x1pt_Bariri_nomask_aave_da_c131213.nc" created withe same tool mkmapdata.sh worked out OK:"* matrix dimensions rows x cols : 259200 x 259200
* number of non-zero elements: 97125
domain_clean: cleaning 259200
Successfully made soil color classes"I really don't understand what was the problem. I recreated all the mapping files again and it still does not help.In the beginging, I was suspecting that the default SCRIPgrid file "SCRIPgrid_3x3_USGS_c120912.nc" for making the map "map_3x3min_USGS_to_1x1pt_Bariri_nomask_aave_da_c131213.nc" was not in a "3x3min" resolution but should be "3x3". Then I checked the default namelist (/lnd/clm/bld/namelist_files/namelist_defaults_clm4_5_tools.xml) , and it indeed considers this USGS SCRIPgrid file as 3x3min. So I guess this should not be a mistake. Right?  Has anyone got the same problem?Really need help. Thanks in advance.Fan ------------------------------------------------------------./mksurfdata.pl -res usrspec -usr_gname 1x1pt_Bariri -usr_gdate 131213 -y 2000 -usr_mapdir $CSMDATA/lnd/clm2/mappingdata/maps/1x1pt_Bariri -pft_frc 95,5 -pft_idx 4,0 -soil_cly 31 -soil_snd 30 -soil_col 18 -dinlc $CSMDATA
CSMDATA is /usr/users/yfan1/cesm_input
resolution: 1x1pt_Bariri rcp=-999.9 sim_year = 2000
namelist: namelist
&clmexp
nglcec = 0
mksrf_fgrid = '/usr/users/yfan1/cesm_input/lnd/clm2/mappingdata/maps/1x1pt_Bariri/map_0.5x0.5_MODIS_to_1x1pt_Bariri_nomask_aave_da_c131213.nc'......./mksurfdata_map < namelist
Attempting to initialize control settings .....
Attempting to create surface boundary data .....
------------------------------------------------------------mksrf_fgrid = /usr/users/yfan1/cesm_input/lnd/clm2/mappingdata/maps/1x1pt_Bariri/map_0.5x0.5_MODIS_to_1x1pt_Bariri_nomask_aave_da_c131213.nc 
mksrf_gridtype = global
Output ALL data in file as 64-bit
Set wetland to 0% over land
Replace soil clay % for all points with: 31.0000000000000
Replace soil sand % for all points with: 30.0000000000000
Set PFT fraction to : 95.0000000000000 5.00000000000000
With PFT index : 4 0
Output ALL data in file as 100% vegetated
calling domain_read
finished domain_read
fsurdat is 2d lat/lon grid
nlon= 720 nlat= 360
(OPNFIL): Successfully opened file surfdata_1x1pt_Bariri_simyr2000_c131214.log
on unit= 99
mksrf_gridtype = global

Attempting to make PFTs .....
Creating surface datasets with the standard # of PFTs = 16
domain_clean WARN: clean domain unecessary
Successfully made PFTs

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
domain_clean: cleaning 25920000
Successfully made %lake

Attempting to make %wetland .....
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDMASK
domain_clean: cleaning 259200
Successfully made %wetland

Attempting to make %glacier .....
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDMASK
Open glacier file:
/usr/users/yfan1/cesm_input/lnd/clm2/rawdata/mksrf_glacier.060929.nc
domain_clean: cleaning 259200
Successfully made %glacier

Attempting to make %sand and %clay .....
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDMASK
Open soil texture file:
/usr/users/yfan1/cesm_input/lnd/clm2/rawdata/mksrf_soitex.10level.c010119.nc
domain_clean: cleaning 9331200
Successfully made %sand and %clayAttempting to make soil color classes .....
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDMASK
Open soil color file:
/usr/users/yfan1/cesm_input/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-20
05.c090630/mksrf_soilcol_global_c090324.nc
nsoicol = 20
(gridmap_map_read) reading mapping matrix data...
(gridmap_map_read) * file name : /usr/users/yfan1/cesm_input/lnd/clm2/mappingdata/maps/1x1pt_Bariri/map_0.5x0.5_MODIS_to_1x1pt_Bariri_nomask_aave_da_c131213.nc
* matrix dimensions rows x cols : 259200 x 259200
* number of non-zero elements: 97125
domain_clean: cleaning 259200
Successfully made soil color classes

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:
/usr/users/yfan1/cesm_input/lnd/clm2/rawdata/mksrf_fmax_3x3min_USGS_c120911.nc
(gridmap_map_read) reading mapping matrix data...
(gridmap_map_read) * file name : /usr/users/yfan1/cesm_input/lnd/clm2/mappingdata/maps/1x1pt_Bariri/map_3x3min_USGS_to_1x1pt_Bariri_nomask_aave_da_c131213.nc
* matrix dimensions rows x cols : 4147200 x 1
* number of non-zero elements: 4
domain_checksame ERROR: output domain size and gridmap destination size are not the same size
domain size = 259200
map dst size= 1
abort:
ERROR in mksurfdata_map: 34304

 

[1057398146@qq_com]

Excese me ,I have met similar question to you,I may ask do you have solve the problem?

 

[1057398146@qq_com]

Excese me ,I have met similar question to you,I may ask do you have solve the problem?

 

[yfan1@uni-goettingen_de]

Hi Ruby,i fixed the problem but it was a long time ago. I cannot really remember the details but the issue was "output domain size and gridmap destination size are not the same size', so please check the input data and output setting and make sure that they have the same grid size. If you used old dataset from another point simulation, that might be the problem.sorry i cannot remember more,Fan

 

[yfan1@uni-goettingen_de]

Hi Ruby,i fixed the problem but it was a long time ago. I cannot really remember the details but the issue was "output domain size and gridmap destination size are not the same size', so please check the input data and output setting and make sure that they have the same grid size. If you used old dataset from another point simulation, that might be the problem.sorry i cannot remember more,Fan

 

[1057398146@qq_com]

Dear Fan :Thank you very much for your quick reply.I'm sorry about I don't konw the specific meaning about what you said about the input data and output setting.The error occours in making lake parameters:(gridmap_map_read) * file name                  : /vol6/home/cmip_fang/cesm1_2_1/models/lnd/clm/tools/shared/mkmapdata/map_3x3min_MODIS_to_1x1_maduo_nomask_aave_da_c141217.nc * matrix dimensions rows x cols :    25920000  x           1 * number of non-zero elements:           10 domain_checksame ERROR: input domain lon and gridmap lon not the same at ni =     18019947  domain  lon=   0.000000000000000E+000  gridmap lon=    97.3250122070312Does that mean my input data is map_3x3min_MODIS_to_1x1_maduo_nomask_aave_da_c141217.nc?But what's the outputdata ?Can I setting it?I really what to know which data it refers to about the so-called domain lon? I'm looking forward to your answer.Thank you very much!

 

[1057398146@qq_com]

Dear Fan :Thank you very much for your quick reply.I'm sorry about I don't konw the specific meaning about what you said about the input data and output setting.The error occours in making lake parameters:(gridmap_map_read) * file name                  : /vol6/home/cmip_fang/cesm1_2_1/models/lnd/clm/tools/shared/mkmapdata/map_3x3min_MODIS_to_1x1_maduo_nomask_aave_da_c141217.nc * matrix dimensions rows x cols :    25920000  x           1 * number of non-zero elements:           10 domain_checksame ERROR: input domain lon and gridmap lon not the same at ni =     18019947  domain  lon=   0.000000000000000E+000  gridmap lon=    97.3250122070312Does that mean my input data is map_3x3min_MODIS_to_1x1_maduo_nomask_aave_da_c141217.nc?But what's the outputdata ?Can I setting it?I really what to know which data it refers to about the so-called domain lon? I'm looking forward to your answer.Thank you very much!

 

[yangk2010@lzu_edu_cn]

Ruby,Hello !When I created my own surface dataset for CLM4.5, a problem  confused  me .I used mkmapdata.sh script to create mkmapdata, "1km-merge-10min_HYDRO1K-merge-nomask"  could not be made, because it could not known"UGRID_1km-merge-10min_HYDRO1K-merge-nomask_c130402.nc" ,like this :UGRID_1km-merge-10min_HYDRO1K-merge-nomask_c130402.nc: NetCDF: Unknown file formatIf you can solve  this problem, could you give me a helpThank you in advance. 

 

[yangk2010@lzu_edu_cn]

Ruby,Hello !When I created my own surface dataset for CLM4.5, a problem  confused  me .I used mkmapdata.sh script to create mkmapdata, "1km-merge-10min_HYDRO1K-merge-nomask"  could not be made, because it could not known"UGRID_1km-merge-10min_HYDRO1K-merge-nomask_c130402.nc" ,like this :UGRID_1km-merge-10min_HYDRO1K-merge-nomask_c130402.nc: NetCDF: Unknown file formatIf you can solve  this problem, could you give me a helpThank you in advance.