Generating ESMFMesh files for WRF-CTSM

[gunnartl]

Following the documentation on using CTSM with WRF, in section 3.4.1.7 Set ctsm runtime options, we are told to modify lilac_in in the following way:
atm_mesh_filename = '/glade/scratch/negins/wrf_ctsm_files/wrf2ctsm_land_conus_ESMFMesh_c20201110.nc'
lnd_mesh_filename = '/glade/scratch/negins/wrf_ctsm_files/wrf2ctsm_land_conus_ESMFMesh_c20201110.nc'

Could you please let us know how these files are generated?

[negin513]

@gunnartl Would you like access to these files to run the CONUS domain or are you interested in creating the mesh file for another WRF domain?
I've put the mesh file (netcdf) for the CONUS domain in dropbox so you can download it:
https://www.dropbox.com/s/hpwmk96vhjla6a3/wrf2ctsm_land_conus_ESMFMesh_c20201110.nc?dl=0

[gunnartl]

I am interested in creating the mesh file for Scandinavia. But i can use these files for some test runs. Thank you!

[mvertens]

@gunartl - if you can generate a file in SCRIP format - then there is an ESMF offline tool that will convert it to a mesh.
http://earthsystemmodeling.org/docs/nightly/develop/ESMF_refdoc/node3.html#SECTION03040000000000000000

A simple script to help you do this is below:

#!/bin/bash
#PBS -A P93300606
#PBS -N mesh
#PBS -j oe
#PBS -q share
#PBS -l walltime=03:00:00
#PBS -l select=1:mpiprocs=36
#PBS -o log.out

export MPI_USE_ARRAY=false

# load modules
module purge
module load ncarenv/1.3
module load intel/19.0.2
module load mpt/2.19
module load netcdf-mpi/4.7.1
module load pnetcdf/1.11.0
module load ncarcompilers/0.5.0
module use /glade/work/turuncu/PROGS/modulefiles/esmfpkgs/intel/19.0.2
module load esmf-8.0.0-ncdfio-mpt-O
module load nco

file_i=<full pathname of input SCRIP file>
file_o=<full pathname of output mesh file>
mpiexec_mpt `hostname` -np 8 ESMF_Scrip2Unstruct $file_i $file_o  0 ESMF

[irismuz]

Thank you @mvertens. Following your suggestion, we were able to generate the required mesh files.

In the WRF-CTSM running phase, we receive an error that indicates the inconsistency in the land-sea mask (two CPUs illustrate that this issue refers to the two tiles). Kindly see the error message below.

However, we have found no mismatch in the land-sea mask between geo_em.d01.nc and surfdata_1x1_hist_16pfts_Irrig_CMIP6_simyr2000_c210611.nc.
@negin513 could you please provide us with your input?