[ncl-talk] Please help me, I can not turn off contour line label bar
Tao Lu
hakufu.asano at gmail.com
Sat Feb 18 22:55:35 MST 2017
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>
*****************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170219/7d277fab/attachment.html
More information about the ncl-talk
mailing list