aboutsummaryrefslogtreecommitdiff
path: root/src/cal/get_tle/check_tle.pl
diff options
context:
space:
mode:
authorJoseph Hunkeler <jhunkeler@gmail.com>2015-03-07 23:42:46 -0500
committerJoseph Hunkeler <jhunkeler@gmail.com>2015-03-07 23:42:46 -0500
commit727498350ee97618f5a23f93addaef43e7d950af (patch)
treeea77abba309c71e438bbf3fd5aa600c86116e7f9 /src/cal/get_tle/check_tle.pl
parentece6c047eacb388ed4d36ce1eb18cc96b47a047e (diff)
downloadcalfuse-727498350ee97618f5a23f93addaef43e7d950af.tar.gz
Massive fixes. Adding the rest of the code.
Diffstat (limited to 'src/cal/get_tle/check_tle.pl')
-rwxr-xr-xsrc/cal/get_tle/check_tle.pl100
1 files changed, 0 insertions, 100 deletions
diff --git a/src/cal/get_tle/check_tle.pl b/src/cal/get_tle/check_tle.pl
deleted file mode 100755
index 483692e..0000000
--- a/src/cal/get_tle/check_tle.pl
+++ /dev/null
@@ -1,100 +0,0 @@
-#!/usr/local/bin/perl
-use FileHandle;
-
-# ***************************************************
-# add_tle.pl
-#
-# This Perl module will read in the latest five orbital elements
-# from the file five.tle (which was created by get_tle.pl) and
-# add any new orbital elements to the file FUSE.TLE.
-#
-# Author: Ed Murphy
-#
-# History: Written July 27, 1999
-#
-# ***************************************************
-
-# Define the file names. old_maintle_filename is not actually opened,
-# but is used in a system call at the end of the program.
-
-$maintle_filename = "/data1/fuse/calfuse/calfiles/FUSE.TLE";
-$outtle_filename = "check_tle.dat";
-
-# Open the files.
-open (TLE_MAINFILE, "<$maintle_filename") || die "Cannot open $maintle_filename";
-open (TLE_OUT, ">$outtle_filename") || die "Cannot open $outtle_filename";
-
-# Read the first two lines from the maintle file to get the date of the
-# most recent set of TLEs. Save these lines for later output.
-while ($line1 = <TLE_MAINFILE>) {
- $line2 = <TLE_MAINFILE>;
- $line3 = <TLE_MAINFILE>;
-
- &parse_lines;
-
- &calculate_a0;
- &calculate_orbs;
-
- printf TLE_OUT "%7.3f %6.2f %6.2f %7.2f %6.2f %7.4f %6.2f %6.2f %8.6f\n",$day, $apogee, $perigee, $semiax, $raan, $incl, $mean_anom, $arg_perig, $eccen;
-
-}
-
- close (TLE_MAINFILE);
- close (TLE_OUT);
-
-### end of Perl script
-
-
-sub parse_lines {
-$year = substr($line2, 18, 2);
-$doy = substr($line2, 20, 3);
-$dayfrac = substr($line2, 23, 9);
-$day=$doy.$dayfrac;
-$epoch = $year.$doy.$dayfrac;
-$incl = substr($line3, 9, 8);
-$raan = substr($line3, 17,8);
-$eccen = substr($line3, 26, 7);
-$eccen = "0." . $eccen;
-$mean_mot = substr($line3, 52, 11);
-$mean_anom = substr($line3, 43, 8);
-$arg_perig = substr($line3, 34, 8);
-
-
-if ($year > 50) {
- $year += 1900;
-} else {
- $year += 2000;
-}
-
-$dayfrac = "0".$dayfrac;
-
-}
-
-sub calculate_a0 {
-
-$mu = 3.986005E5;
-$pi = 3.14159265358979;
-$radian = 0.0174532925200;
-
-$mm=2.0*$pi*$mean_mot/(24.0*60.0);
-$ke=0.74366916E-1;
-$k2=5.413080E-4;
-$a1=($ke/$mm)**(2.0/3.0);
-$inrad = $incl * $radian;
-
-$d1=3.0/2.0*$k2/($a1*$a1)*(3.0*cos($inrad)*cos($inrad)-1.0)/((1.0-$eccen*$eccen)**(3.0/2.0));
-$a0=$a1*(1.0-$d1/3.0-$d1*$d1-134.0/81.0*$d1*$d1*$d1);
-$d0=3.0/2.0*$k2/($a0*$a0)*(3.0*cos($inrad)*cos($inrad)-1.0)/((1.0-$eccen*$eccen)**(3.0/2.0));
-$semiax=$a0/(1-$d0)*6378.135;
-
-}
-
-sub calculate_orbs {
-
-$period = 60.0*24.0/($mean_mot/(1+$d0));
-$cax=$eccen*$semiax;
-$perigee=$semiax-$cax-6378.1;
-$apogee=$semiax+$cax-6378.1;
-
-}
-