;
; after running gain loss. try to make an image
;
ii=where(gIar.lambdacm eq 6,cnt)
xx=giar[ii].x0f
yy=giar[ii].y0f
xoutL=findgen(1000)*1 - 500.
youtL=xoutL
triangulate,xx,yy,tr
;
imgloss=griddata(xx,yy,giar[ii].gainloss,xout=xoutL,yout=youtL,$
                     /inverse_distance,triangles=tr,/grid)
imglossdb=alog10(imgloss)*10.

iimage,imglossdb,xoutL,youtL
;

hor,-500,500
ver,-500,500
plot,xx,yy,psym=2
;
; try gridding in x then y
;
nstep=5
nn=100
xN=findgen(nn)/nn * 1000. - 500
xx=fltarr(nx,nx)
yy=fltarr(ny,ny)
ll=fltarr(ny,ny)
ii=where(giar.lambdacm eq 6)
g=giar[ii]
;
; first interpolate across x for each y
;
nstep=5
y1=fltarr(nn,nstep)
g1=fltarr(nn,nstep)
y0=-400.
step=200.
for i=0,nstep-1 do begin &$
	y=y0 + i*step &$
	ii=where(g.y0f eq y,cnt) &$
	x1[*,i]=xn &$
	y1[*,i]=y  &$
	g1[*,i]=interpol(g[ii].gainloss,g[ii].x0f,xN,/quad) &$
endfor
;
; now interpolate y at each of the new x
;:q

xx=fltarr(nn,nn)
yy=fltarr(nn,nn)
gl=fltarr(nn,nn)
for i=0,nn-1 do begin &$
	xx[*,i]=xn &$
	yy[i,*]=xn &$
	gl[i,*]=interpol(g1[i,*],y1[i,*],xN,/quad) &$ &$
endfor
r=sqrt(xx^2  + yy^2)
ii=where(r lt 500)
xoutL=(findgen(100)/100 - .5 )*1000.
youtL=xoutL
searchellipse=20.
missingdata=0.
;
triangulate,xx[ii],yy[ii],tr
imgloss=griddata(xx[ii],yy[ii],gl[ii],xout=xoutL,yout=youtL,$
                     /inverse_distance,triangles=tr,/grid,$
					missing=missingdata,search_Ellipse=searchellipse)
ii=where(imgloss ne 0.)
imglossdb=imgloss
imglossdb[ii]=alog10(imgloss[ii])*10.
minv=min(imglossdb[ii])
ii=where(imgloss eq 0.)
eps=.5
imglossdb[ii]=minv -eps
iimage,imglossdb,xoutL,youtL