subroutine write_rfquad

use bmad
use parameters_bmad

implicit none
integer lun
real(rp) t, scale, scalex,scaley, deltat/1.e-8/
!
lun=lunget()
open (unit=lun, file = trim(directory)//'/'//'rfquad.dat')
 print '(3(a,es12.4))',' horizontal  amp = ',rf_quad(1)%amp_h,' freq = ', rf_quad(1)%freq_h,' phase = ', rf_quad(1)%phi_h
 print '(3(a,es12.4))',' vertical  amp = ',rf_quad(1)%amp_v,' freq = ', rf_quad(1)%freq_v,' phase = ', rf_quad(1)%phi_v
 write(lun,'(4a12)')'time','x', 'y', 'sum'
t=0
do while(t < 25.e-6)
 scalex = rf_quad(1)%amp_h*sin(rf_quad(1)%freq_h * t + rf_quad(1)%phi_h)
 scaley = rf_quad(1)%amp_v*sin(rf_quad(1)%freq_v * t + rf_quad(1)%phi_v)
 scale = rf_quad(1)%amp_h*sin(rf_quad(1)%freq_h * t + rf_quad(1)%phi_h) + rf_quad(1)%amp_v*sin(rf_quad(1)%freq_v * t + rf_quad(1)%phi_v)
 write(lun,'(4es12.4)')t,scalex, scaley, scale
 t = t+deltat
end do
close(lun)
end subroutine