<div dir="ltr"><div>Hello,</div><div><br></div><div>I'm trying to plot 3-D points on an isosurface, however I can't seem to get them to show up, either with or without the isosurface being plotted.  I've attach (I think) then relevant portion of the code below, and the full ncl script.  Any thoughts?</div><div><br></div><div>Thanks,</div><div><br></div><div>Mike</div><div><br></div><br>   px =  new(11, float)<br>   py =  new(11, float)<br>   pz =  new(11, float)<br><br>   mtri = 15000<br>   rtri = new((/mtri,10/), float)<br><br>   do i=0,10<br>   px(i)= 1650 +10.1*i<br>   py(i)= 100 +10.1*i<br>   pz(i)= 30 + 2.1*i<br>   end do<br><br>  wks  = gsn_open_wks("X11","gsn_contour")    ; open a ps file<br><br>  ntri = 0<br>  dotsize = 0.02<div><br>  do i = 0, 10<br>    tdmtri(-5, (/px(i), py(i), pz(i)/), dotsize, \<br>           rtri, ntri, 4, (/1649.,50.,30./),(/1799.,300.,90./))<br>  end do<br><br>  rtwk = new((/2,mtri/),float)</div><div><br>  itwk = tdotri(rtri, ntri, rtwk, 0)<br>  tddtri(wks,rtri, ntri, itwk)<br><br>  tdgrds(wks,(/1649., 50., 30./), (/1799, 300., 90./), (/10., 10., 10./),11,1)<br>   tdez3d(wks,xi(1648:1798),yi(50:300),zi(10:30),data2(1648:1798,50:300,10:30),.0003,2.,-20.,50.,-3)<br><br>  tdgrds(wks,(/1649., 50., 30./), (/1799, 300., 90./), (/10., 10., 10./),11,0)<div><br></div><div><br></div><div><br></div><div><br></div><div>-- <br></div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div style="font-size:12.8px">Michael S. Buban, Ph.D.</div><div style="font-size:12.8px">Research Associate</div><div style="font-size:12.8px">NOAA/ARL Atmospheric Turbulence and Diffusion Division</div><div style="font-size:12.8px">and Oak Ridge Associated Universities</div><div style="font-size:12.8px">456 South Illinois Avenue</div><div style="font-size:12.8px">Oak Ridge TN 37830</div></div></div></div></div>
</div></div>