From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001 From: Joe Hunkeler Date: Tue, 11 Aug 2015 16:51:37 -0400 Subject: Repatch (from linux) of OSX IRAF --- noao/astutil/astradius.dat | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) create mode 100644 noao/astutil/astradius.dat (limited to 'noao/astutil/astradius.dat') diff --git a/noao/astutil/astradius.dat b/noao/astutil/astradius.dat new file mode 100644 index 00000000..62b9c6a5 --- /dev/null +++ b/noao/astutil/astradius.dat @@ -0,0 +1,27 @@ +# Print images which are within a given radius in the sky. + +# Get parameters. +racenter = clget ("astradius.racenter") +deccenter = clget ("astradius.deccenter") +epcenter = clget ("astradius.epcenter") +radius = clget ("astradius.radius") +ra = imget(clget("keywpars.ra")) +dec = imget(clget("keywpars.dec")) + +epoch = imget(clget("keywpars.epoch")) +if (str(epoch) == "" || real(epoch) == 0.) + date = imget(clget("keywpars.date_obs")) + ut = imget(clget("keywpars.ut")) + epoch = epoch (date, ut) +endif + +# Precess image coordinates to center epoch and compute separation. +radec = precess (ra, dec, epoch, epcenter) +ra1 = ra_precess (ra, dec, epoch, epcenter) +dec1 = dec_precess (ra, dec, epoch, epcenter) +sep = arcsep (racenter, deccenter, ra1, dec1) + +# Print result if within radius. +if (sep < real (radius)) + printf ("%-15s %s\n", $I, imget ("title")) +endif -- cgit