Attempting to create surface boundary data ..... ------------------------------------------------------------ Grid_name: Cols0000486 Input rawdata files and corresponding meshes PFTs from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.25x0.25.LUH2.histsimyr1850-2015.c20230226/mksrf_landuse_ctsm52_histLUH2_2000.c20230226.nc mesh for pft /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc percent lake from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/lake_area/mksurf_lake_0.05x0.05_hist_clm5_hydrolakes_2000.cdf5.c20220325.nc mesh for percent lake /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc lake depth from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_LakePnDepth_3x3min_simyr2017_MODISgrid.cdf5.c200305.nc mesh for lake depth /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc inland wetland from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_lanwat.050425.nc mesh for wetland /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc soil texture mapunits from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_soil_mapunits_5x5min_WISE.c220330.nc soil texture (sand/clay, orgc) lookup: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_soil_lookup.10level.WISE.c220330.nc mesh for soil texture /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_5x5min_nomask_cdf5_c200129.nc soil color from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.25x0.25.LUH2.histsimyr1850-2015.c20230226/mksrf_soilcolor_ctsm52_histLUH2_2005.c20230226.nc mesh for soil color /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc fmax from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_fmax_0.125x0.125_c200220.nc mesh for fmax /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.125x0.125_nomask_cdf5_c200129.nc soil depth from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksf_soilthk_5x5min_ORNL-Soil_simyr1900-2015_c170630.nc mesh for soil depth /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_5x5min_nomask_cdf5_c200129.nc VOC emission factors from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_vocef_0.5x0.5_simyr2000.c110531.nc mesh for VOC pct emis /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc gdp from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_gdp_0.5x0.5_AVHRR_simyr2000.c130228.nc mesh for gdp /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc peat from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_peatf_0.5x0.5_AVHRR_simyr2000.c130228.nc mesh for peatlands /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc harvest from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.25x0.25.LUH2.histsimyr1850-2015.c20230226/mksrf_landuse_ctsm52_histLUH2_2000.c20230226.nc mesh for harvest /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc topography statistics from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_topostats_1km-merge-10min_HYDRO1K-merge-nomask_simyr2000.c130402.nc mesh for topography stats /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UGRID_1km-merge-10min_HYDRO1K-merge-nomask_cdf5_c130402.nc glaciers from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_glacier_3x3min_simyr2000.c20210708.nc with: 10 glacier elevation classes mesh for glacier /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc glacier region ID from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_GlacierRegion_10x10min_nomask_cd5_c220131.nc mesh for glacier region /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_10x10min_nomask_cdf5_c200129.nc urban from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/gao_oneill_urban/historical/urban_properties_GaoOneil_05deg_ThreeClass_2000_cdf5_c20220910.nc mesh for urban /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc urban topography from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_topo.10min.cdf5.c220201.nc mesh for urban topography /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_10x10min_nomask_cdf5_c200129.nc mesh for lai/sai /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc mesh for ag fire pk month /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc abm from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_abm_0.5x0.5_AVHRR_simyr2000.c130201.nc VIC parameters from: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_vic_0.9x1.25_GRDC_simyr2000.c130307.nc mesh for VIC parameters /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.9x1.25_nomask_cdf5_c200129.nc ================================================================================ Model grid configuration variables mksrf_fgrid_mesh = /glade/p/cesm/cseg/inputdata/share/meshes/ne3pg3_ESMFmesh_c221214_cdf5.nc nlon= 486 nlat= 1 Grid_size: 486 ================================================================================ Output configuration variables output file is 1d format Output ALL data in file as 64-bit Output VIC fields Only output ESMF PET files if fatal errors happen in ESMF npes = 64 grid size = 486 Opening file and reading pio_iotype from txt file with the same name (mkpftInit) num_natpft = 14 (mkpftInit) natpft_lb = 0 (mkpftInit) natpft_ub = 14 (mkpftInit) num_cft = 64 (mkpftInit) cft_lb = 15 (mkpftInit) cft_ub = 78 Creating surface datasets with extra types for crops; total pfts = 78 ================================================================================ opening output file surfdata_ne3np4.pg3_hist_2000_78pfts_c240130.nc (mkpio_wopen) creating file surfdata_ne3np4.pg3_hist_2000_78pfts_c240130.nc ================================================================================ Attempting to create model lats and lons from model mesh ..... Successfully made model lats and lons mksrfdata writing out model grid mksrfdata writing out LONGXY mksrfdata writing out LATIXY Initialization is complete, going on to process types of input files Attempting to make LAIs/SAIs/heights ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.25x0.25.LUH2.histsimyr1850-2015.c20230226/mksrf_pftlaihgt_ctsm52_histLUH2_2005.c20230226.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc WARNING: mklai(): parameter numpft+1 = 79 does not equal input dataset numpft = 16 This inconsistency used to stop the program. Now we allow it because crop pfts 17-last are assumed to never use satellite lai data. LAI Output for month 1 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 2.486 2.486 2 1.766 1.766 3 0.332 0.332 4 4.737 4.737 5 3.807 3.807 6 1.676 1.676 7 1.217 1.217 8 0.898 0.898 9 0.462 0.462 10 0.474 0.474 11 0.183 0.183 12 0.356 0.356 13 0.797 0.797 14 0.933 0.933 15 1.499 1.499 Successfully made LAIs/SAIs/heights for month 1 LAI Output for month 2 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 2.392 2.392 2 1.626 1.626 3 0.305 0.305 4 4.727 4.727 5 3.911 3.911 6 1.695 1.695 7 1.268 1.268 8 0.911 0.911 9 0.484 0.484 10 0.521 0.521 11 0.188 0.188 12 0.372 0.372 13 0.813 0.813 14 0.962 0.962 15 1.554 1.554 Successfully made LAIs/SAIs/heights for month 2 LAI Output for month 3 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 2.357 2.357 2 1.441 1.441 3 0.241 0.241 4 4.778 4.778 5 3.986 3.986 6 1.815 1.815 7 1.252 1.252 8 0.714 0.714 9 0.607 0.607 10 0.587 0.587 11 0.177 0.177 12 0.341 0.341 13 0.865 0.865 14 1.043 1.043 15 1.646 1.646 Successfully made LAIs/SAIs/heights for month 3 LAI Output for month 4 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 2.479 2.479 2 1.368 1.368 3 0.228 0.228 4 4.879 4.879 5 4.179 4.179 6 2.024 2.024 7 1.473 1.473 8 0.187 0.187 9 0.741 0.741 10 0.652 0.652 11 0.163 0.163 12 0.315 0.315 13 0.972 0.972 14 1.180 1.180 15 1.818 1.818 Successfully made LAIs/SAIs/heights for month 4 LAI Output for month 5 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 3.008 3.008 2 1.672 1.672 3 0.990 0.990 4 4.841 4.841 5 4.206 4.206 6 2.307 2.307 7 2.315 2.315 8 0.829 0.829 9 0.768 0.768 10 0.740 0.740 11 0.369 0.369 12 0.471 0.471 13 1.194 1.194 14 1.330 1.330 15 2.116 2.116 Successfully made LAIs/SAIs/heights for month 5 LAI Output for month 6 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 3.712 3.712 2 2.833 2.833 3 2.242 2.242 4 4.781 4.781 5 4.054 4.054 6 2.565 2.565 7 2.863 2.863 8 2.466 2.466 9 0.633 0.633 10 0.887 0.887 11 0.945 0.945 12 0.962 0.962 13 1.506 1.506 14 1.415 1.415 15 2.552 2.552 Successfully made LAIs/SAIs/heights for month 6 LAI Output for month 7 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 3.847 3.847 2 3.400 3.400 3 2.776 2.776 4 4.888 4.888 5 4.156 4.156 6 2.835 2.835 7 2.898 2.898 8 3.025 3.025 9 0.537 0.537 10 0.982 0.982 11 1.393 1.393 12 1.360 1.360 13 1.636 1.636 14 1.522 1.522 15 2.814 2.814 Successfully made LAIs/SAIs/heights for month 7 LAI Output for month 8 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 3.638 3.638 2 3.066 3.066 3 2.286 2.286 4 5.089 5.089 5 4.273 4.273 6 2.974 2.974 7 2.769 2.769 8 2.374 2.374 9 0.504 0.504 10 0.916 0.916 11 1.095 1.095 12 1.056 1.056 13 1.452 1.452 14 1.531 1.531 15 2.587 2.587 Successfully made LAIs/SAIs/heights for month 8 LAI Output for month 9 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 3.305 3.305 2 2.644 2.644 3 0.540 0.540 4 5.220 5.220 5 4.253 4.253 6 3.087 3.087 7 2.234 2.234 8 0.555 0.555 9 0.462 0.462 10 0.687 0.687 11 0.521 0.521 12 0.531 0.531 13 1.103 1.103 14 1.364 1.364 15 2.050 2.050 Successfully made LAIs/SAIs/heights for month 9 LAI Output for month 10 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 3.017 3.017 2 2.345 2.345 3 0.166 0.166 4 5.155 5.155 5 4.140 4.140 6 2.909 2.909 7 1.491 1.491 8 0.064 0.064 9 0.430 0.430 10 0.541 0.541 11 0.312 0.312 12 0.355 0.355 13 0.923 0.923 14 1.159 1.159 15 1.710 1.710 Successfully made LAIs/SAIs/heights for month 10 LAI Output for month 11 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 2.797 2.797 2 2.112 2.112 3 0.276 0.276 4 5.050 5.050 5 3.994 3.994 6 2.442 2.442 7 1.181 1.181 8 0.308 0.308 9 0.480 0.480 10 0.478 0.478 11 0.195 0.195 12 0.284 0.284 13 0.855 0.855 14 1.046 1.046 15 1.569 1.569 Successfully made LAIs/SAIs/heights for month 11 LAI Output for month 12 ...................................................................... PFT input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... 0 0.000 0.000 1 2.642 2.642 2 1.918 1.918 3 0.332 0.332 4 4.840 4.840 5 3.849 3.849 6 1.907 1.907 7 1.181 1.181 8 0.363 0.363 9 0.512 0.512 10 0.456 0.456 11 0.173 0.173 12 0.305 0.305 13 0.810 0.810 14 0.971 0.971 15 1.498 1.498 Successfully made LAIs/SAIs/heights for month 12 ================================================================================ Attempting to make PFTs ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.25x0.25.LUH2.histsimyr1850-2015.c20230226/mksrf_landuse_ctsm52_histLUH2_2000.c20230226.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc dimid, length= 1 1440 dimid, length= 2 720 dimid, length= 3 64 num_cft = 64 ...................................................................... PFTs Output plant type output grid area 10**6 km**2 ...................................................................... not vegetated 23.032 needleleaf evergreen temperate tree 3.790 needleleaf evergreen boreal tree 5.426 needleleaf deciduous boreal tree 2.005 broadleaf evergreen tropical tree 11.441 broadleaf evergreen temperate tree 0.846 broadleaf deciduous tropical tree 1.716 broadleaf deciduous temperate tree 2.162 broadleaf deciduous boreal tree 0.901 broadleaf evergreen shrub 0.023 broadleaf deciduous temperate shrub 3.737 broadleaf deciduous boreal shrub 3.901 c3 arctic grass 6.509 c3 non-arctic grass 13.810 c4 grass 16.811 c3_crop 0.000 c3_irrigated 0.000 temperate_corn 0.688 irrigated_temperate_corn 0.119 spring_wheat 1.686 irrigated_spring_wheat 0.349 winter_wheat 0.000 irrigated_winter_wheat 0.000 temperate_soybean 0.498 irrigated_temperate_soybean 0.047 barley 0.446 irrigated_barley 0.039 winter_barley 0.000 irrigated_winter_barley 0.000 rye 0.097 irrigated_rye 0.001 winter_rye 0.000 irrigated_winter_rye 0.000 cassava 0.115 irrigated_cassava 0.003 citrus 0.044 irrigated citrus 0.016 cocoa 0.050 irrigated_cocoa 0.001 coffee 0.060 irrigated_coffee 0.005 cotton 0.167 irrigated_cotton 0.102 datepalm 0.004 irrigated_datepalm 0.005 foddergrass 3.094 irrigated_foddergrass 0.272 grapes 0.037 irrigated_grapes 0.014 groundnuts 0.128 irrigated_groundnuts 0.023 millet 0.300 irrigated_millet 0.033 oilpalm 0.051 irrigated_oilpalm 0.003 potatoes 0.168 irrigated_potatoes 0.017 pulses 0.549 irrigated_pulses 0.111 rapeseed 0.197 irrigated_rapeseed 0.041 rice 0.791 irrigated_rice 0.459 sorghum 0.338 irrigated_sorghum 0.040 sugarbeet 0.044 irrigated_sugarbeet 0.006 sugarcane 0.096 irrigated_sugarcane 0.037 sunflower 0.160 irrigated_sunflower 0.012 miscanthus 0.000 irrigated_miscanthus 0.000 switchgrass 0.000 irrigated_switchgrass 0.000 tropical_corn 0.472 irrigated_tropical_corn 0.050 tropical_soybean 0.197 irrigated_tropical_soybean 0.025 Successfully made PFTs mksrfdata writing land fraction from pft dataset ================================================================================ Attempting to make harvest fields ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.25x0.25.LUH2.histsimyr1850-2015.c20230226/mksrf_landuse_ctsm52_histLUH2_2000.c20230226.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc reading: HARVEST_VH1 mkharvest writing out 1d CONST_HARVEST_VH1 ...................................................................... diagnostics for CONST_HARVEST_VH1 surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... CONST_HARVES 17.605 17.605 all surface 510.100 510.100 reading: HARVEST_VH2 mkharvest writing out 1d CONST_HARVEST_VH2 ...................................................................... diagnostics for CONST_HARVEST_VH2 surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... CONST_HARVES 7.348 7.348 all surface 510.100 510.100 reading: HARVEST_SH1 mkharvest writing out 1d CONST_HARVEST_SH1 ...................................................................... diagnostics for CONST_HARVEST_SH1 surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... CONST_HARVES 9.319 9.319 all surface 510.100 510.100 reading: HARVEST_SH2 mkharvest writing out 1d CONST_HARVEST_SH2 ...................................................................... diagnostics for CONST_HARVEST_SH2 surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... CONST_HARVES 13.619 13.619 all surface 510.100 510.100 reading: HARVEST_SH3 mkharvest writing out 1d CONST_HARVEST_SH3 ...................................................................... diagnostics for CONST_HARVEST_SH3 surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... CONST_HARVES 20.124 20.124 all surface 510.100 510.100 reading: GRAZING mkharvest writing out 1d CONST_GRAZING ...................................................................... diagnostics for CONST_GRAZING surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... CONST_GRAZIN -15678.115 -15678.115 all surface 510.100 510.100 reading: FERTNITRO_CFT mkharvest writing out 2d CONST_FERTNITRO_CFT reading: UNREPRESENTED_PFT_LULCC mkharvest writing out 2d UNREPRESENTED_PFT_LULCC reading: UNREPRESENTED_CFT_LULCC mkharvest writing out 2d UNREPRESENTED_CFT_LULCC Successfully made harvest and grazing ================================================================================ Attempting to make %lake Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/lake_area/mksurf_lake_0.05x0.05_hist_clm5_hydrolakes_2000.cdf5.c20220325.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc Attempting to make %lake ..... ...................................................................... diagnostics for pct lake surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... pct lake 2.916 2.171 all surface 510.100 510.100 Successfully made %lake ================================================================================ Attempting to make lake depth Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_LakePnDepth_3x3min_simyr2017_MODISgrid.cdf5.c200305.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc Attempting to make lake depth ..... mklakdep writing out lakedepth ...................................................................... output diagnostics for lake depth with units = m global input sum= 11.048 global output sum 11.048 ...................................................................... Successfully made lake depth Attempting to make %wetland ..... mkwetlnd setting PCT_WETLAND to zero ================================================================================ Attempting to make %glacier ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_glacier_3x3min_simyr2000.c20210708.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc ...................................................................... diagnostics for pct glacier surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... pct glacier 15.991 12.217 all surface 510.100 510.100 Successfully made %glacier ================================================================================ Attempting to make glacier region ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_GlacierRegion_10x10min_nomask_cd5_c220131.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_10x10min_nomask_cdf5_c200129.nc mkglacierregion writing out glacier_region ............................................................ index input grid area output grid area 10**6 km**2 10**6 km**2 ............................................................ 0 0.000 0.000 1 1.904 2.269 2 2.329 1.963 3 34.170 34.170 Successfully made glacier region ================================================================================ Attempting to make %sand, %clay, orgc, cfrag, bulk, phaq ..... Input mapunit file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_soil_mapunits_5x5min_WISE.c220330.nc Input lookup table file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_soil_lookup.10level.WISE.c220330.nc Input mesh/grid file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_5x5min_nomask_cdf5_c200129.nc reading mesh_i directly in mksoiltex Reading in mapunit data in mksoiltex Setting mask in mesh where mapunit data is 0 mksoiltex before route handle creation mksoiltex after route handle creation mksoiltex before call to dynamic mask set creation mksoiltex after call to dynamic mask set creation mksoiltex mksoiltex writing out mapunits WARNING: assigning sand_o = -4 to 99% WARNING: assigning other sand_o < 0 to 43% WARNING: assigning clay_o = -4 to 1% WARNING: assigning other clay_o < 0 to 18% WARNING: assigning orgc_o = -4 to 1 WARNING: assigning other orgc_o < 0 to 0 WARNING: same warnings for cfrag_o as for orgc_o WARNING: assigning bulk_o < 0 to 1.5 WARNING: assigning phaq_o < 0 to 7 mksoiltex writing out soil percent sand mksoiltex writing out soil percent clay mksoiltex writing out soil organic matter mksoiltex writing out soil organic carbon content mksoiltex writing out coarse fragments in soil mksoiltex writing out soil bulk density mksoiltex writing out soil pH measured in H2O Successfully made %sand, %clay, orgc, cfrag, bulk, phaq ..... ================================================================================ Attempting to make soil color classes ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.25x0.25.LUH2.histsimyr1850-2015.c20230226/mksrf_soilcolor_ctsm52_histLUH2_2005.c20230226.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.25x0.25_nomask_cdf5_c200129.nc mksoilcol writing out soil color mksoilcol writing out mksoil_color ............................................................ index input grid area output grid area 10**6 km**2 10**6 km**2 ............................................................ 0 0.000 0.000 1 1.826 2.198 2 1.647 6.445 3 1.674 2.033 4 1.766 3.075 5 1.707 1.112 6 1.571 0.000 7 1.587 0.000 8 1.916 0.949 9 2.404 1.065 10 20.101 44.911 11 4.652 3.074 12 7.168 4.767 13 10.450 11.748 14 14.405 33.599 15 20.822 216.461 16 23.782 64.184 17 13.615 21.481 18 9.997 23.065 19 9.341 37.152 20 6.368 32.780 Successfully made soil color classes ================================================================================ Attempting to make %fmax ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_fmax_0.125x0.125_c200220.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.125x0.125_nomask_cdf5_c200129.nc ...................................................................... diagnostics for Max Fractional Sataturated Area surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... Max Fraction 0.629 0.629 all surface 510.100 510.100 mksoilfmax writing out soil fmax (maximum fraction saturated area) Successfully made %fmax ================================================================================ Attempting to make GDP..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_gdp_0.5x0.5_AVHRR_simyr2000.c130228.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc mkgdp writing out gdp ...................................................................... output diagnostics for GDP with units = x1000 US$ per capita global input sum= 7.679 global output sum 7.679 ...................................................................... Successfully made GDP ================================================================================ Attempting to make peat ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_peatf_0.5x0.5_AVHRR_simyr2000.c130228.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc mkpeat writing out peatland fraction ...................................................................... diagnostics for Peat surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... Peat 2.157 2.157 all surface 510.100 510.100 Successfully made peat Attempting to make soildepth ..... mksoildepth writing out soil depth ...................................................................... diagnostics for Soildepth surface type input grid area output grid area 10**6 km**2 10**6 km**2 ...................................................................... Soildepth 2048.572 2048.572 all surface 510.100 510.100 Successfully made soildepth ================================================================================ Attempting to make agricultural fire peak month data ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_abm_0.5x0.5_AVHRR_simyr2000.c130201.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc mkagfirepkmon writing abm (agricultural fire peak month) ............................................................ index input grid area output grid area 10**6 km**2 10**6 km**2 ............................................................ 1 2.516 7.498 2 2.653 8.316 3 3.780 17.348 4 5.981 32.346 5 4.796 24.014 6 2.794 12.360 7 4.168 15.449 8 6.095 29.400 9 5.286 17.098 10 4.286 13.518 11 2.686 6.482 12 2.669 7.533 13 115.921 318.737 Successfully made Agricultural fire peak month ================================================================================ Attempting to make %urban ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/gao_oneill_urban/historical/urban_properties_GaoOneil_05deg_ThreeClass_2000_cdf5_c20220910.nc Successfully made %urban ================================================================================ Attempting to make urban region ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/gao_oneill_urban/historical/urban_properties_GaoOneil_05deg_ThreeClass_2000_cdf5_c20220910.nc max urban regions = 33 Successfully made urban region ............................................................ index input grid area output grid area 10**6 km**2 10**6 km**2 ............................................................ 1 1.582 1.295 2 7.809 7.790 3 8.559 8.623 4 6.609 6.483 5 5.680 5.525 6 10.158 11.031 7 0.295 0.287 8 9.436 9.489 9 5.824 5.725 10 0.712 0.746 11 2.181 1.520 12 2.251 2.238 13 3.194 3.329 14 2.548 2.510 15 6.032 6.969 16 9.509 8.724 17 1.240 0.909 18 0.842 0.124 19 1.774 1.756 20 1.281 1.415 21 1.025 1.207 22 17.019 17.180 23 3.087 3.267 24 1.943 1.000 25 1.405 1.923 26 1.296 1.752 27 4.844 4.371 28 1.122 1.509 29 1.801 1.511 30 4.210 3.936 31 5.183 5.267 32 5.043 5.717 33 1.118 1.484 mksrfdata writing out urban region id ================================================================================ Attempting to make urban topo elevation ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_topo.10min.cdf5.c220201.nc ...................................................................... output diagnostics for Urban elev variable with units = m global input sum= 230.498 global output sum 230.498 ...................................................................... Successfully made elevation ================================================================================ Attempting to make Topography statistics..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_topostats_1km-merge-10min_HYDRO1K-merge-nomask_simyr2000.c130402.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UGRID_1km-merge-10min_HYDRO1K-merge-nomask_cdf5_c130402.nc mktopostats creating a routehandle mktopostats finished creating a routehandle ...................................................................... output diagnostics for Topo Std Dev with units = m global input sum= 802.194 global output sum 233.081 ...................................................................... ...................................................................... output diagnostics for Slope with units = degrees global input sum= 1.483 global output sum 1.166 ...................................................................... mktopostats writing topo_stddev mktopostats writing slope Successfully made Topography statistics ================================================================================ Attempting to make VIC parameters..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_vic_0.9x1.25_GRDC_simyr2000.c130307.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.9x1.25_nomask_cdf5_c200129.nc ...................................................................... output diagnostics for VIC b parameter with units = unitless global input sum= 0.099 global output sum 0.099 ...................................................................... ...................................................................... output diagnostics for VIC Ws parameter with units = unitless global input sum= 0.734 global output sum 0.734 ...................................................................... ...................................................................... output diagnostics for VIC Dsmax parameter with units = mm/day global input sum= 26.410 global output sum 26.410 ...................................................................... ...................................................................... output diagnostics for VIC Ds parameter with units = unitless global input sum= 0.100 global output sum 0.100 ...................................................................... mkVICparams writing out VIC parameters Successfully made VIC parameters ================================================================================ Attempting to make VOC emission factors ..... Input data file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_vocef_0.5x0.5_simyr2000.c110531.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_0.5x0.5_nomask_cdf5_c200129.nc mkvocef writing out voc emission factors Successfully made VOC Emission Factors sum over domain of pft 0 5641.64718259635 sum over domain of pft 1 1098.93039578132 sum over domain of pft 2 967.460658856876 sum over domain of pft 3 240.613502926321 sum over domain of pft 4 4623.77333524818 sum over domain of pft 5 600.677310495383 sum over domain of pft 6 457.425024621020 sum over domain of pft 7 565.174601865952 sum over domain of pft 8 236.629880672086 sum over domain of pft 9 7.68017334306267 sum over domain of pft 10 919.791653308026 sum over domain of pft 11 801.958110402672 sum over domain of pft 12 2422.28954091748 sum over domain of pft 13 3664.57142420037 sum over domain of pft 14 5854.57104004200 sum over domain of cft 15 0.000000000000000E+000 sum over domain of cft 16 0.000000000000000E+000 sum over domain of cft 17 94.6395178219310 sum over domain of cft 18 20.5145011447597 sum over domain of cft 19 257.326579764620 sum over domain of cft 20 53.4069777795058 sum over domain of cft 21 0.000000000000000E+000 sum over domain of cft 22 0.000000000000000E+000 sum over domain of cft 23 66.5604728337356 sum over domain of cft 24 9.49097711868583 sum over domain of cft 25 79.3239203630029 sum over domain of cft 26 7.01599800502267 sum over domain of cft 27 0.000000000000000E+000 sum over domain of cft 28 0.000000000000000E+000 sum over domain of cft 29 10.8197692014652 sum over domain of cft 30 0.283397110879167 sum over domain of cft 31 0.000000000000000E+000 sum over domain of cft 32 0.000000000000000E+000 sum over domain of cft 33 32.5637395583574 sum over domain of cft 34 2.22994989044559 sum over domain of cft 35 10.9500785216857 sum over domain of cft 36 4.01051846412920 sum over domain of cft 37 24.7367243828886 sum over domain of cft 38 1.53644494473786 sum over domain of cft 39 29.0912866373941 sum over domain of cft 40 3.82758027919017 sum over domain of cft 41 35.3286513115782 sum over domain of cft 42 18.0042583055540 sum over domain of cft 43 0.585961572467184 sum over domain of cft 44 0.618691628717020 sum over domain of cft 45 708.709784302611 sum over domain of cft 46 79.2485206836792 sum over domain of cft 47 10.0592946997224 sum over domain of cft 48 3.53691258278960 sum over domain of cft 49 37.7304466769359 sum over domain of cft 50 6.76424148262196 sum over domain of cft 51 45.7389748833770 sum over domain of cft 52 5.21554892467355 sum over domain of cft 53 28.5222517586780 sum over domain of cft 54 1.17380739995150 sum over domain of cft 55 24.5642313455394 sum over domain of cft 56 3.34585667988274 sum over domain of cft 57 100.709177586017 sum over domain of cft 58 21.4195207548954 sum over domain of cft 59 31.8878791882411 sum over domain of cft 60 5.89867629665545 sum over domain of cft 61 242.922919418953 sum over domain of cft 62 129.003496854350 sum over domain of cft 63 53.0610681928877 sum over domain of cft 64 8.47216847694177 sum over domain of cft 65 7.35415400225224 sum over domain of cft 66 1.02110102356800 sum over domain of cft 67 54.3285271431544 sum over domain of cft 68 14.8515929536043 sum over domain of cft 69 26.6519627175583 sum over domain of cft 70 2.57039489910371 sum over domain of cft 71 0.000000000000000E+000 sum over domain of cft 72 0.000000000000000E+000 sum over domain of cft 73 0.000000000000000E+000 sum over domain of cft 74 0.000000000000000E+000 sum over domain of cft 75 125.640923308026 sum over domain of cft 76 17.2702091266790 sum over domain of cft 77 25.6976515763794 sum over domain of cft 78 4.59938780280288 mksrfdata writing out percnt urban Attempting to make Urban Parameters ..... Opening input urban parameter file: /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/gao_oneill_urban/historical/urban_properties_GaoOneil_05deg_ThreeClass_2000_cdf5_c20220910.nc Successfully made Urban Parameters mksrfdata writing out PCT_URBAN mksrfdata writing out PCT_GLACIER mksrfdata writing out PCT_LAKE mksrfdata writing out PCT_WETLAND mksrfdata writing out PCT_OCEAN mksrfdata writing PCT_NATVEG mksrfdata writing PCT_CROP mksrfdata writing PCT_NAT_PFT mksrfdata writing PCT_CFT mksrfdata writing LANDFRAC_MKSURFDATA mkglcmecInit writing out GLC_MEC ================================================================================ Attempting to make percent elevation class and mean elevation for glaciers ..... Input file is /glade/p/cesm/cseg/inputdata/lnd/clm2/rawdata/mksrf_glacier_3x3min_simyr2000.c20210708.nc Input mesh file is /glade/p/cesm/cseg/inputdata/lnd/clm2/mappingdata/grids/UNSTRUCTgrid_3x3min_nomask_cdf5_c200129.nc mkglcmec writing out pct_glc_mec mkglcmec writing out topo_glc_mec Successfully made percent elevation class and mean elevation for glaciers Successfully created surface data output file = surfdata_ne3np4.pg3_hist_2000_78pfts_c240130.nc This file contains the land model surface data Successfully ran mksurfdata_esmf