subroutine track_write_backleg_trajectory(lat,co) use bmad use parameters_bmad implicit none type(lat_struct) lat type(coord_struct), allocatable:: co(:) integer lun49, i lun49=lunget() open(unit=lun49, file=trim(directory)//'/'//'track_backleg_inflector.dat') do i=1,lat%branch(0)%n_ele_track print '(7es12.4,1x,a)',lat%ele(i)%s,co(i)%vec(1:6), lat%ele(i)%name write(lun49,'(7es12.4,1x,a)')lat%ele(i)%s,co(i)%vec(1:6), lat%ele(i)%name end do close(unit=lun49) return end