NCL画图加地图

杨雨轩

                                                我想用兰溪的地图脚本，但是用不成功，请各位帮我看看哪里出了问题
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"

begin

a=addfile("/Users/alain/Desktop/2015southwinterrain/numerical_simulation/POD/wrf_3dvar_input_d02_2015-12-08_19_00_00_meizhou_assfield.nc","r")
b=addfile("/Users/alain/Desktop/2015southwinterrain/numerical_simulation/POD/wrf_3dvar_input_d02_2015-12-08_19_00_00_meizhou_bgfield.nc","r")

type="pdf"
wks=gsn_open_wks(type,"d02_POD")

;=========================================================

mpres = True
lats = (/18.0,34.0/)
lons = (/105.0,123.0/)
loc = wrf_user_ll_to_ij(a,lons,lats,True)
x_start = loc(0,0)-1
x_end = loc(0,1)-1
y_start = loc(1,0)-1
y_end = loc(1,1)-1

mpres@ZoomIn = True
mpres@Xstart = x_start
mpres@Xend = x_end
mpres@Ystart = y_start
mpres@Yend = y_end

mpres@mpGeophysicalLineColor = "black"
mpres@mpGeophysicalLineThicknessF = 2.0
mpres@mpNationalLineColor = "black"
mpres@mpNationalLineThicknessF = 2.0
;mpres@mpUSStateLineThicknessF = "black"
;mpres@mpUSStateLineColor = 2.0
mpres@mpProvincialLineColor = "black"
mpres@mpProvincialLineThicknessF = 2.0

mpres@mpGridAndLimbOn = False
mpres@mpDataBaseVersion = "MediumRes"
mpres@mpDataSetName = "Earth..4"

;mpres@mpFillOn = True
;mpres@mpAreaMaskingOn = True
;mpres@mpMaskAreaSpecifiers = (/"China","Tainwan"/)
;mpres@mpOutLineSpecifiers =(/"China:states","Taiwan"/)

;mpres@mpLandFillColor = "white"
;mpres@mpInlandWaterFillColor = "white"
;mpres@mpOceanFillColor = "white"
;mpres@mpOutlineBoundarySets = "NoBoundaries"
;mpres@mpOutlineBoundarySets = (/"China:states"/)

mpres@tmXBOn = True
mpres@tmXTOn = False
mpres@tmYLOn = True
mpres@tmYROn = False
mpres@tmXBMode = "Explicit"
mpres@tmXBValues = (/105,108,111,114,117,120,123/)
mpres@tmXBLabels = (/"105~S~o~N~E","108~S~o~N~E","111~S~o~N~E","114~S~o~N~E",\
                     "117~S~o~N~E","120~S~o~N~E","123~S~o~N~E"/)

mpres@tmYLMode = "Explicit"
mpres@tmYLValues = (/18,22,26,30,34/)
mpres@tmYLLabels = (/"18~S~o~N~N","22~S~o~N~N","26~S~o~N~N","30~S~o~N~N","34~S~o~N~N"/)

;=========================================================

res               = True
res@NoHeaderFooter= True

pltres            = True

;=========================================================

u1 = wrf_user_getvar(a,"ua",-1)
v1 = wrf_user_getvar(a,"va",-1)

u2 = wrf_user_getvar(b,"ua",-1)
v2 = wrf_user_getvar(b,"va",-1)

p = wrf_user_getvar(a,"pressure",-1)

du = u1-u2
dv = v1-v2

pressure = 850
du1 = wrf_user_intrp3d(du,p,"h",pressure,0.,False)
dv1 = wrf_user_intrp3d(dv,p,"h",pressure,0.,False)

;---------------------------------------------------------


lat = fspan(18.0,34.0,393)
lon = fspan(105.0,123.0,432)
lat@units = "degrees_north"
lon@units = "degrees_east"

du1!1 = "lat"
du1!2 = "lon"
du1&lat = lat
du1&lon = lon

dv1!1 = "lat"
dv1!2 = "lon"
dv1&lat = lat
dv1&lon = lon

;=========================================================

vcres = True
vcres@FieldTitle = ""
vcres@cnInfoLabelOn = False
vcres@gsnAddCyclic = False
vcres@vcRefAnnoOrthogonalPosF = -1.0
vcres@vcRefMagnitudeF = 1.5
vcres@vcRefLengthF = 0.045
vcres@vcGlyphStyle = "CurlyVector"
vcres@vcMinDistanceF = 0.015

;=========================================================

printVarSummary(du1)
du1_zoom = du1(0,y_start:y_end,x_start:x_end)
dv1_zoom = dv1(0,y_start:y_end,x_start:x_end)
vector = wrf_vector(a,wks,du1_zoom,dv1_zoom,vcres)
delete(vcres)
plot=wrf_map_overlays(a,wks,(/vector/),pltres,mpres)

config = True
config@china = True
config@river = True
config@province = True
config@nanhai = False
config@diqu = False
chinamap = add_china_map(wks,plot,config)
draw(plot)
frame(wks)

end

                        file:///page1image256file:///page2image256               
                                               

杨雨轩

画出来是单独的两张图

尽头的尽头

gsnDraw = False
gsnFrame = False

杨雨轩

尽头的尽头 发表于 2016-12-7 10:33
gsnDraw = False
gsnFrame = False

请问加在什么位置

weiguang1233

在res=True下面，画图之前加
res@gsnDraw = False
res@gsnFrame = False

兰溪之水

推荐使用http://bbs.06climate.com/forum.php?mod=viewthread&tid=45252

19921219

您好~请问问题解决了吗？

pop_lxh

感谢楼主分享

pop_lxh

杨雨轩 发表于 2016-12-7 10:19
画出来是单独的两张图

感谢楼主分享感谢楼主分享

pop_lxh

杨雨轩 发表于 2016-12-7 10:52
请问加在什么位置

感谢楼主分享