Re: Panel plot errors

Hi Rick,

Thankyou for your suggestions. I've changed the script and now it terminates without any error and output. Please suggest to avoid termination and to generate output panel plot.
Thanks. Ibo

The code 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=0,11
   filename1 = "6104_pk_"+sprinti("%0.2i",i+1)+"_avg.nc"
   print(filename1)
   filename2 = "6104_ind_"+sprinti("%0.2i",i+1)+"_avg.nc"
   print(filename2)
   filename3 = "6104_sec_"+sprinti("%0.2i",i+1)+"_avg.nc"
   print(filename3)
   f1 = addfile(filename1,"r")
   f2 = addfile(filename2,"r")
   f3 = addfile(filename3,"r")
   if (i.eq.0) then
          a1=f1->a
          a2=f2->a
          a3=f3->a
          dim1=dimsizes(a1)
          dim2=dimsizes(a2)
          dim3=dimsizes(a3)
          p=new((/12,dim1/),"float")
          a=new((/12,dim2/),"float")
          c=new((/12,dim3/),"float")
        end if
        p(i,:)=f1->a
        a(i,:)=f2->a
        c(i,:)=f3->a
        delete(f1)
        delete(f2)
        delete(f3)
    end do
   x1 = dim_standardize_Wrap(a1,1)
   x2 = dim_standardize_Wrap(a2,1)
   x3 = dim_standardize_Wrap(a3,1)
   dimx1 = dimsizes(x1)
   dimx2 = dimsizes(x2)
   dimx3 = dimsizes(x3)
   dimmx = max((/dimx1,dimx2,dimx3/))
   time = ispan(1961,2004,1)
   data = new((/3,dimmx/),"float")
   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,:) = x1
   data(1,:) = x2
   data(2,:) = x3
   plot(i) = gsn_csm_xy(wks,time,data,res) ; plot correlation
   end do
   gsn_panel(wks,plot,(/4,3/),res) ; plot correlation
;end

________________________________
From: Rick Brownrigg <brownrig@ucar.edu>
To: Ibo Ze <ibo85@ymail.com>; Ncl Talk <ncl-talk@ucar.edu>
Sent: Mon, May 24, 2010 6:21:16 AM
Subject: Re: Panel plot errors

Hi,

There are at least two problems that I see:

1. NCL's array indexing is zero-based. Your do_loop indexing is 1-based, ranging from 1-12, indexing into
the arrays "p", "d", "c". You either want something like "do i=0,11" or "p(i-1,:)=f1->a", etc.

2. On the line that reads "data = new((/3,dim/),"float")", the variable "dim" has never been created. You have variables "dim1", "dim2", "dim3"; perhaps one of those was intended?

Its possible that fixing these two issues will clear up the remaining errors (?) Hope that helps...
Rick

On Sun, 23 May 2010 01:53:10 -0700 (PDT)
Ibo Ze <ibo85@ymail.com> wrote:
> 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 Mon May 24 06:47:20 2010