program main
          character cjunk
          real junk,old(1000000,3),new(1000000,3),sum(3)
          integer i,j,k

   15     i=0
          open(9,status='unknown',file='olddat.dat')
          read(9,*)cjunk
          read(9,*)cjunk
          read(9,*)cjunk
   20     I=I+1
            read(9,*,end=30) junk,junk,(old(i,j), j=1,3)
            go to 20   
   30     k=i-1                   
          close(9)            
          i=0
          open(9,file='newdat.dat')
          read(9,*)cjunk
          read(9,*)cjunk
          read(9,*)cjunk
          do 40 i=1,k
 40         read(9,*)junk,junk,(new(i,j), j=1,3)
          close(9)
          do 60 j=1,3
             sum(j)=0.0
             do 50 i=1,k
 50             sum(j)=sum(j)+(new(i,j)-old(i,j))**2
 60          sum(j)=sum(j)/real(k)
          write(6,*)'avg errors are ',(sqrt(sum(j)),j=1,3)
          end