<div dir="ltr"><div class="gmail_default" style="font-size:small">​Hi Laura,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Would you be able to provide a clean version of your script​ and the data, as per usual? </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">It seems like what you are doing is correct, and for the record, there is no special weighting that goes on when the zonal means plot is created.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I was going to try to create an example for you, but your data sounds unique enough (because it's MPAS and on multiple levels) that it would be better if I could have your data.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks,</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 class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, May 14, 2018 at 5:40 PM, Laura Fowler <span dir="ltr"><<a href="mailto:laura@ucar.edu" target="_blank">laura@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello:<br>
<br>
I would like to plot latitude versus pressure cross-sections from my<br>
MPAS outputs (the resolution of my MPAS mesh is about 30 km). The<br>
first thing that I did was to use the function<br>
ESMF_regrid_with_weights to regrid each output  to a 0.25 deg x 0.25<br>
deg horizontal resolution (720x1440) for each vertical level. Then, I<br>
used the function dim_avg_n_Wrap to calculate the zonal mean prior to<br>
using gsn_csm_pres_hgt.<br>
<br>
The plot "ZonalMeanTemperature.pdf" show the latitude versus pressure<br>
cross-sections which looks really boxy. In contrast, a 500 hPa map<br>
does not show the same boxiness with res@gsnZonalMean = True (see<br>
Temperature_500hPa.pdf). Other levels do not show any "boxiness"<br>
either.<br>
<br>
I assume that the "boxiness" in ZonalMeanTemperature.pdf results<br>
because I do not weigh the zonal means as a function of the area, but<br>
not sure, since I simply use dim_avg_n_Wrap.<br>
<br>
How can I reduce the boxiness of the pressure versus latitude<br>
cross-sections. Is there a "correction" of the zonal mean built in the<br>
function gsnZonalMean? Is there an ncl function that I can apply to<br>
the zonal mean to help reduce the step change in the temperature as a<br>
function of the latitude? Or is there an other gsn_csm_pres_hgt that<br>
includes a weighting function?<br>
<br>
<br>
<br>
Thanks,<br>
Laura<br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
<br>
<br>
<br>
<br>
-- <br>
!-----------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------------<br>
Laura D. Fowler<br>
Mesoscale and Microscale Meteorology Division (MMM)<br>
National Center for Atmospheric Research<br>
P.O. Box 3000, Boulder CO 80307-3000<br>
<br>
e-mail: <a href="mailto:laura@ucar.edu">laura@ucar.edu</a><br>
phone: 303-497-1628<br>
<br>
!-----------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------------<br>
</font></span><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>