[ncl-talk] some problem during vector panel plot
grace
313695096 at qq.com
Tue Oct 13 01:41:00 MDT 2015
Hi all:
I am trying to make a panel plot,but the Xlabelvalue of the third plote not appear?
and how can I turn off the the reference vector box of the first and second plot and make the reference vector box of the third plot samller?
The text (TEST,OBS) that I add on the plot did not make any difference.
How can I slove those problem?
Or do you guys have some advice?
Any information will be appreciated.
This is my plot:
This is my scirpt:
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/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
;
; The WRF ARW input file.
; This needs to have a ".nc" appended, so just do it.
a = addfile("/public/home/huanglei/wrfdata/wrfout_d03_2015-05-01_12:00:00.nc","r")
ij = wrf_user_ll_to_ij(a, 108.966667,34.433333,True)
print("lon location is: " +ij(0))
print("lat location is: " + ij(1))
; First get the variables we will need
u10_2=new((/31,24/),float,"No_FillValue")
v10_2=new((/31,24/),float,"No_FillValue")
filename1=systemfunc("ls /public/home/huanglei/wrfdata/wrfout_d03_2015-05*.nc")
fin1=addfiles(filename1,"r")
do i=0,30
u10_1=fin1[i]->U10(0:23,ij(1)-1, ij(0)-1)
v10_1=fin1[i]->V10(0:23,ij(1)-1, ij(0)-1)
u10_2 (i,:)= u10_1
v10_2 (i,:)= v10_1
end do
u10=ndtooned(u10_2)
v10=ndtooned(v10_2)
u10_d=u10(0:724)
v10_d=v10(0:724)
;;;;;;;;;;;;;;;get obs data;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
ascii_filename = "/public/home/huanglei/obs.txt"
seismic = asciiread(ascii_filename,(/725,6/),"float")
seismic at _FillValue = -999
wdir = seismic(:,4)
wspd = seismic(:,5)
rad = 4.0*atan(1.0)/180.
u_o = wspd*sin(rad*wdir)
v_o = wspd*cos(rad*wdir)
;;;;;;;;;;;;;;;;;;;;;;;;;data for plot1;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
u10_new1=new((/25,220/),float)
u10_new1=0
u10_new1(8,:)=u10_d(0:219)
u10_new1(15,:)=u_o(0:219)
v10_new1=new((/25,220/),float)
v10_new1=0
v10_new1(8,:)=v10_d(0:219)
v10_new1(15,:)=v_o(0:219)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;data for plot2;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
u10_new2=new((/25,240/),float)
u10_new2=0
u10_new2(8,:)=u10_d(220:459)
u10_new2(15,:)=u_o(220:459)
v10_new2=new((/25,240/),float)
v10_new2=0
v10_new2(8,:)=v10_d(220:459)
v10_new2(15,:)=v_o(220:459)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;data for plot3;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
u10_new3=new((/25,264/),float)
u10_new3=0
u10_new3(8,:)=u10_d(461:724)
u10_new3(15,:)=u_o(461:724)
v10_new3=new((/25,264/),float)
v10_new3=0
v10_new3(8,:)=v10_d(461:724)
v10_new3(15,:)=v_o(461:724)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
wks = gsn_open_wks("pdf","gsn_vector_panel") ; open a Graphic wkStation
plot = new(3,graphic)
res = True ; plot mods desired
res at gsnFrame = False
res at gsnDraw = False ; don't draw
; res at tiMainString = "Basic Vector Plot"
res at vcGlyphStyle = "CurlyVector" ; turn on curley vectors
res at gsnMaximize = True
res at vpHeightF = 0.15
res at vpWidthF = 0.8
;res at vcRefMagnitudeF = 5.0 ; add a reference vector
; res at vcRefLengthF = 0.045 ; what the ref length is
; res at vcRefMagnitudeF = 0.15
; res at vcRefLengthF = 0.001
res at vcRefAnnoArrowLineColor = "black" ; change ref vector color
res at pmLegendWidthF = 0.10 ; Change width and
res at pmLegendHeightF = 0.10 ; height of legend.
; res at lgLabelFontHeightF = 0.02 ; change font height
res at tmYLOn = False ; Don't draw labels for left
res at tmYLLabelsOn = False ; Y axes
res at tmYROn = False ; Don't draw labels for right
res at tmYRLabelsOn = False ; Y axes
;;;;;;;;;;;;;;;;;;;;add text TEST / OBS;;;;;;;
res_text = True
res_text at txFontHeightF = 0.01
; res_text at txJust = "CenterLeft"
m1 = "TEST"
m2 = "OBS"
res0=res
res0 at tmXBMode ="Explicit"
res0 at tmXBValues =(/5,29,53,77,101,125,149,173,197/)
res0 at tmXBLabels = (/"5/2","5/3"," 5/4","5/5", "5/6","5/7","5/8","5/9","5/10"/)
res0 at tmXBLabelFontHeightF = 0.01
res0 at vcRefAnnoOrthogonalPosF = -0.33 ; move ref vector up
res0 at vcRefMagnitudeF = 0.5 ; define vector ref mag
res0 at vcRefLengthF = 0.025 ; define length of vec ref
res0 at vcRefAnnoParallelPosF = 0.999
res0 at vpHeightF = 0.15
res0 at vpWidthF = 0.8
plot(0) = gsn_vector(wks,u10_new1,v10_new1,res0)
gsn_text(wks,plot(0),m1,190,10,res_text)
gsn_text(wks,plot(0),m2,190,35,res_text)
res1=res
res1 at tmXBMode ="Explicit"
; res1 at tmXBValues =(/221,245,269,293,317,341,365,389,413,437/)
res1 at tmXBValues =(/1,25,49,73,97,121,145,169,193,217/)
res1 at tmXBLabels = (/"5/11","5/12","5/13","5/14","5/15","5/16","5/17","5/18","5/19","5/20"/)
res1 at tmXBLabelFontHeightF = 0.01
res1 at vcRefAnnoOrthogonalPosF = -0.33 ; move ref vector up
res1 at vcRefMagnitudeF = 0.5 ; define vector ref mag
res1 at vcRefLengthF = 0.025 ; define length of vec ref
res1 at vcRefAnnoParallelPosF = 0.999
res1 at vpHeightF = 0.15
res1 at vpWidthF = 0.8
plot(1) = gsn_vector(wks,u10_new2,v10_new2,res1)
gsn_text(wks,plot(1),m1,200,10,res_text)
gsn_text(wks,plot(1),m2,200,35,res_text)
res2=res
res2 at tmXBMode ="Explicit"
; res2 at tmXBValues =(/461,485,509,533,557,581,605,629,653,677,701/)
res2 at tmXBValues =(/1,25,49,73,97,121,145,169,193,217,241/)
res2 at tmXBLabels = (/"5/21"," 5/22"," 5/23","5/24","5/25","5/27","5/28","5/29","5/30","5/31"/)
res2 at tmXBLabelFontHeightF = 0.01
res2 at vcRefAnnoOrthogonalPosF = -0.33 ; move ref vector up
res2 at vcRefMagnitudeF = 0.5 ; define vector ref mag
res2 at vcRefLengthF = 0.025 ; define length of vec ref
res2 at vcRefAnnoParallelPosF = 0.999
res2 at vpHeightF = 0.15
res2 at vpWidthF = 0.8
plot(2) = gsn_vector(wks,u10_new3,v10_new3,res2)
gsn_text(wks,plot(2),m1,220,10,res_text)
gsn_text(wks,plot(2),m2,220,35,res_text)
;;;;;;;;;;;;;;;;add obs test mark;;;;;;;;;;;;;;
resP = True ; modify the panel plot
gsn_panel(wks,plot,(/3,1/),resP) ; now draw as one plot
;frame(wks)
end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151013/b3b6133c/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/octet-stream
Size: 283643 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151013/b3b6133c/attachment-0001.obj
More information about the ncl-talk
mailing list