X-Git-Url: http://dxcluster.net/gitweb/gitweb.cgi?a=blobdiff_plain;f=perl%2Fconvkeps.pl;h=9fda8e820dcaa6f0b36a2f64f606594ec6610af8;hb=2a44647e3c42fb1c532412f6d4dc90481e248ba9;hp=65e53f40728779ca146748c878e17b2bbe20c8fc;hpb=46960c1b50c3937c073a65422313925f78007a2b;p=spider.git diff --git a/perl/convkeps.pl b/perl/convkeps.pl index 65e53f40..9fda8e82 100755 --- a/perl/convkeps.pl +++ b/perl/convkeps.pl @@ -59,6 +59,16 @@ my $state = 0; my $name; my $ref; my $line; +my $count = 0; + +my %lookup = ( + 'AO-5' => 'AO-05', + 'AO-6' => 'AO-06', + 'AO-7' => 'AO-07', + 'AO-8' => 'AO-08', + +); + my $f = \*STDIN; while (@ARGV) { @@ -79,21 +89,26 @@ while (@ARGV) { while (<$f>) { ++$line; +# print; chomp; + last if m{^-}; + s/^\s+//; - s/\s+$//; + s/[\s\r]+$//; next unless $_; last if m{^/EX}i; - last if m{^-}; - if ($state == 0 && /^TO ALL/) { + if ($state == 0 && /^Decode/i) { $state = 1; } elsif ($state == 1) { - last if m{^/EX/i}; + last if m{^-}; + next if m{^To\s+all}i; - if (/^\w+/) { - s/\s/-/g; - $name = $_; + if (/^([- \w]+)(?:\s+\[[-+\w]\])?$/) { + my $n = uc $1; + $n =~ s/\s/-/g; + $name = $lookup{$n}; + $name ||= $n; $ref = $keps{$name} = {}; $state = 2; } @@ -107,12 +122,12 @@ while (<$f>) { $ref->{mm2} = genenum($mm2); $ref->{bstar} = genenum($bstar); $ref->{elset} = $elset - 0; -# print "$id $number $epoch $decay $mm2 $bstar $elset\n"; -# print "mm2: $ref->{mm2} bstar: $ref->{bstar}\n"; + #print "$id $number $epoch $decay $mm2 $bstar $elset\n"; + #print "mm2: $ref->{mm2} bstar: $ref->{bstar}\n"; $state = 3; } else { - print "out of order on line $line\n"; + #print "out of order on line $line\n"; undef $ref; delete $keps{$name}; $state = 1; @@ -127,8 +142,9 @@ while (<$f>) { $ref->{argperigee} = $peri - 0; $ref->{raan} = $raan - 0; $ref->{orbit} = $orbit - 0; + $count++; } else { - print "out of order on line $line\n"; + #print "out of order on line $line\n"; delete $keps{$name}; } undef $ref; @@ -136,22 +152,28 @@ while (<$f>) { } } -my $dd = new Data::Dumper([\%keps], [qw(*keps)]); -$dd->Indent(1); -$dd->Quotekeys(0); -open(OUT, ">$fn") or die "$fn $!"; -print OUT "#\n# this file is automatically produced by convkeps.pl\n#\n"; -print OUT "# Last update: ", scalar gmtime, "\n#\n"; -print OUT "\npackage Sun;\n\n"; -print OUT $dd->Dumpxs; -print OUT "1;\n"; -close(OUT); +if ($count) { + my $dd = new Data::Dumper([\%keps], [qw(*keps)]); + $dd->Indent(1); + $dd->Quotekeys(0); + open(OUT, ">$fn") or die "$fn $!"; + print OUT "#\n# this file is automatically produced by convkeps.pl\n#\n"; + print OUT "# Last update: ", scalar gmtime, "\n#\n"; + print OUT "\npackage Sun;\n\n"; + print OUT $dd->Dumpxs; + print OUT "1;\n"; + close(OUT); +} + +print "$count keps converted\n"; +exit($count ? 0 : -1); # convert (+/-)00000-0 to (+/-).00000e-0 sub genenum { my ($sign, $frac, $esign, $exp) = unpack "aa5aa", shift; + $esign = '+' if $esign eq ' '; my $n = $sign . "." . $frac . 'e' . $esign . $exp; return $n - 0; }