diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /pkg/tbtables/tbfrsi.x | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'pkg/tbtables/tbfrsi.x')
| -rw-r--r-- | pkg/tbtables/tbfrsi.x | 70 |
1 files changed, 70 insertions, 0 deletions
diff --git a/pkg/tbtables/tbfrsi.x b/pkg/tbtables/tbfrsi.x new file mode 100644 index 00000000..8137d633 --- /dev/null +++ b/pkg/tbtables/tbfrsi.x @@ -0,0 +1,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 |