说明

利用 ncl util/plotgrids_new.ncl 可以得到namelist.wps空间配置的图片，如下图所示。但会存在两个问题，我们一般做wrf模拟最内层嵌套往往远远小于省界，因此需要叠加市界线；并且ncl自带的省界线并不准确，会导致无法将图直接用在论文上。

示例图片

代码

直接再plotgrids_new.ncl最后面的frame(wks)前面添加。
方法一：使用自己的shp，可随意叠加省界、市界、县界、河流、公路等等，过程更随意，结果更准确。

myres = True
myres@gsnMaskLambertConformal = True
myres@mpProjection = "LambertConformal"
myres@gsnMaximize = True
myres@pmTickMarkDisplayMode = "Always"
polyres = True
polyres@gsLineColor = "black"
shpname = "/data/DATA/shp/bou2_4l.shp"
poly = gsn_add_shapefile_polylines(wks,mp,shpname,polyres)
draw(mp)

注：
mp为原始代码中设定的。
只能使用polyline的shp数据，如果是polygon的shp数据会自动颜色填充，如果设置填充颜色为不透明，就看不到d01 d02这些方框了，如果设置为透明就会看不到省市界线。
对于没有x11的windows上，直接出图或连接的服务器出图，均只能设置为pdf，png出图则下半部分会被黑色覆盖（未知原因），x11在windows上根本就弹不出框，不予考虑。
方法二：直接使用ncl的地图，只能添加到省界，无法添加市界、县界等等，并且存在台湾和藏南等问题，无法在论文中使用。

代码图

结果示例：

ncl省界

shp省界

shp省界，红圈内没有是因为shp数据只有国内的。缺点是，东海岸附近的ncl自带的国界线无法去掉。