Could you share the LUH2 data at 0.05x0.05 resolution?

What version of the code are you using?
CESM2.1.3

Have you made any changes to files in the source tree?
NO

Describe your problem or question:

I used LUH2 data at 0.25×0.25° resolution for creating regional surface files, but I found 0.05°×0.05° resolution LUH2, soilcolor, and LAI data under the directory /cesm/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.05x0.05.LUH2.histsimyr2005.c190116. I would like to update the mksrf_fvegtyp, mksrf_fhrvtyp, mksrf_fsoicol, and mksrf_flai parameters in the surface namelist file and regenerate the regional surface files.

I ran the command ./mksurfdata_map < surfdata_hubei0.05degrees_hist_78pfts_CMIP6_simyr2000_c241012.namelist and encountered the following error.

(cesm_py27) [fengjinping@swarm01 mksurfdata_map]$ ./mksurfdata_map < surfdata_hubei0.05degrees_hist_78pfts_CMIP6_simyr2000_c241012.namelist
Attempting to initialize control settings .....
Attempting to create surface boundary data .....
------------------------------------------------------------
mksrf_fgrid =
../mkmapdata/map_0.25x0.25_MODIS_to_hubei0.05degrees_nomask_aave_da_c241011.nc



mksrf_gridtype = global
Output ALL data in file as 64-bit
Set wetland to 0% over land
In mkpftMod::mkpftInit()...
calling domain_read
finished domain_read
fsurdat is 2d lat/lon grid
nlon= 180 nlat= 120
(OPNFIL): Successfully opened file
surfdata_hubei0.05degrees_hist_78pfts_CMIP6_simyr2000_c241012.log on unit=
99
mksrf_gridtype = global

In mkpftMod::mkpft()...
Attempting to make PFTs .....
Creating surface datasets with extra types for crops; total pfts = 78
domain_read_dims_2d read lon and lat dims from lon/lat
domain_read initialized domain
domain_read read LANDFRAC
domain_read read LANDMASK
Open PFT file:
/home/fengjinping/scratch/cesmlab/lnd/clm2/rawdata/pftcftdynharv.0.05x0.05.LUH2
.histsimyr2005.c190116/mksrf_landuse_clm52deg005_histLUH2_2005.c190119.nc
(gridmap_map_read) reading mapping matrix data...
(gridmap_map_read) * file name : ../mkmapdata/map_0.25x0.25_MODIS_to_hubei0.05degrees_nomask_aave_da_c241011.nc
* matrix dimensions rows x cols : 1036800 x 21600
* number of non-zero elements: 25920
domain_checksame ERROR: input domain size and gridmap source size are not the s
ame size
domain size = 25920000
map src size= 1036800
abort:

I believe the error arises because the mapping files are all at 0.25×0.25° or 0.5×0.5° resolution, which aligns with the 0.25×0.25° LUH2 data. To use 0.05×0.05° LUH2 data, the mapping files must also be at 0.05×0.05° resolution.

If my understanding is correct, could you please share the following:

1. The globalclm5histLUH2namelist.txt file for generating 0.05×0.05° LUH2 data, along with the corresponding GLM timeseries files and PFT/CFT description files? Since I am using CLM5.0.30 (not CLM52), I need to generate surface and landuse files through CLM5.0.30.
2. The 0.05×0.05° resolution mapping files (NetCDF format), such as:
   • mksrf_fgrid
   • map_fpft
   • map_fglacier,
   • map_fglacierregion
   • map_fsoicol,
   • map_furban,
   • map_fmax
   • map_forganic,
   • map_flai,
   • map_fharvest
   • map_flakwat,
   • map_fwetlnd,
   • map_fvocef
   • map_fsoitex,
   • map_furbtopo,
   • map_fgdp
   • map_fpeat,
   • map_fsoildepth,
   • map_fabm,
   • mksrf_fsoitex
     I would greatly appreciate your assistance!

The mksurfdata_map tool that you're using has the option --hirespft as documented elsewhere, including in Forums' threads (e.g. high resolution (0.1 degree) CLM5-BGC-crop regional simulation). This should use the data that you intended.

slevis said:

The mksurfdata_map tool that you're using has the option --hirespft as documented elsewhere, including in Forums' threads (e.g. high resolution (0.1 degree) CLM5-BGC-crop regional simulation). This should use the data that you intended.

Click to expand...

@slevis Thank you for your advice. I successfully created the Surface file.

I would like to proceed with creating the landuse_timeseries file. The four TXT files—landuse_timeseries_hist_78pfts_simyr1850-2015.txt, landuse_timeseries_SSP1-2.6_78pfts_CMIP6_simyr1850-2100.txt, landuse_timeseries_SSP3-7.0_78pfts_CMIP6_simyr1850-2100.txt, and landuse_timeseries_SSP5-8.5_78pfts_CMIP6_simyr1850-2100.txt—all use LUH2 data at 0.25×0.25 resolution. However, in the directory /cesm/inputdata/lnd/clm2/rawdata/pftcftdynharv.0.05x0.05.LUH2.histsimyr2005.c190116, the LUH2 data is at 0.05×0.05 resolution. Could you share the original data used to generate these files or specify their source? I want to reproduce the landuse_timeseries files with LUH2 data at 0.05×0.05 resolution for both the historical (hist) and SSP scenarios.

I contacted @Lawrence via email last week to inquire about how to generate mksrf_landuse_clm52deg005_histLUH2_2005.c190119.nc from LUH2 at a 0.05×0.05 resolution. However, Lawrence has been busy and has not yet responded. Since responses on the CESM forum are relatively fast, I decided to post this inquiry here.


Additionally, I have been studying the clm5landusedatatool for a week, but I still do not fully understand the process of generating LUH2 at a 0.05×0.05 resolution. Therefore, I hope experts here can provide me with some guidance.

Peter Lawrence is the expert in this case.