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

[Rick Brownrigg]

Hi Michael,

I'm not very familiar with the TDPACK functions, but do you need to call
tdinit() to set the viewing parameters?

http://ncl.ucar.edu/Document/Functions/Built-in/tdinit.shtml

Rick


On Tue, Sep 26, 2017 at 9:16 AM, Michael Buban - NOAA Affiliate <
michael.buban at noaa.gov> wrote:

> 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
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170926/728d5a7d/attachment.html>