From b9b9fe6412c41594484bac747e7f3bada344c8f0 Mon Sep 17 00:00:00 2001
From: Joe Taylor <joe@princeton.edu>
Date: Mon, 12 Dec 2022 12:37:39 -0500
Subject: [PATCH] More pruning...

---
 q65w/libm65/CMakeLists.txt |  3 ---
 q65w/libm65/coord.f90      | 40 --------------------------------------
 q65w/libm65/deg2grid.f90   | 30 ----------------------------
 q65w/libm65/grid2k.f90     | 12 ------------
 4 files changed, 85 deletions(-)
 delete mode 100644 q65w/libm65/coord.f90
 delete mode 100644 q65w/libm65/deg2grid.f90
 delete mode 100644 q65w/libm65/grid2k.f90

diff --git a/q65w/libm65/CMakeLists.txt b/q65w/libm65/CMakeLists.txt
index d45eaccff..c13dbb910 100644
--- a/q65w/libm65/CMakeLists.txt
+++ b/q65w/libm65/CMakeLists.txt
@@ -8,10 +8,8 @@ set (libm65_FSRCS
   astrosub.f90
   chkhist.f90
   chkmsg.f90
-  coord.f90
   dcoord.f90
   decode0.f90
-  deg2grid.f90
   dot.f90
   fchisq0.f90
   filbig.f90
@@ -22,7 +20,6 @@ set (libm65_FSRCS
   geocentric.f90
   graycode.f90
   grid2deg.f90
-  grid2k.f90
   indexx.f90
   interleave63.f90
   iqcal.f90
diff --git a/q65w/libm65/coord.f90 b/q65w/libm65/coord.f90
deleted file mode 100644
index 8934b612b..000000000
--- a/q65w/libm65/coord.f90
+++ /dev/null
@@ -1,40 +0,0 @@
-SUBROUTINE COORD(A0,B0,AP,BP,A1,B1,A2,B2)
-
-! Examples:
-! 1. From ha,dec to az,el:
-!      call coord(pi,pio2-lat,0.,lat,ha,dec,az,el)
-! 2. From az,el to ha,dec:
-!      call coord(pi,pio2-lat,0.,lat,az,el,ha,dec)
-! 3. From ra,dec to l,b
-!      call coord(4.635594495,-0.504691042,3.355395488,0.478220215,
-!        ra,dec,l,b)
-! 4. From l,b to ra,dec
-!      call coord(1.705981071d0,-1.050357016d0,2.146800277d0,
-!        0.478220215d0,l,b,ra,dec)
-! 5. From ra,dec to ecliptic latitude (eb) and longitude (el):
-!      call coord(0.d0,0.d0,-pio2,pio2-23.443*pi/180,ra,dec,el,eb)
-! 6. From ecliptic latitude (eb) and longitude (el) to ra,dec:
-!      call coord(0.d0,0.d0,-pio2,pio2-23.443*pi/180,el,eb,ra,dec)
-
-  SB0=sin(B0)
-  CB0=cos(B0)
-  SBP=sin(BP)
-  CBP=cos(BP)
-  SB1=sin(B1)
-  CB1=cos(B1)
-  SB2=SBP*SB1 + CBP*CB1*cos(AP-A1)
-  CB2=SQRT(1.e0-SB2**2)
-  B2=atan(SB2/CB2)
-  SAA=sin(AP-A1)*CB1/CB2
-  CAA=(SB1-SB2*SBP)/(CB2*CBP)
-  CBB=SB0/CBP
-  SBB=sin(AP-A0)*CB0
-  SA2=SAA*CBB-CAA*SBB
-  CA2=CAA*CBB+SAA*SBB
-  TA2O2=0.0 !Shut up compiler warnings. -db
-  IF(CA2.LE.0.e0) TA2O2=(1.e0-CA2)/SA2 
-  IF(CA2.GT.0.e0) TA2O2=SA2/(1.e0+CA2)
-  A2=2.e0*atan(TA2O2)
-  IF(A2.LT.0.e0) A2=A2+6.2831853
-  RETURN
-END SUBROUTINE COORD
diff --git a/q65w/libm65/deg2grid.f90 b/q65w/libm65/deg2grid.f90
deleted file mode 100644
index 9ca3602f8..000000000
--- a/q65w/libm65/deg2grid.f90
+++ /dev/null
@@ -1,30 +0,0 @@
-subroutine deg2grid(dlong0,dlat,grid)
-
-  real dlong                        !West longitude (deg)
-  real dlat                         !Latitude (deg)
-  character grid*6
-
-  dlong=dlong0
-  if(dlong.lt.-180.0) dlong=dlong+360.0
-  if(dlong.gt.180.0) dlong=dlong-360.0
-
-! Convert to units of 5 min of longitude, working east from 180 deg.
-  nlong=60.0*(180.0-dlong)/5.0
-  n1=nlong/240                      !20-degree field
-  n2=(nlong-240*n1)/24              !2 degree square
-  n3=nlong-240*n1-24*n2             !5 minute subsquare
-  grid(1:1)=char(ichar('A')+n1)
-  grid(3:3)=char(ichar('0')+n2)
-  grid(5:5)=char(ichar('a')+n3)
-
-! Convert to units of 2.5 min of latitude, working north from -90 deg.
-  nlat=60.0*(dlat+90)/2.5
-  n1=nlat/240                       !10-degree field
-  n2=(nlat-240*n1)/24               !1 degree square
-  n3=nlat-240*n1-24*n2              !2.5 minuts subsquare
-  grid(2:2)=char(ichar('A')+n1)
-  grid(4:4)=char(ichar('0')+n2)
-  grid(6:6)=char(ichar('a')+n3)
-
-  return
-end subroutine deg2grid
diff --git a/q65w/libm65/grid2k.f90 b/q65w/libm65/grid2k.f90
deleted file mode 100644
index f68b1409e..000000000
--- a/q65w/libm65/grid2k.f90
+++ /dev/null
@@ -1,12 +0,0 @@
-subroutine grid2k(grid,k)
-
-  character*6 grid
-
-  call grid2deg(grid,xlong,xlat)
-  nlong=nint(xlong)
-  nlat=nint(xlat)
-  k=0
-  if(nlat.ge.85) k=5*(nlong+179)/2 + nlat-84
-
-  return
-end subroutine grid2k