X-Git-Url: http://dxcluster.net/gitweb/gitweb.cgi?a=blobdiff_plain;f=cmd%2Fshow%2Fmuf.pl;h=0c96764f9a3ec5a52c472ea1c28f94868c5ccf67;hb=99da9d68a4fee5b30fa90b679a9e015add2b75c9;hp=65cc21b5883a9f58e70ca7338699384046773795;hpb=8ccb679964e6761c209fe2c3ee7a1084295dca97;p=spider.git diff --git a/cmd/show/muf.pl b/cmd/show/muf.pl index 65cc21b5..0c96764f 100644 --- a/cmd/show/muf.pl +++ b/cmd/show/muf.pl @@ -10,8 +10,18 @@ # my ($self, $line) = @_; -my ($prefix, $hr2) = split /\s+/, $line; +my @f = split /\s+/, $line; + +my $prefix = uc shift @f; return (1, $self->msg('e4')) unless $prefix; +my $lp; +my $hr2; + +while (@f) { + my $f = shift @f; + $lp++ if $f =~ /^l/; + $hr2 = $f if $f =~ /^\d+$/; +} $hr2 = 2 if !$hr2 || $hr2 < 2; $hr2 = 24 if $hr2 > 24; @@ -45,6 +55,16 @@ $b1 *= $d2r; $b2 *= $d2r; $d = ($d / $R); +# handle long path +if ($lp) { + $d = $pi2 - $d; + $b1 += $pi; + $b1 -= $pi2 if ($b1 >= $pi2); + $b2 += $pi; + $b2 -= $pi2 if ($b2 >= $pi2); +} + + my ($hr1, $day, $month) = (gmtime($main::systime))[2,3,4]; $month++; my $flux = Geomag::sfi; @@ -56,7 +76,7 @@ $theta=$theta+2.*$pi if( $theta <= -$pi); $theta=$theta-2.*$pi if( $theta >= $pi); my ($lats, $lons); # subsolar coordinates (rad) -my $dB1 = 20; # transmitter output power (dBW) +my $dB1 = 26; # transmitter output power (dBW) my $delay; # path delay (ms) my $psi; # sun zenith angle (rad) @@ -70,7 +90,7 @@ my $beta1; # elevation angle (rad) my $dhop; # hop great-circle distance (rad) my $height; # height of F layer (km) my $time; # time of day (hour) -my $rsens = -123; # RX sensitivity +my $rsens = -128; # RX sensitivity my @freq = qw(1.8 3.5 7.0 10.1 14.0 18.1 21.0 24.9 28.0 50.0); # working frequencies (MHz)