program ecloudmod_test
	use bmad
	implicit none

	! Step values
	real(rp) :: dx, dy, sigx, sigy, x, y, mux, muy
	real(rp) :: kx, ky

	real(rp) :: off = 0
	real(rp) :: maxoff = 40
	real(rp) :: off_step_size = 5e-2

	x = -20
	y = 0
	mux = 0
	muy = 0


	sigx = 1
	sigy = 1

	do while (off<maxoff)
		dx = x + off - mux
		dy = y - muy

		call bbi_kick(dx/sigx, dy/sigy, sigy/sigx, kx, ky)
		write(100,'(6es12.4)') dx, dy, sigx, sigy, kx, ky
		off = off + off_step_size
	end do
end program