diff options
Diffstat (limited to 'src/analysis/cf_wrspec7.c')
| -rw-r--r-- | src/analysis/cf_wrspec7.c | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/src/analysis/cf_wrspec7.c b/src/analysis/cf_wrspec7.c new file mode 100644 index 0000000..979046b --- /dev/null +++ b/src/analysis/cf_wrspec7.c @@ -0,0 +1,58 @@ +/******************************************************************************* + * Johns Hopkins University + * Center For Astrophysical Sciences + * FUSE + ******************************************************************************* + * + * Synopsis: cf_wrspec7(fitsfile *outfits, long npts, float *wave, float *flux, + * float *error, long *counts, float *weights, + * float *bkgd,short *pothole) + * + * Description: Write a file containing a FUSE 2D spectrum. + * + * Arguments: fitsfile *fname Output FITS file structure + * + * Returns: none + * + * History: 12/15/03 bjg Begin work from cf_wrspec4. + * 01/12/04 bjg Call to fits_create_tbl with nrow=1 + * instead of npts + * 03/22/04 bjg Change POTHOLE to QUALITY + * 03/25/04 bjg Moved time stamp writing to beginning + * of routine + * 04/07/05 v1.1 tc Create a multi-rows table rather than + * a unique vector in a cell + * + ******************************************************************************/ + +#include <stdio.h> +#include <stdlib.h> +#include "calfuse.h" + +static char CF_PRGM_ID[] = "cf_wrspec7"; +static char CF_VER_NUM[] = "1.1"; + +void cf_wrspec7(fitsfile *outfits, long npts, float *wave, float *flux, + float *error, long *counts, float *weights,float *bkgd,short *pothole) +{ + char *ttype[] = {"WAVE", "FLUX", "ERROR", "COUNTS", "WEIGHTS", "BKGD", "QUALITY" }; + char *tform[] = {"1E", "1E", "1E", "1J", "1E", "1E", "1I"}; + char *tunit[] = {"ANGSTROMS", "ERG/CM2/S/A", "ERG/CM2/S/A", "COUNTS", "COUNTS", "COUNTS", "UNITLESS"}; + char extname[] = "FUSE 2D Spectrum"; + int tfields = 7, status = 0; + + /* Write time stamp to log file. */ + cf_timestamp(CF_PRGM_ID, CF_VER_NUM, ""); + + fits_create_tbl(outfits, BINARY_TBL, 1, tfields, ttype, tform, tunit, extname, &status); + + /* Write out the data. */ + FITS_write_col(outfits, TFLOAT, 1, 1L, 1L, npts, wave, &status); + FITS_write_col(outfits, TFLOAT, 2, 1L, 1L, npts, flux, &status); + FITS_write_col(outfits, TFLOAT, 3, 1L, 1L, npts, error, &status ); + FITS_write_col(outfits, TLONG, 4, 1L, 1L, npts, counts, &status); + FITS_write_col(outfits, TFLOAT, 5, 1L, 1L, npts, weights, &status); + FITS_write_col(outfits, TFLOAT, 6, 1L, 1L, npts, bkgd, &status); + FITS_write_col(outfits, TSHORT, 7, 1L, 1L, npts, pothole, &status); + +} |