<div dir="ltr"><div>Dear NCL users,<br></div>I have been trying to <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/csvoro.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/csvoro.shtml</a>:<br><pre>begin

;
;  Create input arrays.
;
  ndata = 100
  rlati = new(ndata,float)
  rloni = new(ndata,float)
;
  do i=0,ndata-1
    rlati(i) =  -90. + 180.*<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/rand.shtml"><strong>rand</strong></a>()/32767.
    rloni(i) = -180. + 360.*<a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/rand.shtml"><strong>rand</strong></a>()/32767.
  end do

;
;  Obtain the Voronoi polygon for the first point in the input
;  dataset.  nv(0:numv-1) will contain integers to be used as
;  indices for rlato and rlono to determine the polygonal vertices.
;
  rlato = new(2*ndata,float)
  rlono = new(2*ndata,float)
  alen  = new(2*ndata,float)
  nca   = new(1,integer)
  numv  = new(1,integer)
  nv    = new(ndata,integer)
  <strong>csvoro</strong>(rlati,rloni,0,1,rlato,rlono,alen,nca,numv,nv)

end<br><br></pre><pre>I run this sample many times, but always numv will be &lt;= 15 for ndata = 100 or ndata = 1000.<br>Values of alen will be positive only for first numv index&#39;s.<br><br></pre><pre>What&#39;s wrong?<br><br></pre><pre>Many thanks. <br></pre><br></div>