Initial commit

1 files changed, 115 insertions, 0 deletions

diff --git a/src/fes/cf_fes_cal.c b/src/fes/cf_fes_cal.c
new file mode 100644
index 0000000..0f08bbd
--- /dev/null
+++ b/src/fes/cf_fes_cal.c
@@ -0,0 +1,115 @@
+/*******************************************************************************
+ *	      Johns Hopkins University
+ *	      Center For Astrophysical Sciences
+ *	      FUSE
+ *******************************************************************************
+ *
+ * Synopsis:	int cf_fes_cal(int cal_type, fitsfile *fesfits, float **image, 
+ *                                    int inaxis1, int inaxis2)
+ *
+ * Description: Get bias calibration filename for the raw FES image in
+ *                  the current hdu pointed at by *infits.
+ *              Call cf_fes_get_cal_image() to open the cal file and get
+ *                                          the calibration image.
+ *              Call cf_fes_apply_flat();
+ *              return the modified image to the calling routine.
+ *              4/10/99  MLR removed int inhdu as an argument and removed
+ *                       the call to fits_movabs_hdu.  This requires any
+ *                       user of this subroutine to do these steps before
+ *                       calling cf_fes_flat.
+ *              4/10/99  MLR removed (fitsfile *outfits, int outhdu) from 
+ *                       the argument list.  We will leave writing the image
+ *                       out to the calling routine.
+ *   
+ *
+ * Arguments:	fitsfile   *infits	Input FITS file pointer
+ *              float **image           already read in FES image
+ *              int inaxis1             image size  of the input image
+ *              int inaxis2
+ *
+ * Returns:	none
+ *
+ * History:	07/08/98	gak     calfes_design.070898 design documented
+ *              07/22/98        mlr     started work
+ *              04/10/99        mlt     pulled the actual work out of this
+ *                                      routine and into cf_fes_apply_flat
+ *              04/12/99        mlr     modified the input arguments(see above)
+ *                                      use the now generic cf_fes_read to
+ *                                      get the calibration image.
+ *              04/14/99        mlr     pulled the reading out and put it 
+ *                                      into cf_fes_get_cal_image.c
+ *              04/16/99        mlr     finished modifications to utilize
+ *                                      libcf and FITSIO.h(error handling).
+ *
+ * ToDo:  4/12/99 I have added the naxis1 and naxis2 parameters to the 
+ *                argument list.  I need to use them to compare binning
+ *                factors in the two images.
+ *                I should combine cf_fes_mask, cf_fes_bias, cf_fes_flat
+ *                into one generic cf_fes_do_cal("cf_fes_mask", *fesfits,...etc)
+ *
+ ******************************************************************************/
+
+#include <string.h>
+#include <stdio.h>
+#include <stdlib.h>
+#include "calfuse.h"
+#include "cf_calfes.h"
+
+#define CF_VER_NUM  "1.4"
+
+int cf_fes_cal(int cal_type, fitsfile *fesfits, float **image, 
+                              int inaxis1, int inaxis2)
+{
+	char	buffer[FLEN_CARD];
+        char    calfilename[30];
+        fitsfile *calfile;
+	int 	status; 
+        float   *calimage;
+	char *prog_id, *msg, *keyword, *msg2;
+
+	int (*func)(fitsfile *, float **, float *, int, int);
+
+        status = 0;
+
+        switch (cal_type) {
+           case FES_FLAT:
+		prog_id = "cf_fes_flat";
+		msg = "Begin FES flatfield";
+		msg2 = "Done  FES flatfield";
+		keyword = "FFLT1FCL";
+		func =  cf_fes_apply_flat;
+		break;
+           case FES_MASK:
+		prog_id = "cf_fes_mask";
+		msg = "Begin FES mask";
+		msg2 = "Done  FES mask";
+		keyword = "MASK_FCL";
+		func =  cf_fes_apply_mask;
+		break;
+           case FES_BIAS:
+		prog_id = "cf_fes_bias";
+		msg = "Begin FES bias";
+		msg2 = "Done  FES bias";
+		keyword = "BIAS1FCL";
+		func =  cf_fes_apply_bias;
+		break;
+	}
+
+
+	cf_timestamp(prog_id, CF_VER_NUM, msg);
+
+        FITS_read_key(fesfits, TSTRING, keyword, calfilename, buffer, &status);
+
+        if (cf_fes_get_cal_image(calfilename, &calimage, inaxis1, inaxis2) == -1)
+           return (-1);
+       
+        (*func)(fesfits, image, calimage, inaxis1, inaxis2);
+
+        cf_fes_proc_update(fesfits, prog_id, "COMPLETE");
+
+        free(calimage);
+
+        cf_timestamp(prog_id, CF_VER_NUM, msg2);
+
+        return(status);
+}