Ysurac / FlightAirMap

require/libs/Predict/Predict/PassDetail.php

Indentation +21 added lines, -21 removed lines

@@ -12,26 +12,26 @@
  */
 class Predict_PassDetail
 {
-    public $time;   /*!< time in "jul_utc" */
-    public $pos;    /*!< Raw unprocessed position at time */
-    public $vel;    /*!< Raw unprocessed velocity at time */
-    public $velo;
-    public $az;
-    public $el;
-    public $range;
-    public $range_rate;
-    public $lat;
-    public $lon;
-    public $alt;
-    public $ma;
-    public $phase;
-    public $footprint;
-    public $vis;
-    public $orbit;
+	public $time;   /*!< time in "jul_utc" */
+	public $pos;    /*!< Raw unprocessed position at time */
+	public $vel;    /*!< Raw unprocessed velocity at time */
+	public $velo;
+	public $az;
+	public $el;
+	public $range;
+	public $range_rate;
+	public $lat;
+	public $lon;
+	public $alt;
+	public $ma;
+	public $phase;
+	public $footprint;
+	public $vis;
+	public $orbit;
 
-    public function __construct()
-    {
-        $this->pos = new Predict_Vector();
-        $this->vel = new Predict_Vector();
-    }
+	public function __construct()
+	{
+		$this->pos = new Predict_Vector();
+		$this->vel = new Predict_Vector();
+	}
 }

require/libs/Predict/Predict/ObsSet.php

Indentation +4 added lines, -4 removed lines

@@ -5,8 +5,8 @@
  */
 class Predict_ObsSet
 {
-    public $az         = 0.0;  /*!< Azimuth [deg] */
-    public $el         = 0.0;  /*!< Elevation [deg] */
-    public $range      = 0.0;  /*!< Range [km] */
-    public $range_rate = 0.0;  /*!< Velocity [km/sec] */
+	public $az         = 0.0;  /*!< Azimuth [deg] */
+	public $el         = 0.0;  /*!< Elevation [deg] */
+	public $range      = 0.0;  /*!< Range [km] */
+	public $range_rate = 0.0;  /*!< Velocity [km/sec] */
 }

require/libs/Predict/Predict/TLE.php

Indentation +192 added lines, -192 removed lines

@@ -12,221 +12,221 @@
  */
 class Predict_TLE
 {
-    public $header;     /* Header line of TLE file */
-    public $line1;      /* Line 1 of TLE */
-    public $line2;      /* Line 2 of TLE */
-    public $epoch;      /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */
-    public $epoch_year; /*!< Epoch: year */
-    public $epoch_day;  /*!< Epoch: day of year */
-    public $epoch_fod;  /*!< Epoch: Fraction of day. */
-    public $xndt2o;     /*!< 1. time derivative of mean motion */
-    public $xndd6o;     /*!< 2. time derivative of mean motion */
-    public $bstar;      /*!< Bstar drag coefficient. */
-    public $xincl;      /*!< Inclination */
-    public $xnodeo;     /*!< R.A.A.N. */
-    public $eo;         /*!< Eccentricity */
-    public $omegao;     /*!< argument of perigee */
-    public $xmo;        /*!< mean anomaly */
-    public $xno;        /*!< mean motion */
-
-    public $catnr;      /*!< Catalogue Number.  */
-    public $elset;      /*!< Element Set number. */
-    public $revnum;     /*!< Revolution Number at epoch. */
-
-    public $sat_name;   /*!< Satellite name string. */
-    public $idesg;      /*!< International Designator. */
-    public $status;     /*!< Operational status. */
-
-    /* values needed for squint calculations */
-    public $xincl1;
-    public $xnodeo1;
-    public $omegao1;
-
-
-    /* Converts the strings in a raw two-line element set  */
-    /* to their intended numerical values. No processing   */
-    /* of these values is done, e.g. from deg to rads etc. */
-    /* This is done in the select_ephemeris() function.    */
-    public function __construct($header, $line1, $line2)
-    {
-        if (!$this->Good_Elements($line1, $line2)) {
-            throw new Predict_Exception('Invalid TLE contents');
-        }
-
-        $this->header = $header;
-        $this->line1  = $line1;
-        $this->line2  = $line2;
-
-        /** Decode Card 1 **/
-        /* Satellite's catalogue number */
-        $this->catnr = (int) substr($line1, 2, 5);
-
-        /* International Designator for satellite */
-        $this->idesg = substr($line1, 9, 8);
-
-        /* Epoch time; this is the complete, unconverted epoch. */
-        /* Replace spaces with 0 before casting, as leading spaces are allowed */
-        $this->epoch = (float) str_replace(' ', '0', substr($line1, 18, 14));
-
-        /* Now, convert the epoch time into year, day
+	public $header;     /* Header line of TLE file */
+	public $line1;      /* Line 1 of TLE */
+	public $line2;      /* Line 2 of TLE */
+	public $epoch;      /*!< Epoch Time in NORAD TLE format YYDDD.FFFFFFFF */
+	public $epoch_year; /*!< Epoch: year */
+	public $epoch_day;  /*!< Epoch: day of year */
+	public $epoch_fod;  /*!< Epoch: Fraction of day. */
+	public $xndt2o;     /*!< 1. time derivative of mean motion */
+	public $xndd6o;     /*!< 2. time derivative of mean motion */
+	public $bstar;      /*!< Bstar drag coefficient. */
+	public $xincl;      /*!< Inclination */
+	public $xnodeo;     /*!< R.A.A.N. */
+	public $eo;         /*!< Eccentricity */
+	public $omegao;     /*!< argument of perigee */
+	public $xmo;        /*!< mean anomaly */
+	public $xno;        /*!< mean motion */
+
+	public $catnr;      /*!< Catalogue Number.  */
+	public $elset;      /*!< Element Set number. */
+	public $revnum;     /*!< Revolution Number at epoch. */
+
+	public $sat_name;   /*!< Satellite name string. */
+	public $idesg;      /*!< International Designator. */
+	public $status;     /*!< Operational status. */
+
+	/* values needed for squint calculations */
+	public $xincl1;
+	public $xnodeo1;
+	public $omegao1;
+
+
+	/* Converts the strings in a raw two-line element set  */
+	/* to their intended numerical values. No processing   */
+	/* of these values is done, e.g. from deg to rads etc. */
+	/* This is done in the select_ephemeris() function.    */
+	public function __construct($header, $line1, $line2)
+	{
+		if (!$this->Good_Elements($line1, $line2)) {
+			throw new Predict_Exception('Invalid TLE contents');
+		}
+
+		$this->header = $header;
+		$this->line1  = $line1;
+		$this->line2  = $line2;
+
+		/** Decode Card 1 **/
+		/* Satellite's catalogue number */
+		$this->catnr = (int) substr($line1, 2, 5);
+
+		/* International Designator for satellite */
+		$this->idesg = substr($line1, 9, 8);
+
+		/* Epoch time; this is the complete, unconverted epoch. */
+		/* Replace spaces with 0 before casting, as leading spaces are allowed */
+		$this->epoch = (float) str_replace(' ', '0', substr($line1, 18, 14));
+
+		/* Now, convert the epoch time into year, day
            and fraction of day, according to:
 
            YYDDD.FFFFFFFF
         */
 
-        // Adjust for 2 digit year through 2056
-        $this->epoch_year = (int) substr($line1, 18, 2);
-        if ($this->epoch_year > 56) {
-            $this->epoch_year = $this->epoch_year + 1900;
-        } else {
-            $this->epoch_year = $this->epoch_year + 2000;
-        }
+		// Adjust for 2 digit year through 2056
+		$this->epoch_year = (int) substr($line1, 18, 2);
+		if ($this->epoch_year > 56) {
+			$this->epoch_year = $this->epoch_year + 1900;
+		} else {
+			$this->epoch_year = $this->epoch_year + 2000;
+		}
 
-        /* Epoch day */
-        $this->epoch_day = (int) substr($line1, 20, 3);
+		/* Epoch day */
+		$this->epoch_day = (int) substr($line1, 20, 3);
 
-        /* Epoch fraction of day */
-        $this->epoch_fod = (float) substr($line1, 23, 9);
+		/* Epoch fraction of day */
+		$this->epoch_fod = (float) substr($line1, 23, 9);
 
 
-        /* Satellite's First Time Derivative */
-        $this->xndt2o = (float) substr($line1, 33, 10);
+		/* Satellite's First Time Derivative */
+		$this->xndt2o = (float) substr($line1, 33, 10);
 
-        /* Satellite's Second Time Derivative */
-        $this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2));
+		/* Satellite's Second Time Derivative */
+		$this->xndd6o = (float) (substr($line1, 44, 1) . '.' . substr($line1, 45, 5) . 'E' . substr($line1, 50, 2));
 
-        /* Satellite's bstar drag term
+		/* Satellite's bstar drag term
            FIXME: How about buff[0] ????
         */
-        $this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2));
+		$this->bstar = (float) (substr($line1, 53, 1) . '.' . substr($line1, 54, 5) . 'E' . substr($line1, 59, 2));
 
