diff --git a/q65w/libm65/CMakeLists.txt b/q65w/libm65/CMakeLists.txt
index b7da3953f..e23aa85ad 100644
--- a/q65w/libm65/CMakeLists.txt
+++ b/q65w/libm65/CMakeLists.txt
@@ -6,8 +6,6 @@ set (libm65_FSRCS
   astro.f90
   astro0.f90
   astrosub.f90
-  chkhist.f90
-  chkmsg.f90
   dcoord.f90
   decode0.f90
   dot.f90
@@ -20,8 +18,6 @@ set (libm65_FSRCS
   geocentric.f90
   grid2deg.f90
   indexx.f90
-  iqcal.f90
-  iqfix.f90
   lorentzian.f90
   moon2.f90
   moondop.f90
diff --git a/q65w/libm65/chkhist.f90 b/q65w/libm65/chkhist.f90
deleted file mode 100644
index c814c3ca3..000000000
--- a/q65w/libm65/chkhist.f90
+++ /dev/null
@@ -1,23 +0,0 @@
-subroutine chkhist(mrsym,nmax,ipk)
-
-  integer mrsym(63)
-  integer hist(0:63)
-
-  do i=0,63
-     hist(i)=0
-  enddo
-  do j=1,63
-     i=mrsym(j)
-     hist(i)=hist(i)+1
-  enddo
-
-  nmax=0
-  do i=0,63
-     if(hist(i).gt.nmax) then
-        nmax=hist(i)
-        ipk=i+1
-     endif
-  enddo
-
-  return
-end subroutine chkhist
diff --git a/q65w/libm65/chkmsg.f90 b/q65w/libm65/chkmsg.f90
deleted file mode 100644
index 457a80924..000000000
--- a/q65w/libm65/chkmsg.f90
+++ /dev/null
@@ -1,31 +0,0 @@
-subroutine chkmsg(message,cok,nspecial,flip)
-
-  character message*22,cok*3
-
-  nspecial=0
-  flip=1.0
-  cok="   "
-
-  do i=22,1,-1
-     if(message(i:i).ne.' ') go to 10
-  enddo
-  i=22
-
-10 if(i.ge.11) then
-     if ((message(i-3:i).eq.' OOO') .or. (message(20:22).eq.' OO')) then
-        cok='OOO'
-        flip=-1.0
-        if(message(20:22).eq.' OO') then
-           message=message(1:19)
-        else
-           message=message(1:i-4)
-        endif
-     endif
-  endif
-
-  if(message(1:3).eq.'RO ')  nspecial=2
-  if(message(1:4).eq.'RRR ') nspecial=3
-  if(message(1:3).eq.'73 ')  nspecial=4
-
-  return
-end subroutine chkmsg
diff --git a/q65w/libm65/iqcal.f90 b/q65w/libm65/iqcal.f90
deleted file mode 100644
index c0c4fce2d..000000000
--- a/q65w/libm65/iqcal.f90
+++ /dev/null
@@ -1,30 +0,0 @@
-subroutine iqcal(nn,c,nfft,gain,phase,zsum,ipk,reject)
-
-  complex c(0:nfft-1)
-  complex z,zsum,zave
-
-  if(nn.eq.0) then
-     zsum=0.
-  endif
-  nn=nn+1
-  smax=0.
-  ipk=1
-  do i=1,nfft-1                       !Find strongest signal
-     s=real(c(i))**2 + aimag(c(i))**2
-     if(s.gt.smax) then
-        smax=s
-        ipk=i
-     endif
-  enddo
-  pimage=real(c(nfft-ipk))**2 + aimag(c(nfft-ipk))**2
-  p=smax + pimage
-  z=c(ipk)*c(nfft-ipk)/p              !Synchronous detection of image
-  zsum=zsum+z
-  zave=zsum/nn
-  tmp=sqrt(1.0 - (2.0*real(zave))**2)
-  phase=asin(2.0*aimag(zave)/tmp)    !Estimate phase
-  gain=tmp/(1.0-2.0*real(zave))      !Estimate gain
-  reject=10.0*log10(pimage/smax)
-
-  return
-end subroutine iqcal
diff --git a/q65w/libm65/iqfix.f90 b/q65w/libm65/iqfix.f90
deleted file mode 100644
index 93cace4b1..000000000
--- a/q65w/libm65/iqfix.f90
+++ /dev/null
@@ -1,28 +0,0 @@
-subroutine iqfix(c,nfft,gain,phase)
-
-  complex c(0:nfft-1)
-  complex z,h,u,v
-
-  nh=nfft/2
-  h=gain*cmplx(cos(phase),sin(phase))
-
-  do i=1,nh-1
-     u=c(i)
-     v=c(nfft-i)
-     x=real(u)  + real(v)  - (aimag(u) + aimag(v))*aimag(h) +         &
-          (real(u) - real(v))*real(h)
-     y=aimag(u) - aimag(v) + (aimag(u) + aimag(v))*real(h)  +         &
-          (real(u) - real(v))*aimag(h)
-     c(i)=0.5*cmplx(x,y)
-     z=u
-     u=v
-     v=z
-     x=real(u)  + real(v)  - (aimag(u) + aimag(v))*aimag(h) +         &
-          (real(u) - real(v))*real(h)
-     y=aimag(u) - aimag(v) + (aimag(u) + aimag(v))*real(h)  +         &
-          (real(u) - real(v))*aimag(h)
-     c(nfft-i)=0.5*cmplx(x,y)
-  enddo
-
-  return
-end subroutine iqfix