8 use Mojo::IOLoop::Stream;
9 use Mojo::JSON qw(decode_json encode_json);
13 my $devname = "/dev/davis";
14 my $rain_mult = 0.2; # 0.1 or 0.2 mm or 0.01 inches
23 my $s; # the serial port Mojo::IOLoop::Stream
25 our $last_min = int(time/60)*60;
26 our $last_hour = int(time/3600)*3600;
29 0x0, 0x1021, 0x2042, 0x3063, 0x4084, 0x50a5, 0x60c6, 0x70e7,
30 0x8108, 0x9129, 0xa14a, 0xb16b, 0xc18c, 0xd1ad, 0xe1ce, 0xf1ef,
31 0x1231, 0x210, 0x3273, 0x2252, 0x52b5, 0x4294, 0x72f7, 0x62d6,
32 0x9339, 0x8318, 0xb37b, 0xa35a, 0xd3bd, 0xc39c, 0xf3ff, 0xe3de,
33 0x2462, 0x3443, 0x420, 0x1401, 0x64e6, 0x74c7, 0x44a4, 0x5485,
34 0xa56a, 0xb54b, 0x8528, 0x9509, 0xe5ee, 0xf5cf, 0xc5ac, 0xd58d,
35 0x3653, 0x2672, 0x1611, 0x630, 0x76d7, 0x66f6, 0x5695, 0x46b4,
36 0xb75b, 0xa77a, 0x9719, 0x8738, 0xf7df, 0xe7fe, 0xd79d, 0xc7bc,
37 0x48c4, 0x58e5, 0x6886, 0x78a7, 0x840, 0x1861, 0x2802, 0x3823,
38 0xc9cc, 0xd9ed, 0xe98e, 0xf9af, 0x8948, 0x9969, 0xa90a, 0xb92b,
39 0x5af5, 0x4ad4, 0x7ab7, 0x6a96, 0x1a71, 0xa50, 0x3a33, 0x2a12,
40 0xdbfd, 0xcbdc, 0xfbbf, 0xeb9e, 0x9b79, 0x8b58, 0xbb3b, 0xab1a,
41 0x6ca6, 0x7c87, 0x4ce4, 0x5cc5, 0x2c22, 0x3c03, 0xc60, 0x1c41,
42 0xedae, 0xfd8f, 0xcdec, 0xddcd, 0xad2a, 0xbd0b, 0x8d68, 0x9d49,
43 0x7e97, 0x6eb6, 0x5ed5, 0x4ef4, 0x3e13, 0x2e32, 0x1e51, 0xe70,
44 0xff9f, 0xefbe, 0xdfdd, 0xcffc, 0xbf1b, 0xaf3a, 0x9f59, 0x8f78,
45 0x9188, 0x81a9, 0xb1ca, 0xa1eb, 0xd10c, 0xc12d, 0xf14e, 0xe16f,
46 0x1080, 0xa1, 0x30c2, 0x20e3, 0x5004, 0x4025, 0x7046, 0x6067,
47 0x83b9, 0x9398, 0xa3fb, 0xb3da, 0xc33d, 0xd31c, 0xe37f, 0xf35e,
48 0x2b1, 0x1290, 0x22f3, 0x32d2, 0x4235, 0x5214, 0x6277, 0x7256,
49 0xb5ea, 0xa5cb, 0x95a8, 0x8589, 0xf56e, 0xe54f, 0xd52c, 0xc50d,
50 0x34e2, 0x24c3, 0x14a0, 0x481, 0x7466, 0x6447, 0x5424, 0x4405,
51 0xa7db, 0xb7fa, 0x8799, 0x97b8, 0xe75f, 0xf77e, 0xc71d, 0xd73c,
52 0x26d3, 0x36f2, 0x691, 0x16b0, 0x6657, 0x7676, 0x4615, 0x5634,
53 0xd94c, 0xc96d, 0xf90e, 0xe92f, 0x99c8, 0x89e9, 0xb98a, 0xa9ab,
54 0x5844, 0x4865, 0x7806, 0x6827, 0x18c0, 0x8e1, 0x3882, 0x28a3,
55 0xcb7d, 0xdb5c, 0xeb3f, 0xfb1e, 0x8bf9, 0x9bd8, 0xabbb, 0xbb9a,
56 0x4a75, 0x5a54, 0x6a37, 0x7a16, 0xaf1, 0x1ad0, 0x2ab3, 0x3a92,
57 0xfd2e, 0xed0f, 0xdd6c, 0xcd4d, 0xbdaa, 0xad8b, 0x9de8, 0x8dc9,
58 0x7c26, 0x6c07, 0x5c64, 0x4c45, 0x3ca2, 0x2c83, 0x1ce0, 0xcc1,
59 0xef1f, 0xff3e, 0xcf5d, 0xdf7c, 0xaf9b, 0xbfba, 0x8fd9, 0x9ff8,
60 0x6e17, 0x7e36, 0x4e55, 0x5e74, 0x2e93, 0x3eb2, 0xed1, 0x1ef0
65 $bar_trend{-60} = "Falling Rapidly";
66 $bar_trend{196} = "Falling Rapidly";
67 $bar_trend{-20} = "Falling Slowly";
68 $bar_trend{236} = "Falling Slowly";
69 $bar_trend{0} = "Steady";
70 $bar_trend{20} = "Rising Slowly";
71 $bar_trend{60} = "Rising Rapidly";
75 $SIG{TERM} = $SIG{INT} = sub {++$ending; Mojo::IOLoop->stop;};
84 dbg "*** starting $0";
87 dbg scalar gmtime($last_min);
88 dbg scalar gmtime($last_hour);
90 my $dlog = SMGLog->new("day");
91 $did = Mojo::IOLoop->recurring(1 => sub {$dlog->flushall});
105 $d =~ s/([\%\x00-\x1f\x7f-\xff])/sprintf("%%%02X", ord($1))/eg;
106 dbg "read added '$d' buf lth=" . length $buf if isdbg 'raw';
107 if ($state eq 'waitnl' && $buf =~ /[\cJ\cM]+/) {
108 dbg "Got \\n" if isdbg 'chan';
109 Mojo::IOLoop->remove($tid) if $tid;
112 $s->write("LPS 1 1\n");
113 chgstate("waitloop");
114 } elsif ($state eq "waitloop") {
115 if ($buf =~ /\x06/) {
116 dbg "Got ACK 0x06" if isdbg 'chan';
117 chgstate('waitlooprec');
120 } elsif ($state eq 'waitlooprec') {
121 if (length $buf >= 99) {
122 dbg "got loop record" if isdbg 'chan';
132 dbg "start_loop writing \\n" if isdbg 'chan';
134 $tid = Mojo::IOLoop->recurring(0.6 => sub {dbg "writing \\n" if isdbg 'chan'; $s->write("\n")});
141 dbg "state '$state' -> '$_[0]'" if isdbg 'chan';
148 dbg "do reopen on '$name' ending $ending";
150 $s->close if defined $s;
154 Mojo::IOLoop->start unless Mojo::IOLoop->is_running;
161 my $ob = Serial->new($name, 19200) || die "$name $!\n";
162 dbg "streaming $name fileno(" . fileno($ob) . ")" if isdbg 'chan';
164 my $str = Mojo::IOLoop::Stream->new($ob);
165 $str->on(error=>sub {dbg "serial $_[1]"; do_reopen($name) unless $ending});
166 $str->on(close=>sub {dbg "serial closing"; do_reopen($name) unless $ending});
167 $str->on(timeout=>sub {dbg "serial timeout";});
168 $str->on(read=>sub {on_read(@_)});
171 $rid = Mojo::IOLoop->recurring(2.5 => sub {
172 start_loop() if !$state || $state eq "waitnl";
184 my $loo = substr $blk,0,3;
185 unless ( $loo eq 'LOO') {
186 dbg "Block invalid loo -> $loo" if isdbg 'chan'; return;
193 $h{Pressure_Trend} = unpack("C", substr $blk,3,1);
194 $h{Pressure_Trend_txt} = $bar_trend{$h{'Pressure_Trend'}};
195 $t = unpack("s", substr $blk,7,2) / 1000;
196 $h{Pressure} = sprintf("%.0f",in2mb($t))+0;
198 $t = unpack("s", substr $blk,9,2) / 10;
199 $h{Temp_In} = sprintf("%.1f", f2c($t))+0;
201 $t = unpack("s", substr $blk,12,2) / 10;
202 $h{Temp_Out} = sprintf("%.1f", f2c($t))+0;
204 $t = unpack("C", substr $blk,14,1);
205 $h{Wind} = sprintf("%.1f",mph2mps($t))+0;
206 $h{Dir} = unpack("s", substr $blk,16,2)+0;
208 my $wind = {w => $h{Wind}, d => $h{Dir}};
211 $h{'Humidity_Out'} = unpack("C", substr $blk,33,1)+0;
212 $h{'Humidity_In'} = unpack("C", substr $blk,11,1)+0;
214 $t = unpack("C", substr $blk,43,1)+0;
215 $h{'UV'} = $t unless $t >= 255;
216 $t = unpack("s", substr $blk,44,2)+0; # watt/m**2
217 $h{'Solar'} = $t unless $t >= 32767;
219 $h{'Rain_Rate'} = sprintf("%0.1f",unpack("s", substr $blk,41,2) * $rain_mult)+0;
220 $h{'Rain_Day'} = sprintf("%0.1f", unpack("s", substr $blk,50,2) * $rain_mult)+0;
222 # what sort of packet is it?
224 $t = unpack("C", substr $blk,4,1);
228 $t = unpack("C", substr $blk,18,2);
229 # $h{Wind_Avg_10} = sprintf("%.1f",mph2mps($t/10))+0;
230 $t = unpack("C", substr $blk,20,2);
231 # $h{Wind_Avg_2} = sprintf("%.1f",mph2mps($t/10))+0;
232 $t = unpack("C", substr $blk,22,2);
233 # $h{Wind_Gust_10} = sprintf("%.1f",mph2mps($t/10))+0;
235 # $h{Dir_Avg_10} = unpack("C", substr $blk,24,2)+0;
236 $t = unpack("C", substr $blk,30,2);
237 $h{Dew_Point} = sprintf("%0.1f", f2c($t))+0;
242 $t = unpack("C", substr $blk,15,1);
243 # $h{Wind_Avg_10} = sprintf("%.1f",mph2mps($t))+0;
244 $h{'Dew_Point'} = sprintf("%0.1f", dew_point($h{Temp_Out}, $h{'Humidity_Out'}))+0;
246 $h{'Rain_Month'} = sprintf("%0.1f", unpack("s", substr $blk,52,2) * $rain_mult)+0;
247 $h{'Rain_Year'} = sprintf("%0.1f", unpack("s", substr $blk,54,2) * $rain_mult)+0;
250 $h{'Batt_TX_OK'} = (unpack("C", substr $blk,86,1)+0) ^ 1;
251 $h{'Batt_Console'} = sprintf("%0.2f", unpack("s", substr $blk,87,2) * 0.005859375)+0;
252 $h{'Forecast_Icon'} = unpack("C", substr $blk,89,1);
253 $h{'Forecast_Rule'} = unpack("C", substr $blk,90,1);
255 $h{'Sunrise'} = sprintf( "%04d", unpack("S", substr $blk,91,2) );
256 $h{'Sunrise'} =~ s/(\d{2})(\d{2})/$1:$2/;
257 $h{'Sunset'} = sprintf( "%04d", unpack("S", substr $blk,93,2) );
258 $h{'Sunset'} =~ s/(\d{2})(\d{2})/$1:$2/;
261 my $crc = unpack "%n", substr($blk,97,2);
262 my $crc_calc = CRC_CCITT($blk);
265 my $o = gen_hash_diff($last_reading, \%h);
270 if ($t >= $last_min + 60) {
271 my $a = average(@min);
274 $h{Wind_1m} = sprintf("%0.1f", $a->{w})+0;
275 $h{Dir_1m} = sprintf("%0.0f", $a->{d})+0;
277 $j = encode_json(\%h);
278 $s = qq|{"t":$t,"m":$j}|;
281 if ($t >= $last_hour + 3600) {
282 my $a = average(@hour);
286 $h{Wind_1h} = sprintf("%0.1f", $a->{w})+0;
287 $h{Dir_1h} = sprintf("%0.0f", $a->{d})+0;
289 $j = encode_json(\%h);
290 $s = qq|{"t":$t,"h":$j}|;
295 $j = encode_json($o);
296 $s = qq|{"t":$t,"r":$j}|;
298 dbg "loop rec not changed" if isdbg 'chan';
306 dbg "CRC check failed for LOOP data!";
317 while (my ($k, $v) = each %$now) {
318 if ($last->{$k} ne $now->{$k}) {
323 return $count ? \%o : undef;
331 # Using the simplified approximation for dew point
332 # Accurate to 1 degree C for humidities > 50 %
333 # http://en.wikipedia.org/wiki/Dew_point
335 my $dewpoint = $temp - ((100 - $rh) / 5);
337 # this is the more complete one (which doesn't work)
341 #my $ytrh = log(($rh/100) + ($b * $temp) / ($c + $temp));
342 #my $dewpoint = ($c * $ytrh) / ($b - $ytrh);
349 # Expects packed data...
350 my $data_str = shift @_;
353 my @lst = split //, $data_str;
354 foreach my $data (@lst) {
355 my $data = unpack("c",$data);
358 my $index = $crc >> 8 ^ $data;
359 my $lhs = $crc_table[$index];
360 #print "lhs=$lhs, crc=$crc\n";
361 my $rhs = ($crc << 8) & 0xFFFF;
372 return ($_[0] - 32) * 5/9;
377 return $_[0] * 0.44704;
382 return $_[0] * 33.8637526;
391 while (my ($k, $v) = each %$r) {
396 while (my ($k, $v) = each %out) {