#!/usr/bin/env raku

# SOURCE
# https://www.mathworks.com/matlabcentral/mlc-downloads/downloads/submissions/15285/versions/25/previews/geodetic/cal2jd.m/index.html



comment

=begin comment

Matlab function fix description:

"Y = fix(X)" rounds each element of X to the nearest integer toward zero. 
This operation effectively truncates the numbers in X to integers by 
removing the decimal portion of each number:

  + For positive numbers, the behavior of fix is the same as floor.

  + For negative numbers, the behavior of fix is the same as ceil.

=end comment




sub cal2jd($yr, $mn where {0 < $_ < 13}, $dy) is export(:cal2jd) {
    # CAL2JD  Converts calendar date to Julian date using algorithm
    #   from "Practical Ephemeris Calculations" by Oliver Montenbruck
    #   (Springer-Verlag, 1989). Uses astronomical year for B.C. dates
    #   (2 BC = -1 yr). Non-vectorized version. See also DOY2JD, GPS2JD,
    #   JD2CAL, JD2DOW, JD2DOY, JD2GPS, JD2YR, YR2JD.
    # Version: 2011-11-13
    # Usage:   jd=cal2jd(yr,mn,dy)
    # Input:   yr - calendar year (4-digit including century)
    #          mn - calendar month
    #          dy - calendar day (including fractional day)
    # Output:  jd - jJulian date

    # Copyright (c) 2011, Michael R. Craymer
    # All rights reserved.
    # Email: mike@craymer.com

    if $dy < 1 {
        #if (mn == 2 & dy > 29) | (any(mn == [3 5 9 11]) & dy > 30) | (dy > 31)
        if ($mn == 2 && $dy > 29) || (($mn ~~ /3|5|9|11/) && $dy > 30) || ($dy > 31) {
            die "Invalid input day: '$yr/$mn/$dy'";
            return;
        }
    }

    my ($y, $m, $b);
    if $mn > 2 {
        $y = $yr;
        $m = $mn;
    }
    else {
        $y = $yr - 1;
        $m = $mn + 12;
    }

    my $date1 =  4.5 + 31 * (10 + 12 * 1582);  # Last day of Julian calendar (1582.10.04 Noon)
    my $date2 = 15.5 + 31 * (10 + 12 * 1582);  # First day of Gregorian calendar (1582.10.15 Noon)
    my $date  = $dy  + 31 * ($mn + 12 * $yr);

    if $date <= $date1 {
        $b = -2;
    }
    elsif $date >= $date2 {
        #b = fix(y/400) - fix(y/100);
        $b = floor($y/400) - floor($y/100);
    }
    else {
        die "Dates between October 5 & 15, 1582 do not exist, you entered: '$yr/$mn/$dy'";
        return;
    }

    my $jd;
    if $y > 0 {
        $jd = floor(365.25 * $y) + floor(30.6001 * ($m + 1)) + $b + 1720996.5 + $dy;
    }
    else {
        $jd = floor(365.25 * $y - 0.75) + floor(30.6001 * ($m + 1)) + $b + 1720996.5 + $dy;
    }
}


root

my-jd2cal

MathUtils

React

=begin React :component<HeaderCol> :id<menu> 
=para L<Documentation|/doc/introduction>
=para L<Modules| /mods>
=para L<Examples| file:../examples/index.podlite>
=para L<Download |file:../download/index.podlite>
=para L<About|file:../about/index.podlite>
=para 🔍 K<⌘K>/K<ctrl-K>
=end React


=begin React :component<HeaderCol> :id<footer>  
=begin nested :!nested
=item1 B<Language>
    =item2 L<Get started|file:../getting-started/getting-started.podlite>
    =item2 L<Why Raku?|/doc/language/faq#Why-should-I-learn-Raku-What's-so-great-about-it>
    =item2 L<Try Raku|https://glot.io/new/raku>
    =item2 L<Raku cheat sheet|https://github.com/Raku/mu/blob/master/docs/Perl6/Cheatsheet/cheatsheet.txt>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<Raku on Exercism|https://exercism.org/tracks/raku>
    =item2 L<Wikipedia|https://en.wikipedia.org/wiki/Raku_(programming_language)>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>

=end nested
=begin nested :!nested

=item1 B<L<Documentation|/doc>>
    =item2 L<Getting started, Migration guides from other languages, & Tutorials | file:../doc/introduction.podlite>
    =item2 L<Language References|file:../doc/reference.podlite>
    =item2 L<Type Reference| file:../doc/types.podlite>
    =item2 L<Miscellaneous| file:../doc/miscellaneous.podlite>
    =item2 L<FAQs (Frequently Asked Questions)|/doc/language/faq>
    =item2 L<Community|/doc/language/community>
    =item2 L<The list of all documents|file:../doc/index.podlite>

=item1 B<Resources>
    =item2 L<The Raku Guide|https://raku.guide/>
    =item2 L<Books |https://perl6book.com/>
    =item2 L<Rosetta Code | https://www.raku.org/community/rosettacode>

=end nested
=begin nested :!nested

=item1 B<Resoures>
    =item2 L<Download | file:../download/index.podlite>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>

=item1 B<Learning>
    =item2 L<The Raku Guide|https://raku.guide/>
    =item2 L<Wikibook |https://en.wikibooks.org/wiki/Raku_Programming>
    =item2 L<Books |https://perl6book.com/>
    =item2 L<Rosetta Code | https://www.raku.org/community/rosettacode>
    =item2 L<Learn Raku in Y minutes|https://learnxinyminutes.com/docs/raku/>
    =item2 L<Golfing |https://github.com/AlexDaniel/raku-golf-cheatsheet>
=end nested
=begin nested :!nested

=item1 B<Explore>
    =item2 L<Raku Blog Aggregator|https://planet.raku.org/>
    =item2 L<Rakudo Weekly|https://rakudoweekly.blog/>
    =item2 L<The Weekly Challenge |https://perlweeklychallenge.org/>
    =item2 L<Raku Advent Calendar|https://raku-advent.blog/>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>
    =item2 L<empty|#>

=end nested
=end React


=begin React :component<CookieConsent> :id<CookieConsent> :buttonCaption("Got it!")
=para
This website uses cookies for analytics.
=end React


matlab-cal2jd

DateTime::Julian v1.1.0

Authors

License

Dependencies

Test Dependencies

Provides

Documentation