#!/usr/bin/perl # # translate the HH MM SS and DD MM SS to decimal degrees of # RA and Dec for stars in the AGN standards file # (AJ 122, 2055, 2001) # # updated to read Brian Skiff's trimmed version of the catalog instead # -- MWR 11/15/2001 # while () { if (substr($_, 0, 1) eq "#") { printf "%s", $_; next; } if ($_ =~ /^\s+$/) { next; } $name = substr($_, 0, 22); for ($i = 21; $i > 0; $i--) { if (substr($name, $i, 1) ne " ") { last; } } for ($j = 0; $j < $i; $j++) { if (substr($name, $j, 1) eq " ") { substr($name, $j, 1) = "_"; } } $rah = substr($_, 24, 2); $ram = substr($_, 27, 2); $ras = substr($_, 30, 5); $decd = substr($_, 37, 3); $decm = substr($_, 40, 2); $decs = substr($_, 43, 4); if ($decd =~ /-.*/) { $sign = -1; } else { $sign = 1; } $s = substr($_, 41, 1); if ($s eq " ") { $s = "."; } if (0 == 1) { $u = substr($_, 48, 6); if ($u eq " ") { $u = 99.0; } $sigu = substr($_, 55, 5); if ($sigu eq " ") { $sigu = 99.0; } $b = substr($_, 67, 6); if ($b eq " ") { $b = 99.0; } $sigb = substr($_, 74, 5); if ($sigb eq " ") { $sigb = 99.0; } $v = substr($_, 86, 6); if ($v eq " ") { $v = 99.0; } $sigv = substr($_, 93, 5); if ($sigv eq " ") { $sigv = 99.0; } $r = substr($_, 105, 6); if ($r eq " ") { $r = 99.0; } $sigr = substr($_, 112, 5); if ($sigr eq " ") { $sigr = 99.0; } $i = substr($_, 124, 6); if ($i eq " ") { $i = 99.0; } $sigi = substr($_, 131, 5); if ($sigi eq " ") { $sigi = 99.0; } } else { $b = substr($_, 49, 6); if ($b eq " ") { $b = 99.0; } $sigb = substr($_, 56, 5); if ($sigb eq " ") { $sigb = 99.0; } $v = substr($_, 68, 6); if ($v eq " ") { $v = 99.0; } $sigv = substr($_, 75, 5); if ($sigv eq " ") { $sigv = 99.0; } $r = substr($_, 87, 6); if (($r eq " ") || (length($_) < 87) || ($r + 0 == 0)) { $r = 99.0; } $sigr = substr($_, 94, 5); if (($sigr eq " ") || (length($_) < 87) || ($r + 0 == 0)) { $sigr = 99.0; } $i = substr($_, 106, 6); if (($i eq " ") || (length($_) < 106) || ($i + 0 == 0)) { $i = 99.0; } $sigi = substr($_, 113, 5); if (($sigi eq " ") || (length($_) < 106) || ($i + 0 == 0)) { $sigi = 99.0; } } # calculate the RA and Dec in decimal degrees $ra_deg = ($rah + ($ram/60.0) + ($ras/3600.0))*15.0; if ($decd < 0) { $decd = 0.0 - $decd; } $dec_deg = $sign*($decd + ($decm/60.0) + ($decs/3600.0)); # now print out the information if (0 == 1) { printf "%23s %9.4f %9.4f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f %5.2f \n", $name, $ra_deg, $dec_deg, $u, $sigu, $b, $sigb, $v, $sigv, $r, $sigr, $i, $sigi; } else { printf "%23s %9.4f %9.4f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f \n", $name, $ra_deg, $dec_deg, $b, $sigb, $v, $sigv, $r, $sigr, $i, $sigi; } } exit 0;