in reply to Re^2: Moon phase on historical events
in thread Moon phase on historical events
Likewise and always good to hear from you. I extended your script to attempt to determine what this moon might look like from the points we have defined and their antipodes. Any observation also needs a height which I have set at the height of Boise. I threw in the planets Venus, Saturn, and Jupiter too, as they seem to be part of the big celestial show for the next few days.
I also use Log::Log4perl, because these data would overwhelm STDOUT. The .conf files are from the examples.
My first question is how the azimuth is defined. Astro::Coord is suuuper complicated. Looks like there's even fortran in there. Physicists and mathematicians do things differently in a lot of settings. I can't tell if my results are wrong or untransformed.
The same time is used for 4 different places. (Not that they're gleichzeitig.) I'll put the log and then the source between readmore tags:
2021/10/10 19:11:13 INFO ./3.mp.pl 2021/10/10 19:11:13 INFO 817 meters in the 'Ho on 2021-10-14T03:22:31 +(1634156551 seconds unix-epoch) timezone: +0700 (lat: 43.61, lon: -11 +6.2) 2021/10/10 19:11:13 INFO Moon age: 8.12733852952393 days 2021/10/10 19:11:13 INFO Moon phase: 27.5 % of cycle (birth-to-death) 2021/10/10 19:11:13 INFO Moon's illuminated fraction: 57.9 % of full d +isc 2021/10/10 19:11:13 INFO Object Azimuth Elevation 2021/10/10 19:11:13 INFO Venus 149.440165717993 15.727434443 +4415 2021/10/10 19:11:13 INFO Jupiter 85.3855790615229 -27.009312 +5350183 2021/10/10 19:11:13 INFO Moon 107.660568864375 -18.129253458 +383 2021/10/10 19:11:13 INFO Saturn 99.431806664097 -18.48745351 +68097 2021/10/10 19:11:13 INFO 817 meters on the Scylla on 2021-10-14T03:22: +31 (1634170951 seconds unix-epoch) timezone: Africa/Nairobi (lat: 15. +6325, lon: 38.245833) 2021/10/10 19:11:13 INFO Moon age: 8.30043411740843 days 2021/10/10 19:11:13 INFO Moon phase: 28.1 % of cycle (birth-to-death) 2021/10/10 19:11:13 INFO Moon's illuminated fraction: 59.7 % of full d +isc 2021/10/10 19:11:13 INFO Object Azimuth Elevation 2021/10/10 19:11:13 INFO Venus 168.793996164728 -81.01106162 +58686 2021/10/10 19:11:13 INFO Jupiter 258.080321870357 -14.857584 +1846033 2021/10/10 19:11:13 INFO Moon 251.270738258432 -35.268506488 +5712 2021/10/10 19:11:13 INFO Saturn 256.410482564197 -30.2637811 +706649 2021/10/10 19:11:13 INFO 817 meters in Antipo' 'Ho on 2021-10-14T03:22 +:31 (1634196151 seconds unix-epoch) timezone: -0400 (lat: -43.61, lon +: 63.8) 2021/10/10 19:11:13 INFO Moon age: 8.60197696037111 days 2021/10/10 19:11:13 INFO Moon phase: 29.1 % of cycle (birth-to-death) 2021/10/10 19:11:13 INFO Moon's illuminated fraction: 62.8 % of full d +isc 2021/10/10 19:11:13 INFO Object Azimuth Elevation 2021/10/10 19:11:13 INFO Venus 79.6274128242256 46.927926715 +5178 2021/10/10 19:11:13 INFO Jupiter 127.915109889761 -14.257597 +1913744 2021/10/10 19:11:13 INFO Moon 121.660606406209 0.26418104014 +4833 2021/10/10 19:11:13 INFO Saturn 118.937176589619 -1.57596843 +956774 2021/10/10 19:11:13 INFO 817 meters at Antipodlean Scylla on 2021-10-1 +4T03:22:31 (1634145751 seconds unix-epoch) timezone: +1000 (lat: -15. +6325, lon: -141) 2021/10/10 19:11:13 INFO Moon age: 7.99713693338972 days 2021/10/10 19:11:13 INFO Moon phase: 27.1 % of cycle (birth-to-death) 2021/10/10 19:11:13 INFO Moon's illuminated fraction: 56.5 % of full d +isc 2021/10/10 19:11:13 INFO Object Azimuth Elevation 2021/10/10 19:11:13 INFO Venus 116.80202915781 -4.3309268944 +6216 2021/10/10 19:11:13 INFO Jupiter 185.633623950117 -59.076854 +3475159 2021/10/10 19:11:13 INFO Moon 152.685924835122 -45.404404422 +7405 2021/10/10 19:11:13 INFO Saturn 160.016969111463 -52.8842049 +328507
Source:
#!/usr/bin/perl use v5.030; use warnings; #use Time::Piece; using DateTime instead use Astro::Coord::ECI::Utils 'deg2rad'; use Astro::Coords; use Astro::MoonPhase; use Astro::Telescope; use DateTime; use DateTime::Format::ISO8601; use DateTime::TimeZone; use Data::Dumper; use Log::Log4perl; use Try::Tiny; ###### # Calculate the moon phase of an event in time. # Required is the date (time defaults to 00:00:01) of the event # The moon-phase calculator does not need a location # since moon phase is the same across the planet (more-or-less) # with the usual antipodean side-effects # However, by specifying a timezone or even coordinates of the # event conversion of event's time to UTC is more accurate # (else UTC time is assumed) # author: bliako, modified by aldebaran # date: 2021-10-10 ###### my $debug = 0; my $time_str = "03:22:31"; my $alt = 817; my $date_str = "2021-10-14"; my @events = ( # { # # example event specification # "name" => "a name for reference, required", # "date" => "date of the event in YYYY-MM-DD", # "time" => "optional time as hh:mm:ss, default is 00:00:01", # #specify optional location using one of these for calculating # #epoch which uses UTC time, # #otherwise, UTC timezone will be assumed for event's time: # "timezone" => "standard timezone name, e.g. America/Havana", # "location" => "cuba", # "location" => {lon=>-81.137633, lat=>22.17927}, # }, { "name" => "$alt meters in the 'Ho", "date" => $date_str, "time" => $time_str, "location" => { lon => -116.2, lat => 43.61 }, }, { "name" => "$alt meters on the Scylla", "date" => $date_str, "time" => $time_str, "location" => { lon => 38.245833, lat => 15.6325 }, }, { "name" => "$alt meters in Antipo' 'Ho", "date" => $date_str, "time" => $time_str, "location" => { lon => 63.8, lat => -43.61 }, }, { "name" => "$alt meters at Antipodlean Scylla", "date" => $date_str, "time" => $time_str, "location" => { lon => -141, lat => -15.6325 }, }, ); my $file = '/home/hogan/Documents/hogan/logs/3.log4perl.txt'; # unlink $file or warn "Could not unlink $file: $!"; my $log_conf4 = "/home/hogan/Documents/hogan/logs/conf_files/3.conf"; Log::Log4perl::init($log_conf4); #info my $logger = Log::Log4perl->get_logger(); $logger->info("$0"); ## event loop for my $event (@events) { my $epoch = parse_event($event); $logger->info( event2str($event) ); my ( $MoonPhase, $MoonIllum, $MoonAge, $MoonDist, $MoonAng, $SunDist +, $SunAng ) = phase($epoch); $logger->info("Moon age: $MoonAge days"); $logger->info( "Moon phase: " . sprintf( "%.1f", 100.0 * $MoonPhase ) . " % of cycle (birth-to-death)" ); $logger->info( "Moon's illuminated fraction: ", sprintf( "%.1f", 100.0 * $MoonIllum ) . " % of full disc" ); $logger->info("Object\tAzimuth\t\tElevation"); for my $name (qw/Venus Jupiter Moon Saturn/) { my $planet = Astro::Coords->new( planet => $name ); # use DateTime my $date = DateTime->from_epoch( epoch => $epoch, time_zone => 'UT +C' ); $planet->datetime($date); $planet->telescope( Astro::Telescope->new( Name => $event->{name}, Long => deg2rad( $event->{location}->{lon} ), Lat => deg2rad( $event->{location}->{lat} ), Alt => $alt, #meters above sea level ) ); my $az = $planet->az( format => 'deg' ); my $el = $planet->el( format => 'deg' ); $logger->info("$name\t$az\t\t$el"); } ## end planets loop } ## end event loop sub event2str { my $event = shift; if ( !exists $event->{_is_parsed} ) { warn "event has not been parsed, just dumping it..."; print Dump($event); } my $str = $event->{name} . " on " . $event->{datetime} . " (" . $event->{datetime}->epoch . " seconds unix-epoch)" . " timezone: " . $event->{datetime}->time_zone->name; if ( exists $event->{location} ) { if ( ref( $event->{location} ) eq 'HASH' ) { $str .= " (lat: " . $event->{location}->{lat} . ", lon: " . $event->{location}->{lon} . ")"; } else { $str .= "(" . $event->{location} . ")" } } return $str; } sub parse_event { my $event = shift; if ( !exists $event->{date} ) { die "date field is missing from even +t." } my $datestr = $event->{date}; die "event does not have a 'name' field, please specify one, anythin +g really." unless exists $event->{name}; my $timestr = "00:00:01"; if ( exists $event->{time} ) { $timestr = $event->{time}; print "event2epoch(): setting time to '$timestr' ...\n" if $debug > 0; die "time '$timestr' is not valid, it must be in the form 'hh:mm:s +s'." unless $timestr =~ /^\d{2}:\d{2}:\d{2}$/; } else { $event->{time} = $timestr } my $isostr = $datestr . 'T' . $timestr; my $dt = DateTime::Format::ISO8601->parse_datetime($isostr); die "failed to parse date '$isostr', check date and time fields." unless defined $dt; $event->{datetime} = $dt; my $tzstr = 'UTC'; if ( exists $event->{timezone} ) { $tzstr = $event->{timezone}; print "event2epoch(): found a timezone via 'timezone' field as '$tzstr' (tha +t does not mean it is valid) ...\n" if $debug > 0; } elsif ( exists $event->{location} ) { my $loc = $event->{location}; if ( ( ref($loc) eq '' ) && ( $loc =~ /^[a-zA-Z]$/ ) ) { # we have a location string my @alltzs = DateTime::TimeZone->all_names; my $tzstr; for (@alltzs) { if ( $_ =~ /$loc/i ) { $tzstr = $_; last } } die "event's location can not be converted to a timezone, consider specify +ing the 'timezone' directly or setting 'location' coordinates with: \ +[lat,lon\]." unless $tzstr; print "event2epoch(): setting timezone via 'location' name to '$timestr' ... +\n" if $debug > 0; } elsif ( ( ref($loc) eq 'HASH' ) && ( exists $loc->{lat} ) && ( exists $loc->{lon} ) ) { # we have a [lat,lon] array for location require Geo::Location::TimeZone; my $gltzobj = Geo::Location::TimeZone->new(); $tzstr = $gltzobj->lookup( lat => $loc->{lat}, lon => $loc->{lon +} ); if ( !$tzstr ) { die "timezone lookup from location coordinates lat:" . $loc->{lat} . ", lon:" . $loc->{lon} . " has failed."; } print "event2epoch(): setting timezone via 'location' coordinate +s lat:" . $loc->{lat} . ", lon:" . $loc->{lon} . " ...\n" if $debug > 0; } } if ($tzstr) { print "event2epoch(): deduced timezone to '$tzstr' and setting it +...\n" if $debug > 0; try { $dt->set_time_zone($tzstr) } catch { die "$_\n failed to set the timezone '$tzstr', is it valid?" } } $event->{_is_parsed} = 1; $event->{epoch} = $dt->epoch; return $event->{epoch}; } __END__
One thing I noticed is that these data are very sensitive to altitude of the defined observer. A person at 817 meters above the Scylla would have a much better view of rising and setting planets. Is there any way to figure out what a reasonable guesstimate is of altitude given latitude and longitude?
How long does the moon's transit from Venus to Jupiter take?
How can I determine that a priori? I know that I can bound the value in DateTime with output from the original post. Can one loop over this interval with the methods of DateTime adding the duration of a second and then checking values. A month isn't a lot of seconds for a computer.
|
|---|
| Replies are listed 'Best First'. | |
|---|---|
|
Re^4: Moon phase on historical events
by bliako (Abbot) on Oct 11, 2021 at 11:51 UTC | |
|
Re^4: Moon phase on historical events
by Anonymous Monk on Feb 11, 2022 at 07:26 UTC | |
by bliako (Abbot) on Feb 11, 2022 at 10:29 UTC | |
by Aldebaran (Curate) on Feb 15, 2022 at 18:06 UTC |