From 88fceebb02ec21a1959f0d47ebce5195aa9f69ea Mon Sep 17 00:00:00 2001
From: Pavel Demin <pavel-demin@outlook.com>
Date: Sun, 1 Apr 2018 09:16:54 +0200
Subject: [PATCH] adapt subtractft8 to complex samples

---
 subtractft8.f90 | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/subtractft8.f90 b/subtractft8.f90
index 5dac110..e2a95b1 100644
--- a/subtractft8.f90
+++ b/subtractft8.f90
@@ -11,7 +11,8 @@ subroutine subtractft8(dd,itone,f0,dt)
 
   parameter (NMAX=15*12000,NFRAME=1920*79)
   parameter (NFFT=NMAX,NFILT=1400)
-  real*4  dd(NMAX), window(-NFILT/2:NFILT/2)
+  real*4 window(-NFILT/2:NFILT/2)
+  complex dd(NMAX)
   complex cref,camp,cfilt,cw
   integer itone(79)
   logical first
@@ -53,7 +54,7 @@ subroutine subtractft8(dd,itone,f0,dt)
 ! Subtract the reconstructed signal
   do i=1,nframe
      j=nstart+i-1
-     if(j.ge.1 .and. j.le.NMAX) dd(j)=dd(j)-2*REAL(cfilt(i)*cref(i))
+     if(j.ge.1 .and. j.le.NMAX) dd(j)=dd(j)-cfilt(i)*cref(i)
   enddo
 
   return