X-Git-Url: http://dxcluster.net/gitweb/gitweb.cgi?a=blobdiff_plain;f=perl%2Fconvkeps.pl;h=b591ca6b446fc9f2bc8428d128254f85ef48a5ec;hb=5a53dcc0e3311a7575bbe759f886aa8920a8f825;hp=65e53f40728779ca146748c878e17b2bbe20c8fc;hpb=46960c1b50c3937c073a65422313925f78007a2b;p=spider.git

diff --git a/perl/convkeps.pl b/perl/convkeps.pl
index 65e53f40..b591ca6b 100755
--- a/perl/convkeps.pl
+++ b/perl/convkeps.pl
@@ -59,6 +59,8 @@ my $state = 0;
 my $name;
 my $ref;
 my $line;
+my $count = 0;
+
 my $f = \*STDIN;
 
 while (@ARGV) {
@@ -79,21 +81,23 @@ 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{^-};
 		
-		if (/^\w+/) {
+		if (/^[- \w]+$/) {
 			s/\s/-/g;
-			$name = $_;
+			$name = uc $_;
 			$ref = $keps{$name} = {}; 
 			$state = 2;
 		}
@@ -107,12 +111,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 +131,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 +141,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;
 }