-        /* Element Number */
-        $this->elset = (int) substr($line1, 64, 4);
+		/* Element Number */
+		$this->elset = (int) substr($line1, 64, 4);
 
-        /** Decode Card 2 **/
-        /* Satellite's Orbital Inclination (degrees) */
-        $this->xincl = (float) substr($line2, 8, 8);
+		/** Decode Card 2 **/
+		/* Satellite's Orbital Inclination (degrees) */
+		$this->xincl = (float) substr($line2, 8, 8);
 
-        /* Satellite's RAAN (degrees) */
-        $this->xnodeo = (float) substr($line2, 17, 8);
+		/* Satellite's RAAN (degrees) */
+		$this->xnodeo = (float) substr($line2, 17, 8);
 
-        /* Satellite's Orbital Eccentricity */
-        $this->eo = (float) ('.' . substr($line2, 26, 7));
+		/* Satellite's Orbital Eccentricity */
+		$this->eo = (float) ('.' . substr($line2, 26, 7));
 
-        /* Satellite's Argument of Perigee (degrees) */
-        $this->omegao = (float) substr($line2, 34, 8);
+		/* Satellite's Argument of Perigee (degrees) */
+		$this->omegao = (float) substr($line2, 34, 8);
 
-        /* Satellite's Mean Anomaly of Orbit (degrees) */
-        $this->xmo = (float) substr($line2, 43, 8);
+		/* Satellite's Mean Anomaly of Orbit (degrees) */
+		$this->xmo = (float) substr($line2, 43, 8);
 
-        /* Satellite's Mean Motion (rev/day) */
-        $this->xno = (float) substr($line2, 52, 11);
+		/* Satellite's Mean Motion (rev/day) */
+		$this->xno = (float) substr($line2, 52, 11);
 
-        /* Satellite's Revolution number at epoch */
-        $this->revnum = (float) substr($line2, 63, 5);
-    }
+		/* Satellite's Revolution number at epoch */
+		$this->revnum = (float) substr($line2, 63, 5);
+	}
 
-    /* Calculates the checksum mod 10 of a line from a TLE set and */
-    /* returns true if it compares with checksum in column 68, else false.*/
-    /* tle_set is a character string holding the two lines read    */
-    /* from a text file containing NASA format Keplerian elements. */
-    /* NOTE!!! The stuff about two lines is not quite true.
+	/* Calculates the checksum mod 10 of a line from a TLE set and */
+	/* returns true if it compares with checksum in column 68, else false.*/
+	/* tle_set is a character string holding the two lines read    */
+	/* from a text file containing NASA format Keplerian elements. */
+	/* NOTE!!! The stuff about two lines is not quite true.
        The function assumes that tle_set[0] is the begining
        of the line and that there are 68 elements - see the consumer
     */
-    public function Checksum_Good($tle_set)
-    {
-        if (strlen($tle_set) < 69) {
-            return false;
-        }
-
-        $checksum = 0;
-
-        for ($i = 0; $i < 68; $i++) {
-            if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) {
-                $value = $tle_set[$i] - '0';
-            } else if ($tle_set[$i] == '-' ) {
-                $value = 1;
-            } else {
-                $value = 0;
-            }
-
-            $checksum += $value;
-        }
-
-        $checksum   %= 10;
-        $check_digit = $tle_set[68] - '0';
-
-        return $checksum == $check_digit;
-    }
-
-    /* Carries out various checks on a TLE set to verify its validity */
-    /* $line1 is the first line of the TLE, $line2 is the second line */
-    /* from a text file containing NASA format Keplerian elements. */
-    public function Good_Elements($line1, $line2)
-    {
-        /* Verify checksum of both lines of a TLE set */
-        if (!$this->Checksum_Good($line1) || !$this->Checksum_Good($line2)) {
-            return false;
-        }
-
-        /* Check the line number of each line */
-        if (($line1[0] != '1') || ($line2[0] != '2')) {
-            return false;
-        }
-
-        /* Verify that Satellite Number is same in both lines */
-        if (strncmp($line1[2], $line2[2], 5) != 0) {
-            return false;
-        }
-
-        /* Check that various elements are in the right place */
-        if (($line1[23] != '.') ||
-            ($line1[34] != '.') ||
-            ($line2[11] != '.') ||
-            ($line2[20] != '.') ||
-            ($line2[37] != '.') ||
-            ($line2[46] != '.') ||
-            ($line2[54] != '.') ||
-            (strncmp(substr($line1, 61), ' 0 ', 3) != 0)) {
-
-            return false;
-        }
-
-        return true;
-    }
-
-    /**
-     * A function to allow checksum creation of a line.  This is driven by
-     * the fact that some TLEs from SpaceTrack are missing checksum numbers.
-     * You can use this to create a checksum for a line, but you should
-     * probably have confidence that the TLE data itself is good.  YMMV.
-     *
-     * @throws Predict_Exception if the line is not exactly 68 chars
-     * @return string
-     */
-    static public function createChecksum($line)
-    {
-        if (strlen($line) != 68) {
-            throw Predict_Exception('Invalid line, needs to e 68 chars');
-        }
-
-        $checksum = 0;
-
-        for ($i = 0; $i < 68; $i++) {
-            if (($line[$i] >= '0') && ($line[$i] <= '9')) {
-                $value = (int) $line[$i];
-            } else if ($line[$i] == '-' ) {
-                $value = 1;
-            } else {
-                $value = 0;
-            }
-
-            $checksum += $value;
-        }
-
-        $checksum %= 10;
-
-        return $checksum;
-    }
+	public function Checksum_Good($tle_set)
+	{
+		if (strlen($tle_set) < 69) {
+			return false;
+		}
+
+		$checksum = 0;
+
+		for ($i = 0; $i < 68; $i++) {
+			if (($tle_set[$i] >= '0') && ($tle_set[$i] <= '9')) {
+				$value = $tle_set[$i] - '0';
+			} else if ($tle_set[$i] == '-' ) {
+				$value = 1;
+			} else {
+				$value = 0;
+			}
+
+			$checksum += $value;
+		}
+
+		$checksum   %= 10;
+		$check_digit = $tle_set[68] - '0';
+
+		return $checksum == $check_digit;
+	}
+
+	/* Carries out various checks on a TLE set to verify its validity */
+	/* $line1 is the first line of the TLE, $line2 is the second line */
+	/* from a text file containing NASA format Keplerian elements. */
+	public function Good_Elements($line1, $line2)
+	{
+		/* Verify checksum of both lines of a TLE set */
+		if (!$this->Checksum_Good($line1) || !$this->Checksum_Good($line2)) {
+			return false;
+		}
+
+		/* Check the line number of each line */
+		if (($line1[0] != '1') || ($line2[0] != '2')) {
+			return false;
+		}
+
+		/* Verify that Satellite Number is same in both lines */
+		if (strncmp($line1[2], $line2[2], 5) != 0) {
+			return false;
+		}
+
+		/* Check that various elements are in the right place */
+		if (($line1[23] != '.') ||
+			($line1[34] != '.') ||
+			($line2[11] != '.') ||
+			($line2[20] != '.') ||
+			($line2[37] != '.') ||
+			($line2[46] != '.') ||
+			($line2[54] != '.') ||
+			(strncmp(substr($line1, 61), ' 0 ', 3) != 0)) {
+
+			return false;
+		}
+
+		return true;
+	}
+
+	/**
+	 * A function to allow checksum creation of a line.  This is driven by
+	 * the fact that some TLEs from SpaceTrack are missing checksum numbers.
+	 * You can use this to create a checksum for a line, but you should
+	 * probably have confidence that the TLE data itself is good.  YMMV.
+	 *
+	 * @throws Predict_Exception if the line is not exactly 68 chars
+	 * @return string
+	 */
+	static public function createChecksum($line)
+	{
+		if (strlen($line) != 68) {
+			throw Predict_Exception('Invalid line, needs to e 68 chars');
+		}
+
+		$checksum = 0;
+
+		for ($i = 0; $i < 68; $i++) {
+			if (($line[$i] >= '0') && ($line[$i] <= '9')) {
+				$value = (int) $line[$i];
+			} else if ($line[$i] == '-' ) {
+				$value = 1;
+			} else {
+				$value = 0;
+			}
+
+			$checksum += $value;
+		}
+
+		$checksum %= 10;
+
+		return $checksum;
+	}
 }

