function figure_pos, x, y
;
;
;
xarr = [31, 52, 31, 10, 31]
yarr = [31, 31, 11, 31, 51]
r = sqrt( (xarr-x)^2 + (yarr-y)^2)
dummy = min(r, loc)
;
return, loc
end
;------------------------------------------------------------------------------
;;pro plot_mmap_trans, infil
;
;HISTORY:
;	Written Feb-94 by M.Morrison
;V2.0	21-Mar-94 (MDM)
;V3.0	 6-May-94 (MDM) -
;V3.1	 6-May-94 (MDM) - Correct for laser intensity variation
;V3.2	 6-May-94 (MDM) - Adjusted plotting range
;
;
if (n_elements(infil) eq 0) then begin
    ;infil = '$MDI_CAL_TXT/940312_025326-940312_052339.tr_summary'
    infil = '$MDI_CAL_TXT/tlm940408_222903-tlm940409_183439.tr_summary
    input, 'Enter the transmission summary file', infil, infil
end
;
mat = rd_tfile(infil, 13, 6, /conv)
;
if (n_elements(last_code) eq 0) then last_code = 0
print, 'Enter 0 for left lobe
print, '      1 for center
print, '      2 for right lobe'
print, '     -1 for center (940312 data)
input, 'Enter code', code, last_code
last_code = code
case code of
   -1: begin & ss0 = [0,5,10]			& npos=5	& end
    0: begin & ss0 = [34,35,36]			& npos=1	& end	;ss0 is the pos 0  subscript for a given wavelength
    1: begin & ss0 = [11, 16, 26]		& npos=5	& end
    2: begin & ss0 = [31,32,33]			& npos=1	& end	;long wavelength
endcase
;
ss = [1,0,2]	;indgen(3)
;
wave = mat(3,*)
avg_tr = mat(6, *)
laspow = mat(4,*)
norm_tr = avg_tr / float(laspow>1)
xlas = mat(10, *)
ylas = mat(11, *)
;
n = n_elements(wave)
posarr = intarr(n)
for i=0,n-1 do posarr(i) = figure_pos(xlas(i), ylas(i))
;
;for qplot_fit = 0,1 do begin
for qplot_fit = 1,1 do begin
    ;plot, wave, avg_tr, psym=3, tit=infil, xtit='Wavelength (A)', ytit='Transmission'
    for i=0,npos-1 do begin
	;x = wave(ss*5+i)
	;y = avg_tr(ss*5+i)

	x = wave(ss0+i)
	y = norm_tr(ss0+i)
	x = x(ss)	;reorder
	y = y(ss)
	if (i eq 0) then plot, x, y, psym=3, tit=infil, xtit='Wavelength (A)', ytit='Transmission / Laser Power', $
			xrange=[min(wave(ss0))-.2, max(wave(ss0))+.2], yrange=[0, max(norm_tr(ss0))*1.1]
	dx = (!x.crange(1)-!x.crange(0))
	dy = (!y.crange(1)-!y.crange(0))

	if (qplot_fit eq 0) then oplot, x, y, psym=-(i+1), linestyle=i $
				else oplot, x, y, psym=(i+1)
	str = 'Pos ' + strtrim(i,2)
	;
	if (qplot_fit) then begin
	    coeff = poly_fit(x-6767, y, 2, yfit)
	    xx = findgen(500)/499*dx + !x.crange(0) - 6767
	    yy = poly(xx, coeff)
	    oplot, xx+6767, yy, linestyle=i
	    xx2 = findgen(3000)/2999*3*dx + !x.crange(0)-dx - 6767
	    yy2 = poly(xx2, coeff)
	    amp = max(yy2)
	    ss3 = where(yy2 ge amp/2, nn3)
	    fwhm = xx2(ss3(nn3-1)) - xx2(ss3(0))
	    str = str + ' Center=' + strtrim( -coeff(1)/(2*coeff(2)) + 6767, 2) + $
				' Amp=' + strtrim(amp,2) + ' FWHM=' + strtrim(fwhm,2)
        end

	x0 = dx*.3 + !x.crange(0)
	y0 = dy*.3 + !y.crange(0)
	oplot, [x0], [y0-i*dy*.04], psym=-(i+1)
	xyouts, x0+dx*.05, y0-i*dy*.04, str
    end
    ;
    plottime, 0, 0, 'PLOT_MMAP_TRANS  Ver 3.2'
    plottime
    pause
end
;
end