Re: MJO Phase diagram

Abhik,

As far as I understood you want to make a panel plot. So try work on this
template logic.

__________________________________________________________________
begin

  wks = gsn_open_wks(pltType, pltPath) ; open workstation
  gsn_define_colormap(wks,"radar_1")
  plot=new(10,graphic)

do .......................
......................
..................

end do

 plist=(/0,1,............./)
 resP@gsnPanelRowSpec = True
 gsn_panel(wks,plot(plist),(/3,3,....../),resP)

end
__________________________________________________________________

Since I don't have experience in using "Wheeler-Kiladis diagrams" I cannot
help you further.

Anyway Thanks Dennis Shea and Mary Haley for your comments and supports.

Regards

Gibies George

On 17 August 2012 01:45, Mary Haley <haley@ucar.edu> wrote:

> Abhik,
>
> I do finally see where the dependency comes in on the "n" do loop. Still,
> however, you should not be creating "wks" or "plot" inside the do loop,
> because you will clobber both of these every time.
>
> My suggestion is to move at least these three lines:
>
> wks = gsn_open_wks(pltType, pltPath) ; open workstation
> gsn_define_colormap(wks,"radar_1")
>
> plot = mjo_phase_background(wks, opt) ; generic phase space
> backgrou
>
> to above the line:
>
> do n=0,9
>
> You will also need to move these two lines:
>
> draw(plot)
> frame(wks)
>
> to after the very last "end do" (the one for "do n").
>
> If you still continue to have problems, then you will need to provide your
> data and new version of the script.
>
> --Mary
>
> On Aug 16, 2012, at 1:38 PM, Abhik Santra wrote:
>
> > Hello,
> >
> > I want to draw a single plot for all the 10 years data. Mary
> suggested that I should not use the "do loop". However, I put the call
> workstation only inside the loop and the plot call outside following
> earlier suggestion. But it doesn't solve my problem. The pc1 and pc2
> variables for each year are required to store separately. I am not an
> expert in NCL. So please help me to solve the problem.
> >
> > regards,
> > Abhik
> >
> >
> >
> > ----- Original Message -----
> > From: "Mary Haley" <haley@ucar.edu>
> > To: "Abhik Santra" <abhiksantra@tropmet.res.in>
> > Cc: "ncl-talk" <ncl-talk@ucar.edu>
> > Sent: Thursday, August 16, 2012 8:41:12 PM
> > Subject: Re: MJO Phase diagram
> >
> > Abhik,
> >
> > Gibies is correct that you need to move the initial plot call outside
> the loop, because you are clobbering it for every iteration on "n". Also,
> you put the call to open the workstation inside the do loop, which means
> that's going to get clobbered every time too.
> >
> > This code is hard to debug without being able to run it. I can't tell
> why you have the "do n=0,9" loop. Is there a dependency on "n" somewhere
> inside the loop? You have a bunch of code inside this loop that I think
> should be outside of it. For example, you have:
> >
> > do n=0,9
> > ymdStrt = year+"1016" ; start yyyymmdd
> > ymdLast = year+"0415" ; last
> >
> > pltDir = "./" ; plot directory
> > pltType = "eps"
> > pltName = "mjoclivar" ; yrStrt+"_"+yrLast
> >
> > pltMovie= False ; animation
> > pltTitle= "MJO Phase: 15S-15N: "
> >
> > All of these variables don't change inside the loop, so there's no
> reason for them to be there. Do loops can be expensive time-wise, so you
> want to keep them as minimal as possible.
> >
> > Do you want a separate plot for each value of "n", or are you trying to
> draw things on the same plot for each value of "n"?
> >
> > --Mary
> >
> > On Aug 16, 2012, at 3:57 AM, Abhik Santra wrote:
> >
> >>
> >> Hello,
> >>
> >> Earlier I have tried that also and it was giving similar result.
> However, I have realized that the variable "plot" must be stored inside the
> loop and should be drawn by overlapping. I need some help to do so.
> >> Thank you for your suggestion.
> >>
> >> Regards,
> >> Abhik
> >>
> >>
> >> ----- Original Message -----
> >> From: "gibies george" <gibies@tropmet.res.in>
> >> To: "Abhik Santra" <abhiksantra@tropmet.res.in>
> >> Cc: "ncl-talk" <ncl-talk@ucar.edu>
> >> Sent: Thursday, August 16, 2012 3:14:52 PM
> >> Subject: Re: MJO Phase diagram
> >>
> >> Abhik,
> >>
> >> You can use year loop to read data, but you have to take the plotting
> part outside it, because otherwise it will overwrite on the same plot.
> >>
> >> Regads
> >>
> >> Gibies George
> >>
> >>
> >> On 16 August 2012 12:45, Abhik Santra < abhiksantra@tropmet.res.in >
> wrote:
> >>
> >>
> >> Hi,
> >>
> >> I am trying to plot the MJO Phase Diagram for 10 years (2000-2009) in a
> single panel following the example mjoclivar_15.ncl. So I have introduced a
> do loop to read the multiple years data in the original script. However, it
> is reading all the years successfully but plotting only the last year. I am
> attaching the modified script. Please suggest me what to do?
> >>
> >> With regards,
> >> --
> >> Abhik Santra
> >> CSIR Research Fellow,
> >> Indian Institute of Tropical Meteorology,
> >> Pune - 411008.
> >> India
> >>
> >> _______________________________________________
> >> ncl-talk mailing list
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >>
> >>
> >>
> >> --
> >> _______________________________________________
> >> ncl-talk mailing list
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

-- 
*Gibies George,
CSIR-Senior Research Fellow,*
*Monsoon Group,
Climate and Global Modelling Division,
Indian Institute of Tropical Meteorology,
Dr. Homi Bhabha Road,
NCL (P. O.), Pashan,
Pune 411008, India.
Tel. Off. +912025904306*
*Mobile : +919028523910
               +919881707796
*
*http://sites.google.com/site/gibiesge/*

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Aug 16 23:27:58 2012