开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > OpenGL > 查看源码
deepsky.pl
资源名称：celestia-1.6.0.tar.gz　[点击查看]
上传用户：center1979
上传日期：2022-07-26
资源大小：50633k
文件大小：37k
源码类别：

OpenGL

开发平台：

Visual C++

deepsky.pl：源码内容
  1. #!/usr/bin/perl -w
  2. #
  3. # Author: Fridger.Schrempp@desy.de
  4. #
  5. #use Math::Libm ':all';
  6. use Math::Trig;
  7. use Math::Quaternion qw(slerp);
  8. #
  9. open(DSC, ">deepsky.dsc") || die "Can not create deepsky.dscn";
  10. #
  11. # Boiler plate
  12. #-------------
  13. #
  14. ($ver = "Revision: 1.50 ") =~ s/$//g;
  15. ($me = $0) =~ s/.*///;
  16. ($sec,$min,$hour,$mday,$mon,$year,$wday,$yday,$isdst) = gmtime;
  17. $year += 1900;
  18. $mon += 1;
  19. print DSC "# Revised NGC and IC Catalog, Wolfgang Steinicke, January 6, 2006n";
  20. print DSC "# http://www.ngcic.org/steinicke/default.htmn";
  21. print DSC "#n";
  22. print DSC "# Augmented by galaxies from Local Group volume (V <~ (10 Mpc)^3)n";
  23. print DSC "# Catalog of Neighboring Galaxies (Karachentsev+, 2004)n";
  24. print DSC "# http://cdsweb.u-strasbg.fr/viz-bin/Cat?J/AJ/127/2031n";
  25. print DSC "#n";
  26. print DSC "# Augmented by averages of distance entries fromn";
  27. print DSC "# NED-1D: NASA/IPAC Extragalactic Database of Distancesn";
  28. print DSC "# http://nedwww.ipac.caltech.edu/level5/NED1D/n";
  29. print DSC "#n";
  30. print DSC "# Augmented byn";
  31. print DSC "# NASA/IPAC Extragalactic Database (NED, batch)n";
  32. print DSC "# http://nedwww.ipac.caltech.edu/n";
  33. print DSC "#n";
  34. print DSC "# Augmented byn";
  35. print DSC "# Revised 3rd Reference Catalog of Bright Galaxies (RC3,VII/155)n";
  36. print DSC "# http://cdsweb.u-strasbg.fr/viz-bin/Cat?VII/155n";
  37. print DSC "#n";
  38. print DSC "# Augmented byn";
  39. print DSC "# Mark III Catalog of Galaxy Peculiar Velocities (Willick+ 1997,VII/198)n";
  40. print DSC "# http://cdsweb.u-strasbg.fr/viz-bin/Cat?VII/198#sRM2.18n";
  41. print DSC "#n";
  42. print DSC "# Augmented by distances fromn";
  43. print DSC "# The SBF Survey of Galaxy Distances. IV.n";
  44. print DSC "# SBF Magnitudes, Colors, and Distances, n";
  45. print DSC "# J.L. Tonry et al., Astrophys J 546, 681 (2001)n";
  46. print DSC "#n";
  47. print DSC "# Augmented by distances fromn";
  48. print DSC "# Compilation of "200 Brightest Galaxies"n";
  49. print DSC "# http://www.anzwers.org/free/universe/galax200.htmln";
  50. print DSC "#n";
  51. print DSC "# Augmented by redshifts-> distances fromn";
  52. print DSC "# VII/193 The CfA Redshift Catalogue, Version June 1995  (Huchra+ 1995)n";
  53. print DSC "# http://cdsweb.u-strasbg.fr/viz-bin/Cat?VII/193n";
  54. print DSC "#n";
  55. print DSC "# Using today's Hubble constant = 73.2 [km/sec/Mpc]n";  
  56. print DSC "# (WMAP 2007, 3 years running, legacy archive)n";
  57. print DSC "#n";
  58. print DSC "# Abreviations for various distance methods used:n";
  59. print DSC "# SBF= SBF (Surface Brightness Fluctuations), T-F= Tully-Fischern";
  60. print DSC "# V_cmb = rad. velocity in CMB frame & Hubble law n";
  61. print DSC "# cep=Cepheids, P=photometric, N(G)=planetary nebula (globular cluster) luminosity functionn";
  62. print DSC "#n";
  63. print DSC "# Adapted for Celestia with Perl script: $me $vern";
  64. print DSC "# Processed $year-$mon-$mday $wday $yday $isdst $hour:$min:$sec UTCn";
  65. print DSC "#n";
  66. print DSC "# by Dr. Fridger Schrempp, fridger.schrempp@desy.den";
  67. print DSC "# ------------------------------------------------------ n";
  68. print DSC "nn";
  69. #
  70. # Constants
  71. #----------
  72. #
  73. $H0 = 73.2;       # today's Hubble constant [km/sec/Mpc], WMAP 2007, 3 years
  74. $c = 299792.458;  # Celestia's speed of light [km/sec]
  75. $Omega_matter = 0.27;    # flat Universe
  76. $Omega_lambda = 0.73;
  77. $pc2ly = 3.26167; # parsec-to-lightyear conversion
  78. $epsilon = 23.4392911; # ecliptic angle
  80. # from peak of SBF-survey distribution <=== Maple
  81. $M_abs_ell = -20.133; # absolute (Blue) magnitude for elliptical galaxies
  82. $M_abs_spir = -19.07; # absolute (Blue) magnitude for spiral galaxies
  84. # output of LG member galaxy parameters, set $lgout = 0; to suppress.
  85. $lgout = 0;
  87. $count = 0;
  88. #
  89. # Constant hashes and arrays
  90. # --------------------------
  91. #
  92. # hash with distance methods used, key = $dist_flag
  93. %method =  ("0" => "  # distance uncertain!n",
  94.             "1" => "  # method: SBFn",
  95.             "2" => "  # method: T-Fn",
  96.             "3" => "  # method: V_cmbn",
  97.             "4" => "  # method: B200n",
  98.             "5" => "  # method: ZCATn",
  99.             "6" => "  # method: NED-1D averagen",
  100.             "7" => "  # method: othern",
  101.             "8" => "  # method: NED (batch)n"
  102.             );
  104. # Hubble types of nearby galaxies
  105. %HType =("-6" => "E0",
  106.  "-5" => "E0",
  107.  "-4" => "E0",
  108.  "-3" => "S0",
  109.  "-2" => "S0",
  110.  "-1" => "S0",
  111.   "0" => "S0",
  112.   "1" => "Sa",
  113.   "2" => "Sb",
  114.   "3" => "Sb",
  115.   "4" => "Sc",
  116.   "5" => "Sc",
  117.   "6" => "Sc",
  118.   "7" => "Sc",
  119.   "8" => "Sc",
  120.   "9" => "Irr",
  121.  "90" => "Irr",
  122.  "10" => "Irr",
  123.  "11" => "Irr");
  124.  
  125. #
  126. # add Milky Way (J2000)
  127. # --------------
  128. #
  129. $name = "Milky Way";
  130. $altname = " ";
  131.    
  132. # galactic center direction
  134. $ra  = "17 45 37.224";   #IAU 1959
  135. $dec = "-28 56 10.23";   #IAU 1959
  137. # galactic north pole at (J2000)
  139. $ra_pole  = "12 51 26.282";  #IAU 1959
  140. $dec_pole = "+27 07 42.01";  #IAU 1959
  142. # convert to hours, degrees, respectively
  144. $alpha = &ra2h ($ra);
  145. $delta = &dec2deg ($dec);
  147. $delta_pole = &dec2deg($dec_pole);
  148. $alpha_pole = &ra2h ($ra_pole);
  150. # position angle 
  151. $PA = 122.932 - 90;# - 90;  #IAU 1959
  154. $Mabs = -20.6;
  155. $tp = "SBc";
  156. $dist = 2.772e4;
  158. # Orientation (Axis, Angle)
  159. #
  161. &orientation($alpha_pole,$delta_pole,"MilkyWay",0,$PA);
  162.    
  163. print  DSC "# Local Group volume (< (10 Mpc)^3):n";
  164. print  DSC "# Milky Wayn";
  165. print  DSC "Galaxy "Milky Way"n";
  166. print  DSC "{n";
  167. print  DSC "        Type  "$tp"n";
  168. print  DSC "        CustomTemplate "MilkyWay.png"n";  
  169. printf DSC "        RA        %10.6fn",$alpha;
  170. printf DSC "        Dec       %10.6fn",$delta;
  171. printf DSC "        Distance  %10.4gn",$dist;
  172. printf DSC "        Radius    %10.4gn",50000;
  173. printf DSC "        AbsMag    %10.4gn",$Mabs;
  174. printf DSC "        Axis    [%8.4f %8.4f %8.4f]n",@v;
  175. printf DSC "        Angle    %8.4fn",$angle;
  176. print  DSC "        InfoURL  "http://www.seds.org/messier/more/mw.html"n";  
  177. print  DSC "}nn";
  179. if ($lgout){
  180. $dkpc = $dist/($pc2ly * 1e3);
  181. $rar  = "17 45 37.2";
  182. $decr = "-28 56 10"; 
  183. printf "%17s %15s; %10s; %9s; %5s; %8.3f;  %6.2fn",$name,$altname,$rar,$decr,$tp,$dkpc,$Mabs;
  184. }
  186. $count++;
  188. #
  189. # read in the averaged distances from the NED1D distance data base
  190. # http://nedwww.ipac.caltech.edu/level5/NED1D/
  191. #
  192. open(NED1D," < catalogs/ned1d.txt")|| die "Can not read ned1d.txtn";
  193. $name_old = "";
  194. $ct = 1;
  195. $n_ned1d = 0;
  196. while (<NED1D>){
  197. next if (/^Contents/);
  198. next if (/^s+/);
  199. next if(/^Galaxy/);
  200. next if(/^s*(kpc)/);
  201. next if(/^s*(1)/);
  202. ($name,$dist) = split (/ss+/);
  203. ($name1,$namealt) = split(/,/,$name);
  204. if ($namealt){
  205.      $name = $name1;    
  206. }
  208. $dist = &clean($dist);
  209. $name = &clean($name);
  211.     # adapt syntax of galaxy names to RC3/NGC/IC conventions
  212. $name = &standardize($name);
  214. # only consider Messier, NGC and IC galaxies in NED1D data base for now
  216. next unless ($name =~ /^Msd+|^NGC|^IC/);    
  218.     # calculate the distance averages for all entries with the same $name
  219. if ($name eq $name_old){
  220. $ct++;
  221. $d = (($ct-1) * $d + $dist)/$ct;
  222. } else {
  223.     $n_ned1d++;
  224.     $d = $dist;
  225. $ct = 1;
  226. }
  227. # convert kpc to ly
  228.     $distance_ned1d{$name} = $d * $pc2ly * 1e3;
  229.     #print "$name  $distance_ned1d{$name}n" if ($name ne $name_old);
  230. #print "$name       $distance_ned1d{$name}n" if ($name =~/^Msd+/ && $name ne $name_old);
  231. $name_old = $name;
  232. }
  233. close(NED1D);
  235. # +++++++++++++++++++++++++++++++++++++++++++++++++++
  236. # read in revised SBF-survey data IV for fast lookup
  237. # +++++++++++++++++++++++++++++++++++++++++++++++++++
  239. open(SBF," < catalogs/sbf.txt")|| die "Can not read sbf.txtn";
  240. $n_sbf=0;
  241. while (<SBF>) {
  242.     next if (/^Pref_name/);
  243.     $name = &ss(1,8);
  244.     $name =~ s/s+//g;
  245.     $name = &standardize($name);
  247. # only consider NGC and IC galaxies in NED1D data base for now
  249. next unless ($name =~ /^NGC|^IC/);    
  250.     
  251.     #
  252.     # SBF-distances [ly]
  253.     # ------------------
  254.     #
  255.     $deltam   = &ss(77,81);
  256.     $dm{$name} = $deltam;
  257.     $dm{$name} =~ s/ //g;
  258.     $sbf_distance{$name} = 10**(($deltam+5)/5)*$pc2ly;
  259.     #print "$name  $sbf_distance{$name}n";
  260.     $n_sbf++;
  261. }
  262. close (SBF);
  264. # ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  265. # read in Tully-Fisher (T-F) & Dn-sigma distance data for fast lookup
  266. # ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  268. # read in first GINnumber ($gin), name and BT and B-V
  269. open(TFE," < catalogs/TF/egalf1.dat") || die "Can not read egalf1.datn";
  270. while (<TFE>) {
  271.     $gin = &ss(1,5);
  272.     $gin =~ s/ //g;
  273.     $name = &standardize (&ss(7,16));
  275. # only consider NGC and IC galaxies in NED1D data base for now
  277. next unless ($name =~ /^NGC|^IC/);    
  278.     
  279.     $name_f1{$gin} = $name;
  280.     #print "$namen";
  281. }
  282. close (TFE);
  284. # next enter the corresponding T-F distances for elliptical galaxies, 
  285. # match via $gin with egalf1.dat
  287. open(TF," < catalogs/TF/egalf2.dat") || die "Can not read egalf2.datn";
  288. $n_tf_e=0;
  289. while (<TF>) {
  290.     $gin = &ss(1,4);
  291.     $gin =~ s/ //g;;
  292.     $name = $name_f1{$gin} if ($name_f1{$gin});
  293.     $distance_TFE = &ss(35,41);
  294.     $distance_TFE = &clean($distance_TFE);
  295.     $distance_TFE{$name} = $distance_TFE/$H0 * 1.0e6 * $pc2ly; # [ly]
  296.     # $v_CMB_TFE{$name} = &ss(32,38) + &ss(40,46);
  297.     #printf "%5d %15s   %10.4gn",$gin,$name, $distance_TFE{$name};
  298.     $n_tf_e++;
  299. }
  300. close(TF);
  302. # next enter the corresponding T-F distances for spiral  galaxies
  304. open(TF," < catalogs/TF/tf-spiral.dat") || die "Can not read tf-spiral.datn";
  305. $n_tf_s=0;
  306. while (<TF>) {
  307.     $pgc = &ss(6,10);
  308.     $pgc =~ s/ //g;
  309.     $name = &standardize(&ss(12,20));
  311. # only consider NGC and IC galaxies in NED1D data base for now
  313. next unless ($name =~ /^NGC|^IC/);    
  314.     
  315.     #$name_TFS{$pgc} = $name;
  316.     $distance_TF = &ss(82,86);
  317.     $distance_TF =~ s/ //g;
  318.     $distance_TFS{$name} = $distance_TF/$H0 * 1.0e6*$pc2ly; # [ly]
  319.     $logr = &ss(66,70);
  320.     $logr =~ s/ //g;
  321.     #printf "PGC %6d %15s   %10.4gn",$pgc, $name, $distance_TFS{$name};
  322.     $n_tf_s++;
  323. }
  324. close(TF);
  326. # ++++++++++++++++++++++++++++++++++++++
  327. # read in data on 200 brightest galaxies
  328. # ++++++++++++++++++++++++++++++++++++++
  330. $n_b200 = 0;
  331. open( B200," < catalogs/200.txt") || die "Can not read 200.txtn";
  332. while (<B200>) {
  333.     next if (/^s+/);
  334.     $name = &ss(1,8);
  335.     $name = &cmpname ($name) if ($name =~ /^N|^I/);;
  336.     $name =~ s/s*$//g;
  337.     $messier = &ss(10,13);
  338.     $messier =~ s/ //g;
  339.     if ($messier){
  340.      $messier =~ s/M(d+)/M $1/;
  341.      $name = $messier;
  342.     }
  343. # only consider Messier, NGC and IC galaxies in NED1D data base for now
  345.     next unless ($name =~ /^Msd+|^NGC|^IC/);    
  346.     
  347.     $distance_B200{$name} = &ss(74,79);
  348.     $distance_B200{$name} =~ s/ //g;
  349.     $distance_B200{$name} = 1e6 * $distance_B200{$name};
  350.     $method_B200{$name} = &ss(83,90);
  351.     $method_B200{$name} =~ s/s+$//g;
  352.     #print "$name   $distance_B200{$name} $method_B200{$name}n";
  353.     $n_b200++;
  354. }
  355. close(B200);
  358. # +++++++++++++++++++++++++++++++
  359. # read in redshift data from zcat
  360. # +++++++++++++++++++++++++++++++
  362. $n_zcat = 0;
  363. open( ZCAT," < catalogs/zcat.txt") || die "Can not read zcat.txtn";
  364. while (<ZCAT>) {
  365.     next if (/^s+/);
  366.     $name = &clean(&ss(1,11));
  367.     $name = &standardize ($name); 
  369. # only consider NGC and IC galaxies in NED1D data base for now
  371.     next unless ($name =~ /^NGC|^IC/);    
  372.     
  373.     $cz = &ss(32,36);
  374.     $cz =~ s/ //g;
  375. # skip if redshift z is missing    
  376.     next if (! $cz);
  377.     $z = $cz/$c;
  378. # distance inclusive relativistic cosmological  corrections    
  379.     #printf "%10s  %10.3f  %10.3fn", $name, $z,  $distance_zcat{$name};
  380.     $distance_zcat{$name} = &dcmr($Omega_matter,$Omega_lambda,$z);
  381.     $n_zcat++;
  382. }
  383. close(ZCAT);
  385. # +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  386. # Read in missing redshift data from NED batch-job for fast lookup
  387. # +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  389. open(NED," < catalogs/nedz8.txt")|| die "Can not read nedz.txtn";
  390. $n_ned = 0;
  391. while (<NED>){
  392.     if (/(s+d)sG/){
  393.      $name = &standardize(&ss(19,28));
  394.      next if ($name =~ /NE*$/); # drop the NGC|IC xxx NED 1/2 objects
  395. $z = &ss(75,80);
  396. $z =~ s/ //g;   
  397.      if ($name =~ /^NGC 4518/){$z = 0.021868;}
  398.      if ($name =~ /^IC 1947/){$z = 0.03851;}
  399.      if ($name =~ /^IC 3658/){$z = 0.00015;}
  400. #printf "%10s %10.4fn",$name,$z if ($z);    
  401.      if ($z && $z > 0){
  402.      $distance_ned{$name} = &dcmr($Omega_matter,$Omega_lambda,$z);
  403.         #printf "%10d   %15s %10.4fn",$n_ned, $name,$z; 
  404.         $n_ned++;
  405.     }    
  406.     } 
  407. }
  409. # +++++++++++++++++++++++++++++++++++++++++++++++++
  410. # Read in revised RC3 data rc3.txt for fast lookup
  411. # +++++++++++++++++++++++++++++++++++++++++++++++++
  413. open(RC3," < catalogs/rc3.txt")|| die "Can not read rc3.txtn";
  414. $n_rc3 = 0;
  415. while (<RC3>) {
  416.     $pgc = &ss(106,116);
  417.     $pgc =~ s/^PGCs*0*(d+)s*$/$1/g;
  419.     #
  420.     # Galaxy name in RC3 as tag
  421.     #--------------------------
  422.     #
  423.     $n1 = &standardize(&ss(63,74));
  424.     
  425.     # These data are used BOTH for 
  426.     #  -- $PA and $Type values for 'catalog of neighboring galaxies' 
  427.     #  -- distances via Hubble law for NGC/IC catalog
  428.     # Hence evaluate all alternative names for now
  429.         
  430.     # Hierarchy: UGC , ESO , MCG , UGCA, CGCG (Zwicky et al.)
  431.     # Read out 2 alternative names:
  432.     
  433.     $altname =  &standardize(&ss(75,89));
  434.     $altname2 = &standardize(&ss(91,104));
  435.     
  436.     # RC3 does not list Messier names! Hence, must translate NGC to Messier names for the 5 M-galaxies used in "catalog of neighboring galaxies"     
  437.     
  438.     $n1 =~ s/NGC 224/M 31/;
  439.     $n1 =~ s/NGC 598/M 33/;
  440.     $n1 =~ s/NGC 3031/M 81/;
  441.     $n1 =~ s/NGC 3034/M 82/;
  442.     $n1 =~ s/NGC 5457/M 101/;
  443.     
  444.     # Read out position angle for the "catalog of neighboring galaxies" below
  445.     
  446.     $PA{$n1} = &ss(186,188); 
  447. $PA{$n1} =~ s/ //g;
  448. # assign to 0 if undefined
  449. $PA{$n1} = 0 if (! $PA{$n1});
  450.     $PA{$altname} = $PA{$n1};
  451.     $PA{$altname2} = $PA{$n1};
  452.     
  453.     # Read out morphological types for the "catalog of neighboring galaxies" below
  454.     
  455.     $type = &ss(118,124);
  456.     $type =~ s/ //g;
  457.     $Type{$n1} = $type;
  458.     $Type{$altname} = $type;
  459.     $Type{$altname2} = $type;
  460.     
  461.     if(0){
  462.      if ($type){
  463.      $tp = &RC3type($type);
  464.      printf "%15s %15s %15s %10s %10s %10sn",$pgc, $n1, $altname, $altname2, $type, $tp;
  465.      }
  466.     }    
  468.     #
  469.     # RC3 distances, [ly]
  470.     #--------------------
  471.     # from Hubble's law, d ~ v_recession/H0,
  472.     # if velocity >~ +80 km/sec (positive!)
  473.     # use weighted mean recession velocity corrected to the reference frame
  474.     # defined by the 3K microwave background radiation (CMB), units [km/sec]
  475.     
  476.     $vrad = &ss(359,363);
  477.     $vrad =~ s/ //g;
  478.     if ($vrad) {
  479.         $n_rc3++ if ($vrad > 80.0);
  480.      $rc3_distance{$n1} = $vrad/$H0 * 1e6 * $pc2ly;
  481.      $v_CMB_RC3{$n1} = $vrad;
  482.     }
  483. }
  484. close (RC3);
  486. # +++++++++++++++++++++++++++++++++++++++++++++++++++
  487. # Read in data from "Catalog of Neighboring Galaxies"
  488. # +++++++++++++++++++++++++++++++++++++++++++++++++++
  490. # names, major diameters, inclinations, absolute B-band magnitudes
  492. open(LOCALV4," < catalogs/neighborgals4.txt")|| die "Can not read neighborgals4.txtn";
  493. $n_local = 0;
  494. $j = 0;
  495. while (<LOCALV4>){
  497. $names = &ss(1,15);
  498. $names = &clean($names);
  499. if ($names =~/,/){
  500. ($name,$name2) = split(/,/,$names);
  501. } else {
  502. $name = $names;
  503. }
  504. $name = &standardize($name);
  505.     
  506.     # it was forgotten to list main IC name first for DDO 226    
  507.     $name =~ s/^DDO 226/IC 1574/;
  509. # skip all NGC/IC and Messier objects in 'catalog of neighboring galaxies'
  510. # this drops 114 objects from the catalog
  512. #if ($name =~ /^Msd+|^NGC|^IC/){print "$j  $namen";$j++;}
  513. next if ($name =~ /^Msd+|^NGC|^IC/); 
  514.     next if ($name =~ /^Milky Way/);
  515. #
  516. # Major diameters  [ly]
  517. # (corrected for galaxy inclination and the galactic extinction)
  518. #
  519. $D = &ss(17,21);
  520. $D =~ s/ //g;
  521. next if (!$D);
  522. $radius{$name} = 0.5 * $D * $pc2ly * 1e3;
  523. #print "$name $radius{$name}n";
  524. #
  525. # Inclinations from face-on [deg]
  526. #
  527. $i{$name} = &ss(23,24);
  528.        
  529. #
  530. # Absolute B-band magnitudes
  531. #
  532. $BMag{$name} = &ss(35,40);
  533. $BMag{$name} =~ s/ //g;
  535. }
  536. close (LOCALV4);
  538. # names, altnames, J2000 coordinates, distances, position angles and types from RC3
  540. open(LOCALV," < catalogs/neighborgals.txt")|| die "Can not read neighborgals.txtn";
  541. $n_local = 0;
  542. $j=0;
  543. $name2="";
  544. while (<LOCALV>){
  545. $names = &ss(1,15);
  546. $names = &clean($names);
  547. if ($names =~/,/){
  548. ($name1,$name2) = split(/,/,$names);
  549. $name2 = &standardize($name2);
  550. } else {
  551. $name1 = $names;
  552.      $name2 = "";
  553. }
  555. $name1 = &standardize($name1);
  556. #print "$name1n" if ($name1 =~/^[A-Z]+[a-z]+s*[a-zA-Z]*s*$/);
  557.     
  558.     # it was forgotten to list main IC name first for DDO 226    
  559.     $name1 =~ s/^DDO 226/IC 1574/;
  561. # skip all NGC/IC and Messier objects in 'catalog of neighboring galaxies'
  562. next if ($name1 =~ /^Msd+|^NGC|^IC/); 
  563.     next if ($name1 =~ /^Milky Way/);
  565. #$DD = &ss(38,42);
  566. #$DD =~ s/ //g;
  567. next if (!$radius{$name1});
  568. #
  569. # Read out J2000 coordinates
  570. #
  571. # RA, DEC fields:
  572. $ra = &clean(&ss(17,26));
  573. $dec= &clean(&ss(28,36));
  575. # convert to hours, degrees, respectively
  576.   
  577.    $alpha = &ra2h ($ra);
  578.    $delta = &dec2deg ($dec);
  580. #
  581. # Implement morphological Type (RC2/RC3)
  582. #
  583. $T = &ss(61,62);
  584. $T =~ s/ //g;
  585. if ($Type{$name1}){
  586. $tp = &RC3type($Type{$name1});
  587. #print "$j $name1 $Type{$name1} $tpn";
  588. #$j++;
  589. } elsif ($Type{$name2} && $name2){
  590. $tp = &RC3type($Type{$name2});
  591. #print "$j $name2 $Type{$name2} $tpn";
  592. #$j++;
  593. } else {
  594. $tp = $HType{$T};
  595. }
  596. if ($name1 =~ /^Ands+IIb/) {$tp = "E3";}
  597. if ($name1 =~ /^Ands+Ib/)  {$tp = "E0";}
  598. if ($name1 =~ /^Antliab/)   {$tp = "E3";}
  599. if ($name1 =~ /^Ands+IIIb/){$tp = "E6";}
  600.     if ($name1 =~ /^Antliab/)   {$tp = "E3";}
  601. if ($name1 =~ /^Cetusb/)    {$tp = "E4";}
  602. if ($name1 =~ /^Sex dSphb/) {$tp = "E4";}
  603. if ($name1 =~ /^Tucanab/)   {$tp = "E5";}
  604. if ($name1 =~ /^Ands+Vb/)  {$tp = "E0";}
  605. if ($name1 =~ /^UMinb/)     {$tp = "E5";}
  606. if ($name1 =~ /^Dracob/)    {$tp = "E3";}
  607. if ($name1 =~ /^Sag dSphb/) {$tp = "E6";}
  608. if ($name1 =~ /^Cass+dSphb/){$tp = "E3";}
  609. if ($name1 =~ /^Pegs+dSphb/){$tp = "E3";}
  610. #
  611. # distance [ly]
  612. #
  613. $dist_local = &ss(88,92)* $pc2ly * 1e6;
  614.     #$distance_local{$name1} = $dist_local;
  615.     #printf "%15s  %15.3fn", $name1,$dist_local;
  616. # method used
  617. $meth = &clean(&ss(95,98));
  619. #
  620. # Implement PA, degrees E of N, from RC3, else PA = 0
  621. # check both names 
  622. if ($PA{$name1}){
  623. $pa = &clean($PA{$name1});
  624. #print "$j $name1 $pan";
  625. #$j++;
  627. } elsif ($PA{$name2} && $name2){
  628.   $pa = &clean($PA{$name2});
  629. #print "$j $name2 $pan";
  630. #$j++;
  631. } elsif ($name1 =~ /^Sag dSphb/){
  633.     $pa = 90.0;
  634. } else {
  635. $pa = 0.0;
  636. }  
  638. #
  639. # Orientation (Axis, Angle)
  640. #
  641. &orientation($alpha,$delta,$tp,$i{$name1},$pa);
  643. #
  644. # recall radius and BMag from other file
  645. #
  647. $Radius = $radius{$name1};
  648. $Mabs = $BMag{$name1};
  651. # NED conform nomenclature:
  652. # ellipticities from http://www.ast.cam.ac.uk/~mike/local_members.html
  653. $name1 =~ s/(^Tucana|^Draco|^Antlia|^Carina|^Orion|^Phoenix|^Cetus)/$1 Dwarf/;
  654. $name1 =~ s/^SagDIG/Sagittarius DIG/;
  655. $name1 =~ s/UMin/UrsaMinor Dwarf/;
  656. $name1 =~ s/ESO 349-31/Sculptor DIG/;
  657. $name2 =~ s/SDIG/ESO 349-31/;
  658. $name1 =~ s/^Sculptors*$/Sculptor dSph/;
  659. $name2 =  "ESO 351-30" if ($name1 =~/Sculptor dSph/);
  660. $name1 =~ s/^Sexb/Sextans/;
  661. $name1 =~ s/^Sagb/Sagittarius/;
  662. $name1 =~ s/^Cass+dSphb/And VII/;
  663. $name2 = "Cassiopeia dSph" if ($name1 =~/And VII/);
  664. $name1 =~ s/^DDO 210/Aquarius Dwarf/;
  665. $name2 =  "DDO 210" if ($name1 =~/Aquarius Dwarf/);
  666. $name1 =~ s/^Pegs+dSph/And VI/;
  667. $name2 = "Pegasus dSph" if ($name1 =~/And VIb/);
  668. $name1 =~ s/^Fornax/Fornax dSph/;
  669. $name1 =~ s/^Pegasus/Pegasus DIG/;
  670.  
  671. $na = "";
  672. if ($name2){$na = ":".$name2;}
  673.   $name1 =~ s/ $//g;
  675. if($lgout){   
  676.      $dkpc = $dist_local /($pc2ly * 1e3);
  677. if ($dkpc < 1400){
  678. printf "%17s %15s; %10s; %9s; %5s; %8.3f;  %6.2fn",$name1,$name2,$ra,$dec,$tp,$dkpc,$Mabs;
  679. }
  680. }
  681. #
  682. # add 333 galaxies of the local group volume, V < (10 Mpc)^3,
  683. # (without Messier, NGC and IC galaxies)
  684. print  DSC "# Local Group volume (< (10 Mpc)^3):n";
  685. print  DSC "Galaxy "$name1$na"n";
  686. print  DSC "{n";
  687. print  DSC "        Type  "$tp"n";
  688. printf DSC "        RA        %10.4fn",$alpha;
  689. printf DSC "        Dec       %10.4fn",$delta;
  690. printf DSC "        Distance  %10.4g",$dist_local;
  691. print  DSC " # method: $methn";
  692. printf DSC "        Radius    %10.4gn",$Radius;
  693. printf DSC "        AbsMag    %10.4gn",$Mabs;
  694. printf DSC "        Axis    [%8.4f %8.4f %8.4f]n",@v;
  695. printf DSC "        Angle    %8.4fn",$angle;
  696. print  DSC "        InfoURL  "http://simbad.u-strasbg.fr/sim-id.pl?Ident=$name1"n";  
  697. print  DSC "}nn";
  698. $count++;
  699. $n_local++;
  700. }
  701. close (LOCALV);
  702. if(!$lgout){
  703. print "n Catalog Input:n---------------n";
  704. printf "%5d Galaxies in Local Group Volume (d < 10 Mpc)n",$n_local;
  705. printf "%5d Averaged distances from NED-1D Distance Data Basen",$n_ned1d;
  706. printf "%5d Surface Brightness Fluctuation (SBF) distance datan",$n_sbf;
  707. printf "%5d Tully-Fisher elliptical galaxy distance datan",$n_tf_e;
  708. printf "%5d Tully-Fisher spiral galaxy distance datan",$n_tf_s;
  709. printf "%5d ZCAT redshift datan",$n_zcat;
  710. printf "%5d RC3 galaxy distance datan",$n_rc3;
  711. printf "%5d NED redshift data (batch)n",$n_ned;
  712. }
  714. $ngcic_flag = 0;        # =1 for NGC, =2 for IC
  715. $c_sbf = 0;
  716. $c_hubble = 0;
  717. $c_tf = 0;
  718. $c_B200 = 0;
  719. $c_zcat = 0;
  720. $c_ned1d = 0;
  721. $c_ned = 0;
  723. # ++++++++++++++++++++++++++++++++++++++++++++++++++++++
  724. # Read out data from revised NGC/IC catalog (Steinecke), 
  725. # merge with above distance data
  726. # +++++++++++++++++++++++++++++++++++++++++++++++++++++
  727. $name = "";
  728. open(NGCIC,"<catalogs/revngcic.txt")|| die "Can not read revngcic.txtn";
  729. while (<NGCIC>) {
  731. # skip header and identify NGC or IC catalogs
  732. if (/^NGC/) {
  733.    $ngcic_flag = 1;
  734.      next;
  735.    } elsif (/^ IC/) {
  736.      $ngcic_flag = 2;
  737.      next;
  738.    } elsif (/^-/) {
  739.      next;
  740.    }
  742. # status ( reading out the S-column)
  743.    $status = &clean(&ss(11,12));
  745. # drop all but galaxies S=1 (drop also  parts of galaxies S=6)
  746.    next if (!$status);
  747.    next if ($status != 1);
  749. # skip extension lines
  750.    next if (&ss(9,10) =~ /2|10/);
  752.     $pgc = &ss(89,94);
  753.     $pgc =~ s/ //g;
  754.     $pgc =~ s/0*(d+)/$1/g;
  755. #
  756. # J2000 coordinates
  757. #
  758.    $ra =  &clean(&ss(21,30));
  759.    $dec = &clean(&ss(32,40));
  761. # Convert coordinates into hours (ra) and degrees (dec):
  763.    $alpha = &ra2h ($ra);
  764.    $delta = &dec2deg ($dec);
  765. #
  766. # Galaxy type, translation to simplified Hubble classes
  767. #
  768.    $Type = &clean(&ss(79,87));
  769.     $HType = &type($Type);
  771. # uncomment as a check on "punch-through" NON-simple Hubble types    
  772. #   next if ($HType =~ /^S[0a-c]$/);
  773. # next if ($HType =~ /^SB[a-c]$/);
  774. # next if ($HType =~ /^E[0-7]$/);
  775. # next if ($HType =~ /^Irr$/);
  777. #   print STDOUT "$name     $Type    $HTypen";
  779. #
  780. # Galaxy names
  781. #
  783. # NGC number
  784.   
  785.    $basename = &clean(&ss(1,5));
  786. #   $basename =~ s/ //g;
  787. #
  788. # Read alternative galaxy names, too
  789. #
  790.    $ID1 = &ss(97,111);
  791. #    $ID1 =~ s/ //g;  
  792.    $ID2 = &ss(113,127);
  793.   $ID2 =~ s/ //g;
  794.   $ID3 = &ss(129,143);
  795.    $ID3 =~ s/ //g;
  796.   
  797.    if    ($ngcic_flag == 1) {$basename = "NGC ".$basename;} 
  798.    elsif ($ngcic_flag == 2) {$basename =  "IC ".$basename;} 
  799.    else                     {next;}
  800.    
  801.     $partflag = &ss(6,6);
  802.    $partflag =~ s/ //g;
  803. # connect the different components (C-column) of an object with a dash  
  804.   
  805.    if ($partflag){
  806.    $name = $basename."-".$partflag;
  807.    #$basename = &standardize($ID1); # so NED (batch) will find it!
  808.     } else {
  809.     $name = $basename;
  810.     }
  812. # NGC 292 = SMC will be taken care of separately
  813.    next if($name =~ /NGC 292/);
  814.   
  815.    $altname1 = "";
  816.    if ($ID1) {
  817. #
  818. # Exchange NGC number by Messier designation from Steinicke's catalog
  819. # find all 40 galaxies from Messier catalog! (checked)
  820. #
  821.    if ( $ID1 =~ /(^Msd{1,3})/) {
  822.      $an1 =  $name;
  823.        $name = $1;
  824.        $basename = $name;
  825.    } else {
  826.      $an1 = &standardize($ID1);
  827.    }
  828.      if ($an1) {
  829.      $altname1 = ":".$an1 ;
  830.      }
  831.    }
  832.    $altname2 = "";
  833.    $an2 = "";
  834.    if ($ID2) {
  835.      $an2      = &standardize($ID2);
  836.        $altname2 = ":".$an2;
  837.    }
  838.    $altname3 = "";
  839.    if ($ID3) {
  840.      $an3      = &standardize($ID3);
  841.      $an3 = "DDO 226" if ($an3 eq "UGCA 9"); 
  842.      $altname3 = ":".$an3;
  843.    }
  844.   
  845.    @nm = ();
  846.    if ($basename) {
  847. #$name =~ s/^(NGC|IC) (d+)[ABCD-]1*$/$1 $2/; 
  848. #$name =~ s/^(?:NGC|IC) d+[ABCD-][234]$//;
  849. #$basename =~ s/-[12]$//;
  850.    #print "$name     $basename" if ($basename);
  851.    #print "          $an1n" if ($an1);
  852.    #print "n" if (!$an1);
  853.        push(@nm,$basename);# if ($basename);
  854.    }
  855.    if ($an1)  {push(@nm,$an1) };
  856.    if ($an2)  {push(@nm,$an2) };
  857.    if ($an3)  {push(@nm,$an3) };
  858. #
  859. # Max/min extensions $D/$d [arcmin] and  inclination $i [degrees]
  860. # ---------------------------------------------------------------
  861. #
  862. # use $D [arcmin] from Steinicke's data
  863. # for all lines, $D entry well-defined
  864. $D = &ss(62,67);
  865.    $D =~ s/ //g;
  866.    $d = &ss(69,73);
  867.    if ($d =~ /[0-9]/) {
  868.      $i = &inclination($d/$D);
  869.     # if $d lacking, equate min size $d to to max size $D.
  870.    } else {
  871.      $d = $D;
  872.      $i = 3.0;
  873.    }
  874. #
  875. # Apparent blue magnitude $Bmag
  876. # ------------------------------------------------
  877. #
  878.    $Bmag = "";                   # $Bmag undefined
  879.    if (&ss(42,47) =~ /[0-9]/) {
  880.      $Bmag = &ss(42,47);
  881.      $Bmag =~ s/ //g;
  882.    }
  883.   
  884. # Distance [ly] from various methods
  885. #-----------------------------------
  886. #
  887.    $distance = "";
  888.    $dist_flag = "";
  889.    # loop over alternate names
  890.    foreach $n (@nm) {
  891.     #
  892.     # first, use the distance-averages from the accurate NED-1D 
  893.     # Extragalactic Distance Database
  894.     # 
  895.     #
  896.      if ($distance_ned1d{$n}) {
  897.          #print "$nn";
  898.        $distance = $distance_ned1d{$n};
  899.        $dist_flag = 6;
  900.        last;    
  901.     #
  902.     # next, exploit SBF-distance data (IV)
  903.     # they are very accurate
  904.     #
  905.      } elsif ($sbf_distance{$n}) {
  906.        $distance = $sbf_distance{$n};
  907.        $dist_flag = 1;
  908.        #print "$nn";
  909.        last;
  910.     #
  911.     # next add in available Tully-Fischer distances
  912.     # for elliptical galaxies
  913.     #
  914.      } elsif ($distance_TFE{$n}) {
  915.        $distance = $distance_TFE{$n};
  916.        $dist_flag = 2;
  917.        #printf "%10s %10.2g %5s %10.2fn",$n,$distance,$HType,$Bmag;
  918.        last;
  919.     #
  920.     # next add in available Tully-Fischer distances
  921.     # for spiral galaxies
  922.     #
  923.     
  924.     # match via PGC numbers
  925.      } elsif ($distance_TFS{$n}) {
  926.        $distance = $distance_TFS{$n};
  927.        #printf "%10s %10.2g %5s %10.2fn",$n,$distance,$HType,$Bmag;
  928.             $dist_flag = 2;
  929.     
  930.     #
  931.     # Next, use distances from "200 brightest galaxies" compilation
  932.     #
  933.     #
  934.     } elsif ($distance_B200{$n}) {
  935.      $distance = $distance_B200{$n};
  936.      $dist_flag = 4;
  937.      $method{4} = "  # method: $method_B200{$n}n";
  938.      #printf "%10s %10s %10.2gn", $name,$method_B200{$name},$distance_B200{$name};
  939.     #
  940.     # next exploit Hubble's law, distance ~ v_recession/H0,
  941.     # if radial velocity >~ 80 [km/sec] (positive!),
  942.     # use weighted mean recession velocity corrected to
  943.     # the reference frame defined by the 3K microwave background
  944.     # radiation (CMB)
  945.     #
  946.     # recession velocities from RC3_rev catalog
  947.     } elsif ($rc3_distance{$n}) {
  948. if ($v_CMB_RC3{$n} > 80.0) {
  949.      $distance = $rc3_distance{$n};
  950.      #print "$pgc $n   $distance  $v_CMB_RC3{$n}n";
  951.      $dist_flag = 3;  
  952.    }  
  953.   
  954.    #
  955.    # Next, use distances from zcat redshift catalog
  956.    #  
  957.         } elsif ($distance_zcat{$n}) {
  958.    if ($distance_zcat{$n} > 0){  
  959.      $distance = $distance_zcat{$n} ;
  960.      #print "$n  $distancen";
  961.      $dist_flag = 5;
  962.    }  
  963.    } elsif ($distance_ned{$n}) {
  964.      $distance = $distance_ned{$n} ;
  965.      #print "$n  $distancen";        
  966.      $dist_flag = 8;
  967.        }
  968.    }
  970.   #
  971.   # Individual corrections from various recent sources
  972.   #
  973.    if($name =~ /M 51/) {
  974.      $distance = $sbf_distance{"NGC 5195"} - 0.5e6;
  975.      $dist_flag = 4;
  976.      $method{4} = "  # method: $method_B200{$name}n";
  977.    }
  978.    if($name =~ /NGC 1569/) {
  979.      $distance = 2.375 * $pc2ly * 1e6;
  980.      $dist_flag = 7;
  981.      $method{7} = "  # method: Pn";
  982.    }
  983.    if($name =~ /M 90/) {
  984.      $distance = $distance_B200{"M 90"};
  985.      $dist_flag = 4;
  986.      $method{4} = "  # method: $method_B200{$name}n";
  987.    }
  988.    if($name =~ /NGC 4438/) {
  989.      $distance = $rc3_distance{"NGC 4435"};
  990.      $dist_flag = 7;
  991.      $method{7} = "  # method: M. Machacek et al. 2003, ApJn";
  992.    }
  993.   
  994.    #
  995.    # last resort for distances: (<=> TF for eta = eta_peak)
  996.    # Simple estimate using /peak/ absolute mags from (absolute) brightness
  997.    # distributions of SBF-galaxy survey, distinguishing "E" and "S"-type 
  998.    # galaxies.
  999.    #
  1000.   
  1001.    # if(!$distance){printf "%10s  %10.3f %4sn",$name,$Bmag,$HType;}
  1002.    #next if (!$distance);
  1003.   
  1004.    if(! $distance){ 
  1005.         
  1006.      $M_abs = $M_abs_ell if ($HType =~ /^E|^Irr/);
  1007.      $M_abs = $M_abs_spir if ($HType =~ /^S/);
  1008.      $distance = 10**(($Bmag - $M_abs * (1 + 0.0 * (1 - 2 * rand())) + 5)/5) * $pc2ly;
  1009.      $dist_flag = 0;
  1010.         #print "$j     $basenamen" if ($basename);
  1011.         $j++;
  1012.      #printf "%10s  Distance = %8.4g ly  HType =%4sn",$name,$distance,$HType;
  1013.    }
  1015.    # Radius [ly]
  1016.    #-------------
  1017.    #
  1018.    $radius = 0.5*deg2rad($D/60.0)*$distance;
  1019.    #
  1020.    # Position angle
  1021.    #---------------
  1022.    #
  1023.    # use values from Steinicke's catalog
  1024.    $PA = &ss(75,77);
  1025.    $PA =~ s/ //g;
  1026.    # assign to 0 if undefined
  1027.    $PA = 0 if (! $PA);
  1028.    #printf "%10s %10.2f %5sn",$name,$PA,$HType if ($HType eq "Sc");
  1030.    #
  1031.    # Orientation (Axis, Angle)
  1032.    #---------------------------
  1033.    #
  1034.    &orientation($alpha,$delta,$HType,$i,$PA);
  1035.    #
  1036.    # absBMag (<== Bmag)
  1037.    #--------------------
  1038.    # use the distance from the galaxy boundary (rather than the center)
  1039.    # matches with Celestia code.
  1040.     #printf "%10s %15.4g %10dn",$name,$distance,$dist_flag if ($name =~ /^Msd+/|$name =~/^NGCs5196/);
  1041.    $Mabs = $Bmag - 5*log(($distance - $radius)/(10*$pc2ly))/log(10);
  1042.     #$DM = $Bmag - $Mabs;
  1043.     #printf "%10.3f  %10.3fn",$Bmag, $DM if ($dist_flag == 2);  
  1044.    # count the various distance methods used
  1045.    $c_sbf++    if ($dist_flag == 1);
  1046.    $c_tf++     if ($dist_flag == 2);
  1047.    $c_hubble++ if ($dist_flag == 3);
  1048.    $c_B200++   if ($dist_flag == 4);
  1049.    $c_zcat++   if ($dist_flag == 5);
  1050.    $c_ned1d++  if ($dist_flag == 6);
  1051.    $c_ned++    if ($dist_flag == 8);
  1052.     
  1053. if($lgout){   
  1054.      $dkpc = $distance /($pc2ly * 1e3);
  1055. if ($dkpc < 1400){
  1056. printf "%17s %15s; %10s; %9s; %5s; %8.3f;  %6.2fn",$name,$an1,$ra,$dec,$HType,$dkpc,$Mabs;
  1057. }
  1058. }
  1060. print  DSC "Galaxy "$name$altname1$altname2$altname3"n";
  1061.    print  DSC "{n";
  1062.    print  DSC "        Type  "$HType"n";
  1063. printf DSC "        RA        %10.4fn",$alpha;
  1064.    printf DSC "        Dec       %10.4fn",$delta;
  1065.    printf DSC "        Distance  %10.4g",$distance;
  1066.    print  DSC "$method{$dist_flag}";
  1067.    printf DSC "        Radius    %10.4gn",$radius;
  1068.    printf DSC "        AbsMag    %10.4gn",$Mabs;
  1069.    printf DSC "        Axis    [%8.4f %8.4f %8.4f]n",@v;
  1070.    printf DSC "        Angle    %8.4fn",$angle;
  1071.    print  DSC "        InfoURL "http://simbad.u-strasbg.fr/sim-id.pl?Ident=$name"n";
  1072.    print  DSC "}nn";
  1073.   
  1074.    $count++;
  1075. }
  1077. # Summary output of used distance data
  1078. $sum = $n_local + $c_ned1d + $c_sbf + $c_tf + $c_hubble + $c_B200 + $c_zcat + $c_ned;
  1079. print  "nSummary Output:n---------------n";
  1080. printf "%5d Galaxies in Totaln",$count;
  1081. printf "%5d with good distance datann",$sum;
  1082. printf "%5d with distances from 'catalog of neighboring galaxies'n",$n_local;
  1083. printf "%5d with NED-1D averaged distancesn",$c_ned1d;
  1084. printf "%5d with Surface Brightness Fluctuation (SBF) distancesn",$c_sbf;
  1085. printf "%5d with Tully-Fisher distancesn",$c_tf;
  1086. printf "%5d with B200 distancesn",$c_B200;
  1087. printf "%5d with distances from Hubble Law (v_CMB, RC3)n",$c_hubble;
  1088. printf "%5d with distances from Hubble Law (ZCAT).n",$c_zcat;
  1089. printf "%5d with distances from Hubble Law (NED, batch).n",$c_ned;
  1091. # like substr($_,first,last), but one-based.
  1092. sub ss {
  1093.    substr ($_, $_[0]-1, $_[1]-$_[0]+1);
  1094. }
  1095. sub clean {
  1096. # squeeze out superfluous spaces
  1097. my($string) = @_;
  1098. $string =~ s/^s*//;
  1099. $string =~ s/s*$//;
  1100. $string =~ s/s+/ /g;
  1101. $string;
  1102. }
  1103. sub ra2h {
  1104.    my($coord) = @_;
  1105.    my $h =  substr($coord,0,2);
  1106.    my $min =  substr($coord,3,2);
  1107.    my $sec = substr($coord,6,3);
  1108.    $h+$min/60+$sec/3600;
  1109. }
  1110. sub dec2deg {
  1111. my($coord) = @_;
  1112.    my $sign = substr($coord,0,1);
  1113.    my $deg  = substr($coord,1,2);
  1114.    my $min  = substr($coord,4,2);
  1115.    my $sec  = substr($coord,7,3);
  1116.    if ($sign eq "+") {
  1117.      $s =1;
  1118.    } else {
  1119.      $s = -1;
  1120.    }
  1121.    $s*($deg+$min/60+$sec/3600);
  1122. }
  1123. sub ra2h_rc3 {
  1124.     my($h,$min,$sec) = @_;
  1125.     $h+$min/60+$sec/3600;
  1126. }
  1127. sub dec2deg_rc3 {
  1128.     my($sign,$deg,$min,$sec) = @_;
  1129.     if ($sign eq "+"){
  1130.      $s =1;
  1131.     } else {
  1132.      $s = -1;
  1133.     }
  1134.     $s*($deg+$min/60+$sec/3600);
  1135. }
  1136. sub standardize {
  1137.    # eliminate prenulls and bring the ned1d galaxy names in sync with RC3
  1138.    # and NGC/ICrevised
  1139.    my( $name ) = @_;
  1140.     $name =~ s/ //g;
  1141. $name =~ s/bI*G//g;
  1142. $name =~ s/b0+//g;
  1143. $name =~ s/([A-Z]+)0+([1-9]+)/$1 $2/g;
  1144.     $name =~ s/(^[A-Z]+)(-*[1-9].*)/$1 $2/;    
  1145.     $name =~ s/^Eb/ESO/;
  1146.     $name =~ s/^Ub/UGC/;
  1147.     $name =~ s/^UAb/UGCA/;
  1148.     $name =~ s/^Db/DDO/;
  1149.     $name =~ s/^Nb/NGC/;
  1150.     $name =~ s/^Ib/IC/;
  1151.     $name =~ s/^ANDIV/And IV/;
  1152.     $name =~ s/(^[A-Z][a-ce-z]+)(dSph)/$1 $2/;
  1153.     $name =~ s/(^[A-Z][a-z]+G*)(I*[VAB]*)/$1 $2/ unless ($name =~ /SagDIG|dSph/);
  1154.     $name =~ s/^ARP/Arp/;
  1155.     $name =~ s/(^[A-Za-z]+)([1-9]+)/$1 $2/;
  1156.     $name =~ s/LGS -3/LGS 3/;
  1157.     $name =~ s/=w+//;
  1158.     $name =~ s/[V]//;
  1159.     #$name =~ s/^Z (d+)/ZWG $1/;
  1160.     $name =~ s/^C (d+)_0*(d+)/CGCG $1-$2/;
  1161.     $name =~ s/^P (d+)/PGC $1/;
  1162.     $name =~ s/^[KK98](d+)/KK $1/;     
  1163.     $name =~ s/Dwarf//;
  1164.     $name =~ s/-dE1//;
  1165.     $name =~ s//39$//;
  1166.     $name =~ s/&$//;
  1167.     $name =~ s/^NPMs1G(.*)/NPM1G $1/;
  1168.     
  1169.     return $name;
  1170. }
  1171. sub cmpname {
  1172.    my( $name ) = @_;
  1173.    my $num ="";
  1174.    if ($name =~ /(^[A-Z]+)s*([0-9-_AB]++*.*d*)/) {
  1175.      my $let = $1;
  1176.      $num = $2;
  1177.      $num =~ s/^-//;
  1178.      $let." ".$num;
  1179.    }
  1180. }
  1181. sub RC3type {
  1182. # Reduce the RC3 morphology scheme to simple Hubble classes
  1183. my($iType) = @_;
  1184.     return "" if (length($iType) != 7); 
  1185.     $oType = substr($iType,1,length($iType) - 3);
  1186.     $oType =~ s/^[IP].*$/Irr/;
  1187.     $oType =~ s/^L.*$/S0/;    
  1188.     $oType =~ s/^RING.*$/SBa/;
  1189.     $oType =~ s/[AX]//;
  1190.     $oType =~ s/^SB.*0/SBa/;
  1191.     $oType =~ s/^E.*(d).*$/E$1/;
  1192.     $oType =~ s/+*?*.*//g;
  1193.     $oType =~ s/^(SB*).*[89]/Irr/;    
  1194.     $oType =~ s/^(SB*).*[567]$/$1c/;
  1195.     $oType =~ s/^(SB*).*[34]$/$1b/;
  1196.     $oType =~ s/^(SB*).*[12]$/$1a/;
  1197.     $oType =~ s/^S$/S0/;
  1198.     $oType =~ s/^E$/E0/;
  1199.     return $oType;
  1200. }
  1201. sub type {
  1202. my($Type) = @_;
  1203.    if ($Type =~ /S[+?0abcdm]?|^[ES]$|SB[?0abcdm]?|d?E[+?0-7]?|Irr|IB?m?|C|RING|Ring|GxyP|P|D|^R/) {
  1204.      $Type =~ s/^E/SB0/S0/;
  1205.      $Type =~ s/^(w)+C/$1/;
  1206.      $Type =~ s//P?R?G?D?b//;
  1207.      $Type =~ s/(w+??)sB?R?M?/$1/;
  1208.      $Type =~ s/^d//;
  1209.      $Type =~ s/(w+)d-(w)/$1$2/;
  1210.      $Type =~ s/^E-?/?S0/S0/;
  1211.      $Type =~ s/^S(B?)[abc]([abc])/S$1$2/;
  1212.      $Type =~ s/SBR???$/SBa/;
  1213.      $Type =~ s/SB0/S0/;
  1214.      $Type =~ s/^E?? ?B?R?$/E0/;
  1215.      $Type =~ s/^S??$/S0/;
  1216.      $Type =~ s/^Sd$/Irr/;
  1217.      $Type =~ s/^S(B?)c?d/S$1c/;
  1218.      $Type =~ s/^IB?m?$/Irr/;
  1219.      $Type =~ s/^Irr B$/Irr/;
  1220.      $Type =~ s/^SB?d?m/Irr/;
  1221.      $Type =~ s/^C/E0/;
  1222.      $Type =~ s/^P??/Irr/;
  1223.      $Type =~ s/GxyP/Irr/;
  1224.      $Type =~ s/^Ring.*[AB]?$/S0/;
  1225.      $Type =~ s/^SB?cm/Irr/;
  1226.      $Type =~ s/S+S/Sa/;
  1227.      $Type =~ s/E+E??/E0/;
  1228.      $Type =~ s/^[cd]?D ?R?/E0/;
  1229.      $Type =~ s/E+*?C???S?/E0/;
  1230.      $Type =~ s/^Rwd$/Irr/;
  1231.      $Type =~ s/^R/S0/;
  1232.      # new for celestia-1.4.1
  1233.      $Type =~ s/S0+S0$/S0/;
  1234.      $Type =~ s/S[0c]?$/S0/;
  1235.      $Type =~ s/Spec$/Irr/;
  1236.      $Type =~ s/c?E0w+$/E0/;
  1237.      $Type =~ s/[34]S$/Sc/;
  1238.      $Type =~ s/5C$/Irr/;
  1239.      $Type =~ s/E0+*3/E0/;
  1240.      $Type =~ s/SB+C?/SBa/;
  1241.     }
  1242. return $Type;
  1243. }
  1244. sub inclination {
  1245.    my( $r ) = @_;
  1246.    # compute inclination angle $i [degrees] from $d/$D extension ratio
  1247.    my $wu = ($r**2 - 0.2**2)/(1.0 - 0.2**2);
  1248.    if ($wu > 1.0e-8) {
  1249.      3.0 + 180/pi*acos(sqrt($wu));
  1250.    } else {
  1251.      90.0;
  1252.    }
  1253. }
  1255. sub orientation {
  1256.    my($ra, $dec, $Type, $i, $PA) = @_;
  1257.    my $decrot = Math::Quaternion::rotation(deg2rad(90 - $dec),1,0,0);
  1258.    my $rarot  = Math::Quaternion::rotation(deg2rad(90 - $ra * 15),0,1,0);
  1259.    my $radecrot = $decrot * $rarot;
  1260.    my $incrot = Math::Quaternion::rotation(deg2rad($i),0,0,1);
  1261. my $pa =  Math::Quaternion::rotation(deg2rad($PA),0,1,0);
  1262.    my $eclipticsrot = Math::Quaternion::rotation(deg2rad($epsilon),1,0,0);
  1263.    #
  1264.    if ($Type =~ /^E|^Milky/) {
  1265.      $rot = $pa * $radecrot * $eclipticsrot;
  1266. } else {
  1267.      $rot = $incrot * $pa *  $radecrot * $eclipticsrot;
  1268.    }
  1269. if ($rot->rotation_angle > pi){
  1270. $rot = $rot->negate;
  1271. }
  1272. $rot   = $rot->normalize;
  1273.    $angle = $rot->rotation_angle * 180/pi;
  1274. @v     = $rot->rotation_axis;
  1275. }
  1277. sub dcmr {
  1278. # calculate the comoving distance appropriate for Hubble's law, 
  1279. # given $z, $Omega_matter=$Wm, $Omega_lambda=$Wlam.
  1280.     my($Wm, $Wlam, $z) = @_;
  1281.     my $Wr = 0.4165/($H0 * $H0); # 3 massless neutrinos to Omega (radiation)
  1282.     my $Wk = 1 - $Wm -$Wr -$Wlam; # flat universe, Omega (total) = 1
  1283.     my $az = 1/($z + 1);
  1284.     my $numpoints = 50; # support points for numerical integration
  1285.     # do integral over $a=1/(1+$z) from $az to 1 in $numpoints steps, 
  1286.     # midpoint rule
  1287.    $dcmr = 0;
  1288.    #print "$Wm $Wlam $Wk  $Wr $azn";
  1289.    for ($i = 0; $i < $numpoints; $i++) {
  1290.     my $a = $az + (1 - $az) * ($i + 0.5)/$numpoints;
  1291.     my $adot = sqrt($Wk + ($Wm/$a) + ($Wr/($a * $a)) + ($Wlam * $a * $a));
  1292.     $dcmr = $dcmr + 1/($a * $adot);
  1293. }
  1294. return $dcmr = $pc2ly * 1e6 * $c/$H0 * (1 - $az) * $dcmr/$numpoints;
  1295. }