void freeAsnInfo (AsnInfo *asn) {
/* Local variables */
int i, j, p; /* loop index */
int numexp;
/* Function definitions */
void initAsnProduct (ProdInfo *, int);
void initAsnInfo (AsnInfo *);
void initAsnSubProd (SubProdInfo *, int);
void initAsnExp (ExpInfo *);
/* If ASN table was read in, free it. */
if (asn->product != NULL) {
for (p = 0; p < asn->numprod; p++) {
for (i=0; i <= asn->numsp; i++) {
for (j=0; j <= asn->spmems[i]; j++) {
initAsnExp (&(asn->product[p].subprod[i].exp[j]));
}
free (asn->product[p].subprod[i].exp);
if (asn->spmems[i] > 0)
numexp = asn->spmems[i];
else
numexp = 1;
initAsnSubProd (&(asn->product[p].subprod[i]), numexp);
}
free (asn->product[p].subprod);
initAsnProduct (&(asn->product[p]), asn->numsp+1);
}
free (asn->product);
}
free (asn->spmems);
initAsnInfo (asn);
}
int CalAcsRun (char *input, int printtime, int save_tmp, int verbose, int debug, const unsigned nThreads, const unsigned cteAlgorithmGen, const char * pcteTabNameFromCmd) {
/* arguments:
char *input i: name of the FITS file/table to be processed
int printtime i: true --> print time stamps at intermediate steps
int *save_tmp i: true --> save temporary files
int verbose i: true --> print info during processing
int debug i: true --> print debugging info during processing
int onecpu i: true --> turn off use of OpenMP during processing
*/
extern int status;
ACSInfo acshdr; /* calibration switches, etc */
AsnInfo asn; /* association table data */
char *acsdth_input; /* Input list for ACSDTH */
int prod;
void PrBegin (char *);
void PrEnd (char *);
void PrFileName (char *, char *);
void TimeStamp (char *, char *);
/* Association table routines */
void initAsnInfo (AsnInfo *);
void freeAsnInfo (AsnInfo *);
int LoadAsn (AsnInfo *);
int ProcessACSCCD (AsnInfo *, ACSInfo *, int *, int, const unsigned nThreads, const unsigned cteAlgorithmGen, const char * pcteTabNameFromCmd);
int ProcessMAMA (AsnInfo *, ACSInfo *, int);
int AcsDth (char *, char *, int, int, int);
char *BuildDthInput (AsnInfo *, int);
void InitDthTrl (char *, char *);
/* Post error handler */
push_hstioerr (errchk);
PrBegin ("CALACS"); /* *** CALACS -- Version ... *** */
if (printtime)
TimeStamp ("CALACS started", "");
/* Determine if input is a single file or an association
table, then populate ASN structure with appropriate
information to control processing.
*/
initAsnInfo(&asn);
if (debug) {
trlmessage ("Initialized Association data ... ");
}
/* Copy Input filename to ASN structure */
strcpy (asn.input, input);
/* Print image name. */
trlmessage ("\n");
PrFileName ("input", asn.input);
/* Set verbose flag... */
asn.verbose = verbose;
asn.debug = debug;
/* LoadAsn will determine whether input is a single file or an
Association table. If a single image, it will look in that images
association table to see what products are associated with it, and
process them accordingly. If it is just a single image as its
own output, it will proceed as a 1 element table.
Based on routines from n_getAsnTable() and n_setup() in CALNICB
*/
if(LoadAsn(&asn)) {
freeAsnInfo (&asn);
return (status);
}
/* Check to see that detector is known, as it could come through with a
value of 0 or UNKNOWN_DETECTOR.
WJH 2Mar99
*/
if (asn.detector == UNKNOWN_DETECTOR || asn.detector == 0) {
trlwarn ("Unknown detector type for observations.");
freeAsnInfo(&asn);
return (status = NOTHING_TO_DO);
}
if (asn.verbose) {
sprintf (MsgText,"CALACS: Detector %s, type %d ",asn.instr, asn.detector);
trlmessage (MsgText);
}
/* Determine what detector we are working with... */
if (asn.detector != MAMA_DETECTOR ) { /* Process CCD data ... */
if (asn.verbose) {
trlmessage ("CALACS: processing a CCD product");
}
if (ProcessACSCCD(&asn, &acshdr, &save_tmp, printtime, nThreads, cteAlgorithmGen, pcteTabNameFromCmd)) {
if (status == NOTHING_TO_DO) {
trlwarn ("No processing desired for CCD data.");
} else {
trlerror ("Couldn't process CCD data");
}
freeAsnInfo(&asn);
return (status);
}
trlmessage("Finished CCD processing...");
} else { /* Process MAMA observations here */
trlmessage("Starting to process MAMA data now...");
if (ProcessMAMA(&asn, &acshdr, printtime)) {
if (status == NOTHING_TO_DO){
trlwarn ("No processing desired for MAMA data.");
} else{
trlerror ("Couldn't process MAMA data");
}
freeAsnInfo(&asn);
return (status);
}
trlmessage("Finished MAMA processing...");
}
/* Add DTH processing here... */
/* For each DTH product... */
if (asn.process == FULL){
if (asn.verbose) {
trlmessage ("CALACS: Building DTH products");
}
acsdth_input = NULL;
for (prod = 0; prod < asn.numprod; prod++) {
/* Create empty DTH product, header only */
/* This uses only one sub-product for the header template,
but later versions should use a function similar to
BuildSumInput to create list of subproducts as inputs...
*/
acsdth_input = BuildDthInput (&asn, prod);
/* We always want to create a final concatenated trailer file for
the entire association whether there is a product or not. So, we
set up the trailer file based on the association file name itself.
*/
/* If desired, we could optionally use the full _drz.tra filename
as the trailer filename, based on the output dither product name.
if (strcmp(asn.product[prod].prodname,"") != 0) {
InitDthTrl (acsdth_input, asn.product[prod].prodname);
} else { */
InitDthTrl(acsdth_input, asn.rootname);
/* End brace for optional dither product name assignment...
} */
/* Check if we have a PROD-DTH specified...*/
if (strcmp(asn.product[prod].prodname, "") != 0) {
if ((asn.dthcorr == PERFORM || asn.dthcorr == DUMMY)) {
if (AcsDth (acsdth_input, asn.product[prod].prodname, asn.dthcorr, printtime, asn.verbose) )
return (status);
/* Pass posid of 0 to indicate a PRODUCT is to be updated */
updateAsnTable(&asn, prod, NOPOSID);
}
} else {
trlwarn ("No DTH product name specified. No product created.");
/*status = ACS_OK; */
}
}
free (acsdth_input);
}
if (asn.verbose) {
trlmessage ("CALACS: Finished processing product ");
}
freeAsnInfo(&asn);
trlmessage ("\n");
PrEnd ("CALACS"); /* *** CALACS complete *** */
if (printtime)
TimeStamp ("CALACS completed", acshdr.rootname);
/* Return the final value of status, which should be ACS_OK if
all went well or some error condition otherwise. */
return (status);
}
