Hi, All
I have a problem when I overlay a gsn_csm_xy2 on a gsn_csm_xy, which
res@gsnXYBarChart = True
Here is my plot:
http://shinx.myweb.hinet.net/TMIN.jpg
The blue line with triangle marker is for y1(left side y-Axis) , and the red
line with dot marker is for y2(right side y-Axis)
As the plot shows, y2 is not mapped properly to its X-Axis. It should be
like y1 does.
By the way, They were fine until I turned res(sky blue data)@gsnXYBarChart =
True
Any idea? Thanks,
Shinx
;===============================================================
; test.ncl
;===============================================================
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
; ========================================================================
begin
;=========================================================================
;=========================================================================
inputfn="TMIN.txt"
modelname="TMIN"
;=======================
; start to read in data
;=======================
dummy = asciiread(inputfn,-1,"string")
data = new((/8,12/),"float")
data(0,:) = stringtofloat(str_get_field(dummy(1::), 1," "))
data(1,:) = stringtofloat(str_get_field(dummy(1::), 2," "))
data(2,:) = stringtofloat(str_get_field(dummy(1::), 3," "))
data(3,:) = stringtofloat(str_get_field(dummy(1::), 4," "))
data(4,:) = stringtofloat(str_get_field(dummy(1::), 5," "))
data(5,:) = stringtofloat(str_get_field(dummy(1::), 6," "))
data(6,:) = stringtofloat(str_get_field(dummy(1::), 7," "))
data(7,:) = stringtofloat(str_get_field(dummy(1::), 8," "))
dataset=new(8,"string")
dataset(1)=str_get_field(dummy(0), 1," ")
dataset(2)=str_get_field(dummy(0), 2," ")
dataset(3)=str_get_field(dummy(0), 3," ")
dataset(4)=str_get_field(dummy(0), 4," ")
dataset(5)=str_get_field(dummy(0), 5," ")
dataset(6)=str_get_field(dummy(0), 6," ")
dataset(7)=str_get_field(dummy(0), 7," ")
x=ispan(1,12,1)
outputfn=modelname
;====================
; PLOT
;====================
wks = gsn_open_wks("ps",outputfn)
res0 = True ; lines with no marker
resH = True ; Bar chart , Histogram
resL = True ; y1
resR = True ;y2
resH@gsnDraw = False
resH@gsnFrame = False
resH@tmXBMode = "Explicit"
resH@tmXBValues = data(0,:)
resH@tmXBLabels = data(0,:)
resH@gsnMaximize = True
resH@gsnPaperOrientation = "portrait"
resR = resH
resH@trYMinF = 0
resH@trYMaxF = 3.5
res0 = resH
resL = resH
resH@tmXTOn = False
resH@tmYROn = False
resH@gsnXYBarChart = True ; turn on bar chart
resH@gsnXYBarChartBarWidth = 0.75 ; change bar widths
resH@gsnXYBarChartColors = "SkyBlue"
res0@xyDashPattern = 0
res0@xyMarkLineMode = "Lines" ; Markers *and* lines
res0@xyLineColors = "black"
res0@tmXTOn = False
res0@tmYROn = False
resR@xyDashPattern = 0
resR@trYMinF = 0
resR@trYMaxF = 1
resR@tiYAxisString = "COR"
resR@tiYAxisFontHeightF = 0.01
resR@xyMarkLineMode = "MarkLines" ; Markers *and* lines
resR@xyMarkers = 16
resR@xyLineColors = "red"
resR@xyMarkerColors = "red"
resR@xyMarkerSizeF = 0.01
resL@xyDashPattern = 0
resL@tiYAxisString = "RMS"
resL@tiYAxisFontHeightF = 0.01
resL@xyMarkLineMode = "MarkLines" ; Markers *and* lines
resL@xyMarkers = 7
resL@xyLineColors = "blue"
resL@xyMarkerColors = "blue"
resL@xyMarkerSizeF = 0.01
plotH=gsn_csm_xy(wks,data(0,:),data(1,:),resH)
plot0=gsn_csm_xy(wks,data(0,:),data(3:6,:),res0)
plotL=gsn_csm_xy2(wks,data(0,:),data(2,:),data(7,:),resL,resR)
overlay(plotH,plot0)
overlay(plotH,plotL)
draw(plotH)
frame(wks)
system("convert -density 300 -trim "+outputfn+".ps "+outputfn+".png")
system("rm "+outputfn+".ps") ;delete .ps file
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 9 00:10:47 2010
This archive was generated by hypermail 2.1.8 : Thu Feb 11 2010 - 09:50:35 MST