; File: sphscan_intens.pro
; Author: Erik Brisson

ntheta = 100
nphi   = 100
nentries = ntheta*nphi

dat = fltarr(4,nentries)

openr, 1, '../data/xyz_intens.dat'
readf, 1, dat
close, 1

print, 'nentries = ', nentries

theta = dat(0,*)
phi = dat(1,*)
r = dat(2,*)
intens = dat(3,*)

triangulate, theta, phi, tri

x = r * sin(phi) * cos(theta)
y = r * sin(phi) * sin(theta)
z = r * cos(phi)
intens_norm = (intens - min(intens)) / (max(intens)-min(intens))
color_ind = (fix(256.0*intens_norm)) < 255

nvert = nentries
ntri = n_elements(tri)/3

print, "nvert = ", nvert
print, "ntri  = ", ntri

openw, 1, 'results/sph.xyzc'

printf, 1, nvert, ntri

for ivert=0, nvert-1 do begin
    printf, 1, x(ivert), y(ivert), z(ivert), color_ind(ivert)
end

for itri=0, ntri-1 do begin
    printf, 1, tri(1,itri), tri(0,itri), tri(2,itri)
endfor

close, 1

end