1 files changed, 37 insertions, 0 deletions
diff --git a/sys/imfort/db/imputd.x b/sys/imfort/db/imputd.x
new file mode 100644
index 00000000..fc633c23
--- /dev/null
+++ b/sys/imfort/db/imputd.x
@@ -0,0 +1,37 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<mach.h>
+
+# IMPUTD -- Put an image header parameter of type double.
+
+procedure imputd (im, key, dval)
+
+pointer	im			# image descriptor
+char	key[ARB]		# parameter to be set
+double	dval			# double precision value
+
+int	junk, i
+pointer	sp, sval
+int	dtoc(), strlen()
+
+begin
+	call smark (sp)
+	call salloc (sval, SZ_FNAME, TY_CHAR)
+
+        # Reduce the precision of the encoded value if necessary to fit in
+        # the FITS value field.  Start with NDIGITS_DP-1 as the precision
+        # estimate NDIGITS_DP is only approximate, and if we make up half a
+        # digit of precision the result can be 1.00000000000000001 instead
+        # of 1.0.
+
+        for (i=NDIGITS_DP-1;  i >= NDIGITS_RP;  i=i-1) {
+	    junk = dtoc (dval, Memc[sval], SZ_FNAME, i, 'g', SZ_FNAME)
+            if (strlen (Memc[sval]) < 20)
+                break
+        }
+
+	# Write the new value to the header.
+	call impstr (im, key, Memc[sval])
+
+	call sfree (sp)
+end