說明

利用 ncl util/plotgrids_new.ncl 可以得到namelist.wps空間配置的圖片，如下圖所示。但會存在兩個問題，我們一般做wrf模擬最內(nèi)層嵌套往往遠(yuǎn)遠(yuǎn)小于省界，因此需要疊加市界線；并且ncl自帶的省界線并不準(zhǔn)確，會導(dǎo)致無法將圖直接用在論文上。

示例圖片

代碼

直接再plotgrids_new.ncl最后面的frame(wks)前面添加。
方法一：使用自己的shp，可隨意疊加省界、市界、縣界、河流、公路等等，過程更隨意，結(jié)果更準(zhǔn)確。

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為原始代碼中設(shè)定的。
只能使用polyline的shp數(shù)據(jù)，如果是polygon的shp數(shù)據(jù)會自動顏色填充，如果設(shè)置填充顏色為不透明，就看不到d01 d02這些方框了，如果設(shè)置為透明就會看不到省市界線。
對于沒有x11的windows上，直接出圖或連接的服務(wù)器出圖，均只能設(shè)置為pdf，png出圖則下半部分會被黑色覆蓋（未知原因），x11在windows上根本就彈不出框，不予考慮。
方法二：直接使用ncl的地圖，只能添加到省界，無法添加市界、縣界等等，并且存在臺灣和藏南等問題，無法在論文中使用。

代碼圖

結(jié)果示例：

ncl省界

shp省界

shp省界，紅圈內(nèi)沒有是因為shp數(shù)據(jù)只有國內(nèi)的。缺點是，東海岸附近的ncl自帶的國界線無法去掉。