subroutine write_injection_line_twiss(ele)
       use bmad
       use parameters_bmad
       implicit none
       type (ele_struct) ele
       integer, save :: lun72
       logical first/.true./

       if(first)then
         lun72 = lunget()
         open(unit=lun72,file=trim(directory)//'/'//'injection_line_twiss.dat')
         write(lun72,'(9a12,1x,a16)')'s', 'beta x','beta y','alpha x','alpha y', 'eta x','etap x','phi x','phi y', 'element'
         first = .false.
       endif
       write(lun72,'(9es12.4,1x,a16)')ele%s, ele%a%beta,ele%b%beta,ele%a%alpha,ele%b%alpha,ele%x%eta,ele%x%etap,ele%a%phi, ele%b%phi,ele%name
       return
       end