/* ACSCCD -- basic CCD image reduction */

# include <stdio.h>
# include <stdlib.h>		/* calloc */
# include <time.h>
# include <string.h>

int status = 0;			/* zero is OK */

# include <c_iraf.h>		/* for c_irafinit */
# include <ximio.h>
# include <hstio.h>

# include "../acs.h"
# include "../acsinfo.h"
# include "../acserr.h"
# include "../acscorr.h"		/* calibration switch names for acsccd */

static void FreeNames (char *, char *, char *, char *);

/* This is the main module for ACSCCD.  It gets the input and output
   file names, calibration switches, and flags, and then calls ACSccd.

   Warren Hack, 1998 June 1:
   	Revised for ACS... It will work on 1 ACS image at a time. Based on 
	original CALSTIS code by Phil Hodge.
*/

int main (int argc, char **argv) {

	char *inlist;		/* input file name */
	char *outlist;		/* output blev file name */
	int switch_on = 0;	/* was any switch specified? */
	int printtime = NO;	/* print time after each step? */
	int verbose = NO;	/* print additional info? */
	int quiet = NO;	/* print additional info? */
	int too_many = 0;	/* too many command-line arguments? */
	int i, j;		/* loop indexes */
	int k;

	IRAFPointer i_imt, o_imt;	/* imt list pointers */
	char *input;		/* name of input science file */
	char *output;		/* name of output file */
	int n_in, n_out;	/* number of files in each list */
	int n;

	/* Input and output suffixes. */
	char isuffix[] = "_raw";
	char osuffix[] = "_blv_tmp";

	/* A structure to pass the calibration switches to ACSCCD */
	CalSwitch ccd_sw;

	/* reference file keywords and names */
	RefFileInfo refnames;
		
	void InitRefFile (RefFileInfo *);
	void FreeRefFile (RefFileInfo *);
	void initSwitch (CalSwitch *);

	int ACSccd (char *, char *, CalSwitch *, RefFileInfo *, int, int);
	int DefSwitch (char *);
	int MkName (char *, char *, char *, char *, char *, int);
	void WhichError (int);
	int CompareNumbers (int, int, char *);


	c_irafinit (argc, argv);

	/* Allocate space for file names. */
	inlist = calloc (ACS_LINE+1, sizeof (char));
	outlist = calloc (ACS_LINE+1, sizeof (char));
	input = calloc (ACS_LINE+1, sizeof (char));
	output = calloc (ACS_LINE+1, sizeof (char));

	if (inlist == NULL || outlist == NULL ||
		input == NULL || output == NULL) {
	    printf ("Can't even begin; out of memory.\n");
	    exit (ERROR_RETURN);
	}
	inlist[0] = '\0';
	outlist[0] = '\0';
	input[0] = '\0';
	output[0] = '\0';
	
	/* Initialize the lists of reference file keywords and names. */
	InitRefFile (&refnames);

	/* Initial values. */
	initSwitch (&ccd_sw);

	for (i = 1;  i < argc;  i++) {

	    if (strcmp (argv[i], "-dqi") == 0) {	/* turn on */
			ccd_sw.dqicorr = PERFORM;
			switch_on = 1;
	    } else if (strcmp (argv[i], "-atod") == 0) {
			ccd_sw.atodcorr = PERFORM;
			switch_on = 1;
	    } else if (strcmp (argv[i], "-blev") == 0) {
			ccd_sw.blevcorr = PERFORM;
			switch_on = 1;
	    } else if (strcmp (argv[i], "-bias") == 0) {
			ccd_sw.biascorr = PERFORM;
			switch_on = 1;
	    } else if (strcmp (argv[i], "-flash") == 0) {
			ccd_sw.flashcorr = PERFORM;
			switch_on = 1;
	    } else if (argv[i][0] == '-') {
		for (j = 1;  argv[i][j] != '\0';  j++) {
		    if (argv[i][j] == 't') {
				printtime = YES;
		    } else if (argv[i][j] == 'v') {
				verbose = YES;
		    } else if (argv[i][j] == 'q') {
				quiet = YES;
		    } else {
				printf (MsgText, "Unrecognized option %s\n", argv[i]);
			    FreeNames (inlist, outlist, input, output);
				exit (1);
		    }
		}
	    } else if (inlist[0] == '\0') {
			strcpy (inlist, argv[i]);
	    } else if (outlist[0] == '\0') {
			strcpy (outlist, argv[i]);
	    } else {
			too_many = 1;
	    }
	}
	if (inlist[0] == '\0' || too_many) {
	    printf ("syntax:  acsccd [-t] [-v] [-q] input output\n");
	    printf ("  command-line switches:\n");
	    printf ("       -dqi  -atod -blev -bias\n");
	    FreeNames (inlist, outlist, input, output);
	    exit (ERROR_RETURN);
	}
	/* Initialize the structure for managing trailer file comments */
	InitTrlBuf ();
	
	/* Copy command-line value for QUIET to structure */
	SetTrlQuietMode(quiet);

	/* Was no calibration switch specified on command line? */
	if (!switch_on) {	/* default values (mostly PERFORM) */
	    ccd_sw.dqicorr  = DefSwitch ("dqicorr");
	    ccd_sw.atodcorr = DefSwitch ("atodcorr");
	    ccd_sw.blevcorr = DefSwitch ("blevcorr");
	    ccd_sw.biascorr = DefSwitch ("biascorr");
	    ccd_sw.flashcorr = DefSwitch ("postflsh");
	}

	/* Expand the templates. */
	i_imt = c_imtopen (inlist);
	o_imt = c_imtopen (outlist);
	n_in = c_imtlen (i_imt);
	n_out = c_imtlen (o_imt);

	/* The number of input and output files must be the same. */
	if (CompareNumbers (n_in, n_out, "output"))
	    status = 1;
	if (status) {
	    FreeNames (inlist, outlist, input, output);
		CloseTrlBuf();
	    exit (ERROR_RETURN);
	}

	/* Loop over the list of input files. */
	for (n = 0;  n < n_in;  n++) {

	    k = c_imtgetim (i_imt, input, ACS_LINE);
		
	    if (n_out > 0)
		    k = c_imtgetim (o_imt, output, ACS_LINE);
	    else
		    output[0] = '\0';

	    if (MkName (input, isuffix, osuffix, "", output, ACS_LINE)) {
			WhichError (status);
			sprintf (MsgText, "Skipping %s", input);
			trlmessage (MsgText);
			continue;
	    }

	    /* Calibrate the current input file. */
	    if (ACSccd (input, output, &ccd_sw, &refnames, printtime, verbose)) {
			sprintf (MsgText, "Error processing %s.", input);
			trlerror (MsgText);
			WhichError (status);
	    }
	}

	/* Close lists of file names, and free name buffers. */
	c_imtclose (i_imt);
	c_imtclose (o_imt);
	CloseTrlBuf();
	FreeRefFile (&refnames);
	FreeNames (inlist, outlist, input, output);

	if (status)
	    exit (ERROR_RETURN);
	else
	    exit (0);
}


static void FreeNames (char *inlist, char *outlist, char *input, char *output) {

	free (output);
	free (input);
	free (outlist);
	free (inlist);
}