<div dir="ltr"><div class="gmail_default" style="font-size:small">The wmbarb and wmbarbmap functions are from an older interface, and they don't have the same "resource" style of setting options like the other plot types do. Everything has to be set through the wmsetp interface.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Please see this URL, which describes the parameters you can set using wmsetp, that will affect wind barbs:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncarg.ucar.edu//supplements/wmap/index.html#HEADING1-140">http://www.ncarg.ucar.edu//supplements/wmap/index.html#HEADING1-140</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">I think the parameter you want is WBS, and not ARS.</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncarg.ucar.edu//supplements/wmap/index.html#WBS">http://www.ncarg.ucar.edu//supplements/wmap/index.html#WBS</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Nov 5, 2016 at 5:33 AM, <a href="mailto:dyjbean@gmail.com">dyjbean@gmail.com</a> <span dir="ltr"><<a href="mailto:dyjbean@gmail.com" target="_blank">dyjbean@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>
<div><span></span>hi,</div><div> i want add several obervation wind vector data on a map, but i cannot find appropriate function to plot.</div><div> the following is my script:</div><div><br></div><div>++++++++++++++++++++++++++++++<wbr>++++++++</div><div><span style="color:rgb(0,0,0);background-color:rgba(0,0,0,0)">load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"<br>load "$GEODIAG_ROOT/geodiag.<wbr>ncl"<br><br>begin<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr>;;;;;;;;;;;;;;;;;;<br>;; uv data<br>bj_csvs=systemfunc("ls uvdata/<wbr>112920_inside_obs.csv")<br><br>lines=asciiread(bj_csvs,-1,"<wbr>string")<br>nlines=dimsizes(lines) - 1<br><br>delim=","<br>field_names=str_split(lines(0)<wbr>,delim)<br>nfields=dimsizes(field_names)<br><br>;<br>fields=new((/nfields,nlines/),<wbr>string)<br>do nf=0,nfields-1<br> fields(nf,:)=str_get_field(<wbr>lines(1:),nf+1,delim)<br>end do<br><br>print(fields(1,:) + " "+<wbr>fields(2,:) + " "+fields(3,:<wbr>)+" "+fields(4,:))<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<wbr>;;;;;;;;;;;;;;;;<br>;; pm2.5 data<br><br>ncol = 3 <wbr> ; number of columns <wbr>is 3£»ÄãµÄÁÐÊý<br> stationdata= readAsciiTable(<wbr>"pm2.5_txt/15112920.txt", <wbr>ncol, "float", 1)<br> print(dimsizes(stationdata))<br> nps =dimsizes(stationdata)<br> npts = nps(0) ; <wbr>Number of points.<br><br> lat = stationdata(:,1) ;<wbr> latitude values<br> lon = stationdata(:,0) ;<wbr> longitude values<br> pm25= stationdata(:,2) ; <wbr>station numbers to appear on <wbr>the map£»Õâ¸öµØ·½Äã¿ÉÒÔдÎÛȾ<wbr>Êý¾Ý<br><br> printMinMax(pm25,0)<br>;-------interpolate-----------<wbr>------------------------------<wbr>------------<br> olon =new(80,"float")<br> olat =new(80,"float")<br> data1=new((/80,80/),"float")<br> do j=0,79<br> olon(j)=113+j*0.1<br> end do<br> do k=0,79<br> olat(k)=36+k*0.1<br> end do<br> <br><br><br>;---------set interpolate res-<wbr>------------------------------<wbr>-------------<br> olon!0 ="lon"<br> olon@long_name="lontitude"<br> olon@units ="degrees-<wbr>east"<br> olon&lon =olon<br><br> olat!0 ="lat"<br> olat@long_name="latitude"<br> olat@units ="degrees-<wbr>north"<br> olat&lat =olat<br><br>;--------Call interpolation <wbr>function----------------------<wbr>----------------<br> ;pm25@_FillValue=999999.<wbr>000000<br> rscan=(/10,5,3,2,1,0.5,0.2,<wbr>0.1/)<br> data1=obj_anal_ic_<wbr>deprecated(lon, lat, pm25, <wbr>olon, olat, rscan, False);<wbr>creanm<br> printMinMax(data1 ,0)<br> data1!0="lat"<br> data1!1="lon"<br> data1&lat=olat<br> data1&lon=olon<br> data1@units=""<br> data1@long_name="pm2.5"<br> ;data1@_FillValue = pm25@_<wbr>FillValue<br> printMinMax(olat,0)<br> printMinMax(olon, 0)<br><br>;; plot<br> type="png"<br> wks = gsn_open_wks(type,"<a href="tel:2015112920" value="+12015112920" target="_blank">201<wbr>5112920</a>") ; Open a <wbr>workstation and<br> colors=(/"White","Black","(/<wbr>.027,.729,.145/)","(/.431,.<wbr>824,.035/)",\<br> "(/.800,.937,.012/)"<wbr>,"(/1.00,.969,0.00/)","(/1.00,<wbr>.871,0.00/)",\<br> "(/1.00,.710,0.00/)"<wbr>,"(/1.00,.549,0.00/)","(/0.00,<wbr>.757,1.00/)",\<br> "(/0.00,0.455,1.00/)<wbr>","(/0.00,.122,1.00/)","(/1.<wbr>00,.275,0.00/)",\<br> "(/1.00,0.059,0.00/)<wbr>"/) <br> gsn_define_colormap(wks,<wbr>colors)<br><br> res = True<br> res@cnFillOn = True<br> res@cnLinesOn = False<br> res@cnLineLabelsOn = False<br> res@gsnDraw <wbr> = False<br> res@gsnFrame <wbr> = False<br> res@gsnMaximize <wbr> = True<br> res@gsnAddCyclic = False<br><br> <br> res@mpLimitMode = "LatLon"<br> res@mpMinLatF = 39.40<br> res@mpMaxLatF = 41.00<br> res@mpMinLonF = 115.35<br> res@mpMaxLonF = 117.55<br><br> res@isShowProvince = False<br> res@isShowBeijingCounty = <wbr>True<br> res@beijingCountyColor = "<wbr>black"<br> res@beijingCountyThickness =<wbr> 1.0<br> res@isShowBeijingRing2 = <wbr>True<br> res@beijingRing2Color = "<wbr>black"<br> res@beijingRing2Thickness = <wbr>1.50<br> res@isShowBeijingRing3 = <wbr>True<br> res@beijingRing3Color = "<wbr>black"<br> res@beijingRing3Thickness = <wbr>1.50<br> res@isShowBeijingRing4 = <wbr>True<br> res@beijingRing4Color = "<wbr>black"<br> res@beijingRing4Thickness = <wbr>1.5<br> res@isShowBeijingRing5 = <wbr>True<br> res@beijingRing5Color = "<wbr>black"<br> res@beijingRing5Thickness = <wbr>1.5<br> res@isShowBeijingRing6 = <wbr>True<br> res@beijingRing6Color = "<wbr>black"<br> res@beijingRing6Thickness = <wbr>1.5<br> setup_china_map(res)<br><br> plot = gsn_csm_contour_map(<wbr>wks,data1,res)<br><br> attach_china_map(wks, plot)<br><br> draw(plot)<br><br> ; wmsetp("vch", 0.05)<br> ; wmsetp("vcc", 1)<br> ; wmsetp("vcw", 3.)<br> ; wmsetp("wbs", 0.05)<br> ; wmsetp("col", 1)<br> ; wmsetp("wbs", 0.06)<br> ; wmsetp("ezf",2)<br> ; wmsetp("wbs",0.12)<br> </span><span style="background-color:rgba(0,0,0,0)"><font color="#ff0000"><b>wmsetp("ars",0.05)<br> wmsetp("arl",1.5)</b></font><br> ; wmvectmap(wks, tofloat(<wbr>fields(2,:)), tofloat(fields(<wbr>1,:)), tofloat(fields(3,:)), <wbr>tofloat(fields(4,:)))<br> ; wmvlbl(wks, 1.0, 0.0)<br> <b><font color="#ff0000">wmbarbmap(wks, tofloat(<wbr>fields(2,:)), tofloat(fields(<wbr>1,:)), tofloat(fields(3,:)), <wbr>tofloat(fields(4,:)))</font></b><br><br> frame(wks)<br><br><br>end<br><br></span></div><div>++++++++++++++++++++++++++++++<wbr>++++++++ </div><div>the red and bold font is for wind barb.</div><div><br></div><div>the figure is as follows:</div><div> <img src="cid:_Foxmail.1@aab1cd9b-dc9f-64d1-7714-3db79824f84c" border="0" style="font-size:10.5pt;line-height:1.5;background-color:window"></div><div>the wind barb are so small on t<span style="background-color:window;font-size:10.5pt;line-height:1.5">his plot, but those sizes are so tiny , which function can be used to enlarge the wind rod for appropriate size.</span></div><div><br></div><div>thanks </div>
<div><br></div><hr style="width:210px;height:1px" color="#b5c4df" size="1" align="left"><span class="HOEnZb"><font color="#888888">
<div><span><div style="MARGIN:10px;FONT-FAMILY:verdana;FONT-SIZE:10pt"><div><a href="mailto:dyjbean@gmail.com" target="_blank">dyjbean@gmail.com</a></div></div></span></div>
</font></span></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>