program test_muon_pulse_shape
 use precision_def
  implicit none
 real(rp) time
 integer pulse, i
 integer unit

do pulse = 1,8 
unit = 20 + pulse
write(unit,'(a10,i10)')'pulse = ',pulse-1
do i=1,10000
  call muon_pulse_shape(pulse, time)
  write(unit, '(i10, es12.4 )')i,time
end do
end do


end