int main (int argc, char **argv) {
int status; /* zero is OK */
char *input, *output; /* file names */
int printtime = 0; /* print time after each step? */
int verbose = 0; /* print additional info? */
int too_many = 0; /* too many command-line arguments? */
int i, j; /* loop indexes */
c_irafinit (argc, argv);
input = calloc (STIS_LINE+1, sizeof (char));
output = calloc (STIS_LINE+1, sizeof (char));
if (input == NULL || output == NULL) {
printf ("ERROR: Can't even begin: out of memory.\n");
exit (ERROR_RETURN);
}
/* Get names of input and output files. */
for (i = 1; i < argc; i++) {
if (argv[i][0] == '-') {
if (strcmp (argv[i], "--version") == 0) {
PrVersion();
exit (0);
}
if (!(strcmp(argv[i],"--gitinfo")))
{
printGitInfo();
exit(0);
}
if (!(strcmp(argv[i],"--help")))
{
printHelp();
exit(0);
}
if (strcmp (argv[i], "-r") == 0) {
PrFullVersion();
exit (0);
}
for (j = 1; argv[i][j] != '\0'; j++) {
if (argv[i][j] == 't') {
printtime = 1;
} else if (argv[i][j] == 'v') {
verbose = 1;
} else {
printf ("ERROR: Unrecognized option %s\n", argv[i]);
printSyntax();
exit (1);
}
}
} else if (input[0] == '\0') {
strcpy (input, argv[i]);
} else if (output[0] == '\0') {
strcpy (output, argv[i]);
} else {
too_many = 1;
}
}
if (input[0] == '\0' || too_many) {
printSyntax();
exit (ERROR_RETURN);
}
if (output[0] == '\0') {
if ((status = MkName (input, "_x2d", "_sx2", output, STIS_LINE)))
exit (status);
}
/* Sum imsets. */
if ((status = CalStis8 (input, output, printtime, verbose))) {
printf ("Error processing %s.\n", input);
WhichError (status);
}
free (input);
free (output);
if (status)
exit (ERROR_RETURN);
else
exit (0);
}
int main (int argc, char **argv) {
int status; /* zero is OK */
char *inlist; /* list of input file names */
char *outlist; /* list of output file names */
int switch_on = 0; /* was any switch specified? */
int sgeocorr = OMIT; /* calibration switches */
int helcorr = OMIT;
int fluxcorr = OMIT;
int statcorr = OMIT;
int err_algorithm = WGT_VARIANCE;
int printtime = 0; /* print time after each step? */
int verbose = 0; /* print additional info? */
int center_target = 0; /* center target in output image? */
int too_many = 0; /* too many command-line arguments? */
int i, j; /* loop indexes */
double blazeshift = NO_VALUE;
IRAFPointer i_imt, o_imt; /* imt list pointers */
char *input; /* name of input science file */
char *output; /* optional name of output file */
int n_in, n_out; /* number of files in each list */
int n;
/* Input and output suffixes. */
char *isuffix[] =
{"_flt", "_crj", "_fwv", "_cwv", "_fwv_tmp", "_cwv_tmp"};
char *osuffix[] =
{"_x2d", "_sx2", "_w2d", "_w2d", "_w2d_tmp", "_w2d_tmp"};
int nsuffix = 6;
/* reference file keywords and names */
RefFileInfo refnames;
c_irafinit (argc, argv);
inlist = calloc (STIS_LINE+1, sizeof (char));
outlist = calloc (STIS_LINE+1, sizeof (char));
input = calloc (STIS_LINE+1, sizeof (char));
output = calloc (STIS_LINE+1, sizeof (char));
if (inlist == NULL || outlist == NULL ||
input == NULL || output == NULL) {
printf ("ERROR: Can't even begin: out of memory.\n");
exit (ERROR_RETURN);
}
/* Get command-line arguments. */
for (i = 1; i < argc; i++) {
if (strcmp (argv[i], "-x2d") == 0 ||
strcmp (argv[i], "-geo") == 0) {
switch_on = 1;
} else if (strcmp (argv[i], "-sgeo") == 0) { /* turn on */
sgeocorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-hel") == 0) {
helcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-flux") == 0) {
fluxcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-stat") == 0) {
statcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-wgt_err") == 0) {
err_algorithm = WGT_ERROR;
} else if (argv[i][0] == '-') {
if (strcmp (argv[i], "--version") == 0) {
PrVersion();
exit (0);
}
if (strcmp (argv[i], "-r") == 0) {
PrFullVersion();
exit (0);
}
for (j = 1; argv[i][j] != '\0'; j++) {
if (argv[i][j] == 't') {
printtime = 1;
} else if (argv[i][j] == 'v') {
verbose = 1;
} else if (argv[i][j] == 'c') {
center_target = 1; /* yes, center target */
} else if (argv[i][j] == 'b') {
blazeshift = (double) atof (argv[++i]);
if (i == argc-1)
break;
} else {
printf ("ERROR: Unrecognized option %s\n", argv[i]);
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: cs7.e [-t] [-v] [-c] [-wgt_err] [-b blazeshift] input output\n");
printf (" command-line switches: -x2d -sgeo -hel -flux -stat\n");
FreeNames (inlist, outlist, input, output);
exit (ERROR_RETURN);
}
/* Was no calibration switch specified on command line? */
if (!switch_on) { /* default values (mostly PERFORM) */
sgeocorr = DefSwitch ("sgeocorr");
helcorr = DefSwitch ("helcorr");
fluxcorr = DefSwitch ("fluxcorr");
statcorr = DefSwitch ("statcorr");
}
/* Initialize the list of reference file keywords and names. */
InitRefFile (&refnames);
/* 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")) {
FreeNames (inlist, outlist, input, output);
exit (ERROR_RETURN);
}
/* Loop over the list of input files. */
for (n = 0; n < n_in; n++) {
j = c_imtgetim (i_imt, input, STIS_LINE);
if (n_out > 0)
j = c_imtgetim (o_imt, output, STIS_LINE);
else
output[0] = '\0';
status = 0;
if ((status = MkOutName (input, isuffix, osuffix, nsuffix,
output, STIS_LINE))) {
WhichError (status);
printf ("Skipping %s\n", input);
continue;
}
/* Calibrate the current input file. */
if ((status = CalStis7 (input, output,
sgeocorr, helcorr, fluxcorr, statcorr,
&refnames, printtime, verbose, center_target,
blazeshift, err_algorithm))) {
printf ("Error processing %s.\n", input);
WhichError (status);
}
}
/* Close lists of file names, and free name buffers. */
c_imtclose (i_imt);
c_imtclose (o_imt);
FreeRefFile (&refnames);
FreeNames (inlist, outlist, input, output);
if (status)
exit (ERROR_RETURN);
else
exit (0);
}
int main (int argc, char **argv) {
int status; /* zero is OK */
char *inlist; /* list of input file names */
char *outlist; /* list of output file names */
char *blevlist; /* list of output blev file names */
int switch_on = 0; /* was any switch specified? */
int printtime = 0; /* print time after each step? */
int verbose = 0; /* print additional info? */
int too_many = 0; /* too many command-line arguments? */
int i, j; /* loop indexes */
int junk;
IRAFPointer i_imt, o_imt, b_imt; /* imt list pointers */
char *input; /* name of input science file */
char *output; /* optional name of output file */
char *outblev; /* optional file for blev values */
int n_in, n_out, n_blev; /* number of files in each list */
int n;
/* Input and output suffixes. */
char *isuffix[] = {"_raw", "_blv_tmp", "_crj_tmp", "_wav"};
char *osuffix[] = {"_flt", "_flt", "_crj", "_fwv"};
int nsuffix = 4;
/* A structure to pass the calibration switches to CalStis1 */
cs1_switch cs1_sw;
/* reference file keywords and names */
RefFileInfo refnames;
c_irafinit (argc, argv);
/* Allocate space for file names. */
inlist = calloc (1, sizeof (char)); /* allocated later */
outlist = calloc (1, sizeof (char));
blevlist = calloc (1, sizeof (char));
input = calloc (STIS_LINE+1, sizeof (char));
output = calloc (STIS_LINE+1, sizeof (char));
outblev = calloc (STIS_LINE+1, sizeof (char));
if (inlist == NULL || outlist == NULL || blevlist == NULL ||
input == NULL || output == NULL || outblev == NULL) {
printf ("ERROR: Can't even begin; out of memory.\n");
exit (ERROR_RETURN);
}
/* Initialize the lists of reference file keywords and names. */
InitRefFile (&refnames);
/* Initial values. */
cs1_sw.dqicorr = OMIT;
cs1_sw.atodcorr = OMIT;
cs1_sw.blevcorr = OMIT;
cs1_sw.doppcorr = OMIT;
cs1_sw.lorscorr = OMIT;
cs1_sw.glincorr = OMIT;
cs1_sw.lflgcorr = OMIT;
cs1_sw.biascorr = OMIT;
cs1_sw.darkcorr = OMIT;
cs1_sw.flatcorr = OMIT;
cs1_sw.shadcorr = OMIT;
cs1_sw.photcorr = OMIT;
cs1_sw.statcorr = OMIT;
strcpy (cs1_sw.darkscale_string, "");
for (i = 1; i < argc; i++) {
if (strcmp (argv[i], "--version") == 0) {
PrVersion();
exit (0);
}
if (strcmp (argv[i], "-r") == 0) {
PrFullVersion();
exit (0);
}
if (strcmp (argv[i], "-dqi") == 0) { /* turn on */
cs1_sw.dqicorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-atod") == 0) {
cs1_sw.atodcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-blev") == 0) {
cs1_sw.blevcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-dopp") == 0) {
cs1_sw.doppcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-lors") == 0) {
cs1_sw.lorscorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-glin") == 0) {
cs1_sw.glincorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-lflg") == 0) {
cs1_sw.lflgcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-bias") == 0) {
cs1_sw.biascorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-dark") == 0) {
cs1_sw.darkcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-flat") == 0) {
cs1_sw.flatcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-shad") == 0) {
cs1_sw.shadcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-phot") == 0) {
cs1_sw.photcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-stat") == 0) {
cs1_sw.statcorr = PERFORM;
switch_on = 1;
} else if (strcmp (argv[i], "-dscl") == 0) {
strcpy (cs1_sw.darkscale_string, argv[++i]);
switch_on = 1;
} else if (argv[i][0] == '-') {
for (j = 1; argv[i][j] != '\0'; j++) {
if (argv[i][j] == 't') {
printtime = 1;
} else if (argv[i][j] == 'v') {
verbose = 1;
} else {
printf ("ERROR: Unrecognized option %s\n", argv[i]);
exit (1);
}
}
} else if (inlist[0] == '\0') {
free (inlist);
if ((inlist = calloc (strlen(argv[i])+1, sizeof(char)))
== NULL) {
printf ("ERROR: Out of memory.\n");
exit (ERROR_RETURN);
}
strcpy (inlist, argv[i]);
} else if (outlist[0] == '\0') {
free (outlist);
if ((outlist = calloc (strlen(argv[i])+1, sizeof(char)))
== NULL) {
printf ("ERROR: Out of memory.\n");
exit (ERROR_RETURN);
}
strcpy (outlist, argv[i]);
} else if (blevlist[0] == '\0') {
free (blevlist);
if ((blevlist = calloc (strlen(argv[i])+1, sizeof(char)))
== NULL) {
printf ("ERROR: Out of memory.\n");
exit (ERROR_RETURN);
}
strcpy (blevlist, argv[i]);
} else {
too_many = 1;
}
}
if (inlist[0] == '\0' || too_many) {
printf ("syntax: cs1.e [-t] [-v] input output [outblev]\n");
printf (" command-line switches:\n");
printf (" -dqi -atod -blev\n");
printf (" -dopp -lors -glin -lflg\n");
printf (" -bias -dark -flat -shad -phot -stat\n");
FreeNames (inlist, outlist, blevlist, input, output, outblev);
exit (ERROR_RETURN);
}
/* Was no calibration switch specified on command line? */
if (!switch_on) { /* default values (mostly PERFORM) */
cs1_sw.dqicorr = DefSwitch ("dqicorr");
cs1_sw.atodcorr = DefSwitch ("atodcorr");
cs1_sw.blevcorr = DefSwitch ("blevcorr");
cs1_sw.doppcorr = DefSwitch ("doppcorr");
cs1_sw.lorscorr = DefSwitch ("lorscorr");
cs1_sw.glincorr = DefSwitch ("glincorr");
cs1_sw.lflgcorr = DefSwitch ("lflgcorr");
cs1_sw.biascorr = DefSwitch ("biascorr");
cs1_sw.darkcorr = DefSwitch ("darkcorr");
cs1_sw.flatcorr = DefSwitch ("flatcorr");
cs1_sw.shadcorr = DefSwitch ("shadcorr");
cs1_sw.photcorr = DefSwitch ("photcorr");
cs1_sw.statcorr = DefSwitch ("statcorr");
}
/* Expand the templates. */
i_imt = c_imtopen (inlist);
o_imt = c_imtopen (outlist);
b_imt = c_imtopen (blevlist);
n_in = c_imtlen (i_imt);
n_out = c_imtlen (o_imt);
n_blev = c_imtlen (b_imt);
/* The number of input and output files must be the same. */
status = 0;
if (CompareNumbers (n_in, n_out, "output"))
status = ERROR_RETURN;
if (CompareNumbers (n_in, n_blev, "outblev"))
status = ERROR_RETURN;
if (status) {
FreeNames (inlist, outlist, blevlist, input, output, outblev);
exit (ERROR_RETURN);
}
/* Loop over the list of input files. */
for (n = 0; n < n_in; n++) {
junk = c_imtgetim (i_imt, input, STIS_LINE);
if (n_out > 0)
junk = c_imtgetim (o_imt, output, STIS_LINE);
else
output[0] = '\0';
if (n_blev > 0)
junk = c_imtgetim (b_imt, outblev, STIS_LINE);
else
outblev[0] = '\0';
status = 0;
if ((status = MkOutName (input, isuffix, osuffix, nsuffix,
output, STIS_LINE))) {
WhichError (status);
printf ("Skipping %s\n", input);
continue;
}
/* Calibrate the current input file. */
if ((status = CalStis1 (input, output, outblev,
&cs1_sw, &refnames, printtime, verbose))) {
printf ("Error processing %s.\n", input);
WhichError (status);
}
}
/* Close lists of file names, and free name buffers. */
c_imtclose (i_imt);
c_imtclose (o_imt);
c_imtclose (b_imt);
FreeRefFile (&refnames);
FreeNames (inlist, outlist, blevlist, input, output, outblev);
if (status)
exit (ERROR_RETURN);
else
exit (0);
}