From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- noao/digiphot/daophot/daolib/dppadu.x | 36 +++++++++++++++++++++++++++++++++++ 1 file changed, 36 insertions(+) create mode 100644 noao/digiphot/daophot/daolib/dppadu.x (limited to 'noao/digiphot/daophot/daolib/dppadu.x') diff --git a/noao/digiphot/daophot/daolib/dppadu.x b/noao/digiphot/daophot/daolib/dppadu.x new file mode 100644 index 00000000..b1ca7c25 --- /dev/null +++ b/noao/digiphot/daophot/daolib/dppadu.x @@ -0,0 +1,36 @@ +include +include "../lib/daophotdef.h" + +# DP_PADU -- Read the gain value from the image header. + +procedure dp_padu (im, dao) + +pointer im # pointer to IRAF image +pointer dao # pointer to the daophot structure + +pointer sp, key +real padu +real imgetr(), dp_statr() + +begin + call smark (sp) + call salloc (key, SZ_FNAME, TY_CHAR) + call dp_stats (dao, CCDGAIN, Memc[key], SZ_FNAME) + if (Memc[key] == EOS) + padu = dp_statr (dao, PHOTADU) + else { + iferr { + padu = imgetr (im, Memc[key]) + } then { + padu = dp_statr (dao, PHOTADU) + call eprintf ("Warning: Image %s Keyword %s not found.\n") + call pargstr (IM_HDRFILE(im)) + call pargstr (Memc[key]) + } + } + if (IS_INDEFR(padu) || padu <= 0.0) + call dp_setr (dao, PHOTADU, 1.0) + else + call dp_setr (dao, PHOTADU, padu) + call sfree (sp) +end -- cgit