require/libs/Predict/Predict/Geodetic.php

Indentation +4 added lines, -4 removed lines

@@ -6,8 +6,8 @@
  */
 class Predict_Geodetic
 {
-    public $lat; /*!< Lattitude [rad] */
-    public $lon; /*!< Longitude [rad] */
-    public $alt; /*!< Altitude [km] */
-    public $theta;
+	public $lat; /*!< Lattitude [rad] */
+	public $lon; /*!< Longitude [rad] */
+	public $alt; /*!< Altitude [km] */
+	public $theta;
 }

require/libs/Predict/Predict/SGSDPStatic.php

Indentation +28 added lines, -28 removed lines

@@ -8,32 +8,32 @@
  */
 class Predict_SGSDPStatic
 {
-    public $aodp;
-    public $aycof;
-    public $c1;
-    public $c4;
-    public $c5;
-    public $cosio;
-    public $d2;
-    public $d3;
-    public $d4;
-    public $delmo;
-    public $omgcof;
-    public $eta;
-    public $omgdot;
-    public $sinio;
-    public $xnodp;
-    public $sinmo;
-    public $t2cof;
-    public $t3cof;
-    public $t4cof;
-    public $t5cof;
-    public $x1mth2;
-    public $x3thm1;
-    public $x7thm1;
-    public $xmcof;
-    public $xmdot;
-    public $xnodcf;
-    public $xnodot;
-    public $xlcof;
+	public $aodp;
+	public $aycof;
+	public $c1;
+	public $c4;
+	public $c5;
+	public $cosio;
+	public $d2;
+	public $d3;
+	public $d4;
+	public $delmo;
+	public $omgcof;
+	public $eta;
+	public $omgdot;
+	public $sinio;
+	public $xnodp;
+	public $sinmo;
+	public $t2cof;
+	public $t3cof;
+	public $t4cof;
+	public $t5cof;
+	public $x1mth2;
+	public $x3thm1;
+	public $x7thm1;
+	public $xmcof;
+	public $xmdot;
+	public $xnodcf;
+	public $xnodot;
+	public $xlcof;
 }

require/libs/Predict/Predict/QTH.php

Indentation +9 added lines, -9 removed lines

@@ -7,14 +7,14 @@
  */
 class Predict_QTH
 {
-     public $name;   /*!< Name, eg. callsign. */
-     public $loc;    /*!< Location, eg City, Country. */
-     public $desc;   /*!< Short description. */
-     public $lat;    /*!< Latitude in dec. deg. North. */
-     public $lon;    /*!< Longitude in dec. deg. East. */
-     public $alt;    /*!< Altitude above sea level in meters. */
-     public $qra;    /*!< QRA locator */
-     public $wx;     /*!< Weather station code (4 chars). */
+	 public $name;   /*!< Name, eg. callsign. */
+	 public $loc;    /*!< Location, eg City, Country. */
+	 public $desc;   /*!< Short description. */
+	 public $lat;    /*!< Latitude in dec. deg. North. */
+	 public $lon;    /*!< Longitude in dec. deg. East. */
+	 public $alt;    /*!< Altitude above sea level in meters. */
+	 public $qra;    /*!< QRA locator */
+	 public $wx;     /*!< Weather station code (4 chars). */
 
-     public $data;   /*!< Raw data from cfg file. */
+	 public $data;   /*!< Raw data from cfg file. */
 }

require/libs/Predict/Predict/Sat.php

Indentation +296 added lines, -296 removed lines

