/* This function is an interface for calling ACSRej. It sets up the
default parameters, copies printtime and verbose into the par structure,
and calls ACSRej.
*/
static int ACSRej_0 (char *input, char *output, char *mtype, int readnoise_only, int printtime, int verbose) {
extern int status;
clpar par; /* parameters used */
int newpar[MAX_PAR+1]; /* user specifiable parameters */
void rej_reset (clpar *, int []);
int AcsRej (char *, char *, char *, clpar *, int []);
rej_reset (&par, newpar);
par.printtime = printtime;
par.verbose = verbose;
par.readnoise_only = readnoise_only;
status = AcsRej (input, output, mtype, &par, newpar);
return (status);
}
int rej_command (int argc, char **argv, char **input, char *output,
clpar *par, int newpar[])
{
/* arguments
int argc; i: input command-line parameters
char **argv;
char **input; o: input file name or file list
char *output; o: output file name
clpar *par; o: user specified parameters
int newpar[]; o: array of parameters set by the user
*/
extern int status;
int ctoken; /* current command-line token being processed */
/* reset the parameters */
rej_reset (par, newpar);
/* not enough arguments
List all options for user
*/
if (argc < 3) {
syntax_error ("acsrej input output [-t] [-v]");
printf(" [-shadcorr] [-crmask] [-newbias] \n ");
printf(" [-table <filename>] [-scale #]\n");
printf(" [-init (med|min)] [-sky (none|mode)]\n");
printf(" [-sigmas #] [-radius #] [-thresh #]\n");
printf(" [-pdq #]\n");
return(status);
}
/* Get names of input and output files. These are mandatory.
Input list is also limited by tables/c_imt.c */
if ((*input = calloc(strlen(argv[1]) + 1, sizeof(char))) == NULL) {
printf ("Can't even begin; out of memory.\n");
return (status = OUT_OF_MEMORY);
}
strcpy (*input, argv[1]);
strcpy (output, argv[2]);
/* Scan remaining parameters. */
if (argc > 3) {
ctoken = 3;
while (ctoken < argc) {
/* switch must begin with a "-" */
if (argv[ctoken][0] != '-') {
return(syntax_error (argv[ctoken]));
} else {
/* These do not require additional arguments. */
if (strcmp("t", argv[ctoken]+1) == 0) {
par->printtime = 1;
ctoken++;
} else if (strcmp("v", argv[ctoken]+1) == 0) {
par->verbose = 1;
ctoken++;
} else if (strcmp("shadcorr", argv[ctoken]+1) == 0) {
par->shadcorr = 1;
ctoken++;
} else if (strcmp("newbias", argv[ctoken]+1) == 0) {
par->newbias = 1;
ctoken++;
} else if (strcmp("crmask", argv[ctoken]+1) == 0) {
par->mask = 1;
newpar[CRMASK] = 1;
newpar[TOTAL]++;
ctoken++;
/* These require one additional argument, which is
handled by the getArg functions.
*/
} else if (strcmp("table", argv[ctoken]+1) == 0) {
if (getArgT (argv, argc, &ctoken, par->tbname))
return (status = INVALID_VALUE);
} else if (strcmp("scale", argv[ctoken]+1) == 0) {
newpar[TOTAL]++;
newpar[SCALENSE] = 1;
if (getArgR (argv, argc, &ctoken, &par->scalense))
return (status = INVALID_VALUE);
} else if (strcmp("init", argv[ctoken]+1) == 0) {
newpar[TOTAL]++;
newpar[INITGUES] = 1;
if (getArgT (argv, argc, &ctoken, par->initgues))
return (status = INVALID_VALUE);
} else if (strcmp("sky", argv[ctoken]+1) == 0) {
newpar[TOTAL]++;
newpar[SKYSUB] = 1;
if (getArgT (argv, argc, &ctoken, par->sky))
return (status = INVALID_VALUE);
} else if (strcmp("sigmas", argv[ctoken]+1) == 0) {
newpar[TOTAL]++;
newpar[CRSIGMAS] = 1;
if (getArgT (argv, argc, &ctoken, par->sigmas))
return (status = INVALID_VALUE);
} else if (strcmp("radius", argv[ctoken]+1) == 0) {
newpar[TOTAL]++;
newpar[CRRADIUS] = 1;
if (getArgR (argv, argc, &ctoken, &par->radius))
return (status = INVALID_VALUE);
} else if (strcmp("thresh", argv[ctoken]+1) == 0) {
newpar[TOTAL]++;
newpar[CRTHRESH] = 1;
if (getArgR (argv, argc, &ctoken, &par->thresh))
return (status = INVALID_VALUE);
} else if (strcmp("pdq", argv[ctoken]+1) == 0) {
newpar[TOTAL]++;
newpar[BADINPDQ] = 1;
if (getArgS (argv, argc, &ctoken, &par->badinpdq))
return (status = INVALID_VALUE);
/* No match. */
} else
return (syntax_error (argv[ctoken]));
}
}
}
return (0);
}
