1 files changed, 83 insertions, 0 deletions
diff --git a/src/libcf/cf_screen_saa.c b/src/libcf/cf_screen_saa.c
new file mode 100644
index 0000000..c79dc6e
--- /dev/null
+++ b/src/libcf/cf_screen_saa.c
@@ -0,0 +1,83 @@
+/*****************************************************************************
+ *              Johns Hopkins University
+ *              Center For Astrophysical Sciences
+ *              FUSE
+ *****************************************************************************
+ *
+ * Synopsis:    cf_screen_saa(fitsfile *infits, long nseconds,
+ *                            unsigned char *timeline_status, 
+ *			      float *timeline_long, float *timeline_lat);
+ *
+ * Description: Set the screening bit due to SAA boundary.
+ *
+ * Arguments:   fitsfile  *infits       Input FITS file pointer
+ *              long      nseconds      Number of points in the timeline table
+ *              unsigned char  *timeline_status The status flag array
+ *              float     *timeline_long The longitude array
+ *              float     *timeline_lat  The latitude array
+ *
+ * Calls:	saa
+ *
+ * Returns:	0 on success
+ *
+ * History:	10/29/02   1.1   jch    Initial coding
+ * 		12/20/02   1.3   jch    Include reading SAA table here
+ * 		12/20/02   1.4   wvd    Delete call to calfuse.h 
+ *              05/20/03   1.5   rdr    Add call to cf_proc_check
+ *              08/28/03   1.6   bjg    New definition of structure saareg
+ * 		09/02/03   1.7   wvd    Change calfusettag.h to calfuse.h
+ * 		09/10/03   1.8   wvd    Tidy up code.
+ *              02/17/05   1.9   wvd    Place parentheses around assignment
+ *                                      used as truth value.
+ *
+ ****************************************************************************/
+
+#include <stdlib.h>
+
+#include "calfuse.h"
+
+int
+cf_screen_saa(fitsfile *infits, long nseconds, unsigned char *timeline_status,
+	      float *timeline_long, float *timeline_lat)
+{
+    char CF_PRGM_ID[] = "cf_screen_saa";
+    char CF_VER_NUM[] = "1.9";
+
+    char  saa_cal[FLEN_VALUE];
+    int   errflg=0, status=0, hdutype;
+    long  i;
+    saareg saa_bounds;
+    fitsfile *saafits;
+
+    cf_error_init(CF_PRGM_ID, CF_VER_NUM, stderr);
+    cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Begin Processing");
+    if ((errflg = cf_proc_check(infits, CF_PRGM_ID))) return errflg;
+
+    /*
+     *  Read SAA boundaries from the specified calibration file.
+     */
+    FITS_read_key(infits, TSTRING, "SAAC_CAL", saa_cal, NULL, &status);
+    FITS_open_file(&saafits, cf_cal_file(saa_cal), READONLY, &status);
+    FITS_movabs_hdu(saafits, 2, &hdutype, &status);
+
+    saa_bounds.n_points =
+	   cf_read_col(saafits, TFLOAT, "LATITUDE",  (void **) &(saa_bounds.lat));
+    (void) cf_read_col(saafits, TFLOAT, "LONGITUDE", (void **) &(saa_bounds.lon));
+    
+    fits_close_file(saafits, &status);
+
+    /*
+     *  Set the SAA boundary bit if in SAA
+     */
+    for (i = 0; i < nseconds; i++) {
+	if (saa(&saa_bounds, (double)timeline_long[i], (double)timeline_lat[i])) {
+	    timeline_status[i] |= TEMPORAL_SAA; 
+	}
+    }
+    free(saa_bounds.lon);
+    free(saa_bounds.lat);
+
+    cf_proc_update(infits, CF_PRGM_ID, "COMPLETE");
+    cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Done processing");
+    return status;
+}