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
68
69
70
|
include <mach.h> # for SZB_CHAR
include <tbset.h>
include "tbtables.h"
# tbfrsi -- read size info
# This routine reads NAXIS2 (number of rows) and TFIELDS (number of columns)
# from a FITS file and saves the values in the table descriptor.
#
# Phil Hodge, 6-Jul-1995 Subroutine created
# Phil Hodge, 2-Feb-1996 Check whether current HDU is a table.
# Phil Hodge, 7-Jun-1999 Use TB_SUBTYPE instead of TB_HDUTYPE.
# Phil Hodge, 29-Mar-2001 Assign NAXIS1 / 2 to TB_ROWLEN and TB_COLUSED.
procedure tbfrsi (tp)
pointer tp # i: pointer to table descriptor
#--
pointer sp
pointer comment # comment from FITS file
int status # used for fitsio
int keysexist # number of header keywords
int keysadd # space available for new header keywords
int naxis1 # row length in bytes
errchk tbferr
begin
status = 0
call smark (sp)
call salloc (comment, SZ_LINE, TY_CHAR)
if (TB_SUBTYPE(tp) == TBL_SUBTYPE_BINTABLE ||
TB_SUBTYPE(tp) == TBL_SUBTYPE_ASCII) {
call fsgkyj (TB_FILE(tp), "NAXIS1",
naxis1, Memc[comment], status)
if (status != 0)
call tbferr (status)
call fsgkyj (TB_FILE(tp), "NAXIS2",
TB_NROWS(tp), Memc[comment], status)
if (status != 0)
call tbferr (status)
call fsgkyj (TB_FILE(tp), "TFIELDS",
TB_NCOLS(tp), Memc[comment], status)
if (status != 0)
call tbferr (status)
} else {
# The current extension (or primary HDU) is not a table.
naxis1 = 0
TB_NROWS(tp) = 0
TB_NCOLS(tp) = 0
}
call fsghsp (TB_FILE(tp), keysexist, keysadd, status)
if (status != 0)
call tbferr (status)
TB_ROWLEN(tp) = (naxis1 + SZB_CHAR-1) / SZB_CHAR
TB_COLUSED(tp) = TB_ROWLEN(tp)
TB_MAXCOLS(tp) = TB_NCOLS(tp)
TB_NPAR(tp) = keysexist
TB_MAXPAR(tp) = keysexist + keysadd
TB_BOD(tp) = 0
call sfree (sp)
end
|
