<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Matt,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Check out the PyNGL gallery and see the irregular.py example:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.pyngl.ucar.edu/Examples/gallery.shtml">http://www.pyngl.ucar.edu/Examples/gallery.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default" style="font-size:small">I think you need to set:</div><div class="gmail_default" style="font-size:small"><pre style="color:rgb(0,0,0);font-weight:normal;text-decoration:none;white-space:pre-wrap">resources.nglYAxisType = "LinearAxis"</pre>which may seem counterintuitive, but what this says is to linearize the Y axis such that the irregular values on this axis are linearly spaced. The example above graphically shows what this looks like.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 18, 2018 at 6:09 PM, Fearon, Dr. Matthew, Contractor, Code 7533 via pyngl-talk <span dir="ltr"><<a href="mailto:pyngl-talk@ucar.edu" target="_blank">pyngl-talk@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple">
<div class="m_-6837964633838654641WordSection1">
<p class="MsoNormal">Dear PyNgl users,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I’m trying to plot a vertical cross section with regional model data on the native terrain-following vertical coordinate, in my case, that coordinate is sigma height. The levels are *not* equally spaced. However, when I plot the data, the
vertical spacing between the levels is by default assumed to be equally spaced, at least it appears that way. Please see attached cross section plot, and in particular, the values on the vertical axis. How can I plot this properly such that the levels are
not assumed to be equally spaced? Below is my plotting code where I use sfYArray to define my vertical axis/increment. I thought this approach with sfYArray would accommodate the unequal spacing, but it does not. I’ve explored a bit trYAxisType = IrregularAxis,
but have been unsuccessful in getting it work.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Also, here are my vertical levels stored in the variable “yl” used in the plotting routine.<u></u><u></u></p>
<p class="MsoNormal">[5.0, 15.0, 26.0, 38.0, 50.0, 62.0, 75.0, 89.0, 104.0, 120.0, 137.0, 155.0, 174.0, 195.0, 218.0, 242.0, 267.0, 294.0, 323.0, 353.0, 385.0, 419.0, 455.0, 494.0, 535.0, 578.0, 624.0, 673.0, 725.0, 780.0, 838.0, 900.0, 966.0, 1036.0, 1111.0,
1195.0, 1288.0, 1391.0, 1506.0, 1634.0, 1776.0, 1934.0, 2110.0, 2306.0, 2525.0, 2769.0, 3040.0, 3342.0, 3678.0, 4052.0, 4469.0, 4933.0, 5449.0, 6023.0, 6662.0, 7373.0, 8164.0, 9045.0, 10026.0, 11118.0, 12333.0, 13685.0, 15190.0, 16865.0, 18729.0, 20804.0,
23114.0, 25685.0, 28546.0, 31730.0]<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thanks for the help,<u></u><u></u></p>
<p class="MsoNormal">Matt <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">def map_vslice(pname,plabel,<wbr>cvarone,cvartwo,u,v,xv,xl,yv,<wbr>yl,yd):<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> wks_type = "png"<u></u><u></u></p>
<p class="MsoNormal"> wks = Ngl.open_wks(wks_type, pname)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> #cmap = Ngl.read_colormap_file("<wbr>WhiteBlueGreenYellowRed")<u></u><u></u></p>
<p class="MsoNormal"> #cmap = Ngl.read_colormap_file("MPL_<wbr>PuOr")<u></u><u></u></p>
<p class="MsoNormal"> #cmap = Ngl.read_colormap_file("MPL_<wbr>PuBu")<u></u><u></u></p>
<p class="MsoNormal"> cmap = Ngl.read_colormap_file("MPL_<wbr>Greens")<u></u><u></u></p>
<p class="MsoNormal"> #cmap = cmap[::-1,:]<u></u><u></u></p>
<p class="MsoNormal"> #cmap = Ngl.read_colormap_file("<wbr>cmocean_tempo")<u></u><u></u></p>
<p class="MsoNormal"> #cmap = Ngl.read_colormap_file("<wbr>WhiteYellowOrangeRed")<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> cfres = Ngl.Resources()<u></u><u></u></p>
<p class="MsoNormal"> cfres.nglFrame = False<u></u><u></u></p>
<p class="MsoNormal"> cfres.nglDraw = False<u></u><u></u></p>
<p class="MsoNormal"> clres = Ngl.Resources()<u></u><u></u></p>
<p class="MsoNormal"> clres.nglFrame = False<u></u><u></u></p>
<p class="MsoNormal"> clres.nglDraw = False<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> ylimit = max(yl) #)-12<u></u><u></u></p>
<p class="MsoNormal"> cfres.sfXArray = xv<u></u><u></u></p>
<p class="MsoNormal"> cfres.sfYArray = yl<u></u><u></u></p>
<p class="MsoNormal"> clres.sfXArray = xv<u></u><u></u></p>
<p class="MsoNormal"> clres.sfYArray = yl<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> #For Y axis<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLMode = "Explicit"<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLLabelFont = "times-roman"<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLValues = yl[::5]<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLLabels = yl[::5]<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLMinorOn = True<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLMinorValues = yl<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> #For X axis<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmXBMode = "Explicit"<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmXBLabelFont = "times-roman"<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmXBValues = xv[::20]<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmXBLabels = xl[::20]<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmXBMinorOn = True<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmXBMinorValues = xv[::5]<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> cfres.tmXBLabelFontHeightF = 0.006<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLLabelFontHeightF = 0.006<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> clres.cnFillOn = False<u></u><u></u></p>
<p class="MsoNormal"> clres.cnLinesOn = True<u></u><u></u></p>
<p class="MsoNormal"> clres.cnLineColor = "Black"<u></u><u></u></p>
<p class="MsoNormal"> clres.cnLineThicknessF = 3<u></u><u></u></p>
<p class="MsoNormal"> clres.cnLineLabelsOn = True<u></u><u></u></p>
<p class="MsoNormal"> clres.cnLineLabelFontHeightF= 0.007<u></u><u></u></p>
<p class="MsoNormal"> clres.cnLineLabelPerimOn = False<u></u><u></u></p>
<p class="MsoNormal"> clres.<wbr>cnLineLabelBackgroundColor = -1<u></u><u></u></p>
<p class="MsoNormal"> clres.cnInfoLabelOn = False<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> cfres.tmXBLabelFontHeightF = 0.007<u></u><u></u></p>
<p class="MsoNormal"> cfres.tmYLLabelFontHeightF = 0.007<u></u><u></u></p>
<p class="MsoNormal"> cfres.tiYAxisString = "Sigma-Height/m"<u></u><u></u></p>
<p class="MsoNormal"> cfres.tiXAxisString = "Latitude/deg"<u></u><u></u></p>
<p class="MsoNormal"> cfres.tiXAxisFontHeightF = 0.012<u></u><u></u></p>
<p class="MsoNormal"> cfres.tiYAxisFontHeightF = 0.012<u></u><u></u></p>
<p class="MsoNormal"> cfres.tiMainOn = True<u></u><u></u></p>
<p class="MsoNormal"> cfres.tiMainFontHeightF = 0.013<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> cfres.cnFillOn = True<u></u><u></u></p>
<p class="MsoNormal"> cfres.cnLinesOn = False<u></u><u></u></p>
<p class="MsoNormal"> cfres.cnLineLabelsOn = False<u></u><u></u></p>
<p class="MsoNormal"> cfres.cnFillPalette = cmap<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> cfres.lbAutoManage = True #False<u></u><u></u></p>
<p class="MsoNormal"> cfres.lbLabelStride = 1<u></u><u></u></p>
<p class="MsoNormal"> cfres.pmLabelBarWidthF = 0.02<u></u><u></u></p>
<p class="MsoNormal"> cfres.lbLabelFontHeightF = 0.016<u></u><u></u></p>
<p class="MsoNormal"> cfres.lbPerimThicknessF = 2<u></u><u></u></p>
<p class="MsoNormal"> cfres.pmLabelBarOrthogonalPosF = 0.002<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> plot = Ngl.contour(wks,cvarone,cfres)<u></u><u></u></p>
<p class="MsoNormal"> oplot= Ngl.contour(wks,cvartwo,clres)<u></u><u></u></p>
<p class="MsoNormal"> Ngl.overlay(plot,oplot)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"> Ngl.draw(plot)<u></u><u></u></p>
<p class="MsoNormal"> Ngl.frame(wks)<u></u><u></u></p>
<p class="MsoNormal"> Ngl.destroy(wks)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<br>______________________________<wbr>_________________<br>
pyngl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/pyngl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/pyngl-talk</a><br>
<br></blockquote></div><br></div>