mksurfdata_map netcdf error from mkharvest_init rcode (CESM 1.2.2, CLM4_0)

[matthias_demuzere@ees_kuleuven_be]

Hi all, I am currently creating my own surfdata set for a high-resolution regional domain centered over Belgium. I managed to created all the mapdata, and now I am iterating over all input files using the mksurfdata_map routine. This is currently done for -phys clm4_0 from cesm1_2_2.The namelist I use for this is attached.Now, after a while the script fails with the following message:calling mklai Attempting to make LAIs/SAIs/heights ..... domain_read read lon and lat dims domain_read initialized domain domain_read read LANDMASK Open LAI file:  /scratch/leuven/300/vsc30061/cesm/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x 0.5.simyr1850-2005.c090630/mksrf_lai_global_c090506.nc(gridmap_map_read) reading mapping matrix data...(gridmap_map_read) * file name                  : ../../shared/mkmapdata/map_0.5x0.5_MODIS_to_0.025x0.025_nomask_aave_da_c160323.nc * matrix dimensions rows x cols :      259200  x       48400 * number of non-zero elements:        41163 Successfully made LAIs/SAIs/heights for month            1 Successfully made LAIs/SAIs/heights for month            2 Successfully made LAIs/SAIs/heights for month            3 Successfully made LAIs/SAIs/heights for month            4 Successfully made LAIs/SAIs/heights for month            5 Successfully made LAIs/SAIs/heights for month            6 Successfully made LAIs/SAIs/heights for month            7 Successfully made LAIs/SAIs/heights for month            8 Successfully made LAIs/SAIs/heights for month            9 Successfully made LAIs/SAIs/heights for month           10 Successfully made LAIs/SAIs/heights for month           11 Successfully made LAIs/SAIs/heights for month           12  domain_clean: cleaning       259200------------------------------------------------------------  land model surface data set successfully created for grid of size        48400 creating dynamic land use dataset netcdf error from mkharvest_init rcode =          -49  error =  NetCDF: Variable not found       All listed files are available, including the pftdyn_hist_simyr1850-2005.txt file. Here, all links are adjusted. Since I only use year 2000, I only have this file locally. I did try to replace pftdyn_hist_simyr1850-2005.txt with pftdyn_hist_simyr2000.txt containing only the line for year 2000, but that also did not work. In addition, in this namelist, I have:mksrf_fvegtyp  = '/scratch/leuven/300/vsc30061/cesm/inputdata/lnd/clm2/rawdata/pftlanduse.3minx3min.simyr2000.c110913/mksrf_landuse_rc2000_c110913.nc'But this file does not contain any info on harvest. On the other hand, cesm/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_landuse_rc2000_c090630.nc does contain this information. Unfortunately, using this file directly gives a domain checksame error.Any ideas on how to solve this issue?

 

[oleson]

Unfortunately, I found out that there is no wood harvest data (or even dynamic pft data) at the 3min resolution. The data supplied by George Hurtt is at 0.5 degrees.  So it looks like you would have to use the 0.5 degree landcover change and harvest data and use a mapping file to go from 0.5 to your 0.025x0.025 resolution if you want to do transient landcover change experiments.

 

[matthias_demuzere@ees_kuleuven_be]

Hi Keith, Thanks for this suggestion.I am not really interested in doing transient landcover change experiments. I just want to create the surfdata file for my high-resolution domain using cesm1.2.2 and clm4.0. I did manage to get such a file for clm4.5, but with clm4.0 the routine complains about mkharvest.  I am not sure what the reason could be here?In attach I add the namelists that worked for clm4.5, and the one that now does not work for clm4.0.

 

[oleson]

I can't seem to find your attached namelists...But, If you don't want transient pfts, then just set the fdyn items to blank: mksrf_fdynuse  = ''
 fdyndat        = ''

I tried to create a 0.01X0.01 surface dataset with cesm1_2_2 and the CLM4 mksurfdata_map tool with the high resolution pft data and successfully got a year 2000 dataset, which is what I think you are trying to do.