[ncl-talk] panel_33.ncl

Tue Nov 20 10:29:18 MST 2018

Hi Laura,

This is an excellent question.

It is a bit of black magic trying to get the plots to panel correctly.
gsn_panel could really stand an overhaul to make it work better in this
regard.

Part of the issue is that gsn_panel always wants to expand the plots to fit
in the given area. This does not include the title or the side labelbars,
just the plots.

The "given area" is what's defined by gsnPanelBottom, gsnPanelTop,
gsnPanelLeft, and gsnPanelRight, which by default form a square with values
0, 1, 0, and 1.

If you are plotting a set of 3x4 plots and another set of 1x3 plots, then
you need to set up the two panel areas for the two sets of plots to have
the correct ratio as the plots you have. For example, for a 3x4 plot and a
1x4 plot, you want to have 3x as much height for the 3x4 plot as the 1x4
plot.  If you had a 2x5 and 3x5 set of plots, then the 2x5 plots need to
take up 2/5 of the total height, and the 3x5 plots 3/5 of the height.

It gets more complicated because when you add in a title or a labelbar,
then you need to leave space for those.

I will try to come up with an example for you today. This is something I've
been meaning to do.

--Mary

On Tue, Nov 20, 2018 at 6:46 AM Laura Fowler <laura at ucar.edu> wrote:

> Hi Mary:
>
> Many thanks for taking the time to rework the script. Using
> panel_33_mod.ncl, I was able to modify my script to produce the plot I
> needed (see attached totPREC.png and test.ncl in
> /gpfs/fs1/scratch/laura/MPASprojects2018/MPAS.cu_shallow/novemberPUBLICATION/precipitation_figs).
> I need a little more work to align DEC 2015 and JUN 2016 to their
> respective columns).
>
> I further modified test.ncl (test2.ncl in the same directory) to remove
> one row from the top panel so that I can have 2 panels of two rows with one
> color bar for each panel. I do not understand the logic on how to modify
> res at gsnPanelBottom to make this work on a more systematic basic. It seems
> like there is a lot of "guessing" in doing this?
>
> Thanks,
> Laura
>
>
> On Sat, Nov 17, 2018 at 1:59 PM Mary Haley <haley at ucar.edu> wrote:
>
>> Hi Laura,
>>
>> Try this script. I cleaned up the logic of the panelling to base the size
>> of the difference plots on the area used for the top set of plots.
>>
>> I included a subset of the data file so you can run the modified script.
>>
>> I'll include this on the panel page later. I'll also clean up the
>> panel_33.ncl to be more user friendly to modify.
>>
>> --Mary
>>
>>
>> On Sat, Nov 17, 2018 at 4:29 AM Rashed Mahmood <rashidcomsis at gmail.com>
>> wrote:
>>
>>> I think the first issue is related to setting the bottom limit for the
>>> first set of plots:
>>>
>>>     resP at gsnPanelBottom    = 0.33              ; or something similar
>>>
>>> The second issue is likely related to maximize_output(wks,True) that would try to use maximum available space. May be you can try commenting this.
>>>
>>> Cheers,
>>>
>>> Rashed
>>>
>>>
>>>
>>>
>>>
>>> On Fri, Nov 16, 2018 at 1:16 PM Laura Fowler <laura at ucar.edu> wrote:
>>>
>>>> Hello:
>>>>
>>>> I am trying to produce a panel plot that looks like panel_33.ncl but
>>>> instead of having 2 rows for the actual fields, I have 3 rows, and instead
>>>> of having 1 row for the difference fields, I have 2 rows. I have two
>>>> issues: The first issue is that I am running out of room on the page to
>>>> plot the last 2 rows and I do not know what to do next? The second issue
>>>> that the top 3 rows are plotted in a landscape mode although I set:
>>>>
>>>> resP at gsnPaperOrientation = "portrait"
>>>>
>>>>
>>>>
>>>> Here is the script that I used after I created my 10 individual plots:
>>>>
>>>>
>>>> ;--- create panel plot:
>>>>
>>>> resP = True
>>>>
>>>> resP at gsnDraw             = False
>>>>
>>>> resP at gsnFrame            = False
>>>>
>>>> resP at gsnPanelLabelBar    = True
>>>>
>>>> resP at gsnPaperOrientation = "portrait"
>>>>
>>>>
>>>> resP at lbOrientation      = "vertical"
>>>>
>>>> resP at lbLabelFontHeightF = 0.01
>>>>
>>>> ;resP at pmLabelBarWidthF  = 0.075
>>>>
>>>> ;resP at pmLabelBarHeightF = 0.40
>>>>
>>>>
>>>>
>>>> panelid1 = gsn_panel_return(wks,plots(0:5),(/3,2/),resP)
>>>>
>>>> bb = NhlGetBB(panelid1)
>>>>
>>>> top    = max(bb(:,0))
>>>>
>>>> bottom = min(bb(:,1))
>>>>
>>>> height = (top-bottom)/3.
>>>>
>>>> print("--- top    = "+top)
>>>>
>>>> print("--- bottom = "+bottom)
>>>>
>>>> print("--- height = "+height)
>>>>
>>>> print(" ")
>>>>
>>>>
>>>> ;resP at gsnPanelTop   = bottom-0.01    ; lower 1/3 plus some white space
>>>>
>>>> resP at gsnPanelTop    = bottom
>>>>
>>>> resP at gsnPanelBottom = resP at gsnPanelTop - height
>>>>
>>>> print("--- gsnPanelTop = "+resP at gsnPanelTop)
>>>>
>>>> print("--- gsnPanelBot = "+resP at gsnPanelBottom)
>>>>
>>>>
>>>> ;panelid2 = gsn_panel_return(wks,plots(6:9),(/2,2/),resP)
>>>>
>>>> maximize_output(wks,True)
>>>>
>>>>
>>>> ;resP at gsnDraw             = True
>>>>
>>>> ;resP at gsnFrame            = False
>>>>
>>>> ;gsn_panel(wks,plots(0:5),(/3,2/),resP)
>>>>
>>>> ;frame(wks)
>>>>
>>>>
>>>>
>>>> When print gsnPanelTop and gsnPanelBottom, I get:
>>>>
>>>>
>>>> (0) --- end create individual plot:
>>>>
>>>> (0)
>>>>
>>>> (0) --- top    = 0.889811
>>>>
>>>> (0) --- bottom = 0.110189
>>>>
>>>> (0) --- height = 0.259874
>>>>
>>>> (0)
>>>>
>>>> (0) --- gsnPanelTop = 0.110189
>>>>
>>>> (0) --- gsnPanelBot = -0.149686
>>>>
>>>>
>>>>
>>>> I am lost as what to do next. Help and suggestions would be greatly
>>>> appreciated.
>>>>
>>>>
>>>>
>>>> Thanks,
>>>>
>>>> Laura
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> --
>>>>
>>>> !-------------------------------------------------------------------------------------------------------------
>>>> Laura D. Fowler
>>>>
>>>> Mesoscale and Microscale Meteorology Division (MMM)
>>>> National Center for Atmospheric Research
>>>> P.O. Box 3000, Boulder CO 80307-3000
>>>>
>>>> e-mail: laura at ucar.edu
>>>> phone: 303-497-1628
>>>>
>>>>
>>>> !-------------------------------------------------------------------------------------------------------------
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>
> --
>
> !-------------------------------------------------------------------------------------------------------------
> Laura D. Fowler
>
> Mesoscale and Microscale Meteorology Division (MMM)
> National Center for Atmospheric Research
> P.O. Box 3000, Boulder CO 80307-3000
>
> e-mail: laura at ucar.edu
> phone: 303-497-1628
>
>
> !-------------------------------------------------------------------------------------------------------------
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181120/b6ffec01/attachment.html>