diff --git a/map65/libm65/ccf65.f90 b/map65/libm65/ccf65.f90
index f82d0fd06..07e829920 100644
--- a/map65/libm65/ccf65.f90
+++ b/map65/libm65/ccf65.f90
@@ -62,7 +62,7 @@ subroutine ccf65(ss,nhsym,ssmax,sync1,ipol1,jpz,dt1,flipk,      &
      call pctile(s,nhsym-1,50,base)
      s(1:nhsym-1)=s(1:nhsym-1)-base
      s(nhsym:NFFT)=0.
-     call four2a(s,NFFT,1,-1,0)                !Real-to-complex FFT
+     call four2a(cs,NFFT,1,-1,0)                !Real-to-complex FFT
      do i=0,NH
         cs2(i)=cs(i)*conjg(cpr2(i))            !Mult by complex FFT of pr2
         cs(i)=cs(i)*conjg(cpr(i))              !Mult by complex FFT of pr