cols=[200,160,140,255] set_plot,'ps' device,xoffset=1,yoffset=23,filename='6733ind.ps' ,/times, $ ;device,xoffset=1,yoffset=23,filename='6733.65ind.ps' ,/times, $ /landscape, xsize=18,ysize=18 ;device,xoffset=1,yoffset=5,filename='figure7a.port.ps' ,/times, $ ;device,xoffset=1,yoffset=4,filename='figure7b.port.ps' ,/times, $ ;xsize=18,ysize=18 openr,1,'outfile.6733.ind' contlevs=[0,211,236,252] titstr='6733 '+string("305B)+' -22S' openr,3,'bestfit.22' figlab='(a)' ;openr,1,'outfile.6733.65.ind' ;contlevs=[0,193,218,233] ;titstr='6733 '+string("305B)+' -65S' ;openr,3,'bestfit.65' ;figlab='(b)' points=64 taucld=fltarr(points) taeros=fltarr(points) chisq=fltarr(points) sqrtpts=sqrt(points) x=fltarr(sqrtpts,sqrtpts) y=fltarr(sqrtpts,sqrtpts) z=fltarr(sqrtpts,sqrtpts) for i=0,points-1 do begin readf,1,a,b,c,d taucld(i)=a taeros(i)=b*20 chisq(i)=c endfor close,1 k=-1 for j=0,sqrtpts-1 do for i=0,sqrtpts-1 do begin k=k+1 x(i,j)=taucld(k) y(i,j)=taeros(k) z(i,j)=chisq(k) endfor contour,z,x,y,levels=contlevs,/closed,/fill,c_colors=cols, $ ystyle=1,yrange=[0,17],xtitle='!7s!3!Dcld!N',ytitle='f!Dhze!N', $ title=titstr,charsize=2.0 oplot,x,y,psym=6 print,min(chisq,h),h,max(chisq,h2),h2 print,'Best taeros= ',taeros(h) print,'Worst taeros= ',taeros(h2) readf,3,a,b plots,a,b,psym=7,symsize=3,thick=10 close,3 xyouts,/normal,0.25,0.95,figlab,charsize=2 device,/close set_plot,'x' end