loop.pl

   1 #!/usr/bin/perl
   2 use strict;
   3
   4 use v5.10.1;
   5
   6 use DBI;
   7 use Serial;
   8 use Mojo::IOLoop;
   9 use Mojo::IOLoop::Stream;
  10
  11 my $devname = "/dev/davis";
  12 my $tid;
  13 my $rid;
  14 my $count;
  15 my $state = "ready";
  16 my $buf;
  17
  18 our     @crc_table = (
  19                                   0x0, 0x1021, 0x2042, 0x3063, 0x4084, 0x50a5, 0x60c6, 0x70e7,
  20                                   0x8108, 0x9129, 0xa14a, 0xb16b, 0xc18c, 0xd1ad, 0xe1ce, 0xf1ef,
  21                                   0x1231, 0x210, 0x3273, 0x2252, 0x52b5, 0x4294, 0x72f7, 0x62d6,
  22                                   0x9339, 0x8318, 0xb37b, 0xa35a, 0xd3bd, 0xc39c, 0xf3ff, 0xe3de,
  23                                   0x2462, 0x3443, 0x420, 0x1401, 0x64e6, 0x74c7, 0x44a4, 0x5485,
  24                                   0xa56a, 0xb54b, 0x8528, 0x9509, 0xe5ee, 0xf5cf, 0xc5ac, 0xd58d,
  25                                   0x3653, 0x2672, 0x1611, 0x630, 0x76d7, 0x66f6, 0x5695, 0x46b4,
  26                                   0xb75b, 0xa77a, 0x9719, 0x8738, 0xf7df, 0xe7fe, 0xd79d, 0xc7bc,
  27                                   0x48c4, 0x58e5, 0x6886, 0x78a7, 0x840, 0x1861, 0x2802, 0x3823,
  28                                   0xc9cc, 0xd9ed, 0xe98e, 0xf9af, 0x8948, 0x9969, 0xa90a, 0xb92b,
  29                                   0x5af5, 0x4ad4, 0x7ab7, 0x6a96, 0x1a71, 0xa50, 0x3a33, 0x2a12,
  30                                   0xdbfd, 0xcbdc, 0xfbbf, 0xeb9e, 0x9b79, 0x8b58, 0xbb3b, 0xab1a,
  31                                   0x6ca6, 0x7c87, 0x4ce4, 0x5cc5, 0x2c22, 0x3c03, 0xc60, 0x1c41,
  32                                   0xedae, 0xfd8f, 0xcdec, 0xddcd, 0xad2a, 0xbd0b, 0x8d68, 0x9d49,
  33                                   0x7e97, 0x6eb6, 0x5ed5, 0x4ef4, 0x3e13, 0x2e32, 0x1e51, 0xe70,
  34                                   0xff9f, 0xefbe, 0xdfdd, 0xcffc, 0xbf1b, 0xaf3a, 0x9f59, 0x8f78,
  35                                   0x9188, 0x81a9, 0xb1ca, 0xa1eb, 0xd10c, 0xc12d, 0xf14e, 0xe16f,
  36                                   0x1080, 0xa1, 0x30c2, 0x20e3, 0x5004, 0x4025, 0x7046, 0x6067,
  37                                   0x83b9, 0x9398, 0xa3fb, 0xb3da, 0xc33d, 0xd31c, 0xe37f, 0xf35e,
  38                                   0x2b1, 0x1290, 0x22f3, 0x32d2, 0x4235, 0x5214, 0x6277, 0x7256,
  39                                   0xb5ea, 0xa5cb, 0x95a8, 0x8589, 0xf56e, 0xe54f, 0xd52c, 0xc50d,
  40                                   0x34e2, 0x24c3, 0x14a0, 0x481, 0x7466, 0x6447, 0x5424, 0x4405,
  41                                   0xa7db, 0xb7fa, 0x8799, 0x97b8, 0xe75f, 0xf77e, 0xc71d, 0xd73c,
  42                                   0x26d3, 0x36f2, 0x691, 0x16b0, 0x6657, 0x7676, 0x4615, 0x5634,
  43                                   0xd94c, 0xc96d, 0xf90e, 0xe92f, 0x99c8, 0x89e9, 0xb98a, 0xa9ab,
  44                                   0x5844, 0x4865, 0x7806, 0x6827, 0x18c0, 0x8e1, 0x3882, 0x28a3,
  45                                   0xcb7d, 0xdb5c, 0xeb3f, 0xfb1e, 0x8bf9, 0x9bd8, 0xabbb, 0xbb9a,
  46                                   0x4a75, 0x5a54, 0x6a37, 0x7a16, 0xaf1, 0x1ad0, 0x2ab3, 0x3a92,
  47                                   0xfd2e, 0xed0f, 0xdd6c, 0xcd4d, 0xbdaa, 0xad8b, 0x9de8, 0x8dc9,
  48                                   0x7c26, 0x6c07, 0x5c64, 0x4c45, 0x3ca2, 0x2c83, 0x1ce0, 0xcc1,
  49                                   0xef1f, 0xff3e, 0xcf5d, 0xdf7c, 0xaf9b, 0xbfba, 0x8fd9, 0x9ff8,
  50                                   0x6e17, 0x7e36, 0x4e55, 0x5e74, 0x2e93, 0x3eb2, 0xed1, 0x1ef0
  51                                  );
  52
  53 our %bar_trend;
  54
  55 $bar_trend{-60} = "Falling Rapidly";
  56 $bar_trend{196} = "Falling Rapidly";
  57 $bar_trend{-20} = "Falling Slowly";
  58 $bar_trend{236} = "Falling Slowly";
  59 $bar_trend{0} = "Steady";
  60 $bar_trend{20} = "Rising Slowly";
  61 $bar_trend{60} = "Rising Rapidly";
  62
  63 #$SIG{TERM} = $SIG{INT} = sub {Mojo::IOLoop->reset if Mojo::IOLoop->is_running && !$DB::VERSION};
  64
  65 my $s = do_open($devname);
  66 start_loop();
  67
  68 Mojo::IOLoop->start unless Mojo::IOLoop->is_running;
  69
  70 exit 0;
  71
  72 sub on_read
  73 {
  74         my ($str, $d) = @_;
  75         $buf .= $d;
  76         $d =~ s/([\%\x00-\x1f\x7f-\xff])/sprintf("%%%02X", ord($1))/eg;
  77 #       say "read added '$d' buf lth=" . length $buf;
  78         if ($state eq 'waitnl' && $buf =~ /[\cJ\cM]+/) {
  79                 undef $tid;
  80                 undef $buf;
  81                 $s->write("LOOP 1\n");
  82                 chgstate("waitloop");
  83         } elsif ($state eq "waitloop") {
  84                 if ($buf =~ /\x06/) {
  85                         chgstate('waitlooprec');
  86                         undef $buf;
  87                 }
  88         } elsif ($state eq 'waitlooprec') {
  89                 if (length $buf >= 99) {
  90                         say "got loop record\n";
  91                         chgstate('');
  92                         process($buf);
  93                         undef $buf;
  94                 }
  95         }
  96 }
  97
  98 sub start_loop
  99 {
 100         say "writing \\n";
 101         $s->write("\n");
 102         $tid = Mojo::IOLoop->timer(0.6 => sub {say "writing \\n"; $s->write("\n")});
 103         chgstate("waitnl");
 104 }
 105
 106 sub chgstate
 107 {
 108         say "state '$state' -> '$_[0]'";
 109         $state = $_[0];
 110 }
 111
 112 sub do_open
 113 {
 114         my $name = shift;
 115         my $ob = Serial->new($name, 19200) || die "$name $!\n";
 116         say "streaming $name fileno(", fileno($ob), ")";
 117         my $str = Mojo::IOLoop::Stream->new($ob);
 118         $str->on(error=>sub {say "serial $_[1]"; undef $s; Mojo::IOLoop->reset;});
 119         $str->on(close=>sub {say "serial closing"; undef $s;  Mojo::IOLoop->reset;});
 120         $str->on(timeout=>sub {say "serial timeout";});
 121         $str->on(read=>sub {on_read(@_)});
 122         $str->start;
 123
 124         $rid = Mojo::IOLoop->recurring(2.5 => sub {start_loop() if !$state || $state eq "waitnl";});
 125
 126         return $str;
 127 }
 128
 129 sub process
 130 {
 131         my $blk = shift;
 132         my $loo =  substr $blk,0,3;
 133         unless ( $loo eq 'LOO') {
 134                 say "Block invalid loo -> $loo"; return;
 135         }
 136
 137         my $t;
 138         my %hsh;
 139
 140
 141         #$hsh{'next_rec'}     = unpack("s", substr $blk,5,2);
 142
 143         $hsh{'Barometric_Trend'}    = unpack("C", substr $blk,3,1);
 144         $hsh{'Barometric_Trend_txt'} = $bar_trend{$hsh{'Barometric_Trend'}};
 145         $t = unpack("s", substr $blk,7,2) / 1000;
 146         #       $hsh{'Barometric_Press_hg'}          = $t;
 147         $hsh{'Barometric_Press_mb'} = sprintf("%.2f",$t*33.8637526);
 148
 149
 150         $t = unpack("s", substr $blk,9,2) / 10;
 151         #       $hsh{'Air_Temp_Inside_f'}      = $t;
 152         $hsh{'Air_Temp_Inside_c'} = sprintf("%.1f",($t - 32) * 5/9);
 153         my $tf  = unpack("s", substr $blk,12,2) / 10;
 154         #       $hsh{'Air_Temp_Outside_f'}     = $tf;
 155         $hsh{'Air_Temp_Outside_c'}  = sprintf("%.1f",($tf - 32) * 5/9);
 156
 157         $hsh{'Wind_Speed_mph'}   = unpack("C", substr $blk,14,1);
 158         #       $hsh{'Wind_Speed_mps'} = sprintf("%.1f",$hsh{'Wind_Speed_mph'}*0.44704);
 159         $hsh{'Wind_Speed_10min_Avg_mph'} = unpack("C", substr $blk,15,1);
 160         #       $hsh{'Wind_Speed_10min_Avg_mps'} = sprintf("%.1f",$hsh{'Wind_Speed_10min_Avg_mph'}*0.44704);
 161         $hsh{'Wind_Dir'}     = unpack("s", substr $blk,16,2);
 162
 163
 164         $hsh{'Humidity_Outside'} = unpack("C", substr $blk,33,1);
 165         $hsh{'Humidity_Inside'}  = unpack("C", substr $blk,11,1);
 166         $hsh{'Dew_Point'}  = dew_point($tf, $hsh{'Humidity_Outside'});
 167
 168         #       $hsh{'UV'}         = unpack("C", substr $blk,43,1);
 169         #       $hsh{'Solar'}  = unpack("s", substr $blk,44,2); # watt/m**2
 170
 171         $hsh{'Rain_Rate'}  = (unpack("s", substr $blk,41,2) / 100) * 25.4; # Inches per hr converted to mm
 172         $hsh{'Rain_Storm'} = (unpack("s", substr $blk,46,2) / 100) * 25.4; # Inches per storm
 173         #$hsh{'Storm_Date'} = unpack("s", substr $blk,48,2);  # Need to parse data (not sure what this is)
 174         $hsh{'Rain_Day'}   = (unpack("s", substr $blk,50,2)/100) * 25.4;
 175         $hsh{'Rain_Month'}  = (unpack("s", substr $blk,52,2)/100) * 25.4;
 176         $hsh{'Rain_Year'}  = (unpack("s", substr $blk,54,2)/100) * 25.4;
 177
 178         $hsh{'ET_Day'}   = unpack("s", substr $blk,56,2)/1000;
 179         $hsh{'ET_Month'}  = unpack("s", substr $blk,58,2)/100;
 180         $hsh{'ET_Year'}  = unpack("s", substr $blk,60,2)/100;
 181
 182         #$hsh{'Alarms_Inside'}  = unpack("b8", substr $blk,70,1);
 183         #$hsh{'Alarms_Rain'}  = unpack("b8", substr $blk,70,1);
 184         #$hsh{'Alarms_Outside'}  = unpack("b8", substr $blk,70,1);
 185
 186         $hsh{'Batt_Transmitter'}  = unpack("C", substr $blk,86,1); #  * 0.005859375
 187         $hsh{'Batt_Console'}  = unpack("s", substr $blk,87,2) * 0.005859375;
 188
 189         $hsh{'Forecast_Icon'}  = unpack("C", substr $blk,89,1);
 190         $hsh{'Forecast_Rule'}  = unpack("C", substr $blk,90,1);
 191
 192         $hsh{'Sunrise'}  = sprintf( "%04d", unpack("S", substr $blk,91,2) );
 193         $hsh{'Sunrise'}  =~ s/(\d{2})(\d{2})/$1:$2/;
 194         $hsh{'Sunset'}   = sprintf( "%04d", unpack("S", substr $blk,93,2) );
 195         $hsh{'Sunset'}  =~ s/(\d{2})(\d{2})/$1:$2/;
 196
 197         #my $nl  =  ord substr $blk,95,1;
 198         #my $cr  =  ord substr $blk,96,1;
 199
 200         my $crc = unpack "%n", substr($blk,97,2);
 201         my $crc_calc = CRC_CCITT($blk);
 202
 203         if ($crc_calc==0) {
 204                 say "inside: $hsh{Air_Temp_Inside_c} degC $hsh{Humidity_Inside}\% outside: $hsh{Air_Temp_Outside_c} degC $hsh{Humidity_Outside}\% wind: $hsh{Wind_Speed_mph} $hsh{Wind_Dir} deg $hsh{Barometric_Press_mb} mB";
 205
 206         } else {
 207                 print "CRC check failed for LOOP data!\n";
 208                 return 1;
 209         }
 210         #delete @hsh{'crc', 'crc_calc', 'next_rec'};
 211         #delete($hsh{crc})||die"cant delete crc";
 212         #delete($hsh{crc_calc})||die"cant delete crc_calc";
 213         #delete($hsh{next_rec})||die"cant delete next_rec";
 214 }
 215
 216 sub dew_point
 217 {
 218         my $temp = shift @_;
 219         my $rh   = shift @_;
 220
 221         #  Using the simplified approximation for dew point
 222         #  Accurate to 1 degree C for humidities > 50 %
 223         #  http://en.wikipedia.org/wiki/Dew_point
 224
 225         my $dew_point = $temp - ( (100 - $rh)/5 );
 226
 227         return $dew_point;
 228 }
 229
 230 sub CRC_CCITT
 231 {
 232     # Expects packed data...
 233     my $data_str = shift @_;
 234
 235         my $crc = 0;
 236         my @lst = split //, $data_str;
 237         foreach my $data (@lst) {
 238                 my $data = unpack("c",$data);
 239
 240                 my $crc_prev = $crc;
 241                 my $index = $crc >> 8 ^ $data;
 242                 my $lhs = $crc_table[$index];
 243                 #print "lhs=$lhs, crc=$crc\n";
 244                 my $rhs = ($crc << 8) & 0xFFFF;
 245                 $crc = $lhs ^ $rhs;
 246
 247
 248         }
 249
 250         return $crc;
 251 }
 252