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/obsutil/src/sptime/mkcircle.cl | 16 ++++++++++++++++ 1 file changed, 16 insertions(+) create mode 100644 noao/obsutil/src/sptime/mkcircle.cl (limited to 'noao/obsutil/src/sptime/mkcircle.cl') diff --git a/noao/obsutil/src/sptime/mkcircle.cl b/noao/obsutil/src/sptime/mkcircle.cl new file mode 100644 index 00000000..2fc61210 --- /dev/null +++ b/noao/obsutil/src/sptime/mkcircle.cl @@ -0,0 +1,16 @@ +# MKCIRCLE -- Fraction of Gaussian profile going through a circular aperture +# of specified diameter in units of the profile FWHM. + +real d, logd + +printf ("## CIRCLE -- Fraction of Gaussian profile going through a") +printf (" cicular\n## aperture as a function of the diameter in units") +printf (" of FWHM.\n\n") + +printf ("%4.2f %5.3f\n", 0, 0) +for (logd=-1; logd<=0.6; logd=logd+0.1) { + d = 10.**logd + z = d * 0.8325546111577 + z = 1 - exp (-(z * z)) + printf ("%4.2f %5.3f\n", d, z) +} -- cgit