[ncl-talk] attach column of plots
Noelia otero
noeli1680 at gmail.com
Wed Jun 22 09:09:09 MDT 2016
Hi Guido,
Yes, I have different plot per seasons..the thing is that each plot is a
panel plot, I mean plotD contains 10 plot, and so on..since I have 10
models. Then, I wanted to plot all seasons for all plots..that's why I
wanted to do that..maybe there is another simpler way..?
Thanks!
2016-06-22 15:51 GMT+02:00 Guido Cioni <guidocioni at gmail.com>:
> Hi Noelia,
> I assume that your plots are referred to different season, so why do you
> want to attach them? If your primary goal is to panel your 4 plots in a 2x2
> matrix then I would suggest to use the gsn_panel functionality instead.
> You can easily achieve that by doing something like
>
> gsn_panel(wks,(/plotD,plotM,plotJ,plotS/),(/2,2/),True)
>
> Maybe I misunderstood the question, let me know.
> cheers
>
>
> Guido Cioni
> http://guidocioni.altervista.org
>
> On 22 Jun 2016, at 13:51, Noelia otero <noeli1680 at gmail.com> wrote:
>
> Hi all,
>
> I would like to ask for some help to create a plot attaching 4 column of
> plots. I saw the example panel_10.ncl , but it is not clear to me how to do
> it.
>
>
> First I am plotting 4 different column of plots:
>
> plotD=new(10, graphic)
> plotM=new(10, graphic)
> plotJ=new(10, graphic)
> plotS=new(10, graphic)
>
>
> res = True
> res at gsnMaximize = False ; Maximize in frame
> res at gsnDraw = False ; Don't draw plots
> res at gsnFrame = False ; Don't advance frame
> res at cnFillOn = True ; Turn on contour fill
> res at cnInfoLabelOn = False ; Turn off info label
> res at cnLineLabelsOn = False ; Turn off line labels
> res at cnLinesOn = False ; Turn off contour lines
> res at cnLevelSelectionMode = "ExplicitLevels"
> res at lbLabelBarOn = False
> res at gsnAddCyclic = False
>
> res at mpMinLatF = min(lat)
> res at mpMaxLatF = max(lat)
> res at mpMinLonF = min(lon)
> res at mpMaxLonF = max(lon)
>
>
> :get the values
> do i = 0,nplots-1 ;nplots =10
>
> cwtM = systemfunc("ls " + dirM2(i) + "freq*nc")
> do k = 0,3 ;loop seasons
>
> fileM = addfile(cwtM(k), "r")
> names = getfilevarnames(fileM)
> model = fileM->$names(j+3)$(0,:,:)
> ;print(names(j))
> ;print(max(model))
>
> if(k.eq.0) then ; plotting seasons
> plotD(i)= gsn_csm_contour_map(wks,model,res)
> else if(k.eq.1) then
> plotJ(i)= gsn_csm_contour_map(wks,model,res)
> else if(k.eq.2) then
> plotM(i)= gsn_csm_contour_map(wks,model,res)
> else if(k.eq.3) then
> plotS(i)= gsn_csm_contour_map(wks,model,res)
> end if
> end if
> end if
> end if
>
>
> end do
> end do
>
> My idea now would be attach all of them to get the 4 columns together. I
> tried to attach (first, with 2), but it is not working:
>
> attach_id1 = gsn_attach_plots(plotD,plotM,True,True)
>
> I am not sure if attaching the plots would work here..or if there would be
> another way to combine different columns of plots.
> I'd really appreciate any suggestion ,
>
> many thanks
>
> Noelia
>
> _______________________________________________
> 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/20160622/8af5132a/attachment.html
More information about the ncl-talk
mailing list