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

Adam Phillips asphilli at ucar.edu
Wed Mar 29 08:11:51 MDT 2017


Hi Najib,
Please send the current version of your script to ncl-talk so all can see and offer assistance.
Adam

> On Mar 29, 2017, at 7:40 AM, Najib Yusuf <najibgal at yahoo.com> wrote:
> 
> 
> Dear all,
> 
> Thank you Adam for the help. I change it as you pointed out but I still got a single plot.I don't understand where the problem might be. Please anyone help.Details in the previous messages.
> 
> Thank you.
> 
> Najib
> 
> 
> On Tuesday, March 28, 2017 8:55 PM, Adam Phillips <asphilli at ucar.edu> wrote:
> 
> 
> 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
> 
> 
> 
> 
> -- 
> Adam Phillips 
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/ asphilli/   303-497-1726
> 
> 
> 
> 
> 
> 
> -- 
> Adam Phillips 
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
> 
> 
> 
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170329/62c8fe15/attachment.html 


More information about the ncl-talk mailing list