c For reading RAVE averaging files:
c There are about 4 different kinds of lines:
c Corr. coeff. for operator   2 =  0.86952
c R-factor for operator   2 w.r.t. operator 1 =  0.54242
c Shape similarity : (   0.917)
c Overall Totals:                       60056     0.268     0.316  3125.235 3042.129  1297.157     0.027  1297.157     0.875


	character*132 line
	character*40 flname,outname,imagefile,rflfile
	character*3 astring(50)
	character*2 number
	character*10 base
	real*4 ccop2(50)

	WRITE(6,*)'PROGRAM TO EXTRACT STATISTICS FROM AVERAGING OUTPUT'
    	WRITE(6,*)' AND WRITE THEM INTO cctstsum.rsl FOR PLOTTING.  970403'
	WRITE(6,*)
	write(6,*)' Enter name of msum output file:'
	read(5,11) flname
11	format(a40)
	open(unit=2,file=flname,status='old',recl=132)
	open (unit=3,file='cctstsum.rsl',status='new')

28	count=0
30	read(2,31,end=200) line
31	format(a132)
c 	write(6,*) line
c         1         2         3         4         5         6         7         8
c12345678901234567890123456789012345678901234567890123456789012345678901234567890
cncscore CCoefficients for myx:
cCoefficients for chk:

32 	if (line(1:7).EQ.'ncscore') goto 40
 	if (line(1:12).EQ.'Coefficients') goto 60
 	if (line(1:8).EQ.'--------') goto 80
c	write(6,*)'xxx',line(1:8),'xxx'
	goto 30

40	count=count+1
c	write(6,*)40
	decode(80,43,line) astring(count)
43	format(26x,A3)
c	write(6,*) astring(count)
	goto 30

60	count=count+1
c	write(6,*)60
	decode(80,63,line) astring(count)
63	format(17x,A3)
c	write(6,*) astring(count)
	goto 30

80	write(3,83)(astring(i),i=1,count)
83	format(20(' __',A3,'__'))
	goto 125

125	continue
	close (unit=2)
	open(unit=2,file=flname,status='old',recl=132)

128	count=0
130	read(2,31,end=200) line
c	write(6,*) line
c Predicted Correlation Coefficient =     0.807797
c Correlation Coefficient =     0.769942
c12345678901234567890123456789012345678901234567890123456789012345678901234567890
c         1         2         3         4         5         6         7         8

132 	if (line(1:31).EQ.' Predicted Correlation Coeffici') goto 140
 	if (line(1:26).EQ.' Correlation Coefficient =') goto 160
 	if (line(1:8).EQ.'--------') goto 180
c	write(6,*)'xxx',line(1:8),'xxx'
	goto 130

140	count=count+1
c	write(6,*)40
	decode(80,143,line) ccop2(count)
143	format(40x,f9.6)
c	write(6,*) ccop2(count)
	goto 130

160	count=count+1
c	write(6,*)60
	decode(80,163,line) ccop2(count)
163	format(30x,f9.6)
c	write(6,*) ccop2(count)
	goto 130

180	write(3,183)(ccop2(i),i=1,count)
183	format(20f8.5)
	goto 128


200	close (unit=2)
c	close (unit=3)
	end