  do n=1,ntot
     do i=1,iA
       read(20,71) xn1(i,n),yn1(i,n),zn1(i,n),iso(i)
    enddo
 enddo

71 format(3(1x,f12.8),1x,i1)