subroutine slowft( adj, add, nyq, t0,dt,nt,tt,  f0,df,  nf,ff)
integer it,ie,     adj, add, nyq,       nt,                nf
complex cexp, cmplx,                    tt(nt),         ff(nf)
real pi2, freq, time, scale,       t0,dt,        f0,df
call adjnull(      adj, add,            tt,2*nt,        ff,2*nf)
			pi2= 2. * 3.14159265;		scale = 1./sqrt( 1.*nt)
			df = (1./dt) / nf
		if( nyq>0)
			f0 =   - .5/dt  
		else
			f0 =   - .5/dt + df/2.
do ie = 1, nf {   freq= f0 + df*(ie-1)
do it = 1, nt {   time= t0 + dt*(it-1)
	if( adj == 0 )
		ff(ie)= ff(ie) + tt(it) * cexp(cmplx(0., pi2*freq*time)) * scale
	else
		tt(it)= tt(it) + ff(ie) * cexp(cmplx(0.,-pi2*freq*time)) * scale
	}}
return; end