[ncl-talk] Multiple time series plot of same parametre but from different ascii files

[Adam Phillips]

Hi Najib,
I see the issue. I did not see from your original ncl-talk posting that
your n do loop is closed rather quickly between lines 26 and 28:
do n= 0,nfiles-1
      nfiles = n
end do

A quick test shows that that syntax will result in nfiles = 0 and then the
do loop concludes. I think what you want is to change your coding from this:

do n= 0,nfiles-1
   nfiles = n
end do
print(files(n))
z1 = asciiread(files(n),(/nday+1/),"string")
printVarSummary(z1)
....
wks = gsn_open_wks("png","AOT_550")
....
plot = gsn_csm_xy(wks, x, AOT_550, res)
draw(plot)
frame(wks)
end

to this:

do n= 0,nfiles-1
  print(files(n))
  z1 = asciiread(files(n),(/nday+1/),"string")
  printVarSummary(z1)

....

  wks = gsn_open_wks("png","AOT_550_"+n)

....

  plot = gsn_csm_xy(wks, x, AOT_550, res)
  draw(plot)
  frame(wks)
end do

If after you implemented the above you are still having problems respond to
the ncl-talk list.
Adam




On Sun, Mar 26, 2017 at 12:55 AM, Najib Yusuf <najibgal at yahoo.com> wrote:

> Dear ncl,
>
> I effected the correction but still getting a single plot.
> Kindly assist.
>
> Thank you.
>
> Najib
>
>
> On Friday, March 24, 2017 6:02 PM, Adam Phillips <asphilli at ucar.edu>
> wrote:
>
>
> Hi Najib,
> I think your script is close. You are redefining the workstation each time
> through your n do loop. Thus, each time through the do loop the output png
> file gets overwritten. Change this line:
> wks = gsn_open_wks("png","AOT_550")
> to this:
> wks = gsn_open_wks("png","AOT_550_"+n)
>
> and I think you will get what you want.
>
> I hope that helps. If not, or if you have further questions, please
> respond to the ncl-talk email list and not to me personally.
> Adam
>
>
> On Fri, Mar 24, 2017 at 3:47 AM, Najib Yusuf <najibgal at yahoo.com> wrote:
>
>
> Dear ncl-talk,
>
> I am trying as I'm still at the initial stage learning ncl going through
> user-guides tutorials etc. I want to plot multiple time series with the
> parametre from different files (AOT_550_2005 to AOT_550_2013). I have
> writing the codes (see attached file) but when I run it I got only one plot
> of AOT_550_2005 not with the remaining years.
>
> Kindly assist on how the remaining years will appear on the plot.I also
> tried looping but the errors are getting much so I live it this way since
> there is no error.
>
> Thank you
>
> Najib
>
> ______________________________ _________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/ mailman/listinfo/ncl-talk
> <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>
>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726 <(303)%20497-1726>
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>
>
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170328/587f78ab/attachment.html