- 积分
- 2648
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-5-22
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本人使用循环语句读取资料,用panel 一页画多图,加上gsn_add_polyline画框后,单独出图均有框,但是多图中仅最后一张图有框,请大神们指点,哪里出了问题?
系统提示:warning:TransformPostDraw: tfPolyDrawList element 0 is invalid,没太搞懂是什么意思
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/contrib/cd_string.ncl"
begin
nx=144
ny=72
n=15 ;except mme and obs
lat = fspan(-88.75,88.75,72)
lon = fspan(1.25,358.75,144)
lat@units = "degrees_north"
lon@units = "degrees_east"
eqm=new((/n/),string)
eqm=(/"ACCESS1-0","bcc-csm1-1","CNRM-CM5","FGOALS-s2","GFDL-ESM2G","GISS-E2-H","inmcm4","IPSL-CM5A-LR","IPSL-CM5A-MR","MIROC5","MIROC-ESM","MPI-ESM-LR","MPI-ESM-MR","MRI-CGCM3","NorESM1-M"/)
fig=new((/n/),string)
fig=(/"(a)","(b)","(c)","(d)","(e)","(f)","(g)","(h)","(i)","(j)","(k)","(l)","(m)","(n)","(o)"/)
wks=gsn_open_wks("png","panel-c4xp")
gsn_define_colormap(wks,"BlueRedGray")
plot=new(15,graphic)
;dum=new(15,graphic)
a=addfile("cmip_obs_cha.nc","r")
arry = new((/72,144/),double)
;diri = "/mnt/d/ncl/qihou/pr/"
diri = "/mnt/d/ncl/qihou/pr/cli/"
obs=addfile(diri+"OBSp.nc","r")
op=obs->pr(:,::-1,:)
op_a=dim_avg_n(op,0)
op_a!0="lat"
op_a!1="lon"
op_a&lat=lat
op_a&lon=lon
arryb = new((/72,144/),double)
gres = True
gres@gsLineColor = "Black";"Gray"
gres@gsLineThicknessF= 2.0
gres@tfPolyDrawOrder = "postDraw" ; this can be used for polylines, polymarkers, or polygons
lon1 = (/107.25,114.24,114.24,107.25,107.25/)
lat1 = (/24.63,24.63,30.29,30.29,24.63/)
do k =0,n-1
print(eqm(k))
arry=a->pr(k,:,:)
arry!0="lat"
arry!1="lon"
arry&lat=lat
arry&lon=lon
arryb=arry/op_a*100.
arryb!0="lat"
arryb!1="lon"
arryb&lat=lat
arryb&lon=lon
res = True
; res@gsnDraw = False ; don't draw yet 不要翻页!!!!
res@gsnFrame = False ; don't advance frame yet
res@cnInfoLabelOn = False
res@cnFillOn =True ; turn on color
res@gsnSpreadColors = True ; use full color map
res@cnLinesOn = False ; no contour lines
res@lbLabelBarOn =False ;关闭每张图的色标
res@cnLevelSelectionMode="ExplicitLevels"
res@cnLevels=(/-80,-50,-25,0,25,50,80/)
res@cnFillColors=(/2,34,66,127,128,187,220,253/)
res@tmXBLabelFontHeightF = 0.014 ; adjust some font heights
res@tmYLLabelFontHeightF = 0.014
res@tiMainFontHeightF = 0.022
res@txFontHeightF = 0.017
res@mpMinLatF = 22.5
res@mpMaxLatF = 32.5
res@mpMinLonF = 105.
res@mpMaxLonF = 115.
res@mpFillOn = False
res@mpOutlineOn = True
;设置map
res@mpDataSetName = "Earth..4"
res@mpDataBaseVersion = "MediumRes"
res@mpOutlineDrawOrder = "PostDraw"
res@mpGridSpacingF = 10
res@mpProvincialLineColor = "Black"
res@mpNationalLineColor = "Black"
res@mpNationalLineThicknessF = 2.0
res@mpOutlineBoundarySets = "AllBoundaries"
res@gsnLeftString=fig(k)+eqm(k)
res@gsnRightString="%"
res@cnInfoLabelOn = False
plot(k)= gsn_csm_contour_map(wks,arryb(:,:),res)
dum=gsn_add_polyline(wks,plot(k),lon1,lat1,gres)
draw(plot(k))
frame(wks)
delete(arry)
delete(dum)
end do
resP=True
resP@gsnPanelLabelBar=True
;resP@gsnPanelSave = True
resP@gsnFrame = False
gsn_panel(wks,plot,(/4,4/),resP)
end
|
-
-
|