plotting bar charts

From: Qiaozhen Mu (qiaozhen AT XXXXXX)
Date: Fri Oct 11 2002 - 08:43:03 MDT


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