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 */
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", "_rac_tmp"};
char *osuffix[] = {"_blv_tmp", "_blc_tmp"};
int nsuffix = 2;
/* A structure to pass the calibration switches to WF3CCD */
CCD_Switch ccd_sw;
/* reference file keywords and names */
RefFileInfo refnames;
void InitRefFile (RefFileInfo *);
void FreeRefFile (RefFileInfo *);
void initCCDSwitches (CCD_Switch *);
int WF3ccd (char *, char *, CCD_Switch *, RefFileInfo *, int, int);
int DefSwitch (char *);
int MkName (char *, char *, char *, char *, char *, int);
void WhichError (int);
int CompareNumbers (int, int, char *);
/*===========================================================================*/
/* Initialize IRAF interface environment */
c_irafinit (argc, argv);
/* Post HSTIO error handler */
push_hstioerr (errchk);
/* Allocate space for file names. */
inlist = calloc (SZ_LINE+1, sizeof (char));
outlist = calloc (SZ_LINE+1, sizeof (char));
input = calloc (SZ_LINE+1, sizeof (char));
output = calloc (SZ_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. */
initCCDSwitches (&ccd_sw);
/* Parse the command-line arguments */
for (i = 1; i < argc; i++) {
if (!(strcmp(argv[i],"--version"))) {
printf("%s\n",WF3_CAL_VER_NUM);
exit(0);
}
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 if (argv[i][j] == 'r'){
printf ("Current version: %s\n", WF3_CAL_VER);
exit(0);
} 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: wf3ccd [-t] [-v] [-q] [-r] 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 ("flshcorr");
ccd_sw.fluxcorr = DefSwitch ("fluxcorr");
ccd_sw.pctecorr = DefSwitch ("pctecorr");
}
/* 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++) {
i = c_imtgetim (i_imt, input, SZ_LINE);
if (n_out > 0) {
i = c_imtgetim (o_imt, output, SZ_LINE);
} else {
output[0] = '\0';
}
if (MkOutName (input, isuffix, osuffix, nsuffix, output, SZ_LINE)) {
WhichError (status);
sprintf (MsgText, "Skipping %s", input);
trlmessage (MsgText);
continue;
}
/* Calibrate the current input file. */
if (WF3ccd (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);
}
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 *);
/* For image header access */
Hdr phdr;
int LoadHdr (char *, Hdr *);
int GetSwitch (Hdr *, char *, int *);
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) {
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 (argv[i][0] == '-') {
**********/
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?
default values (mostly PERFORM) except ATODCORR
if (!switch_on) {*/
ccd_sw.dqicorr = DefSwitch ("dqicorr");
ccd_sw.atodcorr = DefSwitch ("atodcorr");
ccd_sw.blevcorr = DefSwitch ("blevcorr");
ccd_sw.biascorr = DefSwitch ("biascorr");
/*}*/
/* 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';
/* Open input image in order to read its primary header. */
if (LoadHdr (input, &phdr)) {
WhichError (status);
sprintf (MsgText, "Skipping %s", input);
trlmessage (MsgText);
continue;
}
/* Determine osuffix. */
strcpy(osuffix, "_blv_tmp");
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);
}
int main (int argc, char **argv) {
/* Local variables */
char input[CHAR_FNAME_LENGTH+1];
int printtime = NO; /* print time after each step? */
int save_tmp = DUMMY; /* save temporary files? */
int verbose = NO; /* print info during processing? */
int debug = NO; /* print debug statements during processing? */
int quiet = NO; /* suppress STDOUT messages? */
int too_many = NO; /* too many command-line arguments? */
int onecpu = NO; /* suppress openmp usage by using only 1 thread?*/
int i, j; /* loop indexes */
/* Function definitions */
void c_irafinit (int, char **);
int CalWf3Run (char *, int, int, int, int, int);
void WhichError (int);
/* Initialize status to OK and MsgText to null */
status = WF3_OK;
MsgText[0] = '\0';
input[0] = '\0';
/* Initialize IRAF environment */
c_irafinit(argc, argv);
PtrRegister ptrReg;
initPtrRegister(&ptrReg);
/* Command line arguments:
** 0. Check for --version option
** 1. input file name
** 2. print time?
** 3. save intermediate files?
** 4. verbose?
*/
for (i = 1; i < argc; i++) {
if (!(strcmp(argv[i],"--version"))) {
printf("%s\n",WF3_CAL_VER_NUM);
freeOnExit(&ptrReg);
exit(0);
}
if (!(strcmp(argv[i],"--gitinfo")))
{
printGitInfo();
freeOnExit(&ptrReg);
exit(0);
}
if (!(strcmp(argv[i],"--help")))
{
printHelp();
freeOnExit(&ptrReg);
exit(0);
}
if (argv[i][0] == '-') {
for (j = 1; argv[i][j] != '\0'; j++) {
if (argv[i][j] == 't') {
printtime = YES;
} else if (argv[i][j] == 's') {
save_tmp = YES;
} else if (argv[i][j] == 'r'){
printf ("Current version: %s\n", WF3_CAL_VER);
freeOnExit(&ptrReg);
exit(0);
} else if (argv[i][j] == 'v') {
verbose = YES;
} else if (argv[i][j] == 'd') {
debug = NO;
} else if (argv[i][j] == 'q') {
quiet = YES;
} else if (argv[i][j] == '1'){
onecpu = YES;
} else {
printf ("Unrecognized option %s\n", argv[i]);
printSyntax();
freeOnExit(&ptrReg);
exit (ERROR_RETURN);
}
}
} else if (input[0] == '\0') {
strcpy (input, argv[i]);
} else {
too_many = YES;
}
}
if (input[0] == '\0' || too_many) {
printSyntax();
freeOnExit(&ptrReg);
exit (ERROR_RETURN);
}
/* Initialize the structure for managing trailer file comments */
InitTrlBuf ();
addPtr(&ptrReg, &trlbuf, &CloseTrlBuf);
trlGitInfo();
/* Copy command-line value for QUIET to structure */
SetTrlQuietMode (quiet);
/* Call the CALWF3 main program */
if (CalWf3Run (input, printtime, save_tmp, verbose, debug, onecpu)) {
if (status == NOTHING_TO_DO) {
/* If there is just nothing to do,
** as for ACQ images, just quit. */
status = 0;
sprintf (MsgText, "CALWF3 did NOT process %s", input);
trlmessage (MsgText);
freeOnExit(&ptrReg);
exit(0);
} else {
/* Error during processing */
sprintf (MsgText, "CALWF3 processing NOT completed for %s",
input);
trlerror (MsgText);
/* Provide interpretation of error for user */
WhichError (status);
freeOnExit(&ptrReg);
exit (ERROR_RETURN);
}
}
/* Successful completion */
sprintf (MsgText, "CALWF3 completion for %s", input);
trlmessage (MsgText);
freeOnExit(&ptrReg);
/* Exit the program */
exit (0);
}
