Panel plot errors

From: Ibo Ze <ibo85_at_nyahnyahspammersnyahnyah>
Date: Sun May 23 2010 - 02:53:10 MDT

Hi,

I want to plot standardized anomaly for all months for three different regions in a panel plot. Getting following error:

fatal:Subscript out of range, error in subscript #0
fatal:Execute: Error occurred at or near line 45 in file stnd_anomaly.ncl
fatal:Variable (dim) is undefined
fatal:Execute: Error occurred at or near line 59 in file stnd_anomaly.ncl
warning:dim_standardize: 1 rightmost sections of the input array contained all missing values
fatal:Assignment type mismatch, right hand side can't be coerced to type of left hand side
fatal:Execute: Error occurred at or near line 64 in file stnd_anomaly.ncl
warning:dim_standardize: 1 rightmost sections of the input array contained all missing values
fatal:Assignment type mismatch, right hand side can't be coerced to type of left hand side
fatal:Execute: Error occurred at or near line 65 in file stnd_anomaly.ncl
warning:dim_standardize: 1 rightmost sections of the input array contained all missing values
fatal:Assignment type mismatch, right hand side can't be coerced to type of left hand side
fatal:Execute: Error occurred at or near line 66 in file stnd_anomaly.ncl
fatal:Undefined identifier: (data) is undefined, can't continue
fatal:Execute: Error occurred at or near line 88 in file stnd_anomaly.ncl
(0) Error: gsn_panel: all of the plots passed to gsn_panel appear to be invalid

*******************************************************************

Please help to rectify this error.
Thanks. Ibo
The script i am using is as follows:

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/contributed.ncl"
begin
  do i=1,12
   filename1 = "6104_pk_"+sprinti("%0.2i",i)+"_avg.nc"
   print(filename1)
   filename2 = "6104_ind_"+sprinti("%0.2i",i)+"_avg.nc"
   print(filename2)
   filename3 = "6104_sec_"+sprinti("%0.2i",i)+"_avg.nc"
   print(filename3)
   f1 = addfile(filename1,"r")
   f2 = addfile(filename2,"r")
   f3 = addfile(filename3,"r")
   if (i.eq.1) then
      a1 = f1->a
      a2 = f2->a
      a3 = f3->a
      dim1 = dimsizes(a1)
      dim2 = dimsizes(a2)
      dim3 = dimsizes(a3)
      p=new((/12,dim1/),"float")
      d=new((/12,dim2/),"float")
      c=new((/12,dim3/),"float")
   end if
   p(i,:)=f1->a
   d(i,:)=f2->a
   c(i,:)=f3->a
   delete(f1)
   delete(f2)
   delete(f3)
  end do
   time = ispan(1961,2004,1)

   data = new((/3,dim/),"float")
   f1 = dim_standardize_Wrap(p(:,:),1)
   f2 = dim_standardize_Wrap(d(:,:),1)
   f3 = dim_standardize_Wrap(c(:,:),1)

   wks = gsn_open_wks("x11","stanom_pic_06") ; open a ps plot
   plot = new(12,graphic)
   res = True ; make plot mods
   res@gsnDraw = False
   res@gsnFrame = False
   res@xyLineThicknesses = (/2.,2.,2./)
   res@xyLineColors = (/"red","blue","orange"/)
   res@xyDashPatterns = (/0.,0.,0./) ; make all lines solid
   res@trYMinF = -5.0 ; bring bars down to zero
   res@trYMaxF = 5.0 ; adds space on either end
   res@trXMinF = 1960
   res@trXMaxF = 2004
   res@gsnYRefLine = 0. ; reference line
   res@tiXAxisString = "Time(years)" ; x-axis label
   res@tiYAxisString = "Standardized Anomaly" ; x-axis label
   do i=0,11
   data(0,:)=f1(i,:)
   data(1,:)=f2(i,:)
   data(2,:)=f3(i,:)

   plot(i) = gsn_csm_xy(wks,time,data,res) ; plot correlation
   end do
   gsn_panel(wks,plot,(/4,3/),res) ; plot correlation

end

      

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun May 23 02:53:17 2010

This archive was generated by hypermail 2.1.8 : Wed May 26 2010 - 10:39:13 MDT