<div dir="ltr"><div><div>Hi Michael, <br><br></div>I'm not very familiar with the TDPACK functions, but do you need to call tdinit() to set the viewing parameters?<br><br><a href="http://ncl.ucar.edu/Document/Functions/Built-in/tdinit.shtml">http://ncl.ucar.edu/Document/Functions/Built-in/tdinit.shtml</a><br><br></div>Rick<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 26, 2017 at 9:16 AM, Michael Buban - NOAA Affiliate <span dir="ltr"><<a href="mailto:michael.buban@noaa.gov" target="_blank">michael.buban@noaa.gov</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>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_<wbr>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.<wbr>,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(<wbr>50:300),zi(10:30),data2(1648:<wbr>1798,50:300,10:30),.0003,2.,-<wbr>20.,50.,-3)<br><br>  tdgrds(wks,(/1649., 50., 30./), (/1799, 300., 90./), (/10., 10., 10./),11,0)<span class="HOEnZb"><font color="#888888"><div><br></div><div><br></div><div><br></div><div><br></div><div>-- <br></div><div class="m_4012325500900849826gmail_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>
</font></span></div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>