@@ -22,304 +22,304 @@
  */
 class Predict_Sat
 {
-    // Fifth root of a hundred, used for magnitude calculation
-    const POGSONS_RATIO = 2.5118864315096;
-
-    public $name     = null;
-    public $nickname = null;
-    public $website  = null;
-
-    public $tle      = null;   /*!< Keplerian elements */
-    public $flags    = 0;      /*!< Flags for algo ctrl */
-    public $sgps     = null;
-    public $dps      = null;
-    public $deep_arg = null;
-    public $pos      = null;   /*!< Raw position and range */
-    public $vel      = null;   /*!< Raw velocity */
-
-    /*** FIXME: REMOVE */
-    public $bearing = null;   /*!< Az, El, range and vel */
-    public $astro   = null;   /*!< Ra and Decl */
-    /*** END */
-
-    /* time keeping fields */
-    public $jul_epoch = null;
-    public $jul_utc   = null;
-    public $tsince    = null;
-    public $aos       = null;    /*!< Next AOS. */
-    public $los       = null;    /*!< Next LOS */
-
-    public $az         = null;   /*!< Azimuth [deg] */
-    public $el         = null;   /*!< Elevation [deg] */
-    public $range      = null;   /*!< Range [km] */
-    public $range_rate = null;   /*!< Range Rate [km/sec] */
-    public $ra         = null;   /*!< Right Ascension [deg] */
-    public $dec        = null;   /*!< Declination [deg] */
-    public $ssplat     = null;   /*!< SSP latitude [deg] */
-    public $ssplon     = null;   /*!< SSP longitude [deg] */
-    public $alt        = null;   /*!< altitude [km] */
-    public $velo       = null;   /*!< velocity [km/s] */
-    public $ma         = null;   /*!< mean anomaly */
-    public $footprint  = null;   /*!< footprint */
-    public $phase      = null;   /*!< orbit phase */
-    public $meanmo     = null;   /*!< mean motion kept in rev/day */
-    public $orbit      = null;   /*!< orbit number */
-    public $otype      = null;   /*!< orbit type. */
-
-    public function __construct(Predict_TLE $tle)
-    {
-        $headerParts    = explode(' ', $tle->header);
-        $this->name     = $headerParts[0];
-        $this->nickname = $this->name;
-        $this->tle      = $tle;
-        $this->pos      = new Predict_Vector();
-        $this->vel      = new Predict_Vector();
-        $this->sgps     = new Predict_SGSDPStatic();
-        $this->deep_arg = new Predict_DeepArg();
-        $this->dps      = new Predict_DeepStatic();
-
-        $this->select_ephemeris();
-        $this->sat_data_init_sat($this);
-    }
-
-    /* Selects the apropriate ephemeris type to be used */
-    /* for predictions according to the data in the TLE */
-    /* It also processes values in the tle set so that  */
-    /* they are apropriate for the sgp4/sdp4 routines   */
-    public function select_ephemeris()
-    {
-        /* Preprocess tle set */
-        $this->tle->xnodeo *= Predict::de2ra;
-        $this->tle->omegao *= Predict::de2ra;
-        $this->tle->xmo    *= Predict::de2ra;
-        $this->tle->xincl  *= Predict::de2ra;
-        $temp = Predict::twopi / Predict::xmnpda / Predict::xmnpda;
-
-        /* store mean motion before conversion */
-        $this->meanmo       = $this->tle->xno;
-        $this->tle->xno     = $this->tle->xno * $temp * Predict::xmnpda;
-        $this->tle->xndt2o *= $temp;
-        $this->tle->xndd6o  = $this->tle->xndd6o * $temp / Predict::xmnpda;
-        $this->tle->bstar  /= Predict::ae;
-
-        /* Period > 225 minutes is deep space */
-        $dd1 = Predict::xke / $this->tle->xno;
-        $dd2 = Predict::tothrd;
-        $a1 = pow($dd1, $dd2);
-        $r1 = cos($this->tle->xincl);
-        $dd1 = 1.0 - $this->tle->eo * $this->tle->eo;
-        $temp = Predict::ck2 * 1.5 * ($r1 * $r1 * 3.0 - 1.0) / pow($dd1, 1.5);
-        $del1 = $temp / ($a1 * $a1);
-        $ao = $a1 * (1.0 - $del1 * (Predict::tothrd * 0.5 + $del1 *
-                                 ($del1 * 1.654320987654321 + 1.0)));
-        $delo = $temp / ($ao * $ao);
-        $xnodp = $this->tle->xno / ($delo + 1.0);
-
-        /* Select a deep-space/near-earth ephemeris */
-        if (Predict::twopi / $xnodp / Predict::xmnpda >= .15625) {
-            $this->flags |= Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
-        } else {
-            $this->flags &= ~Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
-        }
-    }
-
-    /** Initialise satellite data.
-     *  @param sat The satellite to initialise.
-     *  @param qth Optional QTH info, use (0,0) if NULL.
-     *
-     * This function calculates the satellite data at t = 0, ie. epoch time
-     * The function is called automatically by gtk_sat_data_read_sat.
-     */
-    public function sat_data_init_sat(Predict_Sat $sat, Predict_QTH $qth = null)
-    {
-        $obs_geodetic = new Predict_Geodetic();
-        $obs_set = new Predict_ObsSet();
-        $sat_geodetic = new Predict_Geodetic();
-        /* double jul_utc, age; */
-
-        $jul_utc = Predict_Time::Julian_Date_of_Epoch($sat->tle->epoch); // => tsince = 0.0
-        $sat->jul_epoch = $jul_utc;
-
-        /* initialise observer location */
-        if ($qth != null) {
-            $obs_geodetic->lon = $qth->lon * Predict::de2ra;
-            $obs_geodetic->lat = $qth->lat * Predict::de2ra;
-            $obs_geodetic->alt = $qth->alt / 1000.0;
-            $obs_geodetic->theta = 0;
-        }
-        else {
-            $obs_geodetic->lon = 0.0;
-            $obs_geodetic->lat = 0.0;
-            $obs_geodetic->alt = 0.0;
-            $obs_geodetic->theta = 0;
-        }
-
-        /* execute computations */
-        $sdpsgp = Predict_SGPSDP::getInstance($sat);
-        if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
-            $sdpsgp->SDP4($sat, 0.0);
-        } else {
-            $sdpsgp->SGP4($sat, 0.0);
-        }
-
-        /* scale position and velocity to km and km/sec */
-        Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
-
-        /* get the velocity of the satellite */
-        $sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
-        $sat->velo = $sat->vel->w;
-        Predict_SGPObs::Calculate_Obs($jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
-        Predict_SGPObs::Calculate_LatLonAlt($jul_utc, $sat->pos, $sat_geodetic);
-
-        while ($sat_geodetic->lon < -Predict::pi) {
-            $sat_geodetic->lon += Predict::twopi;
-        }
-
-        while ($sat_geodetic->lon > Predict::pi) {
-            $sat_geodetic->lon -= Predict::twopi;
-        }
-
-        $sat->az = Predict_Math::Degrees($obs_set->az);
-        $sat->el = Predict_Math::Degrees($obs_set->el);
-        $sat->range = $obs_set->range;
-        $sat->range_rate = $obs_set->range_rate;
-        $sat->ssplat = Predict_Math::Degrees($sat_geodetic->lat);
-        $sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
-        $sat->alt = $sat_geodetic->alt;
-        $sat->ma = Predict_Math::Degrees($sat->phase);
-        $sat->ma *= 256.0 / 360.0;
-        $sat->footprint = 2.0 * Predict::xkmper * acos (Predict::xkmper/$sat->pos->w);
-        $age = 0.0;
-        $sat->orbit = floor(($sat->tle->xno * Predict::xmnpda / Predict::twopi +
-                                   $age * $sat->tle->bstar * Predict::ae) * $age +
-                                  $sat->tle->xmo / Predict::twopi) + $sat->tle->revnum - 1;
-
-        /* orbit type */
-        $sat->otype = $sat->get_orbit_type($sat);
-    }
-
-    public function get_orbit_type(Predict_Sat $sat)
-    {
-         $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
-
-         if ($this->geostationary($sat)) {
-              $orbit = Predict_SGPSDP::ORBIT_TYPE_GEO;
-         } else if ($this->decayed($sat)) {
-              $orbit = Predict_SGPSDP::ORBIT_TYPE_DECAYED;
-         } else {
-              $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
-         }
-
-         return $orbit;
-    }
-
-
-    /** Determinte whether satellite is in geostationary orbit.
-     *  @author John A. Magliacane, KD2BD
-     *  @param sat Pointer to satellite data.
-     *  @return TRUE if the satellite appears to be in geostationary orbit,
-     *          FALSE otherwise.
-     *
-     * A satellite is in geostationary orbit if
-     *
-     *     fabs (sat.meanmotion - 1.0027) < 0.0002
-     *
-     * Note: Appearantly, the mean motion can deviate much more from 1.0027 than 0.0002
-     */
-    public function geostationary(Predict_Sat $sat)
-    {
-         if (abs($sat->meanmo - 1.0027) < 0.0002) {
-              return true;
-         } else {
-              return false;
-        }
-    }
-
-
-    /** Determine whether satellite has decayed.
-     *  @author John A. Magliacane, KD2BD
-     *  @author Alexandru Csete, OZ9AEC
-     *  @param sat Pointer to satellite data.
-     *  @return TRUE if the satellite appears to have decayed, FALSE otherwise.
-     *  @bug Modified version of the predict code but it is not tested.
-     *
-     * A satellite is decayed if
-     *
-     *    satepoch + ((16.666666 - sat.meanmo) / (10.0*fabs(sat.drag))) < "now"
-     *
-     */
-    public function decayed(Predict_Sat $sat)
-    {
-        /* tle.xndt2o/(twopi/xmnpda/xmnpda) is the value before converted the
+	// Fifth root of a hundred, used for magnitude calculation
+	const POGSONS_RATIO = 2.5118864315096;
+
+	public $name     = null;
+	public $nickname = null;
+	public $website  = null;
+
+	public $tle      = null;   /*!< Keplerian elements */
+	public $flags    = 0;      /*!< Flags for algo ctrl */
+	public $sgps     = null;
+	public $dps      = null;
+	public $deep_arg = null;
+	public $pos      = null;   /*!< Raw position and range */
+	public $vel      = null;   /*!< Raw velocity */
+
+	/*** FIXME: REMOVE */
+	public $bearing = null;   /*!< Az, El, range and vel */
+	public $astro   = null;   /*!< Ra and Decl */
+	/*** END */
+
+	/* time keeping fields */
+	public $jul_epoch = null;
+	public $jul_utc   = null;
+	public $tsince    = null;
+	public $aos       = null;    /*!< Next AOS. */
+	public $los       = null;    /*!< Next LOS */
+
+	public $az         = null;   /*!< Azimuth [deg] */
+	public $el         = null;   /*!< Elevation [deg] */
+	public $range      = null;   /*!< Range [km] */
+	public $range_rate = null;   /*!< Range Rate [km/sec] */
+	public $ra         = null;   /*!< Right Ascension [deg] */
+	public $dec        = null;   /*!< Declination [deg] */
+	public $ssplat     = null;   /*!< SSP latitude [deg] */
+	public $ssplon     = null;   /*!< SSP longitude [deg] */
+	public $alt        = null;   /*!< altitude [km] */
+	public $velo       = null;   /*!< velocity [km/s] */
+	public $ma         = null;   /*!< mean anomaly */
+	public $footprint  = null;   /*!< footprint */
+	public $phase      = null;   /*!< orbit phase */
+	public $meanmo     = null;   /*!< mean motion kept in rev/day */
+	public $orbit      = null;   /*!< orbit number */
+	public $otype      = null;   /*!< orbit type. */
+
+	public function __construct(Predict_TLE $tle)
+	{
+		$headerParts    = explode(' ', $tle->header);
+		$this->name     = $headerParts[0];
+		$this->nickname = $this->name;
+		$this->tle      = $tle;
+		$this->pos      = new Predict_Vector();
+		$this->vel      = new Predict_Vector();
+		$this->sgps     = new Predict_SGSDPStatic();
+		$this->deep_arg = new Predict_DeepArg();
+		$this->dps      = new Predict_DeepStatic();
+
+		$this->select_ephemeris();
+		$this->sat_data_init_sat($this);
+	}
+
+	/* Selects the apropriate ephemeris type to be used */
+	/* for predictions according to the data in the TLE */
+	/* It also processes values in the tle set so that  */
+	/* they are apropriate for the sgp4/sdp4 routines   */
+	public function select_ephemeris()
+	{
+		/* Preprocess tle set */
+		$this->tle->xnodeo *= Predict::de2ra;
+		$this->tle->omegao *= Predict::de2ra;
+		$this->tle->xmo    *= Predict::de2ra;
+		$this->tle->xincl  *= Predict::de2ra;
+		$temp = Predict::twopi / Predict::xmnpda / Predict::xmnpda;
+
+		/* store mean motion before conversion */
+		$this->meanmo       = $this->tle->xno;
+		$this->tle->xno     = $this->tle->xno * $temp * Predict::xmnpda;
+		$this->tle->xndt2o *= $temp;
+		$this->tle->xndd6o  = $this->tle->xndd6o * $temp / Predict::xmnpda;
+		$this->tle->bstar  /= Predict::ae;
+
+		/* Period > 225 minutes is deep space */
+		$dd1 = Predict::xke / $this->tle->xno;
+		$dd2 = Predict::tothrd;
+		$a1 = pow($dd1, $dd2);
+		$r1 = cos($this->tle->xincl);
+		$dd1 = 1.0 - $this->tle->eo * $this->tle->eo;
+		$temp = Predict::ck2 * 1.5 * ($r1 * $r1 * 3.0 - 1.0) / pow($dd1, 1.5);
+		$del1 = $temp / ($a1 * $a1);
+		$ao = $a1 * (1.0 - $del1 * (Predict::tothrd * 0.5 + $del1 *
+								 ($del1 * 1.654320987654321 + 1.0)));
+		$delo = $temp / ($ao * $ao);
+		$xnodp = $this->tle->xno / ($delo + 1.0);
+
+		/* Select a deep-space/near-earth ephemeris */
+		if (Predict::twopi / $xnodp / Predict::xmnpda >= .15625) {
+			$this->flags |= Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
+		} else {
+			$this->flags &= ~Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG;
+		}
+	}
+
+	/** Initialise satellite data.
+	 *  @param sat The satellite to initialise.
+	 *  @param qth Optional QTH info, use (0,0) if NULL.
+	 *
+	 * This function calculates the satellite data at t = 0, ie. epoch time
+	 * The function is called automatically by gtk_sat_data_read_sat.
+	 */
+	public function sat_data_init_sat(Predict_Sat $sat, Predict_QTH $qth = null)
+	{
+		$obs_geodetic = new Predict_Geodetic();
+		$obs_set = new Predict_ObsSet();
+		$sat_geodetic = new Predict_Geodetic();
+		/* double jul_utc, age; */
+
+		$jul_utc = Predict_Time::Julian_Date_of_Epoch($sat->tle->epoch); // => tsince = 0.0
+		$sat->jul_epoch = $jul_utc;
+
+		/* initialise observer location */
+		if ($qth != null) {
+			$obs_geodetic->lon = $qth->lon * Predict::de2ra;
+			$obs_geodetic->lat = $qth->lat * Predict::de2ra;
+			$obs_geodetic->alt = $qth->alt / 1000.0;
+			$obs_geodetic->theta = 0;
+		}
+		else {
+			$obs_geodetic->lon = 0.0;
+			$obs_geodetic->lat = 0.0;
+			$obs_geodetic->alt = 0.0;
+			$obs_geodetic->theta = 0;
+		}
+
+		/* execute computations */
+		$sdpsgp = Predict_SGPSDP::getInstance($sat);
+		if ($sat->flags & Predict_SGPSDP::DEEP_SPACE_EPHEM_FLAG) {
+			$sdpsgp->SDP4($sat, 0.0);
+		} else {
+			$sdpsgp->SGP4($sat, 0.0);
+		}
+
+		/* scale position and velocity to km and km/sec */
+		Predict_Math::Convert_Sat_State($sat->pos, $sat->vel);
+
+		/* get the velocity of the satellite */
+		$sat->vel->w = sqrt($sat->vel->x * $sat->vel->x + $sat->vel->y * $sat->vel->y + $sat->vel->z * $sat->vel->z);
+		$sat->velo = $sat->vel->w;
+		Predict_SGPObs::Calculate_Obs($jul_utc, $sat->pos, $sat->vel, $obs_geodetic, $obs_set);
+		Predict_SGPObs::Calculate_LatLonAlt($jul_utc, $sat->pos, $sat_geodetic);
+
+		while ($sat_geodetic->lon < -Predict::pi) {
+			$sat_geodetic->lon += Predict::twopi;
+		}
+
+		while ($sat_geodetic->lon > Predict::pi) {
+			$sat_geodetic->lon -= Predict::twopi;
+		}
+
+		$sat->az = Predict_Math::Degrees($obs_set->az);
+		$sat->el = Predict_Math::Degrees($obs_set->el);
+		$sat->range = $obs_set->range;
+		$sat->range_rate = $obs_set->range_rate;
+		$sat->ssplat = Predict_Math::Degrees($sat_geodetic->lat);
+		$sat->ssplon = Predict_Math::Degrees($sat_geodetic->lon);
+		$sat->alt = $sat_geodetic->alt;
+		$sat->ma = Predict_Math::Degrees($sat->phase);
+		$sat->ma *= 256.0 / 360.0;
+		$sat->footprint = 2.0 * Predict::xkmper * acos (Predict::xkmper/$sat->pos->w);
+		$age = 0.0;
+		$sat->orbit = floor(($sat->tle->xno * Predict::xmnpda / Predict::twopi +
+								   $age * $sat->tle->bstar * Predict::ae) * $age +
+								  $sat->tle->xmo / Predict::twopi) + $sat->tle->revnum - 1;
+
+		/* orbit type */
+		$sat->otype = $sat->get_orbit_type($sat);
+	}
+
+	public function get_orbit_type(Predict_Sat $sat)
+	{
+		 $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
+
+		 if ($this->geostationary($sat)) {
+			  $orbit = Predict_SGPSDP::ORBIT_TYPE_GEO;
+		 } else if ($this->decayed($sat)) {
+			  $orbit = Predict_SGPSDP::ORBIT_TYPE_DECAYED;
+		 } else {
+			  $orbit = Predict_SGPSDP::ORBIT_TYPE_UNKNOWN;
+		 }
+
+		 return $orbit;
+	}
+
+
+	/** Determinte whether satellite is in geostationary orbit.
+	 *  @author John A. Magliacane, KD2BD
+	 *  @param sat Pointer to satellite data.
+	 *  @return TRUE if the satellite appears to be in geostationary orbit,
+	 *          FALSE otherwise.
+	 *
+	 * A satellite is in geostationary orbit if
+	 *
+	 *     fabs (sat.meanmotion - 1.0027) < 0.0002
+	 *
+	 * Note: Appearantly, the mean motion can deviate much more from 1.0027 than 0.0002
+	 */
+	public function geostationary(Predict_Sat $sat)
+	{
+		 if (abs($sat->meanmo - 1.0027) < 0.0002) {
+			  return true;
+		 } else {
+			  return false;
+		}
+	}
+
+
+	/** Determine whether satellite has decayed.
+	 *  @author John A. Magliacane, KD2BD
+	 *  @author Alexandru Csete, OZ9AEC
+	 *  @param sat Pointer to satellite data.
+	 *  @return TRUE if the satellite appears to have decayed, FALSE otherwise.
+	 *  @bug Modified version of the predict code but it is not tested.
+	 *
+	 * A satellite is decayed if
+	 *
+	 *    satepoch + ((16.666666 - sat.meanmo) / (10.0*fabs(sat.drag))) < "now"
+	 *
+	 */
+	public function decayed(Predict_Sat $sat)
+	{
+		/* tle.xndt2o/(twopi/xmnpda/xmnpda) is the value before converted the
            value matches up with the value in predict 2.2.3 */
-        /*** FIXME decayed is treated as a static quantity.
+		/*** FIXME decayed is treated as a static quantity.
              It is time dependent. Also sat->jul_utc is often zero
              when this function is called
         ***/
-        if ((10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda))) == 0) {
-    		return true;
-    	} elseif ($sat->jul_epoch + ((16.666666 - $sat->meanmo) /
-                               (10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda)))) < $sat->jul_utc) {
-              return true;
-        } else {
-              return false;
-        }
-    }
-
-    /**
-     * Experimental attempt at calculating apparent magnitude.  Known intrinsic
-     * magnitudes are listed inside the function for now.
-     *
-     * @param float       $time The daynum the satellite is calculated for
-     * @param Predict_QTH $qth  The observer location
-     *
-     * @return null on failure, float otherwise
-     */
-    public function calculateApparentMagnitude($time, Predict_QTH $qth)
-    {
-        // Recorded intrinsic magnitudes and their respective
-        // illumination and distance from heavens-above.com
-        static $intrinsicMagnitudes = array(
-            '25544' => array(
-                'mag'      => -1.3,
-                'illum'    => .5,
-                'distance' => 1000,
-            )
-        );
-
-        // Return null if we don't have a record of the intrinsic mag
-        if (!isset($intrinsicMagnitudes[$this->tle->catnr])) {
-            return null;
-        }
-        $imag = $intrinsicMagnitudes[$this->tle->catnr];
-
-        // Convert the observer's geodetic info to radians and km so
-        // we can compare vectors
-        $observerGeo      = new Predict_Geodetic();
-        $observerGeo->lat = Predict_Math::Radians($qth->lat);
-        $observerGeo->lon = Predict_Math::Radians($qth->lon);
-        $observerGeo->alt = $qth->alt * 1000;
-
-        // Now determine the sun and observer positions
-        $observerPos      = new Predict_Vector();
-        $observerVel      = new Predict_Vector();
-        $solarVector      = new Predict_Vector();
-        Predict_Solar::Calculate_Solar_Position($time, $solarVector);
-        Predict_SGPObs::Calculate_User_PosVel($time, $observerGeo, $observerPos, $observerVel);
-
-        // Determine the solar phase and and thus the percent illumination
-        $observerSatPos = new Predict_Vector();
-        Predict_Math::Vec_Sub($this->pos, $observerPos, $observerSatPos);
-        $phaseAngle = Predict_Math::Degrees(Predict_Math::Angle($solarVector, $observerSatPos));
-        $illum      = $phaseAngle / 180;
-
-        $illuminationChange            = $illum / $imag['illum'];
-        $inverseSquareOfDistanceChange = pow(($imag['distance'] / $this->range), 2);
-        $changeInMagnitude             = log(
-            $illuminationChange * $inverseSquareOfDistanceChange,
-            self::POGSONS_RATIO
-        );
-
-        return $imag['mag'] - $changeInMagnitude;
-    }
+		if ((10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda))) == 0) {
+			return true;
+		} elseif ($sat->jul_epoch + ((16.666666 - $sat->meanmo) /
+							   (10.0 * abs($sat->tle->xndt2o / (Predict::twopi / Predict::xmnpda / Predict::xmnpda)))) < $sat->jul_utc) {
+			  return true;
+		} else {
+			  return false;
+		}
+	}
+
+	/**
+	 * Experimental attempt at calculating apparent magnitude.  Known intrinsic
+	 * magnitudes are listed inside the function for now.
+	 *
+	 * @param float       $time The daynum the satellite is calculated for
+	 * @param Predict_QTH $qth  The observer location
+	 *
+	 * @return null on failure, float otherwise
+	 */
+	public function calculateApparentMagnitude($time, Predict_QTH $qth)
+	{
+		// Recorded intrinsic magnitudes and their respective
+		// illumination and distance from heavens-above.com
+		static $intrinsicMagnitudes = array(
+			'25544' => array(
+				'mag'      => -1.3,
+				'illum'    => .5,
+				'distance' => 1000,
+			)
+		);
+
+		// Return null if we don't have a record of the intrinsic mag
+		if (!isset($intrinsicMagnitudes[$this->tle->catnr])) {
+			return null;
+		}
+		$imag = $intrinsicMagnitudes[$this->tle->catnr];
+
+		// Convert the observer's geodetic info to radians and km so
+		// we can compare vectors
+		$observerGeo      = new Predict_Geodetic();
+		$observerGeo->lat = Predict_Math::Radians($qth->lat);
+		$observerGeo->lon = Predict_Math::Radians($qth->lon);
+		$observerGeo->alt = $qth->alt * 1000;
+
+		// Now determine the sun and observer positions
+		$observerPos      = new Predict_Vector();
+		$observerVel      = new Predict_Vector();
+		$solarVector      = new Predict_Vector();
+		Predict_Solar::Calculate_Solar_Position($time, $solarVector);
+		Predict_SGPObs::Calculate_User_PosVel($time, $observerGeo, $observerPos, $observerVel);
+
+		// Determine the solar phase and and thus the percent illumination
+		$observerSatPos = new Predict_Vector();
+		Predict_Math::Vec_Sub($this->pos, $observerPos, $observerSatPos);
+		$phaseAngle = Predict_Math::Degrees(Predict_Math::Angle($solarVector, $observerSatPos));
+		$illum      = $phaseAngle / 180;
+
+		$illuminationChange            = $illum / $imag['illum'];
+		$inverseSquareOfDistanceChange = pow(($imag['distance'] / $this->range), 2);
+		$changeInMagnitude             = log(
+			$illuminationChange * $inverseSquareOfDistanceChange,
+			self::POGSONS_RATIO
+		);
+
+		return $imag['mag'] - $changeInMagnitude;
+	}
 }

