Problems of generating surfdata using mksurfdatapl tool (from cesm1_0_3) on yellowstone


Code:
Hi all,<br /><br />I would like to run experiment using CESM (cesm1_0_3) with RCP8.5 T31 and I understand that T31 is not supported in the default configuration. <br />So I have to generate my own surface data set in order to get the model run.<br /><br />I followed the user guide to create the surface data set<br />(<a class="text" style="text-decoration: underline;" href="https://webmail.psu.edu/webmail/retrieve.cgi?mailbox=inbox&start_num=0&limit=50&sort=1&display=4&timestamp=20160113150322&mid=638817756%2e37405%2e1452621131507%2eJavaMail%2etomcat%2dextraview%40extraview%2dtomcat%2eucar%2eedu">http://www.cesm.ucar.edu/models/cesm1.0/clm/models/lnd/clm/doc/UsersGuide/x8689.html</a>)<br />under my directory:<br />/glade/scratch/cytsai/CESM_Source/cesm1_0_3/models/lnd/clm/tools/mksurfdata<br /><br />using ./mksurfdata.pl -dinlc /glade/p/cesmdata/cseg/inputdata -years<br />1850-2100 -rcp 8.5 -r 48x96 -nomv -pftlc<br />/glade/p/cesmdata/cseg/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630<br /><br /><br />However, it seemed that the tool cannot give out appropriate "<br />pftdyn_rcp8.5_simyr1850-2100.txt" file (you can see at:/glade/scratch/cytsai/CESM_Source/cesm1_0_3/models/lnd/clm/tools/mksurfdata/pftdyn_rcp8.5_simyr1850-<br />2100.txt), and had the error messages like this:<br /><br />land model surface data set successfully created for   3.8 by   3.8 grid<br /><br /> (OPNFIL): Successfully opened file pftdyn_rcp8.5_simyr1850-2100.txt on<br />unit=<br /><br />           98<br /><br /> input pft dynamic dataset is<br /><br /><br />/glade/p/cesmdata/cseg/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf<br />_landuse_rc1850_c090630.n<br /><br />  year is          1850<br /><br /> Attempting to make PFTs .....<br /><br /> Creating surface datasets with the standard # of PFTs =           16<br /><br /><br /><br /> netcdf error from read_domain rcode =             2  error =<br /><br /> No such file or directory<br /><br /><br />sh: line 1: 26738 Aborted                 (core dumped) ./mksurfdata <<br />namelist<br /><br />ERROR in mksurfdata: 34304<br /><br /><br /><br />So, as in the file "pftdyn_rcp8.5_simyr1850-2100.txt", it seemed that the<br />length of the path of data are too long so that all the ".nc" are missing<br />"c" and that's probably why the mksurfdata tool cannot find the file.<br /><br />I do not know how to fix this because the length of inputdata path on<br />Yellowstone is long....<br /><br /><br />Also, I tried to use the mksurfdata tool in cesm1_2_2 (clm4_5), and it worked.<br /><br />cd<br />/glade/scratch/cytsai/CESM_Source/cesm1_2_2/models/lnd/clm/tools/clm4_5/mksurfdata_map<br /><br />using<br /><br />./mksurfdata.pl -dinlc /glade/p/cesmdata/cseg/inputdata -years 1850-2100<br />-rcp 8.5 -r 48x96 -nomv -pftlc<br />/glade/p/cesmdata/cseg/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630<br /><br />And there's no problem with the length in .txt file...<br /><br />However, I cannot use the "surfdata_48x96_rcp8.5_simyr1850_c160111.nc" I<br />generated from cesm1_2_2 because when I run cesm1_0_3, it couldn't read the<br />file... I'm wondering this is because I used newer version of CESM to<br />generate the surface data so I would like to ask if anyone can help me<br />generate surface data set using older version of CESM?<br /><br /><br />Any help would be appreciated! Thank you
 
Hi,Im getting a similar error. while creating my own surfacedata using ./mksurfdata.pl -res 0.9x1.25 -crop -year 2000 for CESM1.2.1 usib=ng CLM4.5.on running Im getting the following error:resolution: 0.9x1.25 rcp=-999.9 sim_year = 2000
namelist: namelist
&clmexp
 nglcec           = 0
 mksrf_fgrid      = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_landuse_to_0.9x1.25_aave_da_110307.nc'
 map_fpft         = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_landuse_to_0.9x1.25_aave_da_110307.nc'
 map_fglacier     = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_3x3min_GLOBE-Gardner_to_0.9x1.25_nomask_aave_da_c120923.nc'
 map_fsoicol      = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_landuse_to_0.9x1.25_aave_da_110307.nc'
 map_furban       = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_3x3min_LandScan2004_to_0.9x1.25_nomask_aave_da_c120522.nc'
 map_fmax         = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_3x3min_USGS_to_0.9x1.25_nomask_aave_da_c120926.nc'
 map_forganic     = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_5x5min_ISRIC-WISE_to_0.9x1.25_nomask_aave_da_c120525.nc'
 map_flai         = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_landuse_to_0.9x1.25_aave_da_110307.nc'
 map_fharvest     = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_landuse_to_0.9x1.25_aave_da_110307.nc'
 map_flakwat      = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_3x3min_MODIS_to_0.9x1.25_nomask_aave_da_c120523.nc'
 map_fwetlnd      = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_lanwat_to_0.9x1.25_aave_da_110307.nc'
 map_fvocef       = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_lanwat_to_0.9x1.25_aave_da_110307.nc'
 map_fsoitex      = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_5minx5min_soitex_to_0.9x1.25_aave_da_110722.nc'
 map_furbtopo     = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_10minx10min_topo_to_0.9x1.25_aave_da_110630.nc'
 map_flndtopo     = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_10minx10min_topo_to_0.9x1.25_aave_da_110630.nc'
 map_fgdp         = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_lanwat_to_0.9x1.25_aave_da_110307.nc'
 map_fpeat        = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_lanwat_to_0.9x1.25_aave_da_110307.nc'
 map_fabm         = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_lanwat_to_0.9x1.25_aave_da_110307.nc'
 map_ftopostats   = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_1km-merge-10min_HYDRO1K-merge-nomask_to_0.9x1.25_nomask_aave_da_c130405.nc'
 map_fvic         = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.9x1.25_GRDC_to_0.9x1.25_nomask_aave_da_c130308.nc'
 map_fch4         = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_360x720_cruncep_to_0.9x1.25_nomask_aave_da_c130326.nc'
 mksrf_fsoitex    = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_soitex.10level.c010119.nc'
 mksrf_forganic   = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_organic_10level_5x5min_ISRIC-WISE-NCSCD_nlev7_c120830.nc'
 mksrf_flakwat    = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_LakePnDepth_3x3min_simyr2004_c111116.nc'
 mksrf_fwetlnd    = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_lanwat.050425.nc'
 mksrf_fmax       = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_fmax_3x3min_USGS_c120911.nc'
 mksrf_fglacier   = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_glacier_3x3min_simyr2000.c120926.nc'
 mksrf_fvocef     = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_vocef_0.5x0.5_simyr2000.c110531.nc'
 mksrf_furbtopo   = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_topo.10min.c080912.nc'
 mksrf_flndtopo   = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/topodata_10min_USGS_071205.nc'
 mksrf_fgdp       = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_gdp_0.5x0.5_AVHRR_simyr2000.c130228.nc'
 mksrf_fpeat      = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_peatf_0.5x0.5_AVHRR_simyr2000.c130228.nc'
 mksrf_fabm       = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_abm_0.5x0.5_AVHRR_simyr2000.c130201.nc'
 mksrf_ftopostats = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_topostats_1km-merge-10min_HYDRO1K-merge-nomask_simyr2000.c130402.nc'
 mksrf_fvic       = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_vic_0.9x1.25_GRDC_simyr2000.c130307.nc'
 mksrf_fch4       = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_ch4inversion_360x720_cruncep_simyr2000.c130322.nc'
 outnc_double   = .true.
 all_urban      = .false.
 no_inlandwet   = .true.
 mksrf_furban   = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/mksrf_urban_0.05x0.05_simyr2000.c120621.nc'
 mksrf_fvegtyp  = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_24pftirrNT_landuse_rc2000_c120712.nc'
 mksrf_fsoicol  = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_soilcol_global_c090324.nc'
 mksrf_flai     = '/scratch/cas/phd/asz122525/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_lai24pftnontrop_c120712.nc'
 fsurdat        = 'surfdata_0.9x1.25_simyr2000_c160505.nc'
 fsurlog        = 'surfdata_0.9x1.25_simyr2000_c160505.log'    
 mksrf_fdynuse  = ' '
 fdyndat        = ' '
 numpft = 24
/
./mksurfdata_map < namelist
 Attempting to initialize control settings .....
 Attempting to create surface boundary data .....
------------------------------------------------------------
 mksrf_fgrid =
 /scratch/cas/phd/asz122525/inputdata/lnd/clm2/mappingdata/maps/0.9x1.25/map_0.5x0.5_landuse_to_0.9x1.25_aave_da_110307.nc                                                                                                                                      
 mksrf_gridtype = global
 Output ALL data in file as 64-bit
 Set wetland to 0% over land
 netcdf error from mkurbanInit rcode =           -46  error =
 NetCDF: Invalid dimension ID or name                                           
sh: line 1: 25215 Aborted                 (core dumped) ./mksurfdata_map < namelist
ERROR in mksurfdata_map: 34304

  Please let me know if you were able to rectify the error.Any help will be appreciated.Thank you in advance!
 
Back
Top