<div dir="ltr">Hello,<div><br></div><div>I occurred a problem when I plot contour.</div><div><br></div><div>Even though I turn off contour line label bar the graph also shows it.</div><div>" opts_rh = opts_xy</div><div> opts_rh@cnLineLabelsOn = False </div><div> opts_rh@cnFillOn = True</div><div>"</div><div><br></div><div>Is there anyone know why?</div><div><br></div><div>Thank you,</div><div>Tao</div><div><br></div><div><br></div><div>Here is my codes:</div><div>###########################################################################################</div><div><div>begin</div><div>;</div><div>;---Read wrf file </div><div> ;a = addfile("../../wrfout_d01_2015-09-08_12:00:00","r") ; small domain ouput 5</div><div> a = addfile("../wrfout_d01_2015-09-08_12:00:00_big2","r") ; big domain ouput 10</div><div><br></div><div>; read need data</div><div> xlat = wrf_user_getvar(a, "XLAT",0)</div><div> xlon = wrf_user_getvar(a, "XLONG",0)</div><div> ter = wrf_user_getvar(a, "HGT",0)</div><div><br></div><div>;--Set up workstation</div><div> type = "x11"</div><div>; type = "pdf"</div><div>; type = "ps"</div><div>; type = "ncgm"</div><div> type@wkWidth<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= 3000 ; Increase size for a slightly</div><div> type@wkHeight<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= 3000</div><div><br></div><div>; Output folder</div><div> folder="./clean_vertical_wind"</div><div> mkdir = systemfunc ("mkdir -p " + folder)</div><div> wks = gsn_open_wks(type,folder+"/plt_CrossSection_3")</div><div><br></div><div><br></div><div>;---Set some basic resources</div><div> res = True</div><div> res@gsnMaximize =True</div><div> res@gsnDraw = False</div><div> res@gsnFrame = False</div><div><br></div><div><br></div><div> </div><div>;------------------------------------------------------------------------------ </div><div>;---Set bounding box</div><div>; C->C' WE</div><div> lats = (/ 36.59845, 36.59845 /)</div><div> lons = (/ 138.50, 141.0 /)</div><div><br></div><div>; X-axis label</div><div> xoption=0 </div><div><br></div><div>;0 means lon X-axis, 1 means lat X-axis</div><div> if (xoption .eq. 0) then</div><div> X_desc = "longitude" </div><div> xvalue=xlon</div><div> px = 0.2</div><div> end if</div><div><br></div><div> if (xoption .eq. 1) then</div><div> X_desc = "latitude"</div><div> xvalue = xlat</div><div> px = 0.5</div><div> end if</div><div><br></div><div> </div><div>;---Set zoomin indices</div><div>; loc(0,;) is west-east (x) ; loc(1,:) is south-north (y)</div><div>; subtract one since we want to use it as an index in NCL</div><div> loc<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= wrf_user_ll_to_ij(a, lons, lats, True)</div><div> loc<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= loc - 1</div><div> x_start<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= loc(0,0)</div><div> x_end<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= loc(0,1)</div><div> y_start<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= loc(1,0)</div><div> y_end<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= loc(1,1)</div><div> </div><div> </div><div>;---Read time</div><div><br></div><div> times = wrf_user_getvar(a,"times",-1) ; get times in the file</div><div> ntimes = dimsizes(times) ; number of times in the file</div><div> FirstTime = True</div><div><br></div><div> mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for the file</div><div> nd = dimsizes(mdims)</div><div> </div><div> </div><div>;--------------------------------------------------------------------------</div><div>; Loop of time to do plot</div><div>;---------------------------------------------------------------</div><div><br></div><div> ;do it = 356,356;ntimes-1,18 ; for output5 when need 0909_1740 (UTC)</div><div> do it = 178,178;ntimes-1,18 ; for output10 when need 0909_1740 (UTC)</div><div> </div><div><br></div><div> print("Working on time: " + times(it) )</div><div><br></div><div> rh = wrf_user_getvar(a,"rh",it) ; relative humidity</div><div> z = wrf_user_getvar(a, "z",it) ; grid point height</div><div> u = wrf_user_getvar(a,"ua",it) ; u in m/s</div><div> w = wrf_user_getvar(a,"wa",it) ; w in m/s</div><div><br></div><div><br></div><div> if ( FirstTime ) then ; get height info for labels</div><div> zmin = 0.</div><div> zmax = 12;max(z)/1000.</div><div> nz = floattoint(zmax + 1)</div><div> FirstTime = False</div><div> end if</div><div><br></div><div>;---------------------------------------------------------------</div><div><br></div><div>; Plot a cross session that run from point A to point B</div><div><br></div><div> plane = new(4,float)</div><div> plane = (/ x_start,x_end, y_start,y_end /) ; start x;y & end x;y point </div><div> opts = True ; start and end points specified</div><div><br></div><div> rh_plane = wrf_user_intrp3d(rh,z,"v",plane,0.,opts)</div><div> u_plane = wrf_user_intrp3d(u,z,"v",plane,0.,opts)</div><div> w_plane = wrf_user_intrp3d(w,z,"v",plane,0.,opts)</div><div> X_plane = wrf_user_intrp2d(xvalue,plane,0,opts)</div><div><br></div><div><br></div><div> dim = dimsizes(rh_plane) ; Find the data span - for use in labels</div><div> zspan = dim(0)</div><div><br></div><div> ; Options for XY Plots</div><div> ; X-axis lables</div><div> dimsX = dimsizes(X_plane)</div><div> xmin = X_plane(0)</div><div> xmax = X_plane(dimsX(0)-1)</div><div> xspan = dimsX(0)-1</div><div> nx = floattoint( (xmax-xmin)/px + 1)</div><div> print ("dimsX= "+dimsX+" xmin= "+xmin+" xmax= "+xmax+" xspan= "+xspan+ " nx="+nx)</div><div><br></div><div><br></div><div> ;---------------------------------------------------------------</div><div> </div><div> ; Options for XY Plots</div><div> opts_xy = res</div><div> opts_xy@tiXAxisString = X_desc</div><div> opts_xy@tiYAxisString = "Height (km)"</div><div> opts_xy@tmXTOn = False</div><div> opts_xy@tmYROn = False</div><div> opts_xy@tmXBMode = "Explicit"</div><div> if (nx .lt. 0) then</div><div> nx = -nx</div><div> end if</div><div> opts_xy@tmXBValues = fspan(0,xspan,nx) ; Create tick marks</div><div> opts_xy@tmXBLabels = sprintf("%.1f",fspan(xmin,xmax,nx)) ; Create labels</div><div> opts_xy@tmXBLabelFontHeightF = 0.015</div><div> opts_xy@tmYLMode = "Explicit"</div><div> opts_xy@tmYLValues = fspan(0,zspan,nz) ; Create tick marks</div><div> opts_xy@tmYLLabels = sprintf("%.1f",fspan(zmin,zmax,nz)) ; Create labels</div><div> opts_xy@tiXAxisFontHeightF = 0.020</div><div> opts_xy@tiYAxisFontHeightF = 0.020</div><div> opts_xy@tmXBMajorLengthF = 0.02</div><div> opts_xy@tmYLMajorLengthF = 0.02</div><div> opts_xy@tmYLLabelFontHeightF = 0.015</div><div> </div><div><br></div><div> ; Plotting options for RH</div><div> opts_rh = opts_xy</div><div> opts_rh@cnLineLabelsOn = False </div><div> opts_rh@cnFillOn = True</div><div> opts_rh@lbOrientation = "vertical"</div><div> opts_rh@ContourParameters = (/ 10., 90., 10. /)</div><div> opts_rh@cnFillColors = (/"White","White","White", \</div><div> "White","Chartreuse","Green", \</div><div> "Green3","Green4", \</div><div> "ForestGreen","PaleGreen4"/) </div><div> ; Set contour title </div><div> mainstring<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= times(it) + " (UTC)"</div><div> leftstring<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= "Relative humidity %"</div><div> rightstring<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= "Wind (m/s)"</div><div> size<span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= 0.02</div><div><br></div><div> opts_rh@tiMainString <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= mainstring </div><div> opts_rh@gsnLeftString <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= leftstring <span class="gmail-Apple-tab-span" style="white-space:pre"> </span></div><div> opts_rh@gsnRightString <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= rightstring</div><div> opts_rh@tiMainFontHeightF <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= 0.030</div><div> opts_rh@gsnLeftStringFontHeightF <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= size<span class="gmail-Apple-tab-span" style="white-space:pre"> </span></div><div> opts_rh@gsnRightStringFontHeightF <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= size<span class="gmail-Apple-tab-span" style="white-space:pre"> </span></div><div> </div><div> </div><div> </div><div> ; Get the contour info for the rh</div><div> contour_rh = gsn_csm_contour(wks,rh_plane,opts_rh)</div><div> </div><div> </div><div>;------------------------------------------------</div><div>; curly vector plot</div><div>;------------------------------------------------</div><div><br></div><div> vecres = opts_xy ; vector only resources</div><div> vecres@vcGlyphStyle = "LineArrow"</div><div> vecres@vcLineArrowThicknessF = 2</div><div> vecres@vcMinDistanceF = 0.008</div><div> vecres@vcRefLengthF = 0.04</div><div> </div><div> </div><div> </div><div> vecres@gsnLeftString <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= " " <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>; add the gsn titles</div><div> vecres@gsnRightString <span class="gmail-Apple-tab-span" style="white-space:pre"> </span>= " "</div><div> vector = gsn_csm_vector(wks,u_plane,w_plane,vecres)</div><div><br></div><div><br></div><div> ; MAKE PLOTS </div><div> overlay(contour_rh,vector)</div><div> draw(contour_rh)</div><div> frame(wks)</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div><br></div><div> end do ; END OF TIME LOOP</div><div><br></div><div>end</div><div>###########################################################################</div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div style="font-size:12.8px"><span style="font-size:12.8px">******************************************************</span><br></div><div style="font-size:12.8px"><span style="font-size:12.8px">盧 涛 (ル タオ) </span><span style="font-family:arial,helvetica,sans-serif;font-size:small">修士課程2年</span><br style="font-size:12.8px"><p style="font-size:12.8px;margin:0px"><font size="2" face="arial, helvetica, sans-serif">〒<span lang="EN-US">112-8551 </span>東京都文京区春日<span lang="EN-US">1-13-27</span></font></p><p style="font-size:12.8px;margin:0px"><font size="2" face="arial, helvetica, sans-serif">中央大学理工学研究科都市環境学専攻</font></p><p style="font-size:12.8px;margin:0px"><font size="2" face="arial, helvetica, sans-serif">河川・水文<span lang="EN-US"><span lang="EN-US"><span lang="EN-US">研究</span></span><span lang="EN-US"><span lang="EN-US">室</span></span></span></font><span style="font-family:arial,helvetica,sans-serif;font-size:small">（山田正教授）</span></p></div><div style="font-size:12.8px"><p style="font-size:12.8px;margin:0px"><font size="2" face="arial, helvetica, sans-serif"><span lang="EN-US"><br></span></font></p><p style="font-size:12.8px;margin:0px"><font size="2" face="arial, helvetica, sans-serif"><span lang="EN-US">TEL: 03-3817-3406; </span></font><span style="font-family:arial,helvetica,sans-serif;font-size:small">Phone: 070-2188-7509</span></p>Email1: <a href="mailto:mail%3Amet.yamos@gmail.com" style="color:rgb(17,85,204);font-size:12.8px" target="_blank">hakufu.asano@gmail.com</a><br style="font-size:12.8px">Email2: <a href="mailto:mail%3Ayamoto@civil.chuo-u.ac.jp" style="color:rgb(17,85,204);font-size:12.8px" target="_blank">lutao@civil.chuo-u.ac.jp</a><br style="font-size:12.8px"><span style="font-size:12.8px">*******************************************************</span></div><div style="font-size:12.8px"><span style="font-size:small;font-family:arial">******************************</span><span style="font-size:small;font-family:arial">************************</span><br></div><div><span style="font-size:small;font-family:arial">Tao Lu (M. Eng)</span></div><div><span style="font-family:arial,helvetica,sans-serif;font-size:12.8px">Laboratory of River Engineering and Hydrology,</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">Dept. of Civil and Environmental Engineering,</span><br style="font-family:arial;font-size:small"><span style="font-family:arial,helvetica,sans-serif;font-size:12.8px">Graduate School of Science and Engineering,</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">Chuo University</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial"> </span><br style="font-family:arial;font-size:small">1-13-27<span style="font-size:small;font-family:arial">,Kasuga,Bunkyo-ku,Tokyo</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">112-8551, Japan</span><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">TEL: </span>03-3817-1805<font size="2" face="arial, helvetica, sans-serif"><span lang="EN-US">; </span></font><span style="font-family:arial,helvetica,sans-serif;font-size:small">Phone: 070-2188-7509</span><br style="font-family:arial;font-size:small"><span style="font-size:12.8px">Email1: </span><a href="mailto:mail%3Amet.yamos@gmail.com" style="font-size:12.8px;color:rgb(17,85,204)" target="_blank">hakufu.asano@gmail.com</a><br style="font-size:12.8px"><span style="font-size:12.8px">Email2: </span><a href="mailto:mail%3Ayamoto@civil.chuo-u.ac.jp" style="font-size:12.8px;color:rgb(17,85,204)" target="_blank">lutao@civil.chuo-u.ac.jp</a><br style="font-family:arial;font-size:small"><span style="font-size:small;font-family:arial">******************************</span><span style="font-size:small;font-family:arial">***********************</span><span style="font-size:12.8px"><br></span></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div>