[ncl-talk] plotting 3-D points on an isosurface

[Michael Buban - NOAA Affiliate]

Hello,

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?

Thanks,

Mike


   px =  new(11, float)
   py =  new(11, float)
   pz =  new(11, float)

   mtri = 15000
   rtri = new((/mtri,10/), float)

   do i=0,10
   px(i)= 1650 +10.1*i
   py(i)= 100 +10.1*i
   pz(i)= 30 + 2.1*i
   end do

  wks  = gsn_open_wks("X11","gsn_contour")    ; open a ps file

  ntri = 0
  dotsize = 0.02

  do i = 0, 10
    tdmtri(-5, (/px(i), py(i), pz(i)/), dotsize, \
           rtri, ntri, 4, (/1649.,50.,30./),(/1799.,300.,90./))
  end do

  rtwk = new((/2,mtri/),float)

  itwk = tdotri(rtri, ntri, rtwk, 0)
  tddtri(wks,rtri, ntri, itwk)

  tdgrds(wks,(/1649., 50., 30./), (/1799, 300., 90./), (/10., 10.,
10./),11,1)

 tdez3d(wks,xi(1648:1798),yi(50:300),zi(10:30),data2(1648:1798,50:300,10:30),.0003,2.,-20.,50.,-3)

  tdgrds(wks,(/1649., 50., 30./), (/1799, 300., 90./), (/10., 10.,
10./),11,0)




-- 
Michael S. Buban, Ph.D.
Research Associate
NOAA/ARL Atmospheric Turbulence and Diffusion Division
and Oak Ridge Associated Universities
456 South Illinois Avenue
Oak Ridge TN 37830
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170926/4ac5c815/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: hrrr_maria.ncl
Type: application/octet-stream
Size: 7758 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170926/4ac5c815/attachment.obj>