Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

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!
 
Top