program analyze_damping
use bmad
implicit none
type(coord_struct) coord
real(rp) s_total, x(2,1000000),y(2,1000000),s(2,1000000)
integer i,j,jmax,tot_file, ios, lun,n
character*(100) file_name(100)/100*'null'/
character*16 element
namelist /file_list/file_name

  OPEN (UNIT=5, FILE='damping_file_list.dat', STATUS='old', IOSTAT=ios)
print *,'ios = ', ios
  READ (5, NML=file_list, IOSTAT=ios)
!   rewind(unit=5)
!    READ(5, NML=file_list, IOSTAT=ios)
    WRITE(6,NML=file_list)
  print *, 'ios=', ios
  CLOSE(5)

lun=lunget()  
i=1
do while(trim(file_name(i)) /= 'null')
   open(unit=lun, file = trim(file_name(i)))
   print '(i10,1x,a)', i,trim(file_name(i))
   ios=0
   j=0
   do while(ios >= 0)
     j=j+1
     read(lun,*,IOSTAT=ios) n,element,coord%vec, coord%t,coord%r, s_total
     x(i,j) = coord%vec(1)
     y(i,j) = coord%vec(3)
     s(i,j) = s_total
     print '(2i10,3es12.4)',i,j,x(i,j),y(i,j),s(i,j)
     jmax=j
     if(j>=1000000)exit
   end do
  close(unit=lun)
  tot_file = i
  i=i+1
end do ! file list


do j=1,jmax
! write(11,'(10(1x,3es12.4))')x(1,j), y(1,j), s(1,j),x(2,j), y(2,j), s(2,j)
 write(11,'(10(1x,3es12.4))')(x(i,j), y(i,j), s(i,j),i=1,tot_file)
end do
end program