Unfortunately, the mksurfdata_map tools haven't and won't be ported to our new machine Derecho. So we don't have a way to generate mapping and surface dataset files using that tool.
However, we have a new tool which will be released soon, mksurfdata_esmf. There are instructions on how to create and use datasets using this tool here:
Context This started as issue #1773 from which I linked a google doc (now out of date and obsolete) with step-by-step instructions intended initially for a case that @XiulinGao needed to run. To th...
github.com
Basically, you would start by checking out the latest version of ctsm5.2.mksurdata.
Then, using your SCRIP file, create lnd_mesh.nc and mask_mesh.nc files.
Then create your surface dataset.
I'm not going to get into the habit of creating datasets for other people with this tool. However, I'm still learning about this tool so I went through the steps myself to create a 0.01x0.01 regional dataset:
git clone
GitHub - ESCOMP/CTSM: Community Terrestrial Systems Model (includes the Community Land Model of CESM)
cd CTSM
git checkout ctsm5.2.mksurfdata
./manage_externals/checkout_externals
module unload python
module load conda
./py_env_create
conda activate ctsm_py
cd tools/mksurfdata_esmf
/glade/u/apps/derecho/23.09/spack/opt/spack/esmf/8.6.0/cray-mpich/8.1.27/oneapi/2023.2.1/7haa/bin/ESMF_Scrip2Unstruct SCRIPgrid_0.01x0.01_shanxi_nomask_c240315.nc lnd_mesh.nc 0
ncap2 -s 'elementMask(:)=0' lnd_mesh.nc mask_mesh.nc
ncdump -v grid_dims SCRIPgrid_0.01x0.01_shanxi_nomask_c240315.nc (to get nx and ny)
./gen_mksurfdata_namelist --start-year
2005 --end-year
2005 --model-mesh-nx 606 --model-mesh-ny 808 --nocrop --model-mesh lnd_mesh.nc --res 0.01x0.01_shanxi
--hires_pft --hires_soitex
./gen_mksurfdata_jobscript_single --number-of-nodes 24 --tasks-per-node 12 --namelist-file surfdata_0.01x0.01_shanxi_hist_2005_16pfts_c240316.namelist
qsub ./mksurfdata_jobscript_single
I've put the resulting surface dataset and log file on our ftp site:
ftp://ftp.cgd.ucar.edu/pub/oleson/0.01x0.01
But you should try this yourself and see if you can get it to work on your machine.
Note that I selected --hires_pft --hires_soitex, which is the high resolution pft and soils data. These are only available for year 2005 which is why I set --start-year
2005 --end-year
2005. Regardless, your grid is extremely high resolution so a lot of the surface dataset fields are going to look blocky.
Note also that there are new fields on the dataset such that it will only be compatible with newer code. For example, you could create/run your cases using the ctsm5.2.mksurdata checkout itself.