require/libs/Predict/Predict/Vector.php

Indentation +4 added lines, -4 removed lines

@@ -6,8 +6,8 @@
  */
 class Predict_Vector
 {
-    public $x = 0;
-    public $y = 0;
-    public $z = 0;
-    public $w = 0;
+	public $x = 0;
+	public $y = 0;
+	public $z = 0;
+	public $w = 0;
 }

require/libs/Predict/Predict/Time.php

Indentation +188 added lines, -188 removed lines

@@ -30,194 +30,194 @@
  */
 class Predict_Time
 {
-    /* The function Julian_Date_of_Epoch returns the Julian Date of     */
-    /* an epoch specified in the format used in the NORAD two-line      */
-    /* element sets. It has been modified to support dates beyond       */
-    /* the year 1999 assuming that two-digit years in the range 00-56   */
-    /* correspond to 2000-2056. Until the two-line element set format   */
-    /* is changed, it is only valid for dates through 2056 December 31. */
-    public static function Julian_Date_of_Epoch($epoch)
-    {
-        $year = 0;
-
-        /* Modification to support Y2K */
-        /* Valid 1957 through 2056     */
-        $day = self::modf($epoch * 1E-3, $year) * 1E3;
-        if ($year < 57) {
-            $year = $year + 2000;
-        } else {
-            $year = $year + 1900;
-        }
-        /* End modification */
-
-        return self::Julian_Date_of_Year($year) + $day;
-    }
-
-    /* Equivalent to the C modf function */
-    public static function modf($x, &$ipart) {
-        $ipart = (int)$x;
-        return $x - $ipart;
-    }
-
-    /* The function Julian_Date_of_Year calculates the Julian Date  */
-    /* of Day 0.0 of {year}. This function is used to calculate the */
-    /* Julian Date of any date by using Julian_Date_of_Year, DOY,   */
-    /* and Fraction_of_Day. */
-    public static function Julian_Date_of_Year($year)
-    {
-        /* Astronomical Formulae for Calculators, Jean Meeus, */
-        /* pages 23-25. Calculate Julian Date of 0.0 Jan year */
-        $year = $year - 1;
-        $i = (int) ($year / 100);
-        $A = $i;
-        $i = (int) ($A / 4);
-        $B = (int) (2 - $A + $i);
-        $i = (int) (365.25 * $year);
-        $i += (int) (30.6001 * 14);
-        $jdoy = $i + 1720994.5 + $B;
-
-        return $jdoy;
-    }
-
-    /* The function ThetaG calculates the Greenwich Mean Sidereal Time */
-    /* for an epoch specified in the format used in the NORAD two-line */
-    /* element sets. It has now been adapted for dates beyond the year */
-    /* 1999, as described above. The function ThetaG_JD provides the   */
-    /* same calculation except that it is based on an input in the     */
-    /* form of a Julian Date. */
-    public static function ThetaG($epoch, Predict_DeepArg $deep_arg)
-    {
-        /* Reference:  The 1992 Astronomical Almanac, page B6. */
-        // double year,day,UT,jd,TU,GMST,_ThetaG;
-
-        /* Modification to support Y2K */
-        /* Valid 1957 through 2056     */
-        $year = 0;
-        $day = self::modf($epoch * 1E-3, $year) * 1E3;
-
-        if ($year < 57) {
-            $year += 2000;
-        } else {
-            $year += 1900;
-        }
-        /* End modification */
-
-        $UT = fmod($day, $day);
-        $jd = self::Julian_Date_of_Year($year) + $day;
-        $TU = ($jd - 2451545.0) / 36525;
-        $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6));
-        $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday);
-        $deep_arg->ds50 = $jd - 2433281.5 + $UT;
-
-        return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494);
-    }
-
-    /* See the ThetaG doc block above */
-    public static function ThetaG_JD($jd)
-    {
-        /* Reference:  The 1992 Astronomical Almanac, page B6. */
-        $UT   = Predict_Math::Frac($jd + 0.5);
-        $jd   = $jd - $UT;
-        $TU   = ($jd - 2451545.0) / 36525;
-        $GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6));
-        $GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday);
-
-        return Predict::twopi * $GMST / Predict::secday;
-    }
-
-    /**
-     * Read the system clock and return the current Julian day.  From phpPredict
-     *
-     * @return float
-     */
-    public static function get_current_daynum() {
-        // Gets the current decimal day number from microtime
-
-        list($usec, $sec) = explode(' ', microtime());
-        return self::unix2daynum($sec, $usec);
-    }
-
-    /**
-     * Converts a standard unix timestamp and optional
-     * milliseconds to a daynum
-     *
-     * @param int $sec  Seconds from the unix epoch
-     * @param int $usec Optional milliseconds
-     *
-     * @return float
-     */
-    public static function unix2daynum($sec, $usec = 0)
-    {
-        $time = ((($sec + $usec) / 86400.0) - 3651.0);
-        return $time + 2444238.5;
-    }
-
-    /* The function Delta_ET has been added to allow calculations on   */
-    /* the position of the sun.  It provides the difference between UT */
-    /* (approximately the same as UTC) and ET (now referred to as TDT).*/
-    /* This function is based on a least squares fit of data from 1950 */
-    /* to 1991 and will need to be updated periodically. */
-    public static function Delta_ET($year)
-    {
-      /* Values determined using data from 1950-1991 in the 1990
+	/* The function Julian_Date_of_Epoch returns the Julian Date of     */
+	/* an epoch specified in the format used in the NORAD two-line      */
+	/* element sets. It has been modified to support dates beyond       */
+	/* the year 1999 assuming that two-digit years in the range 00-56   */
+	/* correspond to 2000-2056. Until the two-line element set format   */
+	/* is changed, it is only valid for dates through 2056 December 31. */
+	public static function Julian_Date_of_Epoch($epoch)
+	{
+		$year = 0;
+
+		/* Modification to support Y2K */
+		/* Valid 1957 through 2056     */
+		$day = self::modf($epoch * 1E-3, $year) * 1E3;
+		if ($year < 57) {
+			$year = $year + 2000;
+		} else {
+			$year = $year + 1900;
+		}
+		/* End modification */
+
+		return self::Julian_Date_of_Year($year) + $day;
+	}
+
+	/* Equivalent to the C modf function */
+	public static function modf($x, &$ipart) {
+		$ipart = (int)$x;
+		return $x - $ipart;
+	}
+
+	/* The function Julian_Date_of_Year calculates the Julian Date  */
+	/* of Day 0.0 of {year}. This function is used to calculate the */
+	/* Julian Date of any date by using Julian_Date_of_Year, DOY,   */
+	/* and Fraction_of_Day. */
+	public static function Julian_Date_of_Year($year)
+	{
+		/* Astronomical Formulae for Calculators, Jean Meeus, */
+		/* pages 23-25. Calculate Julian Date of 0.0 Jan year */
+		$year = $year - 1;
+		$i = (int) ($year / 100);
+		$A = $i;
+		$i = (int) ($A / 4);
+		$B = (int) (2 - $A + $i);
+		$i = (int) (365.25 * $year);
+		$i += (int) (30.6001 * 14);
+		$jdoy = $i + 1720994.5 + $B;
+
+		return $jdoy;
+	}
+
+	/* The function ThetaG calculates the Greenwich Mean Sidereal Time */
+	/* for an epoch specified in the format used in the NORAD two-line */
+	/* element sets. It has now been adapted for dates beyond the year */
+	/* 1999, as described above. The function ThetaG_JD provides the   */
+	/* same calculation except that it is based on an input in the     */
+	/* form of a Julian Date. */
+	public static function ThetaG($epoch, Predict_DeepArg $deep_arg)
+	{
+		/* Reference:  The 1992 Astronomical Almanac, page B6. */
+		// double year,day,UT,jd,TU,GMST,_ThetaG;
+
+		/* Modification to support Y2K */
+		/* Valid 1957 through 2056     */
+		$year = 0;
+		$day = self::modf($epoch * 1E-3, $year) * 1E3;
+
+		if ($year < 57) {
+			$year += 2000;
+		} else {
+			$year += 1900;
+		}
+		/* End modification */
+
+		$UT = fmod($day, $day);
+		$jd = self::Julian_Date_of_Year($year) + $day;
+		$TU = ($jd - 2451545.0) / 36525;
+		$GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6));
+		$GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday);
+		$deep_arg->ds50 = $jd - 2433281.5 + $UT;
+
+		return Predict_Math::FMod2p(6.3003880987 * $deep_arg->ds50 + 1.72944494);
+	}
+
+	/* See the ThetaG doc block above */
+	public static function ThetaG_JD($jd)
+	{
+		/* Reference:  The 1992 Astronomical Almanac, page B6. */
+		$UT   = Predict_Math::Frac($jd + 0.5);
+		$jd   = $jd - $UT;
+		$TU   = ($jd - 2451545.0) / 36525;
+		$GMST = 24110.54841 + $TU * (8640184.812866 + $TU * (0.093104 - $TU * 6.2E-6));
+		$GMST = Predict_Math::Modulus($GMST + Predict::secday * Predict::omega_E * $UT, Predict::secday);
+
+		return Predict::twopi * $GMST / Predict::secday;
+	}
+
+	/**
+	 * Read the system clock and return the current Julian day.  From phpPredict
+	 *
+	 * @return float
+	 */
+	public static function get_current_daynum() {
+		// Gets the current decimal day number from microtime
+
+		list($usec, $sec) = explode(' ', microtime());
+		return self::unix2daynum($sec, $usec);
+	}
+
+	/**
+	 * Converts a standard unix timestamp and optional
+	 * milliseconds to a daynum
+	 *
+	 * @param int $sec  Seconds from the unix epoch
+	 * @param int $usec Optional milliseconds
+	 *
+	 * @return float
+	 */
+	public static function unix2daynum($sec, $usec = 0)
+	{
+		$time = ((($sec + $usec) / 86400.0) - 3651.0);
+		return $time + 2444238.5;
+	}
+
+	/* The function Delta_ET has been added to allow calculations on   */
+	/* the position of the sun.  It provides the difference between UT */
+	/* (approximately the same as UTC) and ET (now referred to as TDT).*/
+	/* This function is based on a least squares fit of data from 1950 */
+	/* to 1991 and will need to be updated periodically. */
+	public static function Delta_ET($year)
+	{
+	  /* Values determined using data from 1950-1991 in the 1990
          Astronomical Almanac.  See DELTA_ET.WQ1 for details. */
 
-      $delta_et = 26.465 + 0.747622 * ($year - 1950) +
-                 1.886913 * sin(Predict::twopi * ($year - 1975) / 33);
-
-      return $delta_et;
-    }
-
-    /**
-     * Converts a daynum to a unix timestamp.  From phpPredict.
-     *
-     * @param float $dn Julian Daynum
-     *
-     * @return float
-     */
-    public static function daynum2unix($dn) {
-        // Converts a daynum to a UNIX timestamp
-
-        return (86400.0 * ($dn - 2444238.5 + 3651.0));
-    }
-
-    /**
-     * Converts a daynum to a readable time format.
-     *
-     * @param float $dn The julian date
-     * @param string $zone The zone string, defaults to America/Los_Angeles
-     * @param string $format The date() function's format string.  Defaults to m-d-Y H:i:s
-     *
-     * @return string
-     */
-    public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s')
-    {
-        $unix = self::daynum2unix($dn);
-        $date = new DateTime("@" . round($unix));
-        $dateTimezone = new DateTimezone($zone);
-        $date->setTimezone($dateTimezone);
-        return $date->format($format);
-    }
-
-    /**
-     * Returns the unix timestamp of a TLE's epoch
-     *
-     * @param Predict_TLE $tle The TLE object
-     *
-     * @return int
-     */
-    public static function getEpochTimeStamp(Predict_TLE $tle)
-    {
-        $year = $tle->epoch_year;
-        $day  = $tle->epoch_day;
-        $sec  = round(86400 * $tle->epoch_fod);
-
-        $zone = new DateTimeZone('GMT');
-        $date = new DateTime();
-        $date->setTimezone($zone);
-        $date->setDate($year, 1, 1);
-        $date->setTime(0, 0, 0);
-
-        return $date->format('U') + (86400 * $day) + $sec - 86400;
-    }
+	  $delta_et = 26.465 + 0.747622 * ($year - 1950) +
+				 1.886913 * sin(Predict::twopi * ($year - 1975) / 33);
+
+	  return $delta_et;
+	}
+
+	/**
+	 * Converts a daynum to a unix timestamp.  From phpPredict.
+	 *
+	 * @param float $dn Julian Daynum
+	 *
+	 * @return float
+	 */
+	public static function daynum2unix($dn) {
+		// Converts a daynum to a UNIX timestamp
+
+		return (86400.0 * ($dn - 2444238.5 + 3651.0));
+	}
+
+	/**
+	 * Converts a daynum to a readable time format.
+	 *
+	 * @param float $dn The julian date
+	 * @param string $zone The zone string, defaults to America/Los_Angeles
+	 * @param string $format The date() function's format string.  Defaults to m-d-Y H:i:s
+	 *
+	 * @return string
+	 */
+	public static function daynum2readable($dn, $zone = 'America/Los_Angeles', $format = 'm-d-Y H:i:s')
+	{
+		$unix = self::daynum2unix($dn);
+		$date = new DateTime("@" . round($unix));
+		$dateTimezone = new DateTimezone($zone);
+		$date->setTimezone($dateTimezone);
+		return $date->format($format);
+	}
+
+	/**
+	 * Returns the unix timestamp of a TLE's epoch
+	 *
+	 * @param Predict_TLE $tle The TLE object
+	 *
+	 * @return int
+	 */
+	public static function getEpochTimeStamp(Predict_TLE $tle)
+	{
+		$year = $tle->epoch_year;
+		$day  = $tle->epoch_day;
+		$sec  = round(86400 * $tle->epoch_fod);
+
+		$zone = new DateTimeZone('GMT');
+		$date = new DateTime();
+		$date->setTimezone($zone);
+		$date->setDate($year, 1, 1);
+		$date->setTime(0, 0, 0);
+
+		return $date->format('U') + (86400 * $day) + $sec - 86400;
+	}
 }