<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Mark,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">[I'm posting back to ncl-talk here.]</div><div class="gmail_default" style="font-size:small"><br>Thanks for providing the data file. I was able to reproduce the problem, and unfortunately I'm not sure about a work-around at this time.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">That's not entirely true: if I don't reverse the Y axis, then the bug goes away. But, with a pressure/height plot, you want to reverse the Y axis. I tried a few things, like reversing the leftmost dimension of your data (that didn't work), and reversing the Y axis myself and then just using gsn_csm_contour (that didn't work either because then a different contour showed up white).</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I have filed a bug report using your small data file and a modified version of your script. The bug report is NCL-2478 for your reference. I'm sorry I couldn't find a work-around. Hopefully this isn't a difficult bug to fix.</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 class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 19, 2016 at 10:49 AM, Mark Branson <span dir="ltr"><<a href="mailto:mark@atmos.colostate.edu" target="_blank">mark@atmos.colostate.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Hi Mary.</div><div><br></div>Wow, what a cool debugging tool. Attached is the resulting <a href="http://mary_debug.nc" target="_blank">mary_debug.nc</a> file. If need be, I can definitely upload the original data file (<a href="http://uandps.nc" target="_blank">uandps.nc</a>). That file size is around 1GB.<div><br></div><div>Thanks!</div><span class="HOEnZb"><font color="#888888"><div>Mark</div><div><br></div><div></div></font></span></div><br><div style="word-wrap:break-word"><div></div><div><br><div><blockquote type="cite"><div>On Aug 19, 2016, at 10:41 AM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:</div><br><div><div dir="ltr"><div class="gmail_default" style="font-size:small">[Offline for a bit.]</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Yeah, this must point to the coordinate somehow, but this is a little odd.</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 the data? I think if you set this unadvertised resource before you call gsn_csm_contour:</div><div class="gmail_default" style="font-size:small"><p><span>res@gsnDebugWriteFileName = "mary_debug"</span></p><p>Then it should create a "<a href="http://mary_debug.nc/" target="_blank">mary_debug.nc</a>" NetCDF file, with just the "uwnd" data you're plotting, along with the coordinate arrays. This way you don't have to provide me with the full data file if you don't want.</p><p>The "mary_debug,nc" file should probably be small enough to just attach in an email. I think a "mary_debug.ncl" script will also get created, and you can email that to me, although I don't need it since I have your original script.</p><p>Thanks,</p><p>--Mary</p><p><br></p></div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 19, 2016 at 10:34 AM, Mark Branson <span dir="ltr"><<a href="mailto:mark@atmos.colostate.edu" target="_blank">mark@atmos.colostate.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Hi Mary.</div><div><br></div>Here is that. Interestingly, the problem has gone away. So perhaps the coordinate arrays are indeed the source of the problem?<span><font color="#888888"><div><br></div><div>Mark</div><div><br></div><div></div></font></span></div><br><div style="word-wrap:break-word"><div></div><div><br><div><blockquote type="cite"><div>On Aug 19, 2016, at 10:26 AM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:</div><br><div><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Mark,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">That is a little odd. You did good work debugging the problem, so missing values and nans are ruled out.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">To further debug this, can you try plotting a very basic contour plot, without the pressure /height axes, and without any coordinate information and send the image:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px"> plot = gsn_contour (wks,uwind,res) </span><br style="font-size:12.8px"></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px"><br></span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px">The gsn_contour script ignores any coordinate arrays, which I want to make sure are not the source of the problem.</span></div><div class="gmail_default" style="font-size:small"><span style="font-size:12.8px"><br></span></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><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Aug 19, 2016 at 10:07 AM, Mark Branson <span dir="ltr"><<a href="mailto:mark@atmos.colostate.edu" target="_blank">mark@atmos.colostate.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hello.<br><br>I'm running ncl v6.2.1 and here is what I'm doing in the attached script:<br><br>1) read in a long time series of zonal wind and surface pressure data from CAM output.<br>2) interpolate them to constant pressure levels using vinth2p.<br>3) compute long-term monthly means with clmMonTLLL.<br>4) compute the annual mean of the monthly means.<br>5) subtract off the annual mean from the monthly means.<br>6) select one latitude (the equator), and then take the mean of all longitudes.<br>6) make a color-filled height versus time cross section plot.<br><br>I get the plot (see attached), but some of the contours near the top are white. I thought maybe it was because there were some missing or bogus values, so I inserted a check for those with the ismissing and isnan_ieee functions and it seems the data is fine.<br><br>Any ideas?<br><br>Thanks,<br>Mark Branson<br><br><br>begin<br> a = addfile("/pool/mark/spcam_nlev<wbr>51/<a href="http://spcam_nlev51.cam.h0.0003-01.nc/" target="_blank">spcam_nlev51.cam.h0.0003-01<wbr>.nc</a>","r")<br> hyam = a->hyam<br> hybm = a->hybm<br><br> b = addfile("/pool/mark/spcam_nlev<wbr>51/<a href="http://uandps.nc/" target="_blank">uandps.nc</a>","r")<br> u = b->U ;u = [time | 338] x [lev | 53] x [lat | 96] x [lon | 144]<br> ps = b->PS ;ps = [time | 338] x [lat | 96] x [lon | 144]<br> u@_FillValue = default_fillvalue("float")<br><br> if (any(ismissing(u))) then<br> print("### original u array contains some missing data ###")<br> end if<br><br>; interpolation variables<br> plvlM = (/ 30.,50.,70.,100.,150.,200.,250<wbr>.,300.,350.,400.,450.,500.,550<wbr>., \<br> 600.,650.,700.,750.<wbr>,800.,850./)<br> plvlM@units = "mb"<br> nplvlM = dimsizes(plvlM)<br><br> p0 = 1000. ; mb required by vinth2p<br> interp = 2 ; log interpolation<br> extrap = False ; no extrapolation past psfc.<br><br> tmp = vinth2p(u,hyam,hybm,plvlM,ps,i<wbr>nterp,p0,1,extrap)<br> if (typeof(tmp).eq."double") then<br> uwind1 = dble2flt(tmp)<br> else<br> uwind1 = tmp<br> end if<br><br>; compute long-term monthly means<br> uclm = clmMonTLLL(uwind1) ;uclm = [month | 12] x [lev_p | 33] x [lat | 96] x [lon | 144]<br><br>; annual mean of monthly means<br> uann = dim_avg_n_Wrap(uclm,0) ;uann = [lev_p | 33] x [lat | 96] x [lon | 144]<br><br>; subtract off annual mean from monthly means<br> do i = 0,11<br> uclm(i,:,:,:) = (/ uclm(i,:,:,:) - uann(:,:,:) /)<br> end do<br><br> uwind2 = uclm(:,:,47,:) ; select the equator for the latitude<br> uwind3 = dim_avg_Wrap(uwind2) ; mean of all longitudes<br> uwind = uwind3(lev_p|:,month|:) ; swap the array dimensions: month,level --> level,month<br><br> printMinMax(uwind,False)<br><br> if (any(ismissing(uwind))) then<br> print("### uwind array contains some missing data ###")<br> end if<br> if (any(isnan_ieee(uwind))) then<br> print("### uwind array contains some NaNs ###")<br> end if<br><br> uwind@long_name = "SP-CAM, 51 levels, Monthly Means"<br><br> wks = gsn_open_wks ("eps", "test1") ; send graphics to EPS file<br> gsn_define_colormap(wks,"BlRe<wbr>")<br><br> res = True <br> res@cnFillOn = True <br> res@cnLinesOn = True <br> res@lbLabelBarOn = True <br> res@cnInfoLabelOn = False<br><br>; change labels for monthly mean plot<br> res@tmXBMode = "Explicit"<br> res@tmXBLabelFontHeightF = 0.015<br> res@tmXBLabels=(/"Jan","Feb",<wbr>"Mar","Apr","May","Jun","Jul",<wbr>"Aug","Sep", \<br> "Oct","Nov",<wbr>"Dec","Jan"/)<br> res@tmXBValues = ispan(0,12,1)<br><br> plot = gsn_csm_pres_hgt (wks,uwind,res) <br>end<br><br>** OUTPUT FROM THE SCRIPT **<br><br>/Users/mark/mmf/SPCAM5/stratos<wbr>phere_gwdrag/spcam_nlev51> ncl hvstime.ncl<br>Copyright (C) 1995-2014 - All Rights Reserved<br>University Corporation for Atmospheric Research<br>NCAR Command Language Version 6.2.1<br>The use of this software is governed by a License Agreement.<br>See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> f<wbr>or more details.<br>(0)<span style="white-space:pre-wrap"> </span>min=-9.835 max=6.77056<br><br></div><br><div style="word-wrap:break-word"></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
______________________________<wbr>_________________<br>ncl-talk mailing list<br><a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>List instructions, subscriber options, unsubscribe:<br><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br></div></blockquote></div><br></div></div><br></blockquote></div><br></div>
</div></blockquote></div><br></div></div><br></blockquote></div><br></div>