parameter(maxstars=300000)
	dimension data(maxstars),sigmax(maxstars)
	character*6 filenames(9),temp
	data filenames/'rad2.0','rad2.5','rad3.0',
     &  'rad3.5',
     &   'rad4.0','rad4.5','rad5.0','rad5.5',
     &   'rad6.0'/
	ibest=0
	small=1e6
	do 100 i=1,9
	open(unit=1,file=filenames(i),status='old')
	open(unit=2,file='tempone',status='new')
C now read in the data
	ix=0
2	read(1,*,end=5) sigma
	fwhm=2.355*sigma
	if(fwhm.lt.1.5) go to 2
	if(fwhm.gt.7.0) go to 2
	ix=ix+1
	data(ix)=fwhm
	sigmax(ix)=1.
	go to 2
5	call xfit(data,sigmax,ix,0,xmean,sigmam,sigma)
	write(6,6) filenames(i),ix,sigma,xmean
6	format('File:',a10,' Num,sigma,mean=',i5,2f8.3)
	if(sigma.lt.small) then
	small=sigma
	ibest=i
	end if
100	continue
	temp=filenames(ibest)
	write(2,101) temp(4:6)
101	format(a3)
	write(6,7) temp(4:6),small
7	format('Best was:',a3,1x,f6.2)
	stop
	end	
c subroutine xfit.f
c
c source
c   Bevington, page 76.
c
c purpose
c   calculate the mean and estimated errors for a set of data points
c
c usage
c   call xfit (x, sigmax, npts, mode, xmean, sigmam, sigma)
c
c description of parameters
c   x      - array of data points
c   sigmax - array of standard deviations for data points
c   npts   - number of data points
c   mode   - determines method of weighting
c            +1 (instrumental) weight(i) = 1./sigmax(i)**2
c             0 (no weighting) weight(i) = 1.
c            -1 (statistical)  weight(i) = 1.
c   xmean  - weighted mean
c   sigmam - standard deviation of mean
c   sigma  - standard deviation of data
c
c subroutines and function subprograms required
c   none
c
	subroutine xfit (x,sigmax,npts,mode,xmean,sigmam,sigma)
	double precision sum,sumx,weight,free
	dimension x(*),sigmax(*)
c
c accumulate weighted sums
c
11	sum=0.
	sumx=0.
	sigma=0.
	sigmam=0.
20	do 32 i=1,npts
21	if (mode) 22,22,24
22	weight=1.
	goto 31
24	weight=1./sigmax(i)**2
31	sum=sum+weight
32	sumx=sumx+weight*x(i)
c
c evaluate mean and standard deviations
c
41	xmean=sumx/sum
51	do 52 i=1,npts
52	sigma=sigma+(x(i)-xmean)**2
	free=npts-1
54	sigma=dsqrt(sigma/free)
61	if (mode) 62,64,66
62	sigmam=dsqrt(xmean/sum)
	goto 70
64	sigmam=sigma/dsqrt(sum)
	goto 70
66	sigmam=dsqrt(1./sum)
70	return
	end