diff options
Diffstat (limited to 'noao/digiphot/apphot/fitsky/apreadsky.x')
| -rw-r--r-- | noao/digiphot/apphot/fitsky/apreadsky.x | 40 |
1 files changed, 40 insertions, 0 deletions
diff --git a/noao/digiphot/apphot/fitsky/apreadsky.x b/noao/digiphot/apphot/fitsky/apreadsky.x new file mode 100644 index 00000000..3db9a183 --- /dev/null +++ b/noao/digiphot/apphot/fitsky/apreadsky.x @@ -0,0 +1,40 @@ +# AP_READSKY -- Procedure to read sky values, sigma and skew from a file +# The x and y positions, sky mode, sigma, and skew values, number of sky +# pixels and number of rejected sky pixels are assumed to be columns +# 1 to 7 respectively. + +int procedure ap_readsky (fd, x, y, sky_mode, sky_sigma, sky_skew, nsky, + nsky_reject) + +int fd # sky file descriptor +real x, y # center of sky annulus +real sky_mode # sky valye +real sky_sigma # sky sigma +real sky_skew # skew of sky pixels +int nsky # number of sky pixels +int nsky_reject # number of rejected pixesl + +int stat +int fscan(), nscan() + +begin + # Initialize. + sky_mode = INDEFR + sky_sigma = INDEFR + sky_skew = INDEFR + nsky = 0 + nsky_reject = 0 + + # Read in and decode a sky file text line. + stat = fscan (fd) + if (stat == EOF) + return (EOF) + call gargr (x) + call gargr (y) + call gargr (sky_mode) + call gargr (sky_sigma) + call gargr (sky_skew) + call gargi (nsky) + call gargi (nsky_reject) + return (nscan ()) +end |