mirror of
https://github.com/fhem/fhem-mirror.git
synced 2025-01-31 06:39:11 +00:00
First version
git-svn-id: https://svn.fhem.de/fhem/trunk@86 2b470e98-0d58-463d-a4d8-8e2adae1ed80
This commit is contained in:
rudolfkoenig
parent
4f84681f66
commit
e4edfac7c1
@ -341,3 +341,4 @@
|
||||
|
||||
- ==DATE== (4.2)
|
||||
- feature: added archivedir/archivecmd to the the main logfile
|
||||
- feature: 99_Sunrise_EL.pm (does not need any Date modules)
|
||||
|
||||
333
fhem/FHEM/99_SUNRISE_EL.pm
Executable file
333
fhem/FHEM/99_SUNRISE_EL.pm
Executable file
@ -0,0 +1,333 @@
|
||||
##############################################
|
||||
# - Please call sunrise_coord before using this module, else you'll get times
|
||||
# for frankfurt am main (germany). See the "at" entry in commandref.html
|
||||
#
|
||||
# This code is derived from DateTime::Event::Sunrise, version 0.0501.
|
||||
# Simplified and removed further package # dependency (DateTime,
|
||||
# Params::Validate, etc). For comments see the original code.
|
||||
#
|
||||
|
||||
package main;
|
||||
use strict;
|
||||
use warnings;
|
||||
use Math::Trig;
|
||||
|
||||
sub sr($$$$);
|
||||
sub sunrise_rel(@);
|
||||
sub sunset_rel(@);
|
||||
sub sunrise_abs(@);
|
||||
sub sunset_abs(@);
|
||||
sub isday();
|
||||
sub sunrise_coord($$$);
|
||||
sub SUNRISE_Initialize($);
|
||||
|
||||
# See perldoc DateTime::Event::Sunrise for details
|
||||
my $long = "8.686";
|
||||
my $lat = "50.112";
|
||||
my $tz = ""; # will be overwritten
|
||||
my $altit = "-6"; # Civil twilight
|
||||
my $RADEG = ( 180 / pi );
|
||||
my $DEGRAD = ( pi / 180 );
|
||||
my $INV360 = ( 1.0 / 360.0 );
|
||||
|
||||
|
||||
sub
|
||||
SUNRISE_EL_Initialize($)
|
||||
{
|
||||
my ($hash) = @_;
|
||||
}
|
||||
|
||||
##########################
|
||||
# Compute:
|
||||
# rise: 1: event is sunrise (else sunset)
|
||||
# isrel: 1: _relative_ times until the next event (else absolute for today)
|
||||
# seconds: second offset to event
|
||||
# daycheck: if set, then return 1 if the sun is visible, 0 else
|
||||
sub
|
||||
sr($$$$)
|
||||
{
|
||||
my ($rise, $seconds, $isrel, $daycheck) = @_;
|
||||
my $needrise = ($rise || $daycheck) ? 1 : 0;
|
||||
my $needset = (!$rise || $daycheck) ? 1 : 0;
|
||||
|
||||
my $nt = time;
|
||||
my @lt = localtime($nt);
|
||||
my $gmtoff = _calctz($nt,@lt); # in hour
|
||||
my ($rt,$st) = _sr($needrise,$needset, $lt[5]+1900,$lt[4]+1,$lt[3], $gmtoff);
|
||||
|
||||
my $nh = $lt[2] + $lt[1]/60 + $lt[0]/3600;
|
||||
if($daycheck) {
|
||||
return 0 if($nh < $rt || $nh > $st);
|
||||
return 1;
|
||||
}
|
||||
|
||||
$seconds = 0 if(!$seconds);
|
||||
my $sst = ($rise ? $rt : $st);
|
||||
if(!$isrel) {
|
||||
$sst += ($seconds/3600);
|
||||
return h2hms_fmt($sst);
|
||||
}
|
||||
|
||||
$sst += ($seconds/3600);
|
||||
|
||||
my $diff = 0;
|
||||
if($nh > $sst) {
|
||||
$nt += 86400; # Tommorow
|
||||
$diff = 24;
|
||||
@lt = localtime($nt);
|
||||
$gmtoff = _calctz($nt,@lt); # in hour
|
||||
|
||||
($rt,$st) = _sr($needrise,$needset, $lt[5]+1900,$lt[4]+1,$lt[3], $gmtoff);
|
||||
|
||||
$sst = ($rise ? $rt : $st);
|
||||
$sst += ($seconds/3600);
|
||||
}
|
||||
|
||||
$diff = $diff + $sst - $nh;
|
||||
return h2hms_fmt($diff);
|
||||
}
|
||||
|
||||
sub
|
||||
_sr($$$$$$)
|
||||
{
|
||||
my ($needrise, $needset, $y, $m, $dy, $offset) = @_;
|
||||
|
||||
my $d = _days_since_2000_Jan_0($y,$m,$dy) + 0.5 - $long / 360.0;
|
||||
my ( $tmp_rise_1, $tmp_set_1 ) =
|
||||
_sunrise_sunset( $d, $long, $lat, $altit, 15.04107 );
|
||||
|
||||
my ($tmp_rise_2, $tmp_rise_3) = (0,0);
|
||||
|
||||
if($needrise) {
|
||||
$tmp_rise_2 = 9; $tmp_rise_3 = 0;
|
||||
until ( _equal( $tmp_rise_2, $tmp_rise_3, 8 ) ) {
|
||||
|
||||
my $d_sunrise_1 = $d + $tmp_rise_1 / 24.0;
|
||||
( $tmp_rise_2, undef ) =
|
||||
_sunrise_sunset( $d_sunrise_1, $long,
|
||||
$lat, $altit, 15.04107 );
|
||||
$tmp_rise_1 = $tmp_rise_3;
|
||||
my $d_sunrise_2 = $d + $tmp_rise_2 / 24.0;
|
||||
( $tmp_rise_3, undef ) =
|
||||
_sunrise_sunset( $d_sunrise_2, $long,
|
||||
$lat, $altit, 15.04107 );
|
||||
|
||||
}
|
||||
}
|
||||
|
||||
my ($tmp_set_2, $tmp_set_3) = (0,0);
|
||||
if($needset) {
|
||||
$tmp_set_2 = 9; $tmp_set_3 = 0;
|
||||
until ( _equal( $tmp_set_2, $tmp_set_3, 8 ) ) {
|
||||
|
||||
my $d_sunset_1 = $d + $tmp_set_1 / 24.0;
|
||||
( undef, $tmp_set_2 ) =
|
||||
_sunrise_sunset( $d_sunset_1, $long,
|
||||
$lat, $altit, 15.04107 );
|
||||
$tmp_set_1 = $tmp_set_3;
|
||||
my $d_sunset_2 = $d + $tmp_set_2 / 24.0;
|
||||
( undef, $tmp_set_3 ) =
|
||||
_sunrise_sunset( $d_sunset_2, $long,
|
||||
$lat, $altit, 15.04107 );
|
||||
|
||||
}
|
||||
}
|
||||
|
||||
return $tmp_rise_3+$offset, $tmp_set_3+$offset;
|
||||
}
|
||||
|
||||
|
||||
|
||||
sub
|
||||
_sunrise_sunset($$$$$)
|
||||
{
|
||||
my ( $d, $lon, $lat, $altit, $h ) = @_;
|
||||
|
||||
my $sidtime = _revolution( _GMST0($d) + 180.0 + $lon );
|
||||
|
||||
my ( $sRA, $sdec ) = _sun_RA_dec($d);
|
||||
my $tsouth = 12.0 - _rev180( $sidtime - $sRA ) / $h;
|
||||
my $sradius = 0.2666 / $sRA;
|
||||
|
||||
$altit -= $sradius;
|
||||
|
||||
# Compute the diurnal arc that the Sun traverses to reach
|
||||
# the specified altitude altit:
|
||||
|
||||
my $cost =
|
||||
( sind($altit) - sind($lat) * sind($sdec) ) /
|
||||
( cosd($lat) * cosd($sdec) );
|
||||
|
||||
my $t;
|
||||
if ( $cost >= 1.0 ) {
|
||||
$t = 0.0; # Sun always below altit
|
||||
}
|
||||
elsif ( $cost <= -1.0 ) {
|
||||
$t = 12.0; # Sun always above altit
|
||||
}
|
||||
else {
|
||||
$t = acosd($cost) / 15.0; # The diurnal arc, hours
|
||||
}
|
||||
|
||||
# Store rise and set times - in hours UT
|
||||
|
||||
my $hour_rise_ut = $tsouth - $t;
|
||||
my $hour_set_ut = $tsouth + $t;
|
||||
return ( $hour_rise_ut, $hour_set_ut );
|
||||
|
||||
}
|
||||
|
||||
sub
|
||||
_GMST0($)
|
||||
{
|
||||
my ($d) = @_;
|
||||
|
||||
my $sidtim0 =
|
||||
_revolution( ( 180.0 + 356.0470 + 282.9404 ) +
|
||||
( 0.9856002585 + 4.70935E-5 ) * $d );
|
||||
return $sidtim0;
|
||||
|
||||
}
|
||||
|
||||
sub
|
||||
_sunpos($)
|
||||
{
|
||||
my ($d) = @_;
|
||||
|
||||
my $Mean_anomaly_of_sun = _revolution( 356.0470 + 0.9856002585 * $d );
|
||||
my $Mean_longitude_of_perihelion = 282.9404 + 4.70935E-5 * $d;
|
||||
my $Eccentricity_of_Earth_orbit = 0.016709 - 1.151E-9 * $d;
|
||||
|
||||
# Compute true longitude and radius vector
|
||||
my $Eccentric_anomaly =
|
||||
$Mean_anomaly_of_sun + $Eccentricity_of_Earth_orbit * $RADEG *
|
||||
sind($Mean_anomaly_of_sun) *
|
||||
( 1.0 + $Eccentricity_of_Earth_orbit * cosd($Mean_anomaly_of_sun) );
|
||||
|
||||
my $x = cosd($Eccentric_anomaly) - $Eccentricity_of_Earth_orbit;
|
||||
|
||||
my $y =
|
||||
sqrt( 1.0 - $Eccentricity_of_Earth_orbit * $Eccentricity_of_Earth_orbit )
|
||||
* sind($Eccentric_anomaly);
|
||||
|
||||
my $Solar_distance = sqrt( $x * $x + $y * $y ); # Solar distance
|
||||
my $True_anomaly = atan2d( $y, $x ); # True anomaly
|
||||
|
||||
my $True_solar_longitude =
|
||||
$True_anomaly + $Mean_longitude_of_perihelion; # True solar longitude
|
||||
|
||||
if ( $True_solar_longitude >= 360.0 ) {
|
||||
$True_solar_longitude -= 360.0; # Make it 0..360 degrees
|
||||
}
|
||||
|
||||
return ( $Solar_distance, $True_solar_longitude );
|
||||
}
|
||||
|
||||
sub
|
||||
_sun_RA_dec($)
|
||||
{
|
||||
my ($d) = @_;
|
||||
|
||||
my ( $r, $lon ) = _sunpos($d);
|
||||
|
||||
my $x = $r * cosd($lon);
|
||||
my $y = $r * sind($lon);
|
||||
|
||||
my $obl_ecl = 23.4393 - 3.563E-7 * $d;
|
||||
|
||||
my $z = $y * sind($obl_ecl);
|
||||
$y = $y * cosd($obl_ecl);
|
||||
|
||||
my $RA = atan2d( $y, $x );
|
||||
my $dec = atan2d( $z, sqrt( $x * $x + $y * $y ) );
|
||||
|
||||
return ( $RA, $dec );
|
||||
}
|
||||
|
||||
sub
|
||||
_days_since_2000_Jan_0($$$)
|
||||
{
|
||||
my ($y, $m, $d) = @_;
|
||||
my @mn = (31,28,31,30,31,30,31,31,30,31,30,31);
|
||||
|
||||
my $ms = 0;
|
||||
for(my $i = 0; $i < $m-1; $i++) {
|
||||
$ms += $mn[$i];
|
||||
}
|
||||
my $x = ($y-2000)*365.25 + $ms + $d;
|
||||
$x++ if($m > 2 && ($y%4) == 0);
|
||||
return int($x);
|
||||
}
|
||||
|
||||
sub sind($) { sin( ( $_[0] ) * $DEGRAD ); }
|
||||
sub cosd($) { cos( ( $_[0] ) * $DEGRAD ); }
|
||||
sub tand($) { tan( ( $_[0] ) * $DEGRAD ); }
|
||||
sub atand($) { ( $RADEG * atan( $_[0] ) ); }
|
||||
sub asind($) { ( $RADEG * asin( $_[0] ) ); }
|
||||
sub acosd($) { ( $RADEG * acos( $_[0] ) ); }
|
||||
sub atan2d($) { ( $RADEG * atan2( $_[0], $_[1] ) ); }
|
||||
|
||||
sub
|
||||
_revolution($)
|
||||
{
|
||||
my $x = $_[0];
|
||||
return ( $x - 360.0 * int( $x * $INV360 ) );
|
||||
}
|
||||
|
||||
sub
|
||||
_rev180($)
|
||||
{
|
||||
my ($x) = @_;
|
||||
return ( $x - 360.0 * int( $x * $INV360 + 0.5 ) );
|
||||
}
|
||||
|
||||
sub
|
||||
_equal($$$)
|
||||
{
|
||||
my ( $A, $B, $dp ) = @_;
|
||||
return sprintf( "%.${dp}g", $A ) eq sprintf( "%.${dp}g", $B );
|
||||
}
|
||||
|
||||
sub
|
||||
_calctz($@)
|
||||
{
|
||||
my ($nt,@lt) = @_;
|
||||
|
||||
my $off = $lt[2]*3600+$lt[1]*60+$lt[0];
|
||||
$off = 12*3600-$off;
|
||||
$nt += $off; # This is noon, localtime
|
||||
|
||||
my @gt = gmtime($nt);
|
||||
|
||||
return (12-$gt[2]);
|
||||
}
|
||||
|
||||
|
||||
sub
|
||||
h2hms($)
|
||||
{
|
||||
my ($in) = @_;
|
||||
my ($h,$m,$s);
|
||||
$h = int($in);
|
||||
$m = int(60*($in-$h));
|
||||
$s = int(3600*($in-$h)-60*$m);
|
||||
return ($h, $m, $s);
|
||||
}
|
||||
|
||||
sub
|
||||
h2hms_fmt($)
|
||||
{
|
||||
my ($in) = @_;
|
||||
my ($h,$m,$s) = h2hms($in);
|
||||
return sprintf("%02d:%02d:%02d", $h, $m, $s);
|
||||
}
|
||||
|
||||
|
||||
sub sunrise_rel(@) { return sr(1, shift, 1, 0) }
|
||||
sub sunset_rel(@) { return sr(0, shift, 1, 0) }
|
||||
sub sunrise_abs(@) { return sr(1, shift, 0, 0) }
|
||||
sub sunset_abs(@) { return sr(0, shift, 0, 0) }
|
||||
sub isday() { return sr(1, 0, 0, 1) }
|
||||
sub sunrise_coord($$$) { ($long, $lat, $tz) = @_; return undef; }
|
||||
|
||||
1;
|
||||
@ -134,3 +134,7 @@
|
||||
|
||||
- Rudi Wed Aug 29 08:28:34 MEST 2007
|
||||
- archive attributes clarified in the doc
|
||||
|
||||
- Rudi Mon Sep 3 15:47:59 MEST 2007
|
||||
- 99_Sunrise_EL.pm checked in. Replaces 99_Sunrise.pm, and does not need
|
||||
any Date module.
|
||||
|
||||
@ -1,4 +1,4 @@
|
||||
#Thu Aug 23 17:16:14 2007
|
||||
#Sun Sep 2 16:38:48 2007
|
||||
setstate FHZ fhtbuf: 1c
|
||||
setstate FHZ 2006-02-12 14:03:39 fhtbuf 23
|
||||
setstate FHZ 2006-03-26 08:47:36 init2 deadbeefdeadbe
|
||||
@ -68,11 +68,11 @@ setstate wz 2006-03-26 08:52:31 unknown_85 4
|
||||
setstate wz 2006-03-26 08:50:36 wed-from1 06:00
|
||||
setstate wz 2006-03-26 08:50:36 wed-to1 23:00
|
||||
setstate wz 2006-03-26 08:52:31 windowopen-temp 12.0 (Celsius)
|
||||
setstate x 13.66 0.0 30 13.62 15.09 -0.2 2.8
|
||||
setstate x 2007-08-23 16:41:19 Is 0.0
|
||||
setstate x 2007-08-23 16:41:19 Temp 30
|
||||
setstate x 2007-08-23 16:41:19 Vs 13.66
|
||||
setstate x 2007-08-23 16:41:19 maxI 2.8
|
||||
setstate x 2007-08-23 16:41:19 maxV 15.09
|
||||
setstate x 2007-08-23 16:41:19 minI -0.2
|
||||
setstate x 2007-08-23 16:41:19 minV 13.62
|
||||
setstate x 13.66 0.0 30 13.62 15.09 -0.2 2.8
|
||||
setstate x 2007-09-02 10:06:47 Is 0.0
|
||||
setstate x 2007-09-02 10:06:47 Temp 30
|
||||
setstate x 2007-09-02 10:06:47 Vs 13.66
|
||||
setstate x 2007-09-02 10:06:47 maxI 2.8
|
||||
setstate x 2007-09-02 10:06:47 maxV 15.09
|
||||
setstate x 2007-09-02 10:06:47 minI -0.2
|
||||
setstate x 2007-09-02 10:06:47 minV 13.62
|
||||
|
||||
Loading…
Reference in New Issue
Block a user