; File: tri_save_ps.pro
; Author: Erik Brisson


N = 40
d = fltarr(3,N)

openr, 1, 'dat/ex_tri.dat'
readf, 1, d
close, 1

x = d(0,*)
y = d(1,*)
z = d(2,*)

triangulate, x, y, tri
set_plot,'PS'
device, filename='tri_save1.ps'
plot, x, y, psym=1, title="'Triangulation of regular data'
for i=0, n_elements(tri)/3-1 do begin
    t = [tri(*,i), tri(0,i)]
    plots, x(t), y(t)
endfor
device, /close

device, filename='tri_save2.ps'
surface, trigrid(x,y,z,tri)
device, /close

set_plot, 'X'

end