do_fit_gcross, '/umbradata2/optest/pdist616.', idark=80, images=indgen(40)*2+1, /qhr, poutfil='B'

 for i=0,5 do do_fit_gcross, spec_code=i
end