[ncl-talk] Please help me, I can not turn off contour line label bar
Barry Lynn
barry.h.lynn at gmail.com
Sat Feb 18 23:37:04 MST 2017
Hi:
Can you please send a picture/map...
On Sun, Feb 19, 2017 at 7:55 AM, Tao Lu <hakufu.asano at gmail.com> wrote:
> Hello,
>
> I occurred a problem when I plot contour.
>
> Even though I turn off contour line label bar the graph also shows it.
> " opts_rh = opts_xy
> opts_rh at cnLineLabelsOn = False
> opts_rh at cnFillOn = True
> "
>
> Is there anyone know why?
>
> Thank you,
> Tao
>
>
> Here is my codes:
> ############################################################
> ###############################
> begin
> ;
> ;---Read wrf file
> ;a = addfile("../../wrfout_d01_2015-09-08_12:00:00","r") ; small
> domain ouput 5
> a = addfile("../wrfout_d01_2015-09-08_12:00:00_big2","r") ; big domain
> ouput 10
>
> ; read need data
> xlat = wrf_user_getvar(a, "XLAT",0)
> xlon = wrf_user_getvar(a, "XLONG",0)
> ter = wrf_user_getvar(a, "HGT",0)
>
> ;--Set up workstation
> type = "x11"
> ; type = "pdf"
> ; type = "ps"
> ; type = "ncgm"
> type at wkWidth = 3000 ; Increase size for a
> slightly
> type at wkHeight = 3000
>
> ; Output folder
> folder="./clean_vertical_wind"
> mkdir = systemfunc ("mkdir -p " + folder)
> wks = gsn_open_wks(type,folder+"/plt_CrossSection_3")
>
>
> ;---Set some basic resources
> res = True
> res at gsnMaximize =True
> res at gsnDraw = False
> res at gsnFrame = False
>
>
>
> ;-----------------------------------------------------------
> -------------------
> ;---Set bounding box
> ; C->C' WE
> lats = (/ 36.59845, 36.59845 /)
> lons = (/ 138.50, 141.0 /)
>
> ; X-axis label
> xoption=0
>
> ;0 means lon X-axis, 1 means lat X-axis
> if (xoption .eq. 0) then
> X_desc = "longitude"
> xvalue=xlon
> px = 0.2
> end if
>
> if (xoption .eq. 1) then
> X_desc = "latitude"
> xvalue = xlat
> px = 0.5
> end if
>
>
> ;---Set zoomin indices
> ; loc(0,;) is west-east (x) ; loc(1,:) is south-north (y)
> ; subtract one since we want to use it as an index in NCL
> loc = wrf_user_ll_to_ij(a, lons, lats, True)
> loc = loc - 1
> x_start = loc(0,0)
> x_end = loc(0,1)
> y_start = loc(1,0)
> y_end = loc(1,1)
>
>
> ;---Read time
>
> times = wrf_user_getvar(a,"times",-1) ; get times in the file
> ntimes = dimsizes(times) ; number of times in the file
> FirstTime = True
>
> mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for the file
> nd = dimsizes(mdims)
>
>
> ;-----------------------------------------------------------
> ---------------
> ; Loop of time to do plot
> ;---------------------------------------------------------------
>
> ;do it = 356,356;ntimes-1,18 ; for output5 when need
> 0909_1740 (UTC)
> do it = 178,178;ntimes-1,18 ; for output10 when need
> 0909_1740 (UTC)
>
>
> print("Working on time: " + times(it) )
>
> rh = wrf_user_getvar(a,"rh",it) ; relative humidity
> z = wrf_user_getvar(a, "z",it) ; grid point height
> u = wrf_user_getvar(a,"ua",it) ; u in m/s
> w = wrf_user_getvar(a,"wa",it) ; w in m/s
>
>
> if ( FirstTime ) then ; get height info for labels
> zmin = 0.
> zmax = 12;max(z)/1000.
> nz = floattoint(zmax + 1)
> FirstTime = False
> end if
>
> ;---------------------------------------------------------------
>
> ; Plot a cross session that run from point A to point B
>
> plane = new(4,float)
> plane = (/ x_start,x_end, y_start,y_end /) ; start x;y & end x;y
> point
> opts = True ; start and end
> points specified
>
> rh_plane = wrf_user_intrp3d(rh,z,"v",plane,0.,opts)
> u_plane = wrf_user_intrp3d(u,z,"v",plane,0.,opts)
> w_plane = wrf_user_intrp3d(w,z,"v",plane,0.,opts)
> X_plane = wrf_user_intrp2d(xvalue,plane,0,opts)
>
>
> dim = dimsizes(rh_plane) ; Find the data span -
> for use in labels
> zspan = dim(0)
>
> ; Options for XY Plots
> ; X-axis lables
> dimsX = dimsizes(X_plane)
> xmin = X_plane(0)
> xmax = X_plane(dimsX(0)-1)
> xspan = dimsX(0)-1
> nx = floattoint( (xmax-xmin)/px + 1)
> print ("dimsX= "+dimsX+" xmin= "+xmin+" xmax= "+xmax+" xspan=
> "+xspan+ " nx="+nx)
>
>
> ;---------------------------------------------------------------
>
> ; Options for XY Plots
> opts_xy = res
> opts_xy at tiXAxisString = X_desc
> opts_xy at tiYAxisString = "Height (km)"
> opts_xy at tmXTOn = False
> opts_xy at tmYROn = False
> opts_xy at tmXBMode = "Explicit"
> if (nx .lt. 0) then
> nx = -nx
> end if
> opts_xy at tmXBValues = fspan(0,xspan,nx)
> ; Create tick marks
> opts_xy at tmXBLabels = sprintf("%.1f",fspan(xmin,xmax,nx))
> ; Create labels
> opts_xy at tmXBLabelFontHeightF = 0.015
> opts_xy at tmYLMode = "Explicit"
> opts_xy at tmYLValues = fspan(0,zspan,nz)
> ; Create tick marks
> opts_xy at tmYLLabels = sprintf("%.1f",fspan(zmin,zmax,nz))
> ; Create labels
> opts_xy at tiXAxisFontHeightF = 0.020
> opts_xy at tiYAxisFontHeightF = 0.020
> opts_xy at tmXBMajorLengthF = 0.02
> opts_xy at tmYLMajorLengthF = 0.02
> opts_xy at tmYLLabelFontHeightF = 0.015
>
>
> ; Plotting options for RH
> opts_rh = opts_xy
> opts_rh at cnLineLabelsOn = False
> opts_rh at cnFillOn = True
> opts_rh at lbOrientation = "vertical"
> opts_rh at ContourParameters = (/ 10., 90., 10. /)
> opts_rh at cnFillColors = (/"White","White","White", \
> "White","Chartreuse","Green", \
> "Green3","Green4", \
> "ForestGreen","PaleGreen4"/)
> ; Set contour title
> mainstring = times(it) + " (UTC)"
> leftstring = "Relative humidity %"
> rightstring = "Wind (m/s)"
> size = 0.02
>
> opts_rh at tiMainString = mainstring
> opts_rh at gsnLeftString = leftstring
> opts_rh at gsnRightString = rightstring
> opts_rh at tiMainFontHeightF = 0.030
> opts_rh at gsnLeftStringFontHeightF = size
> opts_rh at gsnRightStringFontHeightF = size
>
>
>
> ; Get the contour info for the rh
> contour_rh = gsn_csm_contour(wks,rh_plane,opts_rh)
>
>
> ;------------------------------------------------
> ; curly vector plot
> ;------------------------------------------------
>
> vecres = opts_xy ; vector only
> resources
> vecres at vcGlyphStyle = "LineArrow"
> vecres at vcLineArrowThicknessF = 2
> vecres at vcMinDistanceF = 0.008
> vecres at vcRefLengthF = 0.04
>
>
>
> vecres at gsnLeftString = " " ; add the gsn titles
> vecres at gsnRightString = " "
> vector = gsn_csm_vector(wks,u_plane,w_plane,vecres)
>
>
> ; MAKE PLOTS
> overlay(contour_rh,vector)
> draw(contour_rh)
> frame(wks)
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> end do ; END OF TIME LOOP
>
> end
> ############################################################
> ###############
>
> --
> ******************************************************
> 盧 涛 (ル タオ) 修士課程2年
>
> 〒112-8551 東京都文京区春日1-13-27
>
> 中央大学理工学研究科都市環境学専攻
>
> 河川・水文研究室（山田正教授）
>
>
> TEL: 03-3817-3406; Phone: 070-2188-7509
> Email1: hakufu.asano at gmail.com <mail%3Amet.yamos at gmail.com>
> Email2: lutao at civil.chuo-u.ac.jp <mail%3Ayamoto at civil.chuo-u.ac.jp>
> *******************************************************
> ******************************************************
> Tao Lu (M. Eng)
> Laboratory of River Engineering and Hydrology,
> Dept. of Civil and Environmental Engineering,
> Graduate School of Science and Engineering,
> Chuo University
>
> 1-13-27,Kasuga,Bunkyo-ku,Tokyo
> 112-8551, Japan
> TEL: 03-3817-1805; Phone: 070-2188-7509
> Email1: hakufu.asano at gmail.com <mail%3Amet.yamos at gmail.com>
> Email2: lutao at civil.chuo-u.ac.jp <mail%3Ayamoto at civil.chuo-u.ac.jp>
> *****************************************************
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
--
Barry H. Lynn, Ph.D
Senior Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581
C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
US 914 432 3108
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170219/d885a58f/attachment.html
More information about the ncl-talk
mailing list