diff options
Diffstat (limited to 'src/analysis/cf_coadd.c')
| -rw-r--r-- | src/analysis/cf_coadd.c | 192 |
1 files changed, 192 insertions, 0 deletions
diff --git a/src/analysis/cf_coadd.c b/src/analysis/cf_coadd.c new file mode 100644 index 0000000..eaa5fe2 --- /dev/null +++ b/src/analysis/cf_coadd.c @@ -0,0 +1,192 @@ +/****************************************************************************** + * Johns Hopkins University + * Center For Astrophysical Sciences + * FUSE + ****************************************************************************** + * + * Synopsis: cf_coadd root_name + * + * This program calls the routines that produce the exposure-level + * all files and the observation-level all, ano, and nvo files. + * + * 03/29/04 v0.1 bjg Begin work. + * 04/05/04 v0.2 bjg Remove unused variables + * Include unistd.h + * 06/07/04 v0.3 bjg Updated some keywords + * 08/12/04 v0.4 bjg Keywords QIK_COR and COMB_COR + * now set to COMPLETE. + * 08/26/04 v0.5 bjg Handle case aperture=RFPT + * 03/17/05 v0.6 tc Rename EXPTIME to OBSTIME in all extensions + * 08/22/05 v1.0 wvd Completely re-write for CalFUSE v3.1.7. + * Calls cf_make_all_exp for each exposure and + * cf_make_all_obs.csh only once. + * 08/27/07 v1.1 bot Changed i from int to long to match n + * Added L to corresponding values + * + *****************************************************************************/ + + + +#include <string.h> +#include <stdio.h> +#include <stdlib.h> +#include <unistd.h> +#include "calfuse.h" + +static char CF_PRGM_ID[] = "cf_coadd"; +static char CF_VER_NUM[] = "1.1"; + + +/*********************************************************************/ + +int main(int argc, char *argv[]) +{ + + char asnf_file[120], filename[120], instmode[100], aper_string[100], aper_num[100]; + fitsfile *asnf_fits_ptr, *infits; + long n; + int status = 0; + long i; + int anynull; + FILE * ffile; + + + char command0[256]; + + + char memname[120], *memname_ptr; + char memtype[120], *memtype_ptr; + char ** memname_ptr2; + char ** memtype_ptr2; + + char NULL_STR[] = {'\0'}; + + char file_present; + char expo_num_string[10]; + int expo_num; + + char know=0; + + /* Enter a timestamp into the log. */ + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Begin Processing"); + cf_error_init(CF_PRGM_ID, CF_VER_NUM, stderr); + + + memname_ptr=memname; + memtype_ptr=memtype; + + memname_ptr2=&memname_ptr; + memtype_ptr2=&memtype_ptr; + + + /* Input Parameters */ + if (argc != 2) + cf_if_error("Usage: cf_coadd rootname"); + + /*Open association file*/ + + strcpy(asnf_file, argv[1]); + strcat(asnf_file, "asnf.fit"); + FITS_open_file(&asnf_fits_ptr, asnf_file, READONLY, &status); + FITS_movabs_hdu(asnf_fits_ptr, 2, NULL, &status); + + + FITS_read_key(asnf_fits_ptr, TLONG, "NAXIS2", &n, NULL, &status); + + + for (i=0L;i<n;i++){ + + cf_if_fits_error(fits_read_col_str(asnf_fits_ptr, 1, i+1L, 1,1, NULL_STR, + memname_ptr2, &anynull, &status)); + + cf_if_fits_error(fits_read_col_str(asnf_fits_ptr, 2, i+1L, 1,1, NULL_STR, + memtype_ptr2, &anynull, &status)); + + cf_if_fits_error(fits_read_col_log(asnf_fits_ptr, 3, i+1L, 1, 1, 0, + &file_present, &anynull, &status)); + + if (file_present == FALSE) continue; + if (strncasecmp(memtype,"EXPOSURE",8)) continue; + + strncpy(expo_num_string,&memname[8],3); + expo_num=atoi(expo_num_string); + + if (expo_num==701) continue; + + strcpy(filename,memname); + strcat(filename,"1ahistfraw.fit"); + strcpy(instmode,"hist"); + ffile=fopen(filename,"r"); + if (ffile == NULL) { + strcpy(instmode,"ttag"); + strcpy(filename,memname); + strcat(filename,"1attagfraw.fit"); + ffile=fopen(filename,"r"); + if (ffile == NULL) { + printf("Warning: Expo %s not present\n",expo_num_string); + continue; + } + fclose(ffile); + } + else { + fclose(ffile); + } + + + if (!know) { + + FITS_open_file(&infits, filename, READONLY, &status); + FITS_read_key(infits, TSTRING, "APERTURE", aper_string, NULL, &status); + + if (!strncasecmp(aper_string,"LWRS",4)) strcpy(aper_num,"4"); + if (!strncasecmp(aper_string,"MDRS",4)) strcpy(aper_num,"2"); + if (!strncasecmp(aper_string,"HIRS",4)) strcpy(aper_num,"3"); + if (!strncasecmp(aper_string,"PINH",4)) strcpy(aper_num,"1"); + if (!strncasecmp(aper_string,"RFPT",4)) strcpy(aper_num,"4"); + + FITS_close_file(infits, &status); + know=1; + } + + + strcpy(command0,"cf_make_all_exp "); + strcat(command0,memname); + strcat(command0,"00all"); + strcat(command0,aper_num); + strcat(command0,instmode); + strcat(command0,"fcal.fit "); + strcat(command0,memname); + strcat(command0,"1alif"); + strcat(command0,aper_num); + strcat(command0,instmode); + strcat(command0,"fcal.fit "); + strcat(command0,memname); + strcat(command0,"1blif"); + strcat(command0,aper_num); + strcat(command0,instmode); + strcat(command0,"fcal.fit "); + strcat(command0,memname); + strcat(command0,"2alif"); + strcat(command0,aper_num); + strcat(command0,instmode); + strcat(command0,"fcal.fit "); + strcat(command0,memname); + strcat(command0,"2blif"); + strcat(command0,aper_num); + strcat(command0,instmode); + strcat(command0,"fcal.fit "); + + system(command0); + + } + + FITS_close_file(asnf_fits_ptr,&status); + + sprintf(command0,"cf_make_all_obs.csh %s", asnf_file); + system(command0); + + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, "Done Processing"); + + return(EXIT_SUCCESS); + +} |