PRO gfit_chi2_tour
; 12/5/94 BIG kludge
@labx_common
@gfit_common

m=11
n=11
a_range = 900.0 + findgen(m)*70.0
t_range = 4000.0 + findgen(n)*300.0

print, 'Total chi^2 evaluations for contour:',m,' times ',n

chi2_tour = fltarr(m,n)
for i = 0,m-1 do begin
  for j = 0,n-1 do begin
    afit = a_range(i)
    tfit = t_range(j)
    chi2_tour(i,j) = gfit_chi_2(tfit,afit)
  endfor
print, ' ...',i+1,' times',n,' done ...'
endfor

maxchi2 = max(chi2_tour)
minchi2 = min(chi2_tour)

rangechi2 = maxchi2 - minchi2
tit_str = 'Contours of 8 keV Level I Acceptability'

x_str = 'Grating Line-width (A)' & y_str = 'Grating Thickness (A)'
levels = FLOAT(INDGEN(25)+10)/20.

Contour,chi2_tour,a_range,t_range,LEVELS=levels,TITLE=tit_str,$
XTITLE=x_str,YTITLE=y_str,XRANGE=[min(a_range),max(a_range)],$
YRANGE=[min(t_range),max(t_range)],$
/XSTYLE,/YSTYLE,charsiz=1.2,C_Charsize=1.3

RETURN
END