Astro

 view release on metacpan or  search on metacpan

Astro/Coord.pm  view on Meta::CPAN

my @ad =(+1.245e-3,    -1.580e-3,   -0.659e-3);

my @em = ([+0.9999256782, -0.0111820611, -0.0048579477],
	  [+0.0111820610, +0.9999374784, -0.0000271765],
	  [+0.0048579479, -0.0000271474, +0.9999881997],
	  [-0.000551,	    -0.238565,     +0.435739],
	  [+0.238514,     -0.002667,     -0.008541],
	  [-0.435623,     +0.012254,     +0.002117]);

my @emi = ([+0.9999256795, +0.0111814828, +0.0048590039,
	    -0.00000242389840, -0.00000002710544, -0.00000001177742],
	   [-0.0111814828, +0.9999374849, -0.0000271771,
	    +0.00000002710544, -0.00000242392702, +0.00000000006585],
	   [-0.0048590040, -0.0000271557, +0.9999881946,
	    +0.00000001177742, +0.00000000006585, -0.00000242404995],
	   [-0.000551,     +0.238509,     -0.435614,
	    +0.99990432,       +0.01118145,       +0.00485852],
	   [-0.238560,     -0.002667,     +0.012254,
	    -0.01118145,       +0.99991613,       -0.00002717],
	   [+0.435730,     -0.008541,     +0.002117,
	    -0.00485852,       -0.00002716,       +0.99996684]);

=item B<pol2r>

  ($x, $y, $z) = pol2r($polar1, $polar2);

Astro/Misc.pm  view on Meta::CPAN

F2      7500    0.77    11.00    45.80    33.21    57.29    57.44    44.70
F0      7520    0.78    11.00    45.86    33.27    57.30    57.45    44.71
F0      8000    0.94    11.03    46.78    34.19    57.52    57.67    44.93
A5      8500    1.11    11.06    47.81    35.22    57.74    57.89    45.16
A4      8630    1.16    11.07    48.22    36.63    57.81    57.96    45.23
A3      8840    1.23    11.08    48.79    36.20    57.91    58.06    45.33
A3      9000    1.27    11.09    49.11    36.53    57.97    58.12    45.39
A2      9070    1.29    11.09    49.27    36.69    58.00    58.15    45.42
A1      9320    1.35    11.10    49.77    37.19    58.09    58.24    45.51
A1      9400    1.37    11.10    49.93    37.34    58.12    58.27    45.54
A0      9600    1.43    11.12    50.24    37.65    58.20    58.35    45.62
B9.5   10000    1.55    11.14    50.85    38.26    58.37    58.52    45.78
B9.5   10500    1.69    11.17    51.39    38.81    58.55    58.70    45.96
B9     10700    1.74    11.17    51.62    39.04    58.62    58.77    46.03
B9     11000    1.79    11.18    51.85    39.26    58.69    58.84    46.10
B9     11500    1.88    11.18    52.26    39.67    58.80    58.95    46.22
B8     12000    1.97    11.19    52.62    40.03    58.92    59.07    46.33
B8     12500    2.06    11.20    52.98    40.39    59.03    59.18    46.44
B8     13000    2.13    11.20    53.29    40.71    59.11    59.26    46.53
B7     13600    2.22    11.21    53.60    41.02    59.22    59.37    46.63
B7     14000    2.30    11.22    53.88    41.30    59.30    59.45    46.71

Astro/Time.pm  view on Meta::CPAN

  }

  my $JULIAN_DAY_J2000          = 2451545.0;
  my $JULIAN_DAYS_IN_CENTURY    = 36525.0;
  my $SOLAR_TO_SIDEREAL         = 1.002737909350795;

  my $a=101.0 + 24110.54841/86400.0;
  my $b=8640184.812866/86400.0;
  my $e=0.093104/86400.0;
  my $d=0.0000062/86400.0;
  my $tu = (int($mjd)-($JULIAN_DAY_J2000-2400000.5))/$JULIAN_DAYS_IN_CENTURY;
  my $sidtim = $a + $tu*($b + $tu*($e - $tu*$d));
  $sidtim -= int($sidtim);
  if ($sidtim < 0.0) {$sidtim += 1.0};
  my $gmst = $sidtim + ($mjd - int($mjd) + $dUT1/86400.0)*$SOLAR_TO_SIDEREAL;
  while ($gmst<0.0) {
    $gmst += 1.0;
  }
  while ($gmst>1.0) {
    $gmst -= 1.0;
  }



( run in 1.294 second using v1.01-cache-2.11-cpan-7add2cbd662 )