2016-01-21 12:53:34 -05:00
|
|
|
subroutine makepings(pings,npts,width,sig)
|
|
|
|
|
|
|
|
|
|
real pings(npts)
|
|
|
|
|
real*8 t
|
|
|
|
|
real t0(14)
|
|
|
|
|
|
|
|
|
|
iping0=-999
|
|
|
|
|
dt=1.0/12000.0
|
2016-09-25 16:37:30 -04:00
|
|
|
do i=1,14
|
|
|
|
|
t0(i)=i !Make pings at t=1, 2, ... 14 s.
|
2016-01-21 12:53:34 -05:00
|
|
|
enddo
|
|
|
|
|
amp=sig
|
|
|
|
|
|
|
|
|
|
do i=1,npts
|
|
|
|
|
iping=min(max(1,i/12000),14)
|
2017-06-06 14:57:17 -04:00
|
|
|
t=(i*dt-t0(iping))/width
|
|
|
|
|
if(t.lt.0.d0 .and. t.lt.10.0) then !????
|
2016-01-21 12:53:34 -05:00
|
|
|
fac=0.
|
|
|
|
|
else
|
|
|
|
|
fac=2.718*t*dexp(-t)
|
|
|
|
|
endif
|
|
|
|
|
pings(i)=fac*amp
|
|
|
|
|
enddo
|
|
|
|
|
|
|
|
|
|
return
|
|
|
|
|
end subroutine makepings
|
