]> dxcluster.org Git - dweather.git/blob - loop.pl
prefer rename over copy for loop_data shuffle
[dweather.git] / loop.pl
1 #!/usr/bin/perl
2 use strict;
3
4 use v5.10.1;
5
6 use Mojolicious::Lite;
7 use Serial;
8 use Mojo::IOLoop;
9 use Mojo::IOLoop::Stream;
10 use Mojo::Transaction::WebSocket;
11 #use Mojo::JSON  qw(decode_json encode_json);
12 use JSON;
13 use Debug;
14 use SMGLog;
15 use Math::Round qw(nearest);
16 use File::Copy;
17 use Data::Random qw(rand_chars);
18
19 use constant pi => 3.14159265358979; 
20
21 my $devname = "/dev/davis";
22 my $datafn = ".loop_data";
23 my $dataf;
24
25 my $poll_interval = 2.5;
26 my $rain_mult = 0.2;                            # 0.1 or 0.2 mm or 0.01 inches
27 my $tid;
28 my $rid;
29 our $did;
30 my $nlcount;
31 my $state = "ready";
32 my $buf;
33 my $dbg;
34 my $ser;                                                        # the serial port Mojo::IOLoop::Stream
35 my $last_min_h;
36 my $last_hour_h;
37
38 our $json = JSON->new->canonical(1);
39 our $WS = {};                                   # websocket connections
40
41 our $ld = {};
42
43 our $loop_count;                                # how many LOOPs we have done, used as start indicator
44
45 our     @crc_table = (
46                                   0x0, 0x1021, 0x2042, 0x3063, 0x4084, 0x50a5, 0x60c6, 0x70e7,
47                                   0x8108, 0x9129, 0xa14a, 0xb16b, 0xc18c, 0xd1ad, 0xe1ce, 0xf1ef,
48                                   0x1231, 0x210, 0x3273, 0x2252, 0x52b5, 0x4294, 0x72f7, 0x62d6,
49                                   0x9339, 0x8318, 0xb37b, 0xa35a, 0xd3bd, 0xc39c, 0xf3ff, 0xe3de,
50                                   0x2462, 0x3443, 0x420, 0x1401, 0x64e6, 0x74c7, 0x44a4, 0x5485,
51                                   0xa56a, 0xb54b, 0x8528, 0x9509, 0xe5ee, 0xf5cf, 0xc5ac, 0xd58d,
52                                   0x3653, 0x2672, 0x1611, 0x630, 0x76d7, 0x66f6, 0x5695, 0x46b4,
53                                   0xb75b, 0xa77a, 0x9719, 0x8738, 0xf7df, 0xe7fe, 0xd79d, 0xc7bc,
54                                   0x48c4, 0x58e5, 0x6886, 0x78a7, 0x840, 0x1861, 0x2802, 0x3823,
55                                   0xc9cc, 0xd9ed, 0xe98e, 0xf9af, 0x8948, 0x9969, 0xa90a, 0xb92b,
56                                   0x5af5, 0x4ad4, 0x7ab7, 0x6a96, 0x1a71, 0xa50, 0x3a33, 0x2a12,
57                                   0xdbfd, 0xcbdc, 0xfbbf, 0xeb9e, 0x9b79, 0x8b58, 0xbb3b, 0xab1a,
58                                   0x6ca6, 0x7c87, 0x4ce4, 0x5cc5, 0x2c22, 0x3c03, 0xc60, 0x1c41,
59                                   0xedae, 0xfd8f, 0xcdec, 0xddcd, 0xad2a, 0xbd0b, 0x8d68, 0x9d49,
60                                   0x7e97, 0x6eb6, 0x5ed5, 0x4ef4, 0x3e13, 0x2e32, 0x1e51, 0xe70,
61                                   0xff9f, 0xefbe, 0xdfdd, 0xcffc, 0xbf1b, 0xaf3a, 0x9f59, 0x8f78,
62                                   0x9188, 0x81a9, 0xb1ca, 0xa1eb, 0xd10c, 0xc12d, 0xf14e, 0xe16f,
63                                   0x1080, 0xa1, 0x30c2, 0x20e3, 0x5004, 0x4025, 0x7046, 0x6067,
64                                   0x83b9, 0x9398, 0xa3fb, 0xb3da, 0xc33d, 0xd31c, 0xe37f, 0xf35e,
65                                   0x2b1, 0x1290, 0x22f3, 0x32d2, 0x4235, 0x5214, 0x6277, 0x7256,
66                                   0xb5ea, 0xa5cb, 0x95a8, 0x8589, 0xf56e, 0xe54f, 0xd52c, 0xc50d,
67                                   0x34e2, 0x24c3, 0x14a0, 0x481, 0x7466, 0x6447, 0x5424, 0x4405,
68                                   0xa7db, 0xb7fa, 0x8799, 0x97b8, 0xe75f, 0xf77e, 0xc71d, 0xd73c,
69                                   0x26d3, 0x36f2, 0x691, 0x16b0, 0x6657, 0x7676, 0x4615, 0x5634,
70                                   0xd94c, 0xc96d, 0xf90e, 0xe92f, 0x99c8, 0x89e9, 0xb98a, 0xa9ab,
71                                   0x5844, 0x4865, 0x7806, 0x6827, 0x18c0, 0x8e1, 0x3882, 0x28a3,
72                                   0xcb7d, 0xdb5c, 0xeb3f, 0xfb1e, 0x8bf9, 0x9bd8, 0xabbb, 0xbb9a,
73                                   0x4a75, 0x5a54, 0x6a37, 0x7a16, 0xaf1, 0x1ad0, 0x2ab3, 0x3a92,
74                                   0xfd2e, 0xed0f, 0xdd6c, 0xcd4d, 0xbdaa, 0xad8b, 0x9de8, 0x8dc9,
75                                   0x7c26, 0x6c07, 0x5c64, 0x4c45, 0x3ca2, 0x2c83, 0x1ce0, 0xcc1,
76                                   0xef1f, 0xff3e, 0xcf5d, 0xdf7c, 0xaf9b, 0xbfba, 0x8fd9, 0x9ff8,
77                                   0x6e17, 0x7e36, 0x4e55, 0x5e74, 0x2e93, 0x3eb2, 0xed1, 0x1ef0
78                                  );
79
80 our %bar_trend;
81
82 $bar_trend{-60} = "Falling Rapidly";
83 $bar_trend{196} = "Falling Rapidly";
84 $bar_trend{-20} = "Falling Slowly";
85 $bar_trend{236} = "Falling Slowly";
86 $bar_trend{0} = "Steady";
87 $bar_trend{20} = "Rising Slowly";
88 $bar_trend{60} = "Rising Rapidly";
89
90 our $ending = 0;
91
92 $SIG{TERM} = $SIG{INT} = sub {++$ending; Mojo::IOLoop->stop;};
93 $SIG{HUP} = 'IGNORE';
94
95
96 # WebSocket weather service
97 websocket '/weather' => sub {
98   my $c = shift;
99   my $msg = shift;
100   my $tx = $c->tx;
101
102   # Opened
103   app->log->debug('WebSocket opened.');
104   dbg 'WebSocket opened' if isdbg 'chan';
105   $WS->{$tx} = $tx; 
106
107   # send historical data
108   $c->send($ld->{lasthour_h}) if exists $ld->{lasthour_h};
109   $c->send($ld->{lastmin_h}) if exists $ld->{lastmin_h};
110   
111   # disable timeout
112   $c->inactivity_timeout(3615);
113  
114   # Incoming message
115   $c->on(
116                  message => sub {
117                          my ($c, $msg) = @_;
118                          dbg "websocket: text $msg" if isdbg 'chan';
119                  },
120                  json => sub {
121                          my ($c, $msg) = @_;
122                          dbg "websocket: json $msg" if isdbg 'chan';
123                  }
124   );
125  
126   # Closed
127   $c->on(finish => sub {
128     my ($c, $code, $reason) = @_;
129     app->log->debug("WebSocket closed with status $code.");
130         dbg 'webwocket closed with status $code' if isdbg 'chan';
131         delete $WS->{$tx};
132   });
133 };
134
135 get '/' => {template => 'index'};
136
137
138 dbginit();
139 if (@ARGV) {
140         dbgadd(@ARGV);
141
142 dbgadd('chan');
143
144 dbg '***';
145 dbg "*** starting $0";
146 dbg '***';
147
148 our $dlog = SMGLog->new("day");
149 dbg "before next tick";
150 Mojo::IOLoop->next_tick(sub { loop() });        
151 dbg "before app start";
152 app->start;
153 dbg "after app start";
154
155 write_ld();
156 close $dataf if $dataf;
157
158 # move all the files along one 
159 cycle_loop_data_files();
160
161 dbg '***';
162 dbg "*** ending $0";
163 dbg '***';
164
165 exit 0;
166
167 ##################################################################################
168         
169 sub loop
170 {
171
172         open $dataf, "+>>", $datafn or die "cannot open $datafn $!";
173         $dataf->autoflush(1);
174         
175         read_ld();
176         
177         dbg "last_min: " . scalar gmtime($ld->{last_min});
178         dbg "last_hour: " . scalar gmtime($ld->{last_hour});
179         
180         $did = Mojo::IOLoop->recurring(1 => sub {$dlog->flushall});
181         
182         do_reopen($devname);
183 }
184
185
186
187 sub on_read
188 {
189         my ($ser, $d) = @_;
190         $buf .= $d;
191         $d =~ s/([\%\x00-\x1f\x7f-\xff])/sprintf("%%%02X", ord($1))/eg; 
192         dbg "read added '$d' buf lth=" . length $buf if isdbg 'raw';
193         if ($state eq 'waitnl' && $buf =~ /[\cJ\cM]+/) {
194                 dbg "Got \\n" if isdbg 'state';
195                 Mojo::IOLoop->remove($tid) if $tid;
196                 undef $tid;
197                 undef $buf;
198                 $nlcount = 0;
199                 $ser->write("LPS 1 1\n");
200                 chgstate("waitloop");
201         } elsif ($state eq "waitloop") {
202                 if ($buf =~ /\x06/) {
203                         dbg "Got ACK 0x06" if isdbg 'state';
204                         chgstate('waitlooprec');
205                         undef $buf;
206                 }
207         } elsif ($state eq 'waitlooprec') {
208                 if (length $buf >= 99) {
209                         dbg "got loop record" if isdbg 'chan';
210                         chgstate('');
211                         process($buf);
212                         undef $buf;
213                 }
214         }
215 }
216
217 sub start_loop
218 {
219         $nlcount = 1;
220         dbg "start_loop writing $nlcount \\n" if isdbg 'state';
221         $ser->write("\n");
222         Mojo::IOLoop->remove($tid) if $tid;
223         undef $tid;
224         $tid = Mojo::IOLoop->recurring(0.6 => sub {
225                                                                            if (++$nlcount > 10) {
226                                                                                    dbg "\\n count > 10, closing connection" if isdbg 'chan';
227                                                                                    do_reopen($devname);
228                                                                                    return;
229                                                                            }
230                                                                            dbg "writing $nlcount \\n" if isdbg 'state'; 
231                                                                            $ser->write("\n");
232                                                                    });
233         chgstate("waitnl");
234 }
235
236 sub chgstate
237 {
238         dbg "state '$state' -> '$_[0]'" if isdbg 'state';
239         $state = $_[0];
240 }
241
242 sub do_reopen
243 {
244         my $name = shift;
245         dbg "do reopen on '$name' ending $ending";
246         unless ($ending) {
247                 $ser = do_open($name);
248                 start_loop();
249                 chgstate('');
250                 $nlcount = 0;
251                 Mojo::IOLoop->start unless Mojo::IOLoop->is_running;
252         }
253 }
254
255 sub do_open
256 {
257         my $name = shift;
258
259         $ser->close if $ser;
260         undef $ser;
261
262         my $ob = Serial->new($name, 19200) || die "$name $!\n";
263         dbg "streaming $name fileno(" . fileno($ob) . ")" if isdbg 'chan';
264         
265         my $ser = Mojo::IOLoop::Stream->new($ob);
266         $ser->on(error=>sub {dbg "serial $_[1]"; do_reopen($name) unless $ending});
267         $ser->on(close=>sub {dbg "serial closing"; do_reopen($name) unless $ending});
268         $ser->on(timeout=>sub {dbg "serial timeout";});
269         $ser->on(read=>sub {on_read(@_)});
270         $ser->start;
271
272         Mojo::IOLoop->remove($tid) if $tid;
273         undef $tid;
274         Mojo::IOLoop->remove($rid) if $rid;
275         undef $rid;
276         $rid = Mojo::IOLoop->recurring($poll_interval => sub {
277                                                                            start_loop() if !$state;
278                                                                    });
279         chgstate('');
280         $nlcount = 0;
281         
282         return $ser;
283 }
284
285 my @min;
286 my @hour;
287
288 sub process
289 {
290         my $blk = shift;
291         my $loo =  substr $blk,0,3;
292         unless ( $loo eq 'LOO') {
293                 dbg "Block invalid loo -> $loo" if isdbg 'chan'; return;
294         }
295
296         my $tmp;
297         my $temp;
298         my $rain;
299         my %h;
300
301         my $crc_calc = CRC_CCITT($blk);
302
303         if ($crc_calc==0) {
304
305                 # Common ones
306                 $tmp = unpack("s", substr $blk,7,2) / 1000;
307                 $h{Pressure} = nearest(1, in2mb($tmp));
308
309                 $tmp = unpack("s", substr $blk,9,2) / 10;
310                 $h{Temp_In} = nearest(0.1, f2c($tmp));
311
312                 $temp  = nearest(0.1, f2c(unpack("s", substr $blk,12,2) / 10));
313                 $h{Temp_Out}  = $temp;
314                 if ($temp > 75 || $temp < -75) {
315                         dbg "LOOP Temperature out of range ($temp), record ignored";
316                         return;
317                 }
318
319                 $tmp = unpack("C", substr $blk,14,1);
320                 $h{Wind} = nearest(0.1, mph2mps($tmp));
321                 $h{Dir}     = unpack("s", substr $blk,16,2)+0;
322
323                 my $wind = {w => $h{Wind}, d => $h{Dir}};
324                 $wind = 0 if $wind == 255;
325                 push @{$ld->{wind_min}}, $wind;
326
327                 $tmp = int(unpack("C", substr $blk,33,1)+0);
328                 if ($tmp > 100) {
329                         dbg "LOOP Outside Humidity out of range ($tmp), record ignored";
330                         return;
331                 }
332                 $h{Humidity_Out} = $tmp;
333                 $tmp = int(unpack("C", substr $blk,11,1)+0);
334                 if ($tmp > 100) {
335                         dbg "LOOP Inside Humidity out of range ($tmp), record ignored";
336                         return;
337                 }
338                 $h{Humidity_In}  = $tmp;
339                 
340
341                 $tmp = unpack("C", substr $blk,43,1)+0;
342                 $h{UV} = $tmp unless $tmp >= 255;
343                 $tmp = unpack("s", substr $blk,44,2)+0; # watt/m**2
344                 $h{Solar} = $tmp unless $tmp >= 32767;
345
346                 #       $h{Rain_Rate}  = nearest(0.1,unpack("s", substr $blk,41,2) * $rain_mult);
347                 $rain = $h{Rain_Day}   = nearest(0.1, unpack("s", substr $blk,50,2) * $rain_mult);
348                 my $delta_rain = $h{Rain} = nearest(0.1, ($rain >= $ld->{last_rain} ? $rain - $ld->{last_rain} : $rain)) if $loop_count;
349                 $ld->{last_rain} = $rain;
350
351                 # what sort of packet is it?
352                 my $sort =  unpack("C", substr $blk,4,1);
353                 if ($sort) {
354
355                         # Newer LOOP2 packet
356                         $tmp = unpack("C", substr $blk,18,2);
357                         #               $h{Wind_Avg_10} = nearest(0.1,mph2mps($tmp/10));
358                         $tmp = unpack("C", substr $blk,20,2);
359                         #               $h{Wind_Avg_2} = nearest(0.1,mph2mps($tmp/10));
360                         $tmp = unpack("C", substr $blk,22,2);
361                         #               $h{Wind_Gust_10} = nearest(0.1,mph2mps($tmp/10));
362
363                         #               $h{Dir_Avg_10} = unpack("C", substr $blk,24,2)+0;
364                         $tmp = unpack("C", substr $blk,30,2);
365                         $h{Dew_Point} = nearest(0.1, f2c($tmp));
366
367                 } else {
368
369                         # Older LOOP packet
370                         $tmp = unpack("C", substr $blk,15,1);
371                         #               $h{Wind_Avg_10} = nearest(0.1,mph2mps($tmp));
372                         $h{Dew_Point}  = nearest(0.1, dew_point($h{Temp_Out}, $h{Humidity_Out}));
373                         $h{Rain_Month}  = nearest(0.1, unpack("s", substr $blk,52,2) * $rain_mult);
374                         $h{Rain_Year}  = nearest(0.1, unpack("s", substr $blk,54,2) * $rain_mult);
375                 }
376
377                 my $ts = time;
378                 my $s;
379                 my $dayno = int($ts/86400);
380                 if ($dayno > $ld->{last_day}) {
381                         $ld->{Temp_Out_Max} = $ld->{Temp_Out_Min} = $temp;
382                         $ld->{last_day} = $dayno;
383                 }
384                 $ld->{Temp_Out_Max} = $temp if $temp > $ld->{Temp_Out_Max};
385                 $ld->{Temp_Out_Min} = $temp if $temp < $ld->{Temp_Out_Min};
386
387                 if ($ts >= $ld->{last_hour} + 1800) {
388                         $h{Pressure_Trend}    = unpack("C", substr $blk,3,1);
389                         $h{Pressure_Trend_txt} = $bar_trend{$h{Pressure_Trend}};
390                         $h{Batt_TX_OK}  = (unpack("C", substr $blk,86,1)+0) ^ 1;
391                         $h{Batt_Console}  = nearest(0.01, unpack("s", substr $blk,87,2) * 0.005859375);
392                         $h{Forecast_Icon}  = unpack("C", substr $blk,89,1);
393                         $h{Forecast_Rule}  = unpack("C", substr $blk,90,1);
394                         $h{Sunrise}  = sprintf( "%04d", unpack("S", substr $blk,91,2) );
395                         $h{Sunrise}  =~ s/(\d{2})(\d{2})/$1:$2/;
396                         $h{Sunset}   = sprintf( "%04d", unpack("S", substr $blk,93,2) );
397                         $h{Sunset}  =~ s/(\d{2})(\d{2})/$1:$2/;
398                         $h{Temp_Out_Max} = $ld->{Temp_Out_Max};
399                         $h{Temp_Out_Min} = $ld->{Temp_Out_Min};
400
401                         if ($loop_count) {      # i.e not the first
402                                 my $a = wind_average(scalar @{$ld->{wind_hour}} ? @{$ld->{wind_hour}} : {w => $h{Wind}, d => $h{Dir}});
403
404                                 $h{Wind_1h} = nearest(0.1, $a->{w});
405                                 $h{Dir_1h} = nearest(0.1, $a->{d});
406
407                                 $a = wind_average(@{$ld->{wind_min}});
408                                 $h{Wind_1m} = nearest(0.1, $a->{w});
409                                 $h{Dir_1m} = nearest(1, $a->{d});
410
411                                 ($h{Rain_1m}, $h{Rain_1h}, $h{Rain_24h}) = calc_rain($rain);
412                         }
413                         $ld->{last_rain_min} = $ld->{last_rain_hour} = $rain;
414
415                         $last_hour_h = {%h};
416                         $s = genstr($ts, 'h', \%h);
417                         $ld->{lasthour_h} = $s;
418                         
419                         $ld->{last_hour} = int($ts/1800)*1800;
420                         $ld->{last_min} = int($ts/60)*60;
421                         @{$ld->{wind_hour}} = ();
422                         @{$ld->{wind_min}} = ();
423
424                         write_ld();
425                         
426                 } elsif ($ts >= $ld->{last_min} + 60) {
427                         my $a = wind_average(@{$ld->{wind_min}});
428                         my %save;
429
430                         push @{$ld->{wind_hour}}, $a;
431
432                         if ($loop_count) {      # i.e not the first
433                                 my $rm;
434                                 
435                                 $h{Wind_1m} = nearest(0.1, $a->{w});
436                                 $h{Dir_1m} = nearest(1, $a->{d});
437                                 ($h{Rain_1m}, $h{Rain_1h}, $h{Rain_24h}) = calc_rain($rain);
438                         }
439                         $ld->{last_rain_min} = $rain;
440
441                         $h{Temp_Out_Max} = $ld->{Temp_Out_Max};
442                         $h{Temp_Out_Min} = $ld->{Temp_Out_Min};
443
444                         $last_min_h = {%h};
445                         $s = genstr($ts, 'm', \%h);
446                         $ld->{lastmin_h} = $s;
447                         
448                         $ld->{last_min} = int($ts/60)*60;
449                         @{$ld->{wind_min}} = ();
450                         
451                         write_ld();
452
453                 } else {
454                         my $o = gen_hash_diff($ld->{last_h}, \%h);
455                         if ($o) {
456                                 $s = genstr($ts, 'r', $o);
457                         }
458                         else {
459                                 dbg "loop rec not changed" if isdbg 'chan';
460                         }
461                 }
462                 output_str($s) if $s;
463                 $ld->{last_h} = \%h;
464                 ++$loop_count;
465         } else {
466                 dbg "CRC check failed for LOOP data!";
467         }
468         return;
469 }
470
471 sub genstr
472 {
473         my $ts = shift;
474         my $let = shift;
475         my $h = shift;
476         
477         my $j =  $json->encode($h);
478         my ($sec,$min,$hr) = (gmtime $ts)[0,1,2];
479         my $tm = sprintf "%02d:%02d:%02d", $hr, $min, $sec;
480         
481         return qq|{"tm":"$tm","t":$ts,"$let":$j}|;
482 }
483
484 sub output_str
485 {
486         my $s = shift;
487         dbg $s;
488 #       say $s;
489         $dlog->writenow($s);
490         foreach my $ws (keys $WS) {
491                 my $tx = $WS->{$ws};
492                 if ($tx) {
493                         $tx->send($s);
494                 } else {
495                         delete $WS->{$tx};
496                 }
497         }
498 }
499
500 sub gen_hash_diff
501 {
502         my $last = shift;
503         my $now = shift;
504         my %o;
505         my $count;
506
507         while (my ($k, $v) = each %$now) {
508                 if ($last->{$k} ne $now->{$k}) {
509                         $o{$k} = $v;
510                         ++$count;
511                 }
512         }
513         return $count ? \%o : undef;
514 }
515
516 sub dew_point
517 {
518         my $temp = shift @_;
519         my $rh   = shift @_;
520
521         #  Using the simplified approximation for dew point
522         #  Accurate to 1 degree C for humidities > 50 %
523         #  http://en.wikipedia.org/wiki/Dew_point
524
525         my $dewpoint = $temp - ((100 - $rh) / 5);
526
527         # this is the more complete one (which doesn't work)
528         #my $a = 6.1121;
529         #my $b = 18.678;
530         #my $c = 257.14;
531         #my $ytrh = log(($rh/100) + ($b * $temp) / ($c + $temp));
532         #my $dewpoint = ($c * $ytrh) / ($b - $ytrh); 
533
534         return $dewpoint;
535 }
536
537 sub CRC_CCITT
538 {
539     # Expects packed data...
540     my $data_str = shift @_;
541
542         my $crc = 0;
543         my @lst = split //, $data_str;
544         foreach my $data (@lst) {
545                 my $data = unpack("c",$data);
546
547                 my $crc_prev = $crc;
548                 my $index = $crc >> 8 ^ $data;
549                 my $lhs = $crc_table[$index];
550                 #print "lhs=$lhs, crc=$crc\n";
551                 my $rhs = ($crc << 8) & 0xFFFF;
552                 $crc = $lhs ^ $rhs;
553
554
555         }
556
557         return $crc;
558 }
559
560 sub f2c
561 {
562         return ($_[0] - 32) * 5/9;
563 }
564
565 sub mph2mps
566 {
567         return $_[0] * 0.44704;
568 }
569
570 sub in2mb
571 {
572         return $_[0] * 33.8637526;
573 }
574
575 sub wind_average
576 {
577         my ($sindir, $cosdir, $wind);
578         my $count;
579         
580         foreach my $r (@_) {
581                 $wind += $r->{w};
582                 $sindir += sin(d2r($r->{d})) * $r->{w};
583                 $cosdir += cos(d2r($r->{d})) * $r->{w};
584                 ++$count;
585         }
586
587         my $avhdg = r2d(atan2($sindir, $cosdir));
588         $avhdg += 360 if $avhdg < 0; 
589         return {w => nearest(0.1,$wind / $count), d => nearest(0.1,$avhdg)};
590 }
591
592 # radians to degrees
593 sub r2d
594 {
595     my $n = shift;
596     return ($n / pi) * 180;
597 }
598
599 # degrees to radians
600 sub d2r 
601 {
602     my $n = shift;
603     return ($n / 180) * pi;
604 }
605
606 sub calc_rain
607 {
608         my $rain = shift;
609         
610         $ld->{rain24} ||= [];
611         
612         my $Rain_1h = nearest(0.1, $rain >= $ld->{last_rain_hour} ? $rain - $ld->{last_rain_hour} : $rain); # this is the rate for this hour, so far
613         my $rm = nearest(0.1, $rain >= $ld->{last_rain_min} ? $rain - $ld->{last_rain_min} : $rain);
614         my $Rain_1m = nearest(0.1, $rm);
615         push @{$ld->{rain24}}, $Rain_1m;
616         $ld->{rain_24} += $rm;
617         while (@{$ld->{rain24}} > 24*60) {
618                 $ld->{rain_24} -= shift @{$ld->{rain24}};
619         }
620         my $Rain_24h = nearest(0.1, $ld->{rain_24});
621         return ($Rain_1m, $Rain_1h, $Rain_24h);
622 }
623
624 sub read_ld
625 {
626         return unless $dataf;
627
628         seek $dataf, 0, 0;
629         my $s = <$dataf>;
630         chomp $s;
631         dbg "read loop data: $s" if isdbg 'json';
632         $ld = $json->decode($s) if length $s;
633         
634         # sort out rain stats
635         my $c;
636         if ($ld->{rain24} && ($c = @{$ld->{rain24}}) < 24*60) {
637                 my $diff = 24*60 - $c;
638                 unshift @{$ld->{rain24}}, 0 for 0 .. $diff;  
639         }
640         my $rain;
641         
642         if ($ld->{rain24}) {
643                 $rain += $_ for @{$ld->{rain24}};
644         }
645         
646         $ld->{rain_24} = nearest(0.1, $rain);
647         delete $ld->{hour};
648         delete $ld->{min};
649 }
650
651 sub write_ld
652 {
653         return unless $dataf;
654
655         seek $dataf, 0, 0;
656         truncate $dataf, 0;
657         $ld->{ts} = time;
658         my $s = $json->encode($ld);
659         dbg "write loop data: $s" if isdbg 'json';
660         print $dataf "$s\n";
661 }
662
663 sub cycle_loop_data_files
664 {
665         rename "$datafn.oooo", "$datafn.ooooo";
666         rename "$datafn.ooo", "$datafn.oooo";
667         rename "$datafn.oo", "$datafn.ooo";
668         rename "$datafn.o", "$datafn.oo";
669         copy $datafn, "$datafn.o";
670 }
671
672 __DATA__
673  
674 @@ index.html.ep
675 % my $url = url_for 'weather';
676 <!DOCTYPE html>
677 <html>
678   <head><title>DWeather</title></head>
679   <body>
680         <h1>DWeather</h1>
681
682     <script>
683         var ws;
684
685         function process(key,value) {
686                 var d = document.getElementById(key);
687                 if (d !== null) {
688                         d.innerHTML = value;
689                 }
690         }
691
692     function traverse(o) {
693                 console.log(o);
694                 for (var i in o) {
695                         process(i, o[i]);  
696                         if (o[i] !== null && typeof(o[i])=="object") {
697                                 traverse(o[i]);
698                         }
699                 }
700         }
701
702
703     ws = new WebSocket('<%= $url->to_abs %>');
704     document.body.innerHTML += 'ws connecting to: <%= $url->to_abs %> type_of: ' + typeof(ws) + '<br>';
705     if (typeof(ws) === 'object') {
706                 ws.onmessage = function (event) {
707                         var js = JSON.parse(event.data);
708                         if (js !== null && typeof(js) === 'object') {
709                                 traverse(js);
710                         }
711                 };
712                 ws.onopen = function (event) {
713                         ws.send('WebSocket support works! ♥');
714                 };
715         } else {
716                 document.body.innerHTML += 'Webserver only works with Websocket aware browsers';
717         }
718  
719     </script>
720
721 <br><br>
722 <table border=1 width=80%>
723 <tr>
724 <th>Time:<td><span id="tm"> </span>
725 <th>Sunrise:<td><span id="Sunrise"> </span>
726 <th>Sunset:<td><span id="Sunset"> </span>
727 <th>Console Volts:<td><span id="Batt_Console"> </span>
728 <th>TX Battery OK:<td><span id="Batt_TX_OK"> </span>
729
730 <tr>
731 <th>Pressure:<td><span id="Pressure"> </span>
732 <th>Trend:<td><span id="Pressure_Trend_txt"> </span>
733
734 <tr>
735 <th>Temperature in:<td> <span id="Temp_In"> </span> 
736 <th>Humidity:<td> <span id="Humidity_In"> </span>
737
738 <tr>
739 <th>Temperature out:<td> <span id="Temp_Out"> </span> 
740 <th>Min:<td> <span id="Temp_Out_Min"> </span> 
741 <th>Max:<td> <span id="Temp_Out_Max"> </span> 
742 <th>Humidity:<td> <span id="Humidity_Out"> </span> 
743 <th>Dew Point:<td> <span id="Dew_Point"> </span> 
744
745 <tr>
746 <th>Wind Direction:<td> <span id="Dir"> </span> 
747 <th>Minute Avg:<td> <span id="Dir_1m">  </span> 
748 <th>Speed:<td> <span id="Wind"> </span> 
749 <th>Minute Avg:<td> <span id="Wind_1m">  </span>
750
751 <tr>
752 <th>Rain 30mins:<td> <span id="Rain_1h"> </span>
753 <th>Day:<td> <span id="Rain_Day"> </span>
754 <th>24hrs:<td> <span id="Rain_24h"> </span>
755 <th>Month:<td> <span id="Rain_Month"> </span>
756 <th>Year:<td> <span id="Rain_Year"> </span>
757 </table>
758 </body>
759 </html>