Hi All,
When I use this script to plot the bar charts, the error is
warning:Attempt to reference attribute (gsnXYBarChartBarWidth) which is
undefined
fatal:Subscript out of range, error in subscript #1
fatal:Execute: Error occurred at or near line 7244 in file
$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl
fatal:Execute: Error occurred at or near line 116
I am not sure what's wrong.
Best regards,
Qiaozhen
***************************************
*The University of Texas at Austin *
*Institute for Geophysics *
*4412 Spicewood Springs Rd., Bldg. 600*
*Austin, Texas 78759-8500 *
*phone:(512) 471-0462 *
*fax :(512) 471-0348 *
***************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
fields=(/"CLDTOT","FLNT","FSDS","FSNT", \
"LHFLX", "PRECT","PSL","SHFLX",\
"TREFHT" /)
clim=(/"DJF","MAM","JJA","SON","ANN"/)
inptr = (/\
"alfa3_hsnoc15_","alfa3_ricr5_", \
"alfa3_tau216_", "hsnoc15_tau216_", \
"rhmin86_alfa3_","rhmin86_hsnoc15_", \
"rhmin86_ricr5_","rhmin86_tau216_", \
"ricr5_hsnoc15_","ricr5_tau216_" /)
inptr1 = (/\
"alfa3_","hsnoc15_", \
"ricr5_","ccm86_", \
"tau216_" /)
nfields = dimsizes(fields)
parameter = (/"alfa","hsnoc","ricr","rhmin","tau"/)
ninptr = dimsizes(inptr)
ninptr1 = dimsizes(inptr1)
nclim = dimsizes(clim)
nparameter = dimsizes(parameter)
;************************************************
do infields =0, nfields-1
;************************************************
do j = 0, nclim-1
;************************************************
;************************************************
error = new((/ninptr/), "float")
error!0 ="experiment"
error&experiment= ispan(1,ninptr,1)
;************************************************
error1 = new((/ninptr1/), "float")
error1!0 ="experiment"
error1&experiment= ispan(1,ninptr1,1)
;************************************************
ascii_file = "./nc/"+fields(infields)+"_"+clim(j)+"_to_ccmavg_EOFCOMP_error.asc"
error=asciiread(ascii_file,(/ninptr/),"float")
;************************************************
;************************************************
ascii_file1 = "./nc/"+fields(infields)+"_"+clim(j)+"_indiv_to_ccmavg_EOFCOMP_error.asc"
error1=asciiread(ascii_file1,(/ninptr1/),"float")
;************************************************
error_comp = new((/3*ninptr/), "float")
error_comp = (/0.0/)
error_comp!0 ="all"
error_comp&all= ispan(1,3*ninptr,1)
;************************************************
error_comp(0)=(/error(0)/)
error_comp(1)=(/error1(0)/)
error_comp(2)=(/error1(1)/)
error_comp(3)=(/error(1)/)
error_comp(4)=(/error1(0)/)
error_comp(5)=(/error1(2)/)
error_comp(6)=(/error(2)/)
error_comp(7)=(/error1(0)/)
error_comp(8)=(/error1(4)/)
error_comp(9)=(/error(3)/)
error_comp(10)=(/error1(1)/)
error_comp(11)=(/error1(4)/)
error_comp(12)=(/error(4)/)
error_comp(13)=(/error1(3)/)
error_comp(14)=(/error1(0)/)
error_comp(15)=(/error(5)/)
error_comp(16)=(/error1(3)/)
error_comp(17)=(/error1(1)/)
error_comp(18)=(/error(6)/)
error_comp(19)=(/error1(3)/)
error_comp(20)=(/error1(2)/)
error_comp(21)=(/error(7)/)
error_comp(22)=(/error1(3)/)
error_comp(23)=(/error1(4)/)
error_comp(24)=(/error(8)/)
error_comp(25)=(/error1(2)/)
error_comp(26)=(/error1(1)/)
error_comp(27)=(/error(9)/)
error_comp(28)=(/error1(2)/)
error_comp(29)=(/error1(4)/)
ascii_file = "./nc/"+fields(infields)+"_"+clim(j)+"_to_ccmavg_compare_EOFCOMP_error.asc" ; Create name of ASCII file.
system("/bin/rm -f " + ascii_file) ; Remove ASCII file.
asciiwrite(ascii_file,error_comp) ; Write part of temp to ASCII file.
;************************************************
wks = gsn_open_wks ("ps", "./figure/"+fields(infields)+"_"+clim(j)+"_to_ccmavg_compare_EOFCOMP_error")
gsn_define_colormap(wks,"gui_default") ; choose colormap
res = True ; Use plot options
res@gsnDraw = False ; Do not draw plot
res@gsnFrame = False ; Do not advance frame
res@xyDashPatterns = (/0./) ; choose dash patterns
res@xyLineThicknesses = (/1./) ; choose line thickness
res@xyLineColors = (/"black"/)
res@xyMarkLineModes = "Markers"
res@xyDashPatterns = 0
res@xyMarkers = 16
res@tiYAxisString = fields(infields)
res@tiXAxisString = "combine"
plot = new(1,graphic) ; create graphic array
; res@tiMainString = fields(infields)+"_"+clim(j)+"_to_ccmavg_EOFCOMP_error"
res@gsnCenterString = fields(infields)+"_"+clim(j)+"_to_ccmavg_EOFCOMP_error"
plot = gsn_csm_xy(wks,error_comp&all,error_comp,res)
pres = True ; panel mods desired
pres@gsnMaximize = True ; fill the page
gsn_panel(wks,plot,(/1,1/),pres) ; draw: 2-across, 2-down
;************************************************
delete(error)
delete(error1)
delete(error_comp)
;************************************************
delete(res)
delete(plot)
delete(pres)
delete(wks)
;************************************************
end do
;************************************************
;************************************************
error=new((/ninptr/), "float")
error!0 ="experiment"
error&experiment=ispan(1,ninptr,1)
;************************************************
error1=new((/ninptr1/), "float")
error1!0 ="experiment"
error1&experiment=ispan(1,ninptr1,1)
;************************************************
;************************************************
ascii_file = "./nc/"+fields(infields)+"_to_ccmavg_EOFCOMP_error.asc"
error=asciiread(ascii_file,(/ninptr/),"float")
;************************************************
ascii_file1 = "./nc/"+fields(infields)+"_indiv_to_ccmavg_EOFCOMP_error.asc"
error1=asciiread(ascii_file1,(/ninptr1/),"float")
;************************************************
; error_comp = new((/3*ninptr/), "float")
error_comp = new((/30/), "float")
error_comp = (/0.0/)
error_comp!0 ="all"
error_comp&all= ispan(1,30,1)
;************************************************
error_comp(0)=(/error(0)/)
error_comp(1)=(/error1(0)/)
error_comp(2)=(/error1(1)/)
error_comp(3)=(/error(1)/)
error_comp(4)=(/error1(0)/)
error_comp(5)=(/error1(2)/)
error_comp(6)=(/error(2)/)
error_comp(7)=(/error1(0)/)
error_comp(8)=(/error1(4)/)
error_comp(9)=(/error(3)/)
error_comp(10)=(/error1(1)/)
error_comp(11)=(/error1(4)/)
error_comp(12)=(/error(4)/)
error_comp(13)=(/error1(3)/)
error_comp(14)=(/error1(0)/)
error_comp(15)=(/error(5)/)
error_comp(16)=(/error1(3)/)
error_comp(17)=(/error1(1)/)
error_comp(18)=(/error(6)/)
error_comp(19)=(/error1(3)/)
error_comp(20)=(/error1(2)/)
error_comp(21)=(/error(7)/)
error_comp(22)=(/error1(3)/)
error_comp(23)=(/error1(4)/)
error_comp(24)=(/error(8)/)
error_comp(25)=(/error1(2)/)
error_comp(26)=(/error1(1)/)
error_comp(27)=(/error(9)/)
error_comp(28)=(/error1(2)/)
error_comp(29)=(/error1(4)/)
ascii_file = "./nc/"+fields(infields)+"_to_ccmavg_compare_EOFCOMP_error.asc" ; Create name of ASCII file.
system("/bin/rm -f " + ascii_file) ; Remove ASCII file.
asciiwrite(ascii_file,error_comp) ; Write part of temp to ASCII file.
;************************************************
wks = gsn_open_wks ("ps", "./figure/"+fields(infields)+"_to_ccmavg_compare_EOFCOMP_error")
gsn_define_colormap(wks,"gui_default") ; choose colormap
res = True ; Use plot options
res@gsnDraw = False ; Do not draw plot
res@gsnFrame = False ; Do not advance frame
res@xyDashPatterns = (/0./) ; choose dash patterns
res@xyLineThicknesses = (/1./) ; choose line thickness
res@xyLineColors = (/"black"/)
res@xyMarkLineModes = "Markers"
res@xyDashPatterns = 0
res@xyMarkers = 16
res@tiYAxisString = fields(infields)
res@tiXAxisString = "combine"
plot = new(1,graphic) ; create graphic array
; res@tiMainString = fields(infields)+"_to_ccmavg_EOFCOMP_error"
res@gsnCenterString = fields(infields)+"_to_ccmavg_EOFCOMP_error"
plot = gsn_csm_xy(wks,error_comp&all,error_comp,res)
pres = True ; panel mods desired
pres@gsnMaximize = True ; fill the page
gsn_panel(wks,plot,(/1,1/),pres) ; draw: 2-across, 2-down
;************************************************
delete(error)
delete(error1)
delete(error_comp)
;************************************************
delete(res)
delete(plot)
delete(pres)
delete(wks)
;************************************************
end do
end
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Fri Oct 11 2002 - 08:56:18 MDT