Lines Matching refs:nptf
412 subroutine Filter(Filtered,DataHist,nptf,timestep,cutfreq) argument
416 integer nptf, cutfreq, j, k, m, s, Filtime(nptf) local
417 real*8 DataHist(nptf,numImpSrfs), Window(nptf)
418 real*8 Sinc(nptf), FilterSW(nptf), Filtered(nptf,numImpSrfs)
422 do j=1,nptf
424 Window(j) = 0.42+0.5*cos(2*pi*Filtime(j)/nptf)
425 & +0.08*cos(4*pi*Filtime(j)/nptf)
426 Sinc(j) = sin(pi*Filtime(j)*windK/nptf)/sin(pi*Filtime(j)/nptf)
430 do j=1,nptf
431 FilterSW(j) = Window(nptf+1-j)*Sinc(nptf+1-j) !filter for convolution
435 do m=1,nptf
436 do j=1,nptf
437 s=modulo(m-nptf+j,nptf)
439 s=nptf
442 & +FilterSW(j)*DataHist(s,:)/nptf !filter convolution