3 # The geomagnetic information and calculation module
6 # Copyright (c) 1998 - Dirk Koopman G1TLH
20 use Time::HiRes qw(gettimeofday tv_interval);
24 use vars qw($date $sfi $k $a $r $forecast @allowed @denied $fp $node $from
26 $duplth $dupage $filterdef);
28 $fp = 0; # the DXLog fcb
29 $date = 0; # the unix time of the WWV (notional)
30 $sfi = 0; # the current SFI value
31 $k = 0; # the current K value
32 $a = 0; # the current A value
33 $r = 0; # the current R value
34 $forecast = ""; # the current geomagnetic forecast
35 $node = ""; # originating node
36 $from = ""; # who this came from
37 @allowed = (); # if present only these callsigns are regarded as valid WWV updators
38 @denied = (); # if present ignore any wwv from these callsigns
39 $duplth = 20; # the length of text to use in the deduping
40 $dupage = 12*3600; # the length of time to hold spot dups
42 $dirprefix = "$main::local_data/wwv";
43 $param = "$dirprefix/param";
50 # tag, sort, field, priv, special parser
57 ['origin_dxcc', 'nc', 6],
58 ['origin_itu', 'ni', 7],
59 ['origin_zone', 'nz', 8],
64 $fp = DXLog::new('wwv', 'dat', 'm');
65 do "$param" if -e "$param";
66 # read in existing data
67 my $t0 = [gettimeofday];
68 dbg(sprintf "WWV read in upto %d records into cache", $maxcache);
69 @cache = readfile($main::systime);
70 shift @cache while @cache > $maxcache;
71 dbg(sprintf "WWV read in last %d records into cache in %dmS", scalar @cache, _diffms($t0));
75 # write the current data away
78 my $fh = new IO::File;
79 open $fh, "> $param" or confess "can't open $param $!";
80 print $fh "# Geomagnetic data parameter file last mod:", scalar gmtime, "\n";
81 print $fh "\$date = $date;\n";
82 print $fh "\$sfi = $sfi;\n";
83 print $fh "\$a = $a;\n";
84 print $fh "\$k = $k;\n";
85 print $fh "\$r = $r;\n";
86 print $fh "\$from = '$from';\n";
87 print $fh "\$node = '$node';\n";
88 print $fh "\@denied = qw(", join(' ', @denied), ");\n" if @denied > 0;
89 print $fh "\@allowed = qw(", join(' ', @allowed), ");\n" if @allowed > 0;
93 my $s ="$from^$date^$sfi^$a^$k^$forecast^$node^$r";
94 $fp->writeunix($date, $s);
95 push @cache, [ split /\^/, $s ];
96 shift @cache while @cache > $maxcache;
99 # update WWV info in one go (usually from a PC23)
102 my ($mydate, $mytime, $mysfi, $mya, $myk, $myforecast, $myfrom, $mynode, $myr) = @_;
103 $myfrom =~ s/-\d+$//;
104 if ((@allowed && grep {$_ eq $myfrom} @allowed) ||
105 (@denied && !grep {$_ eq $myfrom} @denied) ||
106 (@allowed == 0 && @denied == 0)) {
108 # my $trydate = cltounix($mydate, sprintf("%02d18Z", $mytime));
109 if ($mydate > $date) {
113 $r = 0 unless abs ($mysfi - $sfi) > 3;
118 $forecast = $myforecast;
128 # add or substract an allowed callsign
133 push @allowed, map {uc $_} @_;
137 @allowed = map {$_ ne uc $c} @allowed;
143 # add or substract a denied callsign
148 push @denied, map {uc $_} @_;
152 @denied = map {$_ ne uc $c} @denied;
158 # accessor routines (when I work how symbolic refs work I might use one of those!)
161 @_ ? $sfi = shift : $sfi ;
166 @_ ? $k = shift : $k ;
171 @_ ? $r = shift : $r ;
176 @_ ? $a = shift : $a ;
181 @_ ? $forecast = shift : $forecast ;
186 # print some items from the log backwards in time
188 # This command outputs a list of n lines starting from line $from to $to
195 my $date = $fp->unixtoj($t);
202 if ($t == $main::systime && ($to <= $maxcache)) {
203 dbg("using wwv cache") if isdbg('wwv');
204 @out = reverse @cache;
205 pop @out while @out > $to;
207 dbg("using wwv file(s))") if isdbg('wwv');
212 for (\$c = \$#in; \$c >= 0; \$c--) {
216 next if \$count < \$from;
218 last if \$count >= \$to; # stop after n
223 $fp->close; # close any open files
225 my $fh = $fp->open($date);
226 for ($count = 0; $count < $to; ) {
231 push @in, [ split '\^' ] if length > 2;
233 eval $eval; # do the search on this file
234 return ("Geomag search error", $@) if $@;
235 last if $count >= $to; # stop after n
237 $fh = $fp->openprev(); # get the next file
246 # the standard log printing interpreting routine.
248 # every line that is printed should call this routine to be actually visualised
250 # Don't really know whether this is the correct place to put this stuff, but where
253 # I get a reference to an array of items
259 my $d = cldate($ref[1]);
260 my ($t) = (gmtime($ref[1]))[2];
262 return sprintf("$d %02d %5d %3d %3d %-37s <%s>", $t, $ref[2], $ref[3], $ref[4], $ref[5], $ref[0]);
266 # read in this month's data
270 my $date = $fp->unixtoj(shift);
271 my $fh = $fp->open($date);
278 push @in, [ split '\^' ] if length > 2;
284 # enter the spot for dup checking and return true if it is already a dup
287 my ($d, $sfi, $k, $a, $text, $call) = @_;
290 return 2 if $d < $main::systime - $dupage;
292 my $dupkey = "W$d|$sfi|$k|$a|$call";
293 return DXDupe::check($dupkey, $main::systime+$dupage);
298 return DXDupe::listdups('W', $dupage, @_);