[ncl-talk] panel_33.ncl
Barry Lynn
barry.h.lynn at gmail.com
Tue Nov 27 12:17:38 MST 2018
Hello:
I just had a similar problem with a lack of space for the label bar.
I plotted 8 panels on a page (all the same size). The label bar was left
on the bottom (with numbers), but there was no label.
I ended up putting the label on the side of the bar.
On Tue, Nov 27, 2018 at 9:10 PM Mary Haley <haley at ucar.edu> wrote:
> Hi Laura,
>
> Since paneling two sets of contour plots each with their own labelbar, is
> a common question, I decided to try to make this easier by creating a
> "panel_two_sets" function.
>
> This function calculates the gsnPanelTop / Bottom / Left / Right values
> for you. It also creates the labelbars and will figure out the proper
> layout of the plots based on whether you have a vertical or horizontal
> orientation.
>
> See panel_15.ncl,. panel_18.ncl, panel_33.ncl, panel_5x2_33.ncl, and
> panel_34.ncl all of which now use this function.
>
> http://www.ncl.ucar.edu/Applications/panel.shtml
>
> The "panel_two_sets" procedure is sort of like gsn_panel, except now you
> give it two sets of plots to panel, each with their own dimensions for
> paneling and their own panel resources. There's a third resource list that
> is solely for setting the gsnPanelXXXX values, if needed
>
> The call to this function would look something like this, if no extra
> resources were set:
>
> ;---Panel the two sets of plots
> dims1 = (/5,3/) ; 5 rows, 3 columns
> dims2 = (/2,3/) ; 2 rows, 3 columns
>
> pres1 = True
> pres2 = True
> pres = True
>
> panel_two_sets(wks,plots,diff_plots,dims1,dims2,pres1,pres2,pres)
>
>
> In the above case, since the number of rows are the same for both plots,
> the function is going put the first set of plots on the left, and the
> second set on the right, each with a vertical labelbar.
>
> One gotcha that I'm struggling with, is getting the two labelbars to line
> up. They will line up just fine if the longest label in both label bars is
> roughly the same size. But, if you have one labelbar with small width
> labels, like "1", "2", and "3", and then the second labelbar has labels
> like "1000", "2000", etc, they will be slightly shifted from each other. I
> think this is because the labelbars are being centered in that space, and
> there's no resource that turns this feature off. I'm going to look into
> this more, because it's been bothering me for years trying to figure out
> how to control these labels!
>
> --Mary
>
>
> On Tue, Nov 20, 2018 at 12:26 PM Laura Fowler <laura at ucar.edu> wrote:
>
>> Hi Mary:
>> Thanks. This is indeed very confusing to me. For instance, in
>> panel_33_mod.ncl, if I comment out the lines (see lines below) but keep
>> panelid1 and maximize_output, then I get the following png file (see
>> attached). I thought it was just going to plot the top six figures but it
>> added an other figure at the bottom.
>>
>> So I am not sure what the combination of gsn_panel_return with
>> maximize_output does. And I did not find documentation on the 2 functions.
>>
>> Cheers,
>> Laura
>>
>>
>> ; bb = NhlGetBB(plots(nyears_cols:))
>>
>> ; bottom = min(bb(:,1))
>>
>> ; top = max(bb(:,0))
>>
>> ; height = (top-bottom)
>>
>> ; print("--- bottom = "+bottom)
>>
>> ; print("--- top = "+top)
>>
>> ; print("--- height = "+height)
>>
>>
>>
>>
>> ;---Panel the difference plots at the bottom
>>
>> ; delete(pres at gsnPanelMainString)
>>
>> ; pres at gsnPanelTop := bottom - 0.01 ; the 0.01 gives a little
>> spaces b/w the two sets of plots
>>
>> ; pres at gsnPanelBottom := pres at gsnPanelTop - height
>>
>> ; pres at lbBoxEndCapStyle = "TriangleBothEnds"
>>
>> ; pres at pmLabelBarHeightF = (pres at gsnPanelTop - pres at gsnPanelBottom) *
>> 0.8
>>
>>
>> ; panelid2 =
>> gsn_panel_return(wks,diff_plots,(/nyears_rows-1,nyears_cols/),pres)
>>
>>
>> ; drawNDCGrid(wks)
>>
>> ; draw(wks)
>>
>> ; frame(wks)
>>
>> On Tue, Nov 20, 2018 at 10:29 AM Mary Haley <haley at ucar.edu> wrote:
>>
>>> 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
>>>>
>>>>
>>>> !-------------------------------------------------------------------------------------------------------------
>>>>
>>>
>>
>> --
>>
>> !-------------------------------------------------------------------------------------------------------------
>> 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
>
--
Barry H. Lynn, Ph.D
Senior Associate Scientist, Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581
C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20181127/085656a4/attachment.html>
More information about the ncl-talk
mailing list