1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
/**************************************************************************
* Johns Hopkins University
* Center for Astrophysical Sciences
* FUSE
*************************************************************************
*
* Synopsis: cf_source_aper(infits, ap)
*
* Description: Procedure to determine which aperture contains the source
*
* Arguments : fitsfile infits Pointer to Integmediate Data File
* int *ap Array containing the two active
* apertures
*
* Returns int srctype Type of source in aperture
*
* HISTORY 12/12/02 v1.1 rdr First delivery
* 12/13/02 v1.3 wvd Changed file name
* 02/24/03 v1.4 peb Changed include file to calfitsio.h
* 04/17/03 v1.5 wvd Change to cf_source_aper throughout.
* 07/11/03 v1.6 rdr Change so that src_type returns
* 1 rather than 3 for extended source
* 08/21/03 v1.7 wvd If APERTURE=RFPT, treat as LWRS.
* 04/09/04 v1.8 bjg Include string.h
* Remove unused variables
*
************************************************************************/
#include "calfuse.h"
#include <string.h>
int
cf_source_aper(fitsfile *infits, int *active_ap )
{
/* char CF_PRGM_ID[] = "cf_source_aper"; */
/* char CF_VER_NUM[] = "1.8"; */
char aper[FLEN_CARD], src[FLEN_CARD];
int status=0, hdutype, src_type;
/* Determine the active aperture */
FITS_movabs_hdu(infits, 1, &hdutype, &status);
FITS_read_key(infits, TSTRING, "APERTURE", aper, NULL, &status);
/* Specify the channels of the active aperture */
if (!strncmp(aper,"HIRS",4)) {
active_ap[0] = 1;
active_ap[1] = 5; }
else if (!strncmp(aper,"MDRS",4)) {
active_ap[0] = 2;
active_ap[1] = 6; }
else if (!strncmp(aper,"LWRS",4) || !strncmp(aper,"RFPT",4)) {
active_ap[0] = 3;
active_ap[1] = 7;
}
else
cf_if_error("APERTURE keyword corrupted.");
cf_verbose(3, "APERTURE = %s, channels are %d and %d",
aper, active_ap[0], active_ap[1]);
FITS_read_key(infits, TSTRING,"SRC_TYPE", src, NULL, &status);
src_type = 0;
if (!strncmp(src,"E",1) || !strncmp(src,"e",1) ) src_type = 1;
cf_verbose(3, "SRC_TYPE = %s, source type = %d", src, src_type);
return src_type;
}
|
