CHARACTER*40 INFILE1,INFILE2,FILENAME
	LOGICAL EFLG
	real dano(4,30000),cc(4),denom(4),mean(4)
c 	read first file into m1
6	WRITE(6,*) 'ENTER NAME OF dano hkl file:'
	READ(5,'(A)') INFILE1
	IF (INFILE1(:5).EQ.'     ') STOP
	  INQUIRE(FILE=INFILE1, NAME=FILENAME,EXIST=EFLG)
	  IF (.NOT.EFLG) GOTO 6
        OPEN (UNIT=2,FILE=INFILE1,STATUS='OLD',READONLY)
	J=1
10	read(2,11,end=500)h,k,l,(dano(i,j),i=1,3)
11	format(3I4,4F10.5)
	J=J+1
	goto 10
C   1   2  86  -3.29907  16.60873  -1.89133
C1234123412341234567890

500	N=J-1
	do 600 J=1,N
	do 600 I=1,3
600	mean(i)=mean(i)+dano(i,j)
	do 610 I=1,3
610	mean(i)=mean(i)/N
	write(6,*) 'mean values:',(mean(i),i=1,3)

	do 700 J=1,N
	do 700 I=1,3
	denom(i)=denom(i)+(dano(i,j)-mean(i))*(dano(i,j)-mean(i))
700	CC(i)=CC(i)+ (dano(i,j)-mean(i))*(dano(1,j)-mean(1))

	do 800 I=1,3
	x= cc(i)/sqrt((denom(i)*denom(1)))
800	write(6,*)'The cc between column',I,'and column 1 is', x
	end