<html>
  <head>
    <meta content="text/html; charset=utf-8" http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    This issue has been solved.  I had NCL v6.1.2 still in the PATH
    environment variable.<br>
    <br>
    Thank you,<br>
    <br>
    Scott<br>
    <div class="moz-cite-prefix">On 11/11/2014 03:10 PM, Rick Brownrigg
      wrote:<br>
    </div>
    <blockquote
cite="mid:CAGKRhbFvFk3WdEYYCVZSJxV3EFaO+HnioriUa0q1nbhpAZGS=A@mail.gmail.com"
      type="cite">
      <div dir="ltr">
        <div>
          <div>
            <div>
              <div>Hi,<br>
                <br>
              </div>
              Is there a typo in this line:<br>
              <br>
                 plres@gsSegments  = f-&gt;segments ;(:,0)<br>
              <br>
            </div>
            I don't think the semicolon should be there.  f-&gt;segments
            should be an N x 2 array.<br>
            <br>
          </div>
          HTH...<br>
        </div>
        Rick<br>
        <br>
      </div>
      <div class="gmail_extra"><br>
        <div class="gmail_quote">On Tue, Nov 11, 2014 at 4:00 PM, Scott
          Capps <span dir="ltr">&lt;<a moz-do-not-send="true"
              href="mailto:scott@allvertum.com" target="_blank">scott@allvertum.com</a>&gt;</span>
          wrote:<br>
          <blockquote class="gmail_quote" style="margin:0 0 0
            .8ex;border-left:1px #ccc solid;padding-left:1ex">Greetings,<br>
            <br>
            I have just upgraded to version 6.2.1
            (ncl_ncarg-6.2.1.Linux_Debian7.6_x86_64_nodap_gcc472.tar.gz)
            and am happy with the many enhancements.  I am especially
            interested in the more efficient shapefile plotting
            abilities.  However, I am getting strange output when trying
            the new shapefile methods.  Attached is a ps file generated
            using the new method (Which has an issue).  I can provide
            the shapefile (Major U.S. highways) as well.<br>
            <br>
                 ; THIS IS THE NEW METHOD WHICH PRODUCES STRANGE OUTPUT:<br>
                 plres             = True           ; resources for
            polylines<br>
                 plres@gsLineColor = "gray25"<br>
                 plres@gsLineThicknessF     = 3<br>
                 f = addfile(ncl_dir+"shapefiles/intrastate_shapefiles/intrstat.shp",
            "r") ; Open shapefile<br>
                 plres@gsSegments  = f-&gt;segments ;(:,0)<br>
                 poly0   = gsn_add_polyline(wks, plot, f-&gt;x, f-&gt;y,
            plres)<br>
            <br>
                 ;<br>
                 ; THIS OLD METHOD WORKS BUT IS VERY SLOW:<br>
                 ; f = addfile(ncl_dir+"shapefiles/intrastate_shapefiles/intrstat.shp",
            "r") ; Open shapefile<br>
                 ; segments = f-&gt;segments<br>
                 ; geometry = f-&gt;geometry<br>
                 ; segsDims = dimsizes(segments)<br>
                 ; geomDims = dimsizes(geometry)<br>
                 ; geom_segIndex = f@geom_segIndex<br>
                 ; geom_numSegs  = f@geom_numSegs<br>
                 ; segs_xyzIndex = f@segs_xyzIndex<br>
                 ; segs_numPnts  = f@segs_numPnts<br>
                 ; numFeatures   = geomDims(0)<br>
                 ; lines = new(segsDims(0),graphic)   ; array to hold
            polylines<br>
            <br>
                 ; plres             = True           ; resources for
            polylines<br>
                 ; plres@gsLineColor = "gray25"<br>
                 ; plres@gsLineThicknessF     = 3<br>
            <br>
                 ; lon    = f-&gt;x<br>
                 ; lat    = f-&gt;y<br>
                 ; segNum = 0       ; Counter for adding polylines<br>
                 ; do i=0, numFeatures-1<br>
                 ;   startSegment = geometry(i, geom_segIndex)<br>
                 ;   numSegments  = geometry(i, geom_numSegs)<br>
                 ;   do seg=startSegment, startSegment+numSegments-1<br>
                 ;     startPT = segments(seg, segs_xyzIndex)<br>
                 ;     endPT   = startPT + segments(seg, segs_numPnts) -
            1<br>
                 ;     lines(segNum) = gsn_add_polyline(wks, plot,
            lon(startPT:endPT),  \<br>
                 ;     lat(startPT:endPT), plres)<br>
                 ;     segNum = segNum + 1<br>
                 ;   end do<br>
                 ; end do<br>
                 ;<br>
            <br>
            The full NCL script is attached.  I can provide sample data
            upon request (<a moz-do-not-send="true"
              href="http://invariant_d01.nc" target="_blank">invariant_d01.nc</a>,
            <a moz-do-not-send="true" href="http://wrftest.nc"
              target="_blank">wrftest.nc</a> and the shapefile).<br>
            <br>
            Thank you,<br>
            <br>
            Scott<br>
            <br>
            <br>
            _______________________________________________<br>
            ncl-talk mailing list<br>
            List instructions, subscriber options, unsubscribe:<br>
            <a moz-do-not-send="true"
              href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk"
              target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
            <br>
          </blockquote>
        </div>
        <br>
      </div>
    </blockquote>
    <br>
  </body>
</html>