view release on metacpan or  search on metacpan

lib/Astro/Montenbruck/Ephemeris/Planet/Mars.pm  view on Meta::CPAN

package Astro::Montenbruck::Ephemeris::Planet::Mars;

use strict;
use warnings;

use base qw/Astro::Montenbruck::Ephemeris::Planet/;

use Math::Trig qw/:pi/;
use Astro::Montenbruck::Ephemeris::Pert qw /pert/;
use Astro::Montenbruck::MathUtils qw /frac ARCS/;
use Astro::Montenbruck::Ephemeris::Planet qw/$MA/;

our $VERSION = 0.01;

sub new {
    my $class = shift;
    $class->SUPER::new( id => $MA );
}

sub heliocentric {
    my ( $self, $t ) = @_;

    # Mean anomalies of planets in [rad]
    my $m2 = pi2 * frac( 0.1382208 + 162.5482542 * $t );
    my $m3 = pi2 * frac( 0.9926208 + 99.9970236 * $t );
    my $m4 = pi2 * frac( 0.0538553 + 53.1662736 * $t );
    my $m5 = pi2 * frac( 0.0548944 + 8.4290611 * $t );
    my $m6 = pi2 * frac( 0.8811167 + 3.3935250 * $t );

    my ( $dl, $dr, $db ) = ( 0, 0, 0 );    # Corrections in longitude ["],
    my $pert_cb = sub { $dl += $_[0]; $dr += $_[1]; $db += $_[2] };

    # Perturbations by Venus
    my $term = pert(
        T        => $t,
        M        => $m4,
        m        => $m2,
        I_min    => 0,
        I_max    => 7,
        i_min    => -2,
        i_max    => -1,
        callback => $pert_cb
    );

    $term->( 0, -1, 0, -0.01, -0.03, 0.10,  -0.04, 0.00,  0.00 );
    $term->( 1, -1, 0, 0.05,  0.10,  -2.08, 0.75,  0.00,  0.00 );
    $term->( 2, -1, 0, -0.25, -0.57, -2.58, 1.18,  0.05,  -0.04 );
    $term->( 2, -2, 0, 0.02,  0.02,  0.13,  -0.14, 0.00,  0.00 );
    $term->( 3, -1, 0, 3.41,  5.38,  1.87,  -1.15, 0.01,  -0.01 );
    $term->( 3, -2, 0, 0.02,  0.02,  0.11,  -0.13, 0.00,  0.00 );
    $term->( 4, -1, 0, 0.32,  0.49,  -1.88, 1.21,  -0.07, 0.07 );
    $term->( 4, -2, 0, 0.03,  0.03,  0.12,  -0.14, 0.00,  0.00 );
    $term->( 5, -1, 0, 0.04,  0.06,  -0.17, 0.11,  -0.01, 0.01 );
    $term->( 5, -2, 0, 0.11,  0.09,  0.35,  -0.43, -0.01, 0.01 );
    $term->( 6, -2, 0, -0.36, -0.28, -0.20, 0.25,  0.00,  0.00 );
    $term->( 7, -2, 0, -0.03, -0.03, 0.11,  -0.13, 0.00,  -0.01 );

    # Keplerian motion and perturbations by the Earth
    $term = pert(
        T        => $t,
        M        => $m4,
        m        => $m3,
        I_min    => -1,
        I_max    => 16,
        i_min    => -9,
        i_max    => 0,
        callback => $pert_cb
    );

    $term->( 1,  0,  0, -5.32, 38481.97, -141856.04, 0.40,  -6321.67, 1876.89 );
    $term->( 1,  0,  1, -1.12, 37.98,    -138.67,    -2.93, 37.28,    117.48 );
    $term->( 1,  0,  2, -0.32, -0.03,    0.12,       -1.19, 1.04,     -0.40 );
    $term->( 2,  0,  0, 28.28, 2285.80,  -6608.37,   0.00,  -589.35,  174.81 );
    $term->( 2,  0,  1, 1.64,  3.37,     -12.93,     0.00,  2.89,     11.10 );
    $term->( 2,  0,  2, 0.00,  0.00,     0.00,       0.00,  0.10,     -0.03 );
    $term->( 3,  0,  0, 5.31,  189.29,   -461.81,    0.00,  -61.98,   18.53 );
    $term->( 3,  0,  1, 0.31,  0.35,     -1.36,      0.00,  0.25,     1.19 );