From 97fac5783043c471f25c07fe74d71f1b27422a69 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Sun, 17 Dec 2023 17:19:23 -0500
Subject: [PATCH] Add file cfom_iq.f90.

---
 qmap/CMakeLists.txt      |  3 +++
 qmap/libqmap/cfom_iq.f90 | 48 ++++++++++++++++++++++++++++++++++++++++
 2 files changed, 51 insertions(+)
 create mode 100644 qmap/libqmap/cfom_iq.f90

diff --git a/qmap/CMakeLists.txt b/qmap/CMakeLists.txt
index d633245e2..4a2a8f81b 100644
--- a/qmap/CMakeLists.txt
+++ b/qmap/CMakeLists.txt
@@ -42,6 +42,9 @@ add_executable (qmap ${qmap_CXXSRCS} ${qmap_CSRCS} ${qmap_GENUISRCS} qmap.rc)
 target_include_directories (qmap PRIVATE ${CMAKE_SOURCE_DIR} ${FFTW3_INCLUDE_DIRS})
 target_link_libraries (qmap wsjt_qt qmap_impl ${FFTW3_LIBRARIES} Qt5::Widgets Qt5::Network Usb::Usb)
 
+add_executable (cfom_iq libqmap/cfom_iq.f90)
+target_link_libraries (cfom_iq qmap_impl)
+
 if (WSJT_CREATE_WINMAIN)
   set_target_properties (qmap PROPERTIES WIN32_EXECUTABLE ON)
 endif (WSJT_CREATE_WINMAIN)
diff --git a/qmap/libqmap/cfom_iq.f90 b/qmap/libqmap/cfom_iq.f90
new file mode 100644
index 000000000..d6f765c1b
--- /dev/null
+++ b/qmap/libqmap/cfom_iq.f90
@@ -0,0 +1,48 @@
+program cfom_iq
+
+  parameter(NMAX=60*96000)
+  integer*2 id2(2,NMAX)
+  complex c,w,wstep
+  real*8 fcenter,uth8,twopi
+  character*6 mygrid
+
+  twopi=8.d0*atan(1.d0)
+  open(10,file='231028_0140.iq',status='old',access='stream')
+  open(12,file='231028_0140.cfom',status='unknown',access='stream')
+
+  mygrid='FN20OG'
+  nyear=2023
+  month=10
+  nday=28
+  uth8=01 + 40.d0/60
+  nfreq=1296
+  call astrosub00(nyear,month,nday,uth8,nfreq,mygrid,ndop0)
+  call astrosub00(nyear,month,nday,uth8+1.d0/60.d0,nfreq,mygrid,ndop1)
+
+  print*,ndop0,ndop1
+  
+  read(10) fcenter,id2(1:2,1:56*96000)
+  id2(1:2,56*96000+1:NMAX)=0
+
+  dop0=0.5*ndop0
+  dop1=0.5*ndop1
+  j=0
+  w=1.0
+  do isec=1,60
+     dop=dop0 + (i-0.5)*(dop1-dop0)/60.
+     dphi=-dop*twopi/96000.0
+     wstep=cmplx(cos(dphi),sin(dphi))
+     do n=1,96000
+        j=j+1
+        x=id2(1,j)
+        y=id2(2,j)
+        w=w*wstep
+        c=100.0*w*cmplx(x,y)
+        id2(1,j)=0.01*real(c)
+        id2(2,j)=0.01*aimag(c)
+     enddo
+  enddo
+
+  write(12) fcenter,id2
+        
+end program cfom_iq