<div dir="ltr"><div class="gmail_default" style="font-size:small">Judy,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm so sorry I let this slip. </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The error is that you have a min latitude that is greater than the max latitude:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div style="font-size:12.8000001907349px"><font face="monospace, monospace">res.mpMaxLatF = -79.</font></div><div style="font-size:12.8000001907349px"><font face="monospace, monospace">res.mpMinLatF = -50.</font></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif"><br></font></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">I modified the script to remove some of the code that processes the data read off the file, because this is an unnecessary step. Also, because the lat/lon arrays have missing values in them, you need to plot this using the "triangular mesh" algorithm under the hood.</font></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif"><br></font></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">Please see the modified script and resultant PNG image which I attached.</font></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif"><br></font></span></div><div style="font-size:12.8000001907349px"><span style="font-size:12.8000001907349px"><font face="arial, helvetica, sans-serif">--Mary </font></span><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 1, 2015 at 10:31 AM, Judy Twedt <span dir="ltr"><<a href="mailto:judytwedt@gmail.com" target="_blank">judytwedt@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Mary, <div> Thanks for letting me know. I appreciate any assistance you can provide. I'm using PyNgl for all of my map plots now, but have been regridding ocean data as a short-term work around. </div><div><br></div><div>best,</div><div><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 1, 2015 at 9:28 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@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 dir="ltr"><div class="gmail_default" style="font-size:small">Hi Judy,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm sorry I haven't gotten back to this. I'm just swamped right now. I will try to look later today.</div><span><font color="#888888"><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></font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Aug 27, 2015 at 12:29 PM, Judy Twedt <span dir="ltr"><<a href="mailto:judytwedt@gmail.com" target="_blank">judytwedt@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Mary and other PyNgl users,<div> My name is Judy Twedt. I've recently started using PyNgl on yellowstone and have persistent problems plotting scalar fields from the popgrid using the Ngl.contour_map function. The lat/lon arrays are 2D, which I have raveled to 1D, but I get the following error message and a plot with unphysical lines which run off of the projection:</div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">warning:MapTransInitialize: latitude min exceeds max: exchanging values</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="arial, helvetica, sans-serif">I made a short script to reproduce the problem which can be found and run at: /glade/p/work/jtwedt/LEdata/simplePOPcontour.py</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">This example plots SH ice area -- I also tried it with ocean temperature (which spans both hemispheres) had have the same error. I'd very much appreciate your help locating the problem here. </font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Thanks in advance for your advice!</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif"> The script reads: </font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><div><font face="monospace, monospace">import numpy as np</font></div><div><font face="monospace, monospace">from netCDF4 import Dataset</font></div><div><font face="monospace, monospace">import <a href="http://numpy.ma" target="_blank">numpy.ma</a> as ma</font></div><div><font face="monospace, monospace">import pdb</font></div><div><span style="font-family:monospace,monospace">import Ngl </span><br></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">#load ice lat,lon arrays</font></div><div><font face="monospace, monospace">ds = Dataset('/glade/p/cesmLE/CESM-CAM5-BGC-LE/ice/proc/tseries/monthly/aice/b.e11.B20TRC5CNBDRD.f09_g16.032. <a href="http://cice.h.aice_sh.192001-200512.nc" target="_blank">cice.h.aice_sh.192001-200512.nc</a>')</font></div><div><font face="monospace, monospace">lon1 = ds.variables['TLON'][:,:]</font></div><div><font face="monospace, monospace">lon = np.ravel(np.ma.filled(lon1, np.nan)) </font></div><div><font face="monospace, monospace">lat1 = ds.variables['TLAT'][:,:]</font></div><div><font face="monospace, monospace">lat = np.ravel(np.ma.filled(lat1, np.nan))</font></div><div><font face="monospace, monospace">ice1= ds.variables['aice'][1,:,:]</font></div><div><font face="monospace, monospace">ice = np.ravel(ma.filled(ice1, np.nan))</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">#pdb.set_trace()</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">wks = Ngl.open_wks('X11', 'scalarPOP')</font></div><div><font face="monospace, monospace">res = Ngl.Resources()</font></div><div><font face="monospace, monospace">res.sfXArray = lon </font></div><div><font face="monospace, monospace">res.sfYArray = lat </font></div><div><font face="monospace, monospace">res.cnFillOn = True</font></div><div><font face="monospace, monospace">res.cnLinesOn = False</font></div><div><font face="monospace, monospace">res.mpProjection = 'Stereographic'</font></div><div><font face="monospace, monospace">res.mpLimitMode = 'LatLon'</font></div><div><font face="monospace, monospace">res.mpCenterLatF = -79.</font></div><div><font face="monospace, monospace">res.mpMaxLatF = -79.</font></div><div><font face="monospace, monospace">res.mpMinLatF = -50.</font></div><div><font face="monospace, monospace"><div>Ngl.contour_map(wks, ice, res)</div><div>Ngl.end()</div><span><font color="#888888"><div><br></div></font></span></font></div><span><font color="#888888"><div><br></div>-- <br><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><img src="http://upload.wikimedia.org/wikipedia/commons/thumb/6/68/Gay_flag.svg/2000px-Gay_flag.svg.png" width="96" height="58"><br></div>Judy Twedt<div><br><div><i>Graduate Student </i></div><div><i style="font-size:12.8000001907349px">Department of Atmospheric Sciences</i><br></div><div><div><div><i>University of Washington</i></div></div></div></div></div></div></div></div></div></div>
</font></span></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><img src="http://upload.wikimedia.org/wikipedia/commons/thumb/6/68/Gay_flag.svg/2000px-Gay_flag.svg.png" width="96" height="58"><br></div>Judy Twedt<div><br><div><i>Graduate Student </i></div><div><i style="font-size:12.8000001907349px">Department of Atmospheric Sciences</i><br></div><div><div><div><i>University of Washington</i></div></div></div></div></div></div></div></div></div></div>
</div>
</div></div></blockquote></div><br></div>