ootk 1.0.7 → 1.5.0

package/lib/{sgp4.js → sgp4/sgp4.js}

@@ -9,29 +9,35 @@
  * in order provide context to the functions and clarify any adjustments made for
  * TypeScript/JavaScript compatibility.
  *
- * @copyright MIT License.
- * Copyright (c) 2020-2021 Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2022 Theodore Kruczek
  *
- * Permission is hereby granted, free of charge, to any person obtaining a copy
- * of this software and associated documentation files (the "Software"), to deal
- * in the Software without restriction, including without limitation the rights
- * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
- * copies of the Software, and to permit persons to whom the Software is
- * furnished to do so, subject to the following conditions:
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
  *
- * The above copyright notice and this permission notice shall be included in all
- * copies or substantial portions of the Software.
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
  *
- * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
- * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
- * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
- * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
- * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
- * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
- * SOFTWARE.
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
-import { DEG2RAD, PI, TAU, temp4, x2o3 } from './utils/constants';
-/*     ----------------------------------------------------------------
+// NOTE: This file is meant to maintain as much of the original format as possible.
+/* eslint-disable no-lonely-if */
+/* eslint-disable prefer-destructuring */
+/* eslint-disable complexity */
+/* eslint-disable max-len */
+/* eslint-disable no-param-reassign */
+/* eslint-disable max-params */
+/* eslint-disable no-shadow */
+/* eslint-disable init-declarations */
+/* eslint-disable max-statements */
+/* eslint-disable max-lines-per-function */
+/* eslint-disable max-lines */
+import { DEG2RAD, PI, TAU, temp4, x2o3 } from '../utils/constants';
+/*
+ *     ----------------------------------------------------------------
  *
  *                               sgp4unit.cpp
  *
@@ -86,76 +92,79 @@ import { DEG2RAD, PI, TAU, temp4, x2o3 } from './utils/constants';
  *                           2nd edition baseline
  *                     80  norad
  *                           original baseline
- *       ----------------------------------------------------------------      */
+ *       ----------------------------------------------------------------
+ */
 var Sgp4 = /** @class */ (function () {
     function Sgp4() {
     }
-    /* -----------------------------------------------------------------------------
-   *
-   *                           procedure dpper
-   *
-   *  this procedure provides deep space long period periodic contributions
-   *    to the mean elements.  by design, these periodics are zero at epoch.
-   *    this used to be dscom which included initialization, but it's really a
-   *    recurring function.
-   *
-   *  author        : david vallado                  719-573-2600   28 jun 2005
-   *
-   *  inputs        :
-   *    e3          -
-   *    ee2         -
-   *    peo         -
-   *    pgho        -
-   *    pho         -
-   *    PInco       -
-   *    plo         -
-   *    se2 , se3 , sgh2, sgh3, sgh4, sh2, sh3, si2, si3, sl2, sl3, sl4 -
-   *    t           -
-   *    xh2, xh3, xi2, xi3, xl2, xl3, xl4 -
-   *    zmol        -
-   *    zmos        -
-   *    ep          - eccentricity                           0.0 - 1.0
-   *    inclo       - inclination - needed for lyddane modification
-   *    nodep       - right ascension of ascending node
-   *    argpp       - argument of perigee
-   *    mp          - mean anomaly
-   *
-   *  outputs       :
-   *    ep          - eccentricity                           0.0 - 1.0
-   *    inclp       - inclination
-   *    nodep        - right ascension of ascending node
-   *    argpp       - argument of perigee
-   *    mp          - mean anomaly
-   *
-   *  locals        :
-   *    alfdp       -
-   *    betdp       -
-   *    cosip  , sinip  , cosop  , sinop  ,
-   *    dalf        -
-   *    dbet        -
-   *    dls         -
-   *    f2, f3      -
-   *    pe          -
-   *    pgh         -
-   *    ph          -
-   *    PInc        -
-   *    pl          -
-   *    sel   , ses   , sghl  , sghs  , shl   , shs   , sil   , sinzf , sis   ,
-   *    sll   , sls
-   *    xls         -
-   *    xnoh        -
-   *    zf          -
-   *    zm          -
-   *
-   *  coupling      :
-   *    none.
-   *
-   *  references    :
-   *    hoots, roehrich, norad spacetrack report #3 1980
-   *    hoots, norad spacetrack report #6 1986
-   *    hoots, schumacher and glover 2004
-   *    vallado, crawford, hujsak, kelso  2006
-   ----------------------------------------------------------------------------*/
+    /*
+     * -----------------------------------------------------------------------------
+     *
+     *                           procedure dpper
+     *
+     *  this procedure provides deep space long period periodic contributions
+     *    to the mean elements.  by design, these periodics are zero at epoch.
+     *    this used to be dscom which included initialization, but it's really a
+     *    recurring function.
+     *
+     *  author        : david vallado                  719-573-2600   28 jun 2005
+     *
+     *  inputs        :
+     *    e3          -
+     *    ee2         -
+     *    peo         -
+     *    pgho        -
+     *    pho         -
+     *    PInco       -
+     *    plo         -
+     *    se2 , se3 , sgh2, sgh3, sgh4, sh2, sh3, si2, si3, sl2, sl3, sl4 -
+     *    t           -
+     *    xh2, xh3, xi2, xi3, xl2, xl3, xl4 -
+     *    zmol        -
+     *    zmos        -
+     *    ep          - eccentricity                           0.0 - 1.0
+     *    inclo       - inclination - needed for lyddane modification
+     *    nodep       - right ascension of ascending node
+     *    argpp       - argument of perigee
+     *    mp          - mean anomaly
+     *
+     *  outputs       :
+     *    ep          - eccentricity                           0.0 - 1.0
+     *    inclp       - inclination
+     *    nodep        - right ascension of ascending node
+     *    argpp       - argument of perigee
+     *    mp          - mean anomaly
+     *
+     *  locals        :
+     *    alfdp       -
+     *    betdp       -
+     *    cosip  , sinip  , cosop  , sinop  ,
+     *    dalf        -
+     *    dbet        -
+     *    dls         -
+     *    f2, f3      -
+     *    pe          -
+     *    pgh         -
+     *    ph          -
+     *    PInc        -
+     *    pl          -
+     *    sel   , ses   , sghl  , sghs  , shl   , shs   , sil   , sinzf , sis   ,
+     *    sll   , sls
+     *    xls         -
+     *    xnoh        -
+     *    zf          -
+     *    zm          -
+     *
+     *  coupling      :
+     *    none.
+     *
+     *  references    :
+     *    hoots, roehrich, norad spacetrack report #3 1980
+     *    hoots, norad spacetrack report #6 1986
+     *    hoots, schumacher and glover 2004
+     *    vallado, crawford, hujsak, kelso  2006
+     * ----------------------------------------------------------------------------
+     */
     Sgp4.dpper = function (satrec, options) {
         var e3 = satrec.e3, ee2 = satrec.ee2, peo = satrec.peo, pgho = satrec.pgho, pho = satrec.pho, PInco = satrec.PInco, plo = satrec.plo, se2 = satrec.se2, se3 = satrec.se3, sgh2 = satrec.sgh2, sgh3 = satrec.sgh3, sgh4 = satrec.sgh4, sh2 = satrec.sh2, sh3 = satrec.sh3, si2 = satrec.si2, si3 = satrec.si3, sl2 = satrec.sl2, sl3 = satrec.sl3, sl4 = satrec.sl4, t = satrec.t, xgh2 = satrec.xgh2, xgh3 = satrec.xgh3, xgh4 = satrec.xgh4, xh2 = satrec.xh2, xh3 = satrec.xh3, xi2 = satrec.xi2, xi3 = satrec.xi3, xl2 = satrec.xl2, xl3 = satrec.xl3, xl4 = satrec.xl4, zmol = satrec.zmol, zmos = satrec.zmos;
         var ep = options.ep, inclp = options.inclp, nodep = options.nodep, argpp = options.argpp, mp = options.mp;
@@ -168,7 +177,7 @@ var Sgp4 = /** @class */ (function () {
         var zel = 0.0549;
         //  --------------- calculate time varying periodics -----------
         var zm = zmos + zns * t;
-        // be sure that the initial call has time set to zero
+        // Be sure that the initial call has time set to zero
         if (init === 'y') {
             zm = zmos;
         }
@@ -210,14 +219,16 @@ var Sgp4 = /** @class */ (function () {
             var sinip = Math.sin(inclp);
             var cosip = Math.cos(inclp);
             /* ----------------- apply periodics directly ------------ */
-            // sgp4fix for lyddane choice
-            // strn3 used original inclination - this is technically feasible
-            // gsfc used perturbed inclination - also technically feasible
-            // probably best to readjust the 0.2 limit value and limit discontinuity
-            // 0.2 rad = 11.45916 deg
-            // use next line for original strn3 approach and original inclination
-            // if (inclo >= 0.2)
-            // use next line for gsfc version and perturbed inclination
+            /*
+             * Sgp4fix for lyddane choice
+             * strn3 used original inclination - this is technically feasible
+             * gsfc used perturbed inclination - also technically feasible
+             * probably best to readjust the 0.2 limit value and limit discontinuity
+             * 0.2 rad = 11.45916 deg
+             * use next line for original strn3 approach and original inclination
+             * if (inclo >= 0.2)
+             * use next line for gsfc version and perturbed inclination
+             */
             if (inclp >= 0.2) {
                 ph /= sinip;
                 pgh -= cosip * ph;
@@ -236,8 +247,10 @@ var Sgp4 = /** @class */ (function () {
                 alfdp += dalf;
                 betdp += dbet;
                 nodep %= TAU;
-                //  sgp4fix for afspc written intrinsic functions
-                //  nodep used without a trigonometric function ahead
+                /*
+                 *  Sgp4fix for afspc written intrinsic functions
+                 *  nodep used without a trigonometric function ahead
+                 */
                 /* istanbul ignore next */
                 if (nodep < 0.0 && opsmode === 'a') {
                     nodep += TAU;
@@ -247,8 +260,10 @@ var Sgp4 = /** @class */ (function () {
                 xls += dls;
                 var xnoh = nodep;
                 nodep = Math.atan2(alfdp, betdp);
-                //  sgp4fix for afspc written intrinsic functions
-                //  nodep used without a trigonometric function ahead
+                /*
+                 *  Sgp4fix for afspc written intrinsic functions
+                 *  nodep used without a trigonometric function ahead
+                 */
                 /* istanbul ignore next */
                 if (nodep < 0.0 && opsmode === 'a') {
                     nodep += TAU;
@@ -265,7 +280,7 @@ var Sgp4 = /** @class */ (function () {
                 mp += pl;
                 argpp = xls - mp - cosip * nodep;
             }
-        } // if init == 'n'
+        } // If init == 'n'
         return {
             ep: ep,
             inclp: inclp,
@@ -273,74 +288,76 @@ var Sgp4 = /** @class */ (function () {
             argpp: argpp,
             mp: mp,
         };
-    }; // dpper
-    /*-----------------------------------------------------------------------------
-   *
-   *                           procedure dscom
-   *
-   *  this procedure provides deep space common items used by both the secular
-   *    and periodics subroutines.  input is provided as shown. this routine
-   *    used to be called dpper, but the functions inside weren't well organized.
-   *
-   *  author        : david vallado                  719-573-2600   28 jun 2005
-   *
-   *  inputs        :
-   *    epoch       -
-   *    ep          - eccentricity
-   *    argpp       - argument of perigee
-   *    tc          -
-   *    inclp       - inclination
-   *    nodep       - right ascension of ascending node
-   *    np          - mean motion
-   *
-   *  outputs       :
-   *    sinim  , cosim  , sinomm , cosomm , snodm  , cnodm
-   *    day         -
-   *    e3          -
-   *    ee2         -
-   *    em          - eccentricity
-   *    emsq        - eccentricity squared
-   *    gam         -
-   *    peo         -
-   *    pgho        -
-   *    pho         -
-   *    PInco       -
-   *    plo         -
-   *    rtemsq      -
-   *    se2, se3         -
-   *    sgh2, sgh3, sgh4        -
-   *    sh2, sh3, si2, si3, sl2, sl3, sl4         -
-   *    s1, s2, s3, s4, s5, s6, s7          -
-   *    ss1, ss2, ss3, ss4, ss5, ss6, ss7, sz1, sz2, sz3         -
-   *    sz11, sz12, sz13, sz21, sz22, sz23, sz31, sz32, sz33        -
-   *    xgh2, xgh3, xgh4, xh2, xh3, xi2, xi3, xl2, xl3, xl4         -
-   *    nm          - mean motion
-   *    z1, z2, z3, z11, z12, z13, z21, z22, z23, z31, z32, z33         -
-   *    zmol        -
-   *    zmos        -
-   *
-   *  locals        :
-   *    a1, a2, a3, a4, a5, a6, a7, a8, a9, a10         -
-   *    betasq      -
-   *    cc          -
-   *    ctem, stem        -
-   *    x1, x2, x3, x4, x5, x6, x7, x8          -
-   *    xnodce      -
-   *    xnoi        -
-   *    zcosg  , zsing  , zcosgl , zsingl , zcosh  , zsinh  , zcoshl , zsinhl ,
-   *    zcosi  , zsini  , zcosil , zsinil ,
-   *    zx          -
-   *    zy          -
-   *
-   *  coupling      :
-   *    none.
-   *
-   *  references    :
-   *    hoots, roehrich, norad spacetrack report #3 1980
-   *    hoots, norad spacetrack report #6 1986
-   *    hoots, schumacher and glover 2004
-   *    vallado, crawford, hujsak, kelso  2006
-   ----------------------------------------------------------------------------*/
+    }; // Dpper
+    /*
+     *-----------------------------------------------------------------------------
+     *
+     *                           procedure dscom
+     *
+     *  this procedure provides deep space common items used by both the secular
+     *    and periodics subroutines.  input is provided as shown. this routine
+     *    used to be called dpper, but the functions inside weren't well organized.
+     *
+     *  author        : david vallado                  719-573-2600   28 jun 2005
+     *
+     *  inputs        :
+     *    epoch       -
+     *    ep          - eccentricity
+     *    argpp       - argument of perigee
+     *    tc          -
+     *    inclp       - inclination
+     *    nodep       - right ascension of ascending node
+     *    np          - mean motion
+     *
+     *  outputs       :
+     *    sinim  , cosim  , sinomm , cosomm , snodm  , cnodm
+     *    day         -
+     *    e3          -
+     *    ee2         -
+     *    em          - eccentricity
+     *    emsq        - eccentricity squared
+     *    gam         -
+     *    peo         -
+     *    pgho        -
+     *    pho         -
+     *    PInco       -
+     *    plo         -
+     *    rtemsq      -
+     *    se2, se3         -
+     *    sgh2, sgh3, sgh4        -
+     *    sh2, sh3, si2, si3, sl2, sl3, sl4         -
+     *    s1, s2, s3, s4, s5, s6, s7          -
+     *    ss1, ss2, ss3, ss4, ss5, ss6, ss7, sz1, sz2, sz3         -
+     *    sz11, sz12, sz13, sz21, sz22, sz23, sz31, sz32, sz33        -
+     *    xgh2, xgh3, xgh4, xh2, xh3, xi2, xi3, xl2, xl3, xl4         -
+     *    nm          - mean motion
+     *    z1, z2, z3, z11, z12, z13, z21, z22, z23, z31, z32, z33         -
+     *    zmol        -
+     *    zmos        -
+     *
+     *  locals        :
+     *    a1, a2, a3, a4, a5, a6, a7, a8, a9, a10         -
+     *    betasq      -
+     *    cc          -
+     *    ctem, stem        -
+     *    x1, x2, x3, x4, x5, x6, x7, x8          -
+     *    xnodce      -
+     *    xnoi        -
+     *    zcosg  , zsing  , zcosgl , zsingl , zcosh  , zsinh  , zcoshl , zsinhl ,
+     *    zcosi  , zsini  , zcosil , zsinil ,
+     *    zx          -
+     *    zy          -
+     *
+     *  coupling      :
+     *    none.
+     *
+     *  references    :
+     *    hoots, roehrich, norad spacetrack report #3 1980
+     *    hoots, norad spacetrack report #6 1986
+     *    hoots, schumacher and glover 2004
+     *    vallado, crawford, hujsak, kelso  2006
+     *----------------------------------------------------------------------------
+     */
     Sgp4.dscom = function (options) {
         var epoch = options.epoch, ep = options.ep, argpp = options.argpp, tc = options.tc, inclp = options.inclp, nodep = options.nodep, np = options.np;
         // -------------------------- constants -------------------------
@@ -355,7 +372,7 @@ var Sgp4 = /** @class */ (function () {
         var zsings = -0.98088458;
         var TAU = 2.0 * Math.PI;
         //  --------------------- local variables ------------------------
-        var s1, s2, s3, s4, s5, s6, s7, ss1, ss2, ss3, ss4, ss5, ss6, ss7, sz1, sz2, sz3, sz11, sz12, sz13, sz21, sz22, sz23, sz31, sz32, sz33, z31, z32, z33, z1, z2, z3, z11, z12, z13, z21, z22, z23;
+        var s1, s2, s3, s4, s5, s6, s7, ss1, ss2, ss3, ss4, ss5, ss6, ss7, sz1, sz11, sz12, sz13, sz2, sz21, sz22, sz23, sz3, sz31, sz32, sz33, z1, z11, z12, z13, z2, z21, z22, z23, z3, z31, z32, z33;
         var nm = np;
         var em = ep;
         var snodm = Math.sin(nodep);
@@ -579,91 +596,94 @@ var Sgp4 = /** @class */ (function () {
             zmol: zmol,
             zmos: zmos,
         };
-    }; // dscom
-    /*-----------------------------------------------------------------------------
-   *
-   *                           procedure dsinit
-   *
-   *  this procedure provides deep space contributions to mean motion dot due
-   *    to geopotential resonance with half day and one day orbits.
-   *
-   *  author        : david vallado                  719-573-2600   28 jun 2005
-   *
-   *  inputs        :
-   *    cosim, sinim-
-   *    emsq        - eccentricity squared
-   *    argpo       - argument of perigee
-   *    s1, s2, s3, s4, s5      -
-   *    ss1, ss2, ss3, ss4, ss5 -
-   *    sz1, sz3, sz11, sz13, sz21, sz23, sz31, sz33 -
-   *    t           - time
-   *    tc          -
-   *    gsto        - greenwich sidereal time                   rad
-   *    mo          - mean anomaly
-   *    mdot        - mean anomaly dot (rate)
-   *    no          - mean motion
-   *    nodeo       - right ascension of ascending node
-   *    nodedot     - right ascension of ascending node dot (rate)
-   *    xPIdot      -
-   *    z1, z3, z11, z13, z21, z23, z31, z33 -
-   *    eccm        - eccentricity
-   *    argpm       - argument of perigee
-   *    inclm       - inclination
-   *    mm          - mean anomaly
-   *    xn          - mean motion
-   *    nodem       - right ascension of ascending node
-   *
-   *  outputs       :
-   *    em          - eccentricity
-   *    argpm       - argument of perigee
-   *    inclm       - inclination
-   *    mm          - mean anomaly
-   *    nm          - mean motion
-   *    nodem       - right ascension of ascending node
-   *    irez        - flag for resonance           0-none, 1-one day, 2-half day
-   *    atime       -
-   *    d2201, d2211, d3210, d3222, d4410, d4422, d5220, d5232, d5421, d5433    -
-   *    dedt        -
-   *    didt        -
-   *    dmdt        -
-   *    dndt        -
-   *    dnodt       -
-   *    domdt       -
-   *    del1, del2, del3        -
-   *    ses  , sghl , sghs , sgs  , shl  , shs  , sis  , sls
-   *    theta       -
-   *    xfact       -
-   *    xlamo       -
-   *    xli         -
-   *    xni
-   *
-   *  locals        :
-   *    ainv2       -
-   *    aonv        -
-   *    cosisq      -
-   *    eoc         -
-   *    f220, f221, f311, f321, f322, f330, f441, f442, f522, f523, f542, f543  -
-   *    g200, g201, g211, g300, g310, g322, g410, g422, g520, g521, g532, g533  -
-   *    sini2       -
-   *    temp        -
-   *    temp1       -
-   *    theta       -
-   *    xno2        -
-   *
-   *  coupling      :
-   *    getgravconst
-   *
-   *  references    :
-   *    hoots, roehrich, norad spacetrack report #3 1980
-   *    hoots, norad spacetrack report #6 1986
-   *    hoots, schumacher and glover 2004
-   *    vallado, crawford, hujsak, kelso  2006
-   ----------------------------------------------------------------------------*/
+    }; // Dscom
+    /*
+     *-----------------------------------------------------------------------------
+     *
+     *                           procedure dsinit
+     *
+     *  this procedure provides deep space contributions to mean motion dot due
+     *    to geopotential resonance with half day and one day orbits.
+     *
+     *  author        : david vallado                  719-573-2600   28 jun 2005
+     *
+     *  inputs        :
+     *    cosim, sinim-
+     *    emsq        - eccentricity squared
+     *    argpo       - argument of perigee
+     *    s1, s2, s3, s4, s5      -
+     *    ss1, ss2, ss3, ss4, ss5 -
+     *    sz1, sz3, sz11, sz13, sz21, sz23, sz31, sz33 -
+     *    t           - time
+     *    tc          -
+     *    gsto        - greenwich sidereal time                   rad
+     *    mo          - mean anomaly
+     *    mdot        - mean anomaly dot (rate)
+     *    no          - mean motion
+     *    nodeo       - right ascension of ascending node
+     *    nodedot     - right ascension of ascending node dot (rate)
+     *    xPIdot      -
+     *    z1, z3, z11, z13, z21, z23, z31, z33 -
+     *    eccm        - eccentricity
+     *    argpm       - argument of perigee
+     *    inclm       - inclination
+     *    mm          - mean anomaly
+     *    xn          - mean motion
+     *    nodem       - right ascension of ascending node
+     *
+     *  outputs       :
+     *    em          - eccentricity
+     *    argpm       - argument of perigee
+     *    inclm       - inclination
+     *    mm          - mean anomaly
+     *    nm          - mean motion
+     *    nodem       - right ascension of ascending node
+     *    irez        - flag for resonance           0-none, 1-one day, 2-half day
+     *    atime       -
+     *    d2201, d2211, d3210, d3222, d4410, d4422, d5220, d5232, d5421, d5433    -
+     *    dedt        -
+     *    didt        -
+     *    dmdt        -
+     *    dndt        -
+     *    dnodt       -
+     *    domdt       -
+     *    del1, del2, del3        -
+     *    ses  , sghl , sghs , sgs  , shl  , shs  , sis  , sls
+     *    theta       -
+     *    xfact       -
+     *    xlamo       -
+     *    xli         -
+     *    xni
+     *
+     *  locals        :
+     *    ainv2       -
+     *    aonv        -
+     *    cosisq      -
+     *    eoc         -
+     *    f220, f221, f311, f321, f322, f330, f441, f442, f522, f523, f542, f543  -
+     *    g200, g201, g211, g300, g310, g322, g410, g422, g520, g521, g532, g533  -
+     *    sini2       -
+     *    temp        -
+     *    temp1       -
+     *    theta       -
+     *    xno2        -
+     *
+     *  coupling      :
+     *    getgravconst
+     *
+     *  references    :
+     *    hoots, roehrich, norad spacetrack report #3 1980
+     *    hoots, norad spacetrack report #6 1986
+     *    hoots, schumacher and glover 2004
+     *    vallado, crawford, hujsak, kelso  2006
+     *----------------------------------------------------------------------------
+     */
     Sgp4.dsinit = function (options) {
-        var 
-        // sgp4fix just send in xke as a constant and eliminate getgravconst call
-        // gravconsttype whichconst,
-        xke = options.xke, cosim = options.cosim, argpo = options.argpo, s1 = options.s1, s2 = options.s2, s3 = options.s3, s4 = options.s4, s5 = options.s5, sinim = options.sinim, ss1 = options.ss1, ss2 = options.ss2, ss3 = options.ss3, ss4 = options.ss4, ss5 = options.ss5, sz1 = options.sz1, sz3 = options.sz3, sz11 = options.sz11, sz13 = options.sz13, sz21 = options.sz21, sz23 = options.sz23, sz31 = options.sz31, sz33 = options.sz33, t = options.t, tc = options.tc, gsto = options.gsto, mo = options.mo, mdot = options.mdot, no = options.no, nodeo = options.nodeo, nodedot = options.nodedot, xPIdot = options.xPIdot, z1 = options.z1, z3 = options.z3, z11 = options.z11, z13 = options.z13, z21 = options.z21, z23 = options.z23, z31 = options.z31, z33 = options.z33, ecco = options.ecco, eccsq = options.eccsq;
+        /*
+         * Sgp4fix just send in xke as a constant and eliminate getgravconst call
+         * gravconsttype whichconst,
+         */
+        var xke = options.xke, cosim = options.cosim, argpo = options.argpo, s1 = options.s1, s2 = options.s2, s3 = options.s3, s4 = options.s4, s5 = options.s5, sinim = options.sinim, ss1 = options.ss1, ss2 = options.ss2, ss3 = options.ss3, ss4 = options.ss4, ss5 = options.ss5, sz1 = options.sz1, sz3 = options.sz3, sz11 = options.sz11, sz13 = options.sz13, sz21 = options.sz21, sz23 = options.sz23, sz31 = options.sz31, sz33 = options.sz33, t = options.t, tc = options.tc, gsto = options.gsto, mo = options.mo, mdot = options.mdot, no = options.no, nodeo = options.nodeo, nodedot = options.nodedot, xPIdot = options.xPIdot, z1 = options.z1, z3 = options.z3, z11 = options.z11, z13 = options.z13, z21 = options.z21, z23 = options.z23, z31 = options.z31, z33 = options.z33, ecco = options.ecco, eccsq = options.eccsq;
         var emsq = options.emsq, em = options.em, argpm = options.argpm, inclm = options.inclm, mm = options.mm, nm = options.nm, nodem = options.nodem, irez = options.irez, atime = options.atime, d2201 = options.d2201, d2211 = options.d2211, d3210 = options.d3210, d3222 = options.d3222, d4410 = options.d4410, d4422 = options.d4422, d5220 = options.d5220, d5232 = options.d5232, d5421 = options.d5421, d5433 = options.d5433, dedt = options.dedt, didt = options.didt, dmdt = options.dmdt, dnodt = options.dnodt, domdt = options.domdt, del1 = options.del1, del2 = options.del2, del3 = options.del3, xfact = options.xfact, xlamo = options.xlamo, xli = options.xli, xni = options.xni;
         /* --------------------- local variables ------------------------ */
         /** Ootk -- Some variables imported from outside the class at the top */
@@ -673,15 +693,19 @@ var Sgp4 = /** @class */ (function () {
         var root22 = 1.7891679e-6;
         var root44 = 7.3636953e-9;
         var root54 = 2.1765803e-9;
-        var rptim = 4.37526908801129966e-3; // equates to 7.29211514668855e-5 rad/sec
+        // IDEA: Any way to fix this?
+        // eslint-disable-next-line @typescript-eslint/no-loss-of-precision
+        var rptim = 4.37526908801129966e-3; // Equates to 7.29211514668855e-5 rad/sec
         var root32 = 3.7393792e-7;
         var root52 = 1.1428639e-7;
         var x2o3 = 2.0 / 3.0;
         var znl = 1.5835218e-4;
         var zns = 1.19459e-5;
-        // sgp4fix identify constants and allow alternate values
-        // just xke is used here so pass it in rather than have multiple calls
-        // getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+        /*
+         * Sgp4fix identify constants and allow alternate values
+         * just xke is used here so pass it in rather than have multiple calls
+         * getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+         */
         // -------------------- deep space initialization ------------
         irez = 0;
         if (nm < 0.0052359877 && nm > 0.0034906585) {
@@ -696,7 +720,7 @@ var Sgp4 = /** @class */ (function () {
         var sls = -zns * ss3 * (sz1 + sz3 - 14.0 - 6.0 * emsq);
         var sghs = ss4 * zns * (sz31 + sz33 - 6.0);
         var shs = -zns * ss2 * (sz21 + sz23);
-        // sgp4fix for 180 deg incl
+        // Sgp4fix for 180 deg incl
         if (inclm < 5.2359877e-2 || inclm > PI - 5.2359877e-2) {
             shs = 0.0;
         }
@@ -710,7 +734,7 @@ var Sgp4 = /** @class */ (function () {
         dmdt = sls - znl * s3 * (z1 + z3 - 14.0 - 6.0 * emsq);
         var sghl = s4 * znl * (z31 + z33 - 6.0);
         var shll = -znl * s2 * (z21 + z23);
-        // sgp4fix for 180 deg incl
+        // Sgp4fix for 180 deg incl
         if (inclm < 5.2359877e-2 || inclm > PI - 5.2359877e-2) {
             shll = 0.0;
         }
@@ -728,14 +752,16 @@ var Sgp4 = /** @class */ (function () {
         argpm += domdt * t;
         nodem += dnodt * t;
         mm += dmdt * t;
-        // sgp4fix for negative inclinations
-        // the following if statement should be commented out
-        // if (inclm < 0.0)
-        // {
-        //   inclm  = -inclm;
-        //   argpm  = argpm - PI;
-        //   nodem = nodem + PI;
-        // }
+        /*
+         * Sgp4fix for negative inclinations
+         * the following if statement should be commented out
+         * if (inclm < 0.0)
+         * {
+         *   inclm  = -inclm;
+         *   argpm  = argpm - PI;
+         *   nodem = nodem + PI;
+         * }
+         */
         // -------------- initialize the resonance terms -------------
         if (irez !== 0) {
             var aonv = Math.pow((nm / xke), x2o3);
@@ -748,7 +774,7 @@ var Sgp4 = /** @class */ (function () {
                 emsq = eccsq;
                 var eoc = em * emsq;
                 var g201 = -0.306 - (em - 0.64) * 0.44;
-                var g211 = void 0, g310 = void 0, g322 = void 0, g410 = void 0, g422 = void 0, g520 = void 0, g533 = void 0, g521 = void 0, g532 = void 0;
+                var g211 = void 0, g310 = void 0, g322 = void 0, g410 = void 0, g422 = void 0, g520 = void 0, g521 = void 0, g532 = void 0, g533 = void 0;
                 if (em <= 0.65) {
                     g211 = 3.616 - 13.247 * em + 16.29 * emsq;
                     g310 = -19.302 + 117.39 * em - 228.419 * emsq + 156.591 * eoc;
@@ -875,79 +901,81 @@ var Sgp4 = /** @class */ (function () {
             xli: xli,
             xni: xni,
         };
-    }; // dsinit
-    /*-----------------------------------------------------------------------------
-   *
-   *                           procedure dspace
-   *
-   *  this procedure provides deep space contributions to mean elements for
-   *    perturbing third body.  these effects have been averaged over one
-   *    revolution of the sun and moon.  for earth resonance effects, the
-   *    effects have been averaged over no revolutions of the satellite.
-   *    (mean motion)
-   *
-   *  author        : david vallado                  719-573-2600   28 jun 2005
-   *
-   *  inputs        :
-   *    d2201, d2211, d3210, d3222, d4410, d4422, d5220, d5232, d5421, d5433 -
-   *    dedt        -
-   *    del1, del2, del3  -
-   *    didt        -
-   *    dmdt        -
-   *    dnodt       -
-   *    domdt       -
-   *    irez        - flag for resonance           0-none, 1-one day, 2-half day
-   *    argpo       - argument of perigee
-   *    argpdot     - argument of perigee dot (rate)
-   *    t           - time
-   *    tc          -
-   *    gsto        - gst
-   *    xfact       -
-   *    xlamo       -
-   *    no          - mean motion
-   *    atime       -
-   *    em          - eccentricity
-   *    ft          -
-   *    argpm       - argument of perigee
-   *    inclm       - inclination
-   *    xli         -
-   *    mm          - mean anomaly
-   *    xni         - mean motion
-   *    nodem       - right ascension of ascending node
-   *
-   *  outputs       :
-   *    atime       -
-   *    em          - eccentricity
-   *    argpm       - argument of perigee
-   *    inclm       - inclination
-   *    xli         -
-   *    mm          - mean anomaly
-   *    xni         -
-   *    nodem       - right ascension of ascending node
-   *    dndt        -
-   *    nm          - mean motion
-   *
-   *  locals        :
-   *    delt        -
-   *    ft          -
-   *    theta       -
-   *    x2li        -
-   *    x2omi       -
-   *    xl          -
-   *    xldot       -
-   *    xnddt       -
-   *    xndt        -
-   *    xomi        -
-   *
-   *  coupling      :
-   *    none        -
-   *
-   *  references    :
-   *    hoots, roehrich, norad spacetrack report #3 1980
-   *    hoots, norad spacetrack report #6 1986
-   *    hoots, schumacher and glover 2004
-   *    vallado, crawford, hujsak, kelso  2006
-   ----------------------------------------------------------------------------*/
+    }; // Dsinit
+    /*
+     *-----------------------------------------------------------------------------
+     *
+     *                           procedure dspace
+     *
+     *  this procedure provides deep space contributions to mean elements for
+     *    perturbing third body.  these effects have been averaged over one
+     *    revolution of the sun and moon.  for earth resonance effects, the
+     *    effects have been averaged over no revolutions of the satellite.
+     *    (mean motion)
+     *
+     *  author        : david vallado                  719-573-2600   28 jun 2005
+     *
+     *  inputs        :
+     *    d2201, d2211, d3210, d3222, d4410, d4422, d5220, d5232, d5421, d5433 -
+     *    dedt        -
+     *    del1, del2, del3  -
+     *    didt        -
+     *    dmdt        -
+     *    dnodt       -
+     *    domdt       -
+     *    irez        - flag for resonance           0-none, 1-one day, 2-half day
+     *    argpo       - argument of perigee
+     *    argpdot     - argument of perigee dot (rate)
+     *    t           - time
+     *    tc          -
+     *    gsto        - gst
+     *    xfact       -
+     *    xlamo       -
+     *    no          - mean motion
+     *    atime       -
+     *    em          - eccentricity
+     *    ft          -
+     *    argpm       - argument of perigee
+     *    inclm       - inclination
+     *    xli         -
+     *    mm          - mean anomaly
+     *    xni         - mean motion
+     *    nodem       - right ascension of ascending node
+     *
+     *  outputs       :
+     *    atime       -
+     *    em          - eccentricity
+     *    argpm       - argument of perigee
+     *    inclm       - inclination
+     *    xli         -
+     *    mm          - mean anomaly
+     *    xni         -
+     *    nodem       - right ascension of ascending node
+     *    dndt        -
+     *    nm          - mean motion
+     *
+     *  locals        :
+     *    delt        -
+     *    ft          -
+     *    theta       -
+     *    x2li        -
+     *    x2omi       -
+     *    xl          -
+     *    xldot       -
+     *    xnddt       -
+     *    xndt        -
+     *    xomi        -
+     *
+     *  coupling      :
+     *    none        -
+     *
+     *  references    :
+     *    hoots, roehrich, norad spacetrack report #3 1980
+     *    hoots, norad spacetrack report #6 1986
+     *    hoots, schumacher and glover 2004
+     *    vallado, crawford, hujsak, kelso  2006
+     *----------------------------------------------------------------------------
+     */
     Sgp4.dspace = function (irez, d2201, d2211, d3210, d3222, d4410, d4422, d5220, d5232, d5421, d5433, dedt, del1, del2, del3, didt, dmdt, dnodt, domdt, argpo, argpdot, t, tc, gsto, xfact, xlamo, no, atime, em, argpm, inclm, xli, mm, xni, nodem, nm) {
         /** Ootk -- Some variables imported from outside the class at the top */
         var fasx2 = 0.13130908;
@@ -958,7 +986,9 @@ var Sgp4 = /** @class */ (function () {
         var g44 = 1.8014998;
         var g52 = 1.050833;
         var g54 = 4.4108898;
-        var rptim = 4.37526908801129966e-3; // equates to 7.29211514668855e-5 rad/sec
+        // IDEA: Any way to fix this?
+        // eslint-disable-next-line @typescript-eslint/no-loss-of-precision
+        var rptim = 4.37526908801129966e-3; // Equates to 7.29211514668855e-5 rad/sec
         var stepp = 720.0;
         var stepn = -720.0;
         var step2 = 259200.0;
@@ -970,29 +1000,33 @@ var Sgp4 = /** @class */ (function () {
         argpm += domdt * t;
         nodem += dnodt * t;
         mm += dmdt * t;
-        // sgp4fix for negative inclinations
-        // the following if statement should be commented out
-        // if (inclm < 0.0)
-        // {
-        //   inclm = -inclm;
-        //   argpm = argpm - PI;
-        //   nodem = nodem + PI;
-        // }
+        /*
+         * Sgp4fix for negative inclinations
+         * the following if statement should be commented out
+         * if (inclm < 0.0)
+         * {
+         *   inclm = -inclm;
+         *   argpm = argpm - PI;
+         *   nodem = nodem + PI;
+         * }
+         */
         /* - update resonances : numerical (euler-maclaurin) integration - */
         /* ------------------------- epoch restart ----------------------  */
-        //   sgp4fix for propagator problems
-        //   the following integration works for negative time steps and periods
-        //   the specific changes are unknown because the original code was so convoluted
-        // sgp4fix take out atime = 0.0 and fix for faster operation
-        // ft = 0.0; /** Ootk -- This has no value */
+        /*
+         *   Sgp4fix for propagator problems
+         *   the following integration works for negative time steps and periods
+         *   the specific changes are unknown because the original code was so convoluted
+         */
+        // Sgp4fix take out atime = 0.0 and fix for faster operation
+        // Ft = 0.0; /** Ootk -- This has no value */
         if (irez !== 0) {
-            //  sgp4fix streamline check
+            //  Sgp4fix streamline check
             if (atime === 0.0 || t * atime <= 0.0 || Math.abs(t) < Math.abs(atime)) {
                 atime = 0.0;
                 xni = no;
                 xli = xlamo;
             }
-            // sgp4fix move check outside loop
+            // Sgp4fix move check outside loop
             var delt = void 0;
             if (t > 0.0) {
                 delt = stepp;
@@ -1000,11 +1034,13 @@ var Sgp4 = /** @class */ (function () {
             else {
                 delt = stepn;
             }
-            var xndt = void 0, xldot = void 0, xnddt = void 0, xomi = void 0, x2omi = void 0, x2li = void 0, ft = void 0;
-            var iretn = 381; // added for do loop
+            var ft = void 0, x2li = void 0, x2omi = void 0, xldot = void 0, xnddt = void 0, xndt = void 0, xomi = void 0;
+            var iretn = 381; // Added for do loop
             while (iretn === 381) {
-                //  ------------------- dot terms calculated -------------
-                //  ----------- near - synchronous resonance terms -------
+                /*
+                 *  ------------------- dot terms calculated -------------
+                 *  ----------- near - synchronous resonance terms -------
+                 */
                 if (irez !== 2) {
                     xndt =
                         del1 * Math.sin(xli - fasx2) + del2 * Math.sin(2.0 * (xli - fasx4)) + del3 * Math.sin(3.0 * (xli - fasx6));
@@ -1045,10 +1081,12 @@ var Sgp4 = /** @class */ (function () {
                             d5433 * Math.cos(-xomi + x2li - g54);
                     xnddt *= xldot;
                 }
-                //  ----------------------- integrator -------------------
-                //  sgp4fix move end checks to end of routine
+                /*
+                 *  ----------------------- integrator -------------------
+                 *  sgp4fix move end checks to end of routine
+                 */
                 if (Math.abs(t - atime) >= stepp) {
-                    // iret = 0; /** Ootk -- This has no value */
+                    // Iret = 0; /** Ootk -- This has no value */
                     iretn = 381;
                 }
                 else {
@@ -1060,7 +1098,7 @@ var Sgp4 = /** @class */ (function () {
                     xni += xndt * delt + xnddt * step2;
                     atime += delt;
                 }
-            } // while iretn = 381
+            } // While iretn = 381
             nm = xni + xndt * ft + xnddt * ft * ft * 0.5;
             var xl = xli + xldot * ft + xndt * ft * ft * 0.5;
             if (irez !== 1) {
@@ -1074,76 +1112,82 @@ var Sgp4 = /** @class */ (function () {
             nm = no + dndt;
         }
         return [em, argpm, inclm, mm, nodem, nm];
-    }; // dspace
-    /*-----------------------------------------------------------------------------
-   *
-   *                           procedure initl
-   *
-   *  this procedure initializes the sgp4 propagator. all the initialization is
-   *    consolidated here instead of having multiple loops inside other routines.
-   *
-   *  author        : david vallado                  719-573-2600   28 jun 2005
-   *
-   *  inputs        :
-   *    satn        - satellite number - not needed, placed in satrec
-   *    xke         - reciprocal of tumin
-   *    j2          - j2 zonal harmonic
-   *    ecco        - eccentricity                           0.0 - 1.0
-   *    epoch       - epoch time in days from jan 0, 1950. 0 hr
-   *    inclo       - inclination of satellite
-   *    no          - mean motion of satellite
-   *
-   *  outputs       :
-   *    ainv        - 1.0 / a
-   *    ao          - semi major axis
-   *    con41       -
-   *    con42       - 1.0 - 5.0 cos(i)
-   *    cosio       - cosine of inclination
-   *    cosio2      - cosio squared
-   *    eccsq       - eccentricity squared
-   *    method      - flag for deep space                    'd', 'n'
-   *    omeosq      - 1.0 - ecco * ecco
-   *    posq        - semi-parameter squared
-   *    rp          - radius of perigee
-   *    rteosq      - square root of (1.0 - ecco*ecco)
-   *    sinio       - sine of inclination
-   *    gsto        - gst at time of observation               rad
-   *    no          - mean motion of satellite
-   *
-   *  locals        :
-   *    ak          -
-   *    d1          -
-   *    del         -
-   *    adel        -
-   *    po          -
-   *
-   *  coupling      :
-   *    getgravconst- no longer used
-   *    gstime      - find greenwich sidereal time from the julian date
-   *
-   *  references    :
-   *    hoots, roehrich, norad spacetrack report #3 1980
-   *    hoots, norad spacetrack report #6 1986
-   *    hoots, schumacher and glover 2004
-   *    vallado, crawford, hujsak, kelso  2006
-   ----------------------------------------------------------------------------*/
+    }; // Dspace
+    /*
+     *-----------------------------------------------------------------------------
+     *
+     *                           procedure initl
+     *
+     *  this procedure initializes the sgp4 propagator. all the initialization is
+     *    consolidated here instead of having multiple loops inside other routines.
+     *
+     *  author        : david vallado                  719-573-2600   28 jun 2005
+     *
+     *  inputs        :
+     *    satn        - satellite number - not needed, placed in satrec
+     *    xke         - reciprocal of tumin
+     *    j2          - j2 zonal harmonic
+     *    ecco        - eccentricity                           0.0 - 1.0
+     *    epoch       - epoch time in days from jan 0, 1950. 0 hr
+     *    inclo       - inclination of satellite
+     *    no          - mean motion of satellite
+     *
+     *  outputs       :
+     *    ainv        - 1.0 / a
+     *    ao          - semi major axis
+     *    con41       -
+     *    con42       - 1.0 - 5.0 cos(i)
+     *    cosio       - cosine of inclination
+     *    cosio2      - cosio squared
+     *    eccsq       - eccentricity squared
+     *    method      - flag for deep space                    'd', 'n'
+     *    omeosq      - 1.0 - ecco * ecco
+     *    posq        - semi-parameter squared
+     *    rp          - radius of perigee
+     *    rteosq      - square root of (1.0 - ecco*ecco)
+     *    sinio       - sine of inclination
+     *    gsto        - gst at time of observation               rad
+     *    no          - mean motion of satellite
+     *
+     *  locals        :
+     *    ak          -
+     *    d1          -
+     *    del         -
+     *    adel        -
+     *    po          -
+     *
+     *  coupling      :
+     *    getgravconst- no longer used
+     *    gstime      - find greenwich sidereal time from the julian date
+     *
+     *  references    :
+     *    hoots, roehrich, norad spacetrack report #3 1980
+     *    hoots, norad spacetrack report #6 1986
+     *    hoots, schumacher and glover 2004
+     *    vallado, crawford, hujsak, kelso  2006
+     *----------------------------------------------------------------------------
+     */
     Sgp4.initl = function (options) {
-        // sgp4fix satn not needed. include in satrec in case needed later
-        // int satn,
-        // sgp4fix just pass in xke and j2
-        // gravconsttype whichconst,
+        /*
+         * Sgp4fix satn not needed. include in satrec in case needed later
+         * int satn,
+         * sgp4fix just pass in xke and j2
+         * gravconsttype whichconst,
+         */
         var opsmode = options.opsmode, ecco = options.ecco, epoch = options.epoch, inclo = options.inclo, xke = options.xke, j2 = options.j2;
         var no = options.no;
         /* --------------------- local variables ------------------------ */
         var PI = Math.PI;
         var TAU = PI * 2;
         var x2o3 = 2.0 / 3.0;
-        // sgp4fix use old way of finding gst
+        // Sgp4fix use old way of finding gst
         /** Ootk -- Some variables imported from outside the class at the top */
-        // ----------------------- earth constants ---------------------
-        // sgp4fix identify constants and allow alternate values
-        // only xke and j2 are used here so pass them in directly
-        // getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 )
+        /*
+         * ----------------------- earth constants ---------------------
+         * sgp4fix identify constants and allow alternate values
+         * only xke and j2 are used here so pass them in directly
+         * getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 )
+         */
         // ------------- calculate auxillary epoch quantities ----------
         var eccsq = ecco * ecco;
         var omeosq = 1.0 - eccsq;
@@ -1156,7 +1200,7 @@ var Sgp4 = /** @class */ (function () {
         var delPrime = d1 / (ak * ak);
         var adel = ak * (1.0 - delPrime * delPrime - delPrime * (1.0 / 3.0 + (134.0 * delPrime * delPrime) / 81.0));
         delPrime = d1 / (adel * adel);
-        no = no / (1.0 + delPrime);
+        no /= 1.0 + delPrime;
         var ao = Math.pow((xke / no), x2o3);
         var sinio = Math.sin(inclo);
         var po = ao * omeosq;
@@ -1166,18 +1210,24 @@ var Sgp4 = /** @class */ (function () {
         var posq = po * po;
         var rp = ao * (1.0 - ecco);
         var method = 'n';
-        //  sgp4fix modern approach to finding sidereal time
+        //  Sgp4fix modern approach to finding sidereal time
         /** Ootk -- Continue allowing AFSPC mode for SGP4 Validation */
         var gsto;
-        if (opsmode == 'a') {
-            //  sgp4fix use old way of finding gst
-            //  count integer number of days from 0 jan 1970
+        if (opsmode === 'a') {
+            /*
+             *  Sgp4fix use old way of finding gst
+             *  count integer number of days from 0 jan 1970
+             */
             var ts70 = epoch - 7305.0;
             var ds70 = Math.floor(ts70 + 1.0e-8);
             var tfrac = ts70 - ds70;
-            //  find greenwich location at epoch
+            //  Find greenwich location at epoch
+            // IDEA: Precision issues in JavaScript...
+            // eslint-disable-next-line @typescript-eslint/no-loss-of-precision
             var c1 = 1.72027916940703639e-2;
+            // eslint-disable-next-line @typescript-eslint/no-loss-of-precision
             var thgr70 = 1.7321343856509374;
+            // eslint-disable-next-line @typescript-eslint/no-loss-of-precision
             var fk5r = 5.07551419432269442e-15;
             var c1p2p = c1 + TAU;
             gsto = (thgr70 + c1 * ds70 + c1p2p * tfrac + ts70 * ts70 * fk5r) % TAU;
@@ -1206,95 +1256,99 @@ var Sgp4 = /** @class */ (function () {
             sinio: sinio,
             gsto: gsto,
         };
-    }; // initl
-    /*-----------------------------------------------------------------------------
-   *
-   *                             procedure sgp4init
-   *
-   *  this procedure initializes variables for sgp4.
-   *
-   *  author        : david vallado                  719-573-2600   28 jun 2005
-   *
-   *  inputs        :
-   *    opsmode     - mode of operation afspc or improved 'a', 'i'
-   *    whichconst  - which set of constants to use  72, 84
-   *    satn        - satellite number
-   *    bstar       - sgp4 type drag coefficient              kg/m2er
-   *    ecco        - eccentricity
-   *    epoch       - epoch time in days from jan 0, 1950. 0 hr
-   *    argpo       - argument of perigee (output if ds)
-   *    inclo       - inclination
-   *    mo          - mean anomaly (output if ds)
-   *    no          - mean motion
-   *    nodeo       - right ascension of ascending node
-   *
-   *  outputs       :
-   *    rec      - common values for subsequent calls
-   *    return code - non-zero on error.
-   *                   1 - mean elements, ecc >= 1.0 or ecc < -0.001 or a < 0.95 er
-   *                   2 - mean motion less than 0.0
-   *                   3 - pert elements, ecc < 0.0  or  ecc > 1.0
-   *                   4 - semi-latus rectum < 0.0
-   *                   5 - epoch elements are sub-orbital
-   *                   6 - satellite has decayed
-   *
-   *  locals        :
-   *    cnodm  , snodm  , cosim  , sinim  , cosomm , sinomm
-   *    cc1sq  , cc2    , cc3
-   *    coef   , coef1
-   *    cosio4      -
-   *    day         -
-   *    dndt        -
-   *    em          - eccentricity
-   *    emsq        - eccentricity squared
-   *    eeta        -
-   *    etasq       -
-   *    gam         -
-   *    argpm       - argument of perigee
-   *    nodem       -
-   *    inclm       - inclination
-   *    mm          - mean anomaly
-   *    nm          - mean motion
-   *    perige      - perigee
-   *    PInvsq      -
-   *    psisq       -
-   *    qzms24      -
-   *    rtemsq      -
-   *    s1, s2, s3, s4, s5, s6, s7          -
-   *    sfour       -
-   *    ss1, ss2, ss3, ss4, ss5, ss6, ss7         -
-   *    sz1, sz2, sz3
-   *    sz11, sz12, sz13, sz21, sz22, sz23, sz31, sz32, sz33        -
-   *    tc          -
-   *    temp        -
-   *    temp1, temp2, temp3       -
-   *    tsi         -
-   *    xPIdot      -
-   *    xhdot1      -
-   *    z1, z2, z3          -
-   *    z11, z12, z13, z21, z22, z23, z31, z32, z33         -
-   *
-   *  coupling      :
-   *    getgravconst-
-   *    initl       -
-   *    dscom       -
-   *    dpper       -
-   *    dsinit      -
-   *    sgp4        -
-   *
-   *  references    :
-   *    hoots, roehrich, norad spacetrack report #3 1980
-   *    hoots, norad spacetrack report #6 1986
-   *    hoots, schumacher and glover 2004
-   *    vallado, crawford, hujsak, kelso  2006
-   ----------------------------------------------------------------------------*/
+    }; // Initl
+    /*
+     *-----------------------------------------------------------------------------
+     *
+     *                             procedure sgp4init
+     *
+     *  this procedure initializes variables for sgp4.
+     *
+     *  author        : david vallado                  719-573-2600   28 jun 2005
+     *
+     *  inputs        :
+     *    opsmode     - mode of operation afspc or improved 'a', 'i'
+     *    whichconst  - which set of constants to use  72, 84
+     *    satn        - satellite number
+     *    bstar       - sgp4 type drag coefficient              kg/m2er
+     *    ecco        - eccentricity
+     *    epoch       - epoch time in days from jan 0, 1950. 0 hr
+     *    argpo       - argument of perigee (output if ds)
+     *    inclo       - inclination
+     *    mo          - mean anomaly (output if ds)
+     *    no          - mean motion
+     *    nodeo       - right ascension of ascending node
+     *
+     *  outputs       :
+     *    rec      - common values for subsequent calls
+     *    return code - non-zero on error.
+     *                   1 - mean elements, ecc >= 1.0 or ecc < -0.001 or a < 0.95 er
+     *                   2 - mean motion less than 0.0
+     *                   3 - pert elements, ecc < 0.0  or  ecc > 1.0
+     *                   4 - semi-latus rectum < 0.0
+     *                   5 - epoch elements are sub-orbital
+     *                   6 - satellite has decayed
+     *
+     *  locals        :
+     *    cnodm  , snodm  , cosim  , sinim  , cosomm , sinomm
+     *    cc1sq  , cc2    , cc3
+     *    coef   , coef1
+     *    cosio4      -
+     *    day         -
+     *    dndt        -
+     *    em          - eccentricity
+     *    emsq        - eccentricity squared
+     *    eeta        -
+     *    etasq       -
+     *    gam         -
+     *    argpm       - argument of perigee
+     *    nodem       -
+     *    inclm       - inclination
+     *    mm          - mean anomaly
+     *    nm          - mean motion
+     *    perige      - perigee
+     *    PInvsq      -
+     *    psisq       -
+     *    qzms24      -
+     *    rtemsq      -
+     *    s1, s2, s3, s4, s5, s6, s7          -
+     *    sfour       -
+     *    ss1, ss2, ss3, ss4, ss5, ss6, ss7         -
+     *    sz1, sz2, sz3
+     *    sz11, sz12, sz13, sz21, sz22, sz23, sz31, sz32, sz33        -
+     *    tc          -
+     *    temp        -
+     *    temp1, temp2, temp3       -
+     *    tsi         -
+     *    xPIdot      -
+     *    xhdot1      -
+     *    z1, z2, z3          -
+     *    z11, z12, z13, z21, z22, z23, z31, z32, z33         -
+     *
+     *  coupling      :
+     *    getgravconst-
+     *    initl       -
+     *    dscom       -
+     *    dpper       -
+     *    dsinit      -
+     *    sgp4        -
+     *
+     *  references    :
+     *    hoots, roehrich, norad spacetrack report #3 1980
+     *    hoots, norad spacetrack report #6 1986
+     *    hoots, schumacher and glover 2004
+     *    vallado, crawford, hujsak, kelso  2006
+     *----------------------------------------------------------------------------
+     */
     Sgp4.sgp4init = function (satrec, options) {
         /* eslint-disable no-param-reassign */
         var _a = options.whichconst, whichconst = _a === void 0 ? 'wgs72' : _a, _b = options.opsmode, opsmode = _b === void 0 ? 'i' : _b, _c = options.satn, satn = _c === void 0 ? satrec.satnum : _c, epoch = options.epoch, xbstar = options.xbstar, xecco = options.xecco, xargpo = options.xargpo, xinclo = options.xinclo, xndot = options.xndot, xnddot = options.xnddot, xmo = options.xmo, xno = options.xno, xnodeo = options.xnodeo;
         /* ------------------------ initialization --------------------- */
-        // sgp4fix divisor for divide by zero check on inclination
-        // the old check used 1.0 + Math.cos(PI-1.0e-9), but then compared it to
-        // 1.5 e-12, so the threshold was changed to 1.5e-12 for consistency
+        /*
+         * Sgp4fix divisor for divide by zero check on inclination
+         * the old check used 1.0 + Math.cos(PI-1.0e-9), but then compared it to
+         * 1.5 e-12, so the threshold was changed to 1.5e-12 for consistency
+         */
         // ----------- set all near earth variables to zero ------------
         satrec.isimp = 0;
         satrec.method = 'n';
@@ -1381,8 +1435,10 @@ var Sgp4 = /** @class */ (function () {
         satrec.xli = 0.0;
         satrec.xni = 0.0;
         /* ------------------------ earth constants ----------------------- */
-        // sgp4fix identify constants and allow alternate values
-        // this is now the only call for the constants
+        /*
+         * Sgp4fix identify constants and allow alternate values
+         * this is now the only call for the constants
+         */
         var gravResults = Sgp4.getgravconst(whichconst);
         satrec.tumin = gravResults.tumin;
         satrec.mus = gravResults.mus;
@@ -1396,11 +1452,12 @@ var Sgp4 = /** @class */ (function () {
         var j4 = gravResults.j4;
         var xke = gravResults.xke;
         var j3oj2 = gravResults.j3oj2;
-        //-------------------------------------------------------------------------
+        // -------------------------------------------------------------------------
         satrec.error = 0;
         satrec.operationmode = opsmode;
-        // new alpha5 or 9-digit number
-        /** Ootk -- Using JS code for string manipulation but same effect
+        // New alpha5 or 9-digit number
+        /**
+         * Ootk -- Using JS code for string manipulation but same effect
          * Ex. A2525 = 102525
          * Ex. Z1234 = 351234
          */
@@ -1411,33 +1468,39 @@ var Sgp4 = /** @class */ (function () {
         else {
             satrec.satnum = satn;
         }
-        // sgp4fix - note the following variables are also passed directly via satrec.
-        // it is possible to streamline the sgp4init call by deleting the "x"
-        // variables, but the user would need to set the satrec.* values first. we
-        // include the additional assignments in case twoline2rv is not used.
+        /*
+         * Sgp4fix - note the following variables are also passed directly via satrec.
+         * it is possible to streamline the sgp4init call by deleting the "x"
+         * variables, but the user would need to set the satrec.* values first. we
+         * include the additional assignments in case twoline2rv is not used.
+         */
         satrec.bstar = xbstar;
-        // sgp4fix allow additional parameters in the struct
+        // Sgp4fix allow additional parameters in the struct
         satrec.ndot = xndot;
         satrec.nddot = xnddot;
         satrec.ecco = xecco;
         satrec.argpo = xargpo;
         satrec.inclo = xinclo;
         satrec.mo = xmo;
-        // sgp4fix rename variables to clarify which mean motion is intended
+        // Sgp4fix rename variables to clarify which mean motion is intended
         satrec.no = xno;
         satrec.nodeo = xnodeo;
-        // single averaged mean elements
-        // satrec.am = satrec.em = satrec.im = satrec.Om = satrec.mm = satrec.nm = 0.0;
-        // ------------------------ earth constants -----------------------
-        // sgp4fix identify constants and allow alternate values
-        // getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+        /*
+         * Single averaged mean elements
+         * satrec.am = satrec.em = satrec.im = satrec.Om = satrec.mm = satrec.nm = 0.0;
+         */
+        /*
+         * ------------------------ earth constants -----------------------
+         * sgp4fix identify constants and allow alternate values
+         * getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+         */
         var ss = 78.0 / satrec.radiusearthkm + 1.0;
-        // sgp4fix use multiply for speed instead of pow
+        // Sgp4fix use multiply for speed instead of pow
         var qzms2ttemp = (120.0 - 78.0) / satrec.radiusearthkm;
         var qzms2t = qzms2ttemp * qzms2ttemp * qzms2ttemp * qzms2ttemp;
         satrec.init = 'y';
         satrec.t = 0.0;
-        // sgp4fix remove satn as it is not needed in initl
+        // Sgp4fix remove satn as it is not needed in initl
         var initlOptions = {
             ecco: satrec.ecco,
             epoch: epoch,
@@ -1453,18 +1516,20 @@ var Sgp4 = /** @class */ (function () {
         satrec.no = initlResult.no;
         satrec.con41 = initlResult.con41;
         satrec.gsto = initlResult.gsto;
-        satrec.a = Math.pow(satrec.no * satrec.tumin, -2.0 / 3.0);
+        satrec.a = Math.pow((satrec.no * satrec.tumin), (-2.0 / 3.0));
         satrec.alta = satrec.a * (1.0 + satrec.ecco) - 1.0;
         satrec.altp = satrec.a * (1.0 - satrec.ecco) - 1.0;
         satrec.error = 0;
-        // sgp4fix remove this check as it is unnecessary
-        // the mrt check in sgp4 handles decaying satellite cases even if the starting
-        // condition is below the surface of te earth
-        // if (rp < 1.0)
-        // {
-        //   printf("// *** satn%d epoch elts sub-orbital ***\n", satn);
-        //   satrec.error = 5;
-        // }
+        /*
+         * Sgp4fix remove this check as it is unnecessary
+         * the mrt check in sgp4 handles decaying satellite cases even if the starting
+         * condition is below the surface of te earth
+         * if (rp < 1.0)
+         * {
+         *   printf("// *** satn%d epoch elts sub-orbital ***\n", satn);
+         *   satrec.error = 5;
+         * }
+         */
         if (omeosq >= 0.0 || satrec.no >= 0.0) {
             satrec.isimp = 0;
             if (rp < 220.0 / satrec.radiusearthkm + 1.0) {
@@ -1479,7 +1544,7 @@ var Sgp4 = /** @class */ (function () {
                 if (perige < 98.0) {
                     sfour = 20.0;
                 }
-                // sgp4fix use multiply for speed instead of pow
+                // Sgp4fix use multiply for speed instead of pow
                 var qzms24temp = (120.0 - sfour) / satrec.radiusearthkm;
                 qzms24 = qzms24temp * qzms24temp * qzms24temp * qzms24temp;
                 sfour = sfour / satrec.radiusearthkm + 1.0;
@@ -1536,7 +1601,7 @@ var Sgp4 = /** @class */ (function () {
             }
             satrec.nodecf = 3.5 * omeosq * xhdot1 * satrec.cc1;
             satrec.t2cof = 1.5 * satrec.cc1;
-            // sgp4fix for divide by zero with xinco = 180 deg
+            // Sgp4fix for divide by zero with xinco = 180 deg
             if (Math.abs(cosio + 1.0) > 1.5e-12) {
                 satrec.xlcof = (-0.25 * j3oj2 * sinio * (3.0 + 5.0 * cosio)) / (1.0 + cosio);
             }
@@ -1544,7 +1609,7 @@ var Sgp4 = /** @class */ (function () {
                 satrec.xlcof = (-0.25 * j3oj2 * sinio * (3.0 + 5.0 * cosio)) / temp4;
             }
             satrec.aycof = -0.5 * j3oj2 * sinio;
-            // sgp4fix use multiply for speed instead of pow
+            // Sgp4fix use multiply for speed instead of pow
             var delmotemp = 1.0 + satrec.eta * Math.cos(satrec.mo);
             satrec.delmo = delmotemp * delmotemp * delmotemp;
             satrec.sinmao = Math.sin(satrec.mo);
@@ -1763,109 +1828,119 @@ var Sgp4 = /** @class */ (function () {
                             12.0 * satrec.cc1 * satrec.d3 +
                             6.0 * satrec.d2 * satrec.d2 +
                             15.0 * cc1sq * (2.0 * satrec.d2 + cc1sq));
-            } // if omeosq = 0 ...
-            /* finally propogate to zero epoch to initialize all others. */
-            // sgp4fix take out check to let satellites process until they are actually below earth surface
-            // if(satrec.error == 0)
+            } // If omeosq = 0 ...
+            /* Finally propogate to zero epoch to initialize all others. */
+            /*
+             * Sgp4fix take out check to let satellites process until they are actually below earth surface
+             * if(satrec.error == 0)
+             */
         }
         Sgp4.propagate(satrec, 0);
         satrec.init = 'n';
-        //sgp4fix return boolean. satrec.error contains any error codes
-        // return satrec; -- no reason to return anything in JS
-    }; // sgp4init
-    /*----------------------------------------------------------------------------
-   *
-   *                             procedure sgp4
-   *
-   *  this procedure is the sgp4 prediction model from space command. this is an
-   *    updated and combined version of sgp4 and sdp4, which were originally
-   *    published separately in spacetrack report //3. this version follows the
-   *    methodology from the aiaa paper (2006) describing the history and
-   *    development of the code.
-   *
-   *  author        : david vallado                  719-573-2600   28 jun 2005
-   *
-   *  inputs        :
-   *    satrec  - initialised structure from sgp4init() call.
-   *    tsince  - time since epoch (minutes)
-   *
-   *  outputs       :
-   *    r           - position vector                     km
-   *    v           - velocity                            km/sec
-   *  return code - non-zero on error.
-   *                   1 - mean elements, ecc >= 1.0 or ecc < -0.001 or a < 0.95 er
-   *                   2 - mean motion less than 0.0
-   *                   3 - pert elements, ecc < 0.0  or  ecc > 1.0
-   *                   4 - semi-latus rectum < 0.0
-   *                   5 - epoch elements are sub-orbital
-   *                   6 - satellite has decayed
-   *
-   *  locals        :
-   *    am          -
-   *    axnl, aynl        -
-   *    betal       -
-   *    cosim   , sinim   , cosomm  , sinomm  , cnod    , snod    , cos2u   ,
-   *    sin2u   , coseo1  , sineo1  , cosi    , sini    , cosip   , sinip   ,
-   *    cosisq  , cossu   , sinsu   , cosu    , sinu
-   *    delm        -
-   *    delomg      -
-   *    dndt        -
-   *    eccm        -
-   *    emsq        -
-   *    ecose       -
-   *    el2         -
-   *    eo1         -
-   *    eccp        -
-   *    esine       -
-   *    argpm       -
-   *    argpp       -
-   *    omgadf      -
-   *    pl          -
-   *    r           -
-   *    rtemsq      -
-   *    rdotl       -
-   *    rl          -
-   *    rvdot       -
-   *    rvdotl      -
-   *    su          -
-   *    t2  , t3   , t4    , tc
-   *    tem5, temp , temp1 , temp2  , tempa  , tempe  , templ
-   *    u   , ux   , uy    , uz     , vx     , vy     , vz
-   *    inclm       - inclination
-   *    mm          - mean anomaly
-   *    nm          - mean motion
-   *    nodem       - right asc of ascending node
-   *    xinc        -
-   *    xincp       -
-   *    xl          -
-   *    xlm         -
-   *    mp          -
-   *    xmdf        -
-   *    xmx         -
-   *    xmy         -
-   *    nodedf      -
-   *    xnode       -
-   *    nodep       -
-   *    np          -
-   *
-   *  coupling      :
-   *    getgravconst-
-   *    dpper
-   *    dspace
-   *
-   *  references    :
-   *    hoots, roehrich, norad spacetrack report //3 1980
-   *    hoots, norad spacetrack report //6 1986
-   *    hoots, schumacher and glover 2004
-   *    vallado, crawford, hujsak, kelso  2006
-   ----------------------------------------------------------------------------*/
+        /*
+         * Sgp4fix return boolean. satrec.error contains any error codes
+         *  return satrec; -- no reason to return anything in JS
+         */
+    }; // Sgp4init
+    /*
+     *----------------------------------------------------------------------------
+     *
+     *                             procedure sgp4
+     *
+     *  this procedure is the sgp4 prediction model from space command. this is an
+     *    updated and combined version of sgp4 and sdp4, which were originally
+     *    published separately in spacetrack report //3. this version follows the
+     *    methodology from the aiaa paper (2006) describing the history and
+     *    development of the code.
+     *
+     *  author        : david vallado                  719-573-2600   28 jun 2005
+     *
+     *  inputs        :
+     *    satrec  - initialised structure from sgp4init() call.
+     *    tsince  - time since epoch (minutes)
+     *
+     *  outputs       :
+     *    r           - position vector                     km
+     *    v           - velocity                            km/sec
+     *  return code - non-zero on error.
+     *                   1 - mean elements, ecc >= 1.0 or ecc < -0.001 or a < 0.95 er
+     *                   2 - mean motion less than 0.0
+     *                   3 - pert elements, ecc < 0.0  or  ecc > 1.0
+     *                   4 - semi-latus rectum < 0.0
+     *                   5 - epoch elements are sub-orbital
+     *                   6 - satellite has decayed
+     *
+     *  locals        :
+     *    am          -
+     *    axnl, aynl        -
+     *    betal       -
+     *    cosim   , sinim   , cosomm  , sinomm  , cnod    , snod    , cos2u   ,
+     *    sin2u   , coseo1  , sineo1  , cosi    , sini    , cosip   , sinip   ,
+     *    cosisq  , cossu   , sinsu   , cosu    , sinu
+     *    delm        -
+     *    delomg      -
+     *    dndt        -
+     *    eccm        -
+     *    emsq        -
+     *    ecose       -
+     *    el2         -
+     *    eo1         -
+     *    eccp        -
+     *    esine       -
+     *    argpm       -
+     *    argpp       -
+     *    omgadf      -
+     *    pl          -
+     *    r           -
+     *    rtemsq      -
+     *    rdotl       -
+     *    rl          -
+     *    rvdot       -
+     *    rvdotl      -
+     *    su          -
+     *    t2  , t3   , t4    , tc
+     *    tem5, temp , temp1 , temp2  , tempa  , tempe  , templ
+     *    u   , ux   , uy    , uz     , vx     , vy     , vz
+     *    inclm       - inclination
+     *    mm          - mean anomaly
+     *    nm          - mean motion
+     *    nodem       - right asc of ascending node
+     *    xinc        -
+     *    xincp       -
+     *    xl          -
+     *    xlm         -
+     *    mp          -
+     *    xmdf        -
+     *    xmx         -
+     *    xmy         -
+     *    nodedf      -
+     *    xnode       -
+     *    nodep       -
+     *    np          -
+     *
+     *  coupling      :
+     *    getgravconst-
+     *    dpper
+     *    dspace
+     *
+     *  references    :
+     *    hoots, roehrich, norad spacetrack report //3 1980
+     *    hoots, norad spacetrack report //6 1986
+     *    hoots, schumacher and glover 2004
+     *    vallado, crawford, hujsak, kelso  2006
+     *----------------------------------------------------------------------------
+     */
     Sgp4.propagate = function (satrec, tsince) {
         /* ------------------ set mathematical constants --------------- */
-        // sgp4fix divisor for divide by zero check on inclination
-        // the old check used 1.0 + cos(PI-1.0e-9), but then compared it to
-        // 1.5 e-12, so the threshold was changed to 1.5e-12 for consistency
-        // sgp4fix identify constants and allow alternate values
-        // getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+        /*
+         * Sgp4fix divisor for divide by zero check on inclination
+         * the old check used 1.0 + cos(PI-1.0e-9), but then compared it to
+         * 1.5 e-12, so the threshold was changed to 1.5e-12 for consistency
+         */
+        /*
+         * Sgp4fix identify constants and allow alternate values
+         * getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+         */
         var xke = satrec.xke, j2 = satrec.j2, j3oj2 = satrec.j3oj2;
         var vkmpersec = (satrec.radiusearthkm * satrec.xke) / 60.0;
         // --------------------- clear sgp4 error flag -----------------
@@ -1884,7 +1959,7 @@ var Sgp4 = /** @class */ (function () {
         var templ = satrec.t2cof * t2;
         if (satrec.isimp !== 1) {
             var delomg = satrec.omgcof * satrec.t;
-            //  sgp4fix use mutliply for speed instead of pow
+            //  Sgp4fix use mutliply for speed instead of pow
             var delmtemp = 1.0 + satrec.eta * Math.cos(xmdf);
             var delm = satrec.xmcof * (delmtemp * delmtemp * delmtemp - satrec.delmo);
             var temp_1 = delomg + delm;
@@ -1903,27 +1978,31 @@ var Sgp4 = /** @class */ (function () {
             var tc = satrec.t;
             var dspaceResult = Sgp4.dspace(satrec.irez, satrec.d2201, satrec.d2211, satrec.d3210, satrec.d3222, satrec.d4410, satrec.d4422, satrec.d5220, satrec.d5232, satrec.d5421, satrec.d5433, satrec.dedt, satrec.del1, satrec.del2, satrec.del3, satrec.didt, satrec.dmdt, satrec.dnodt, satrec.domdt, satrec.argpo, satrec.argpdot, satrec.t, tc, satrec.gsto, satrec.xfact, satrec.xlamo, satrec.no, satrec.atime, em, argpm, inclm, satrec.xli, mm, satrec.xni, nodem, nm);
             em = dspaceResult[0], argpm = dspaceResult[1], inclm = dspaceResult[2], mm = dspaceResult[3], nodem = dspaceResult[4], nm = dspaceResult[5];
-        } // if methjod = d
+        } // If methjod = d
         if (nm <= 0.0) {
-            // printf("// error nm %f\n", nm);
+            // Printf("// error nm %f\n", nm);
             satrec.error = 2;
-            // sgp4fix add return
+            // Sgp4fix add return
             return { position: false, velocity: false };
         }
         var am = Math.pow((xke / nm), x2o3) * tempa * tempa;
         nm = xke / Math.pow(am, 1.5);
         em -= tempe;
-        // fix tolerance for error recognition
-        // sgp4fix am is fixed from the previous nm check
+        /*
+         * Fix tolerance for error recognition
+         * sgp4fix am is fixed from the previous nm check
+         */
         /* istanbul ignore next | This is no longer possible*/
         if (em >= 1.0 || em < -0.001) {
-            // || (am < 0.95)
-            // printf("// error em %f\n", em);
+            /*
+             * || (am < 0.95)
+             * printf("// error em %f\n", em);
+             */
             satrec.error = 1;
-            // sgp4fix to return if there is an error in eccentricity
+            // Sgp4fix to return if there is an error in eccentricity
             return { position: false, velocity: false };
         }
-        //  sgp4fix fix tolerance to avoid a divide by zero
+        //  Sgp4fix fix tolerance to avoid a divide by zero
         if (em < 1.0e-6) {
             em = 1.0e-6;
         }
@@ -1933,15 +2012,17 @@ var Sgp4 = /** @class */ (function () {
         argpm %= TAU;
         xlm %= TAU;
         mm = (xlm - argpm - nodem) % TAU;
-        // sgp4fix recover singly averaged mean elements
-        // TODO: Skip these?
-        // satrec.am = am;
-        // satrec.em = em;
-        // satrec.im = inclm;
-        // satrec.Om = nodem;
-        // satrec.om = argpm;
-        // satrec.mm = mm;
-        // satrec.nm = nm;
+        /*
+         * Sgp4fix recover singly averaged mean elements
+         * TODO: Skip these?
+         * satrec.am = am;
+         * satrec.em = em;
+         * satrec.im = inclm;
+         * satrec.Om = nodem;
+         * satrec.om = argpm;
+         * satrec.mm = mm;
+         * satrec.nm = nm;
+         */
         // ----------------- compute extra mean quantities -------------
         var sinim = Math.sin(inclm);
         var cosim = Math.cos(inclm);
@@ -1973,9 +2054,9 @@ var Sgp4 = /** @class */ (function () {
                 argpp -= PI;
             }
             if (ep < 0.0 || ep > 1.0) {
-                //  printf("// error ep %f\n", ep);
+                //  Printf("// error ep %f\n", ep);
                 satrec.error = 3;
-                //  sgp4fix add return
+                //  Sgp4fix add return
                 return { position: false, velocity: false };
             }
         }
@@ -1984,7 +2065,7 @@ var Sgp4 = /** @class */ (function () {
             sinip = Math.sin(xincp);
             cosip = Math.cos(xincp);
             satrec.aycof = -0.5 * j3oj2 * sinip;
-            //  sgp4fix for divide by zero for xincp = 180 deg
+            //  Sgp4fix for divide by zero for xincp = 180 deg
             /* istanbul ignore next | This is no longer possible*/
             if (Math.abs(cosip + 1.0) > 1.5e-12) {
                 satrec.xlcof = (-0.25 * j3oj2 * sinip * (3.0 + 5.0 * cosip)) / (1.0 + cosip);
@@ -2002,9 +2083,11 @@ var Sgp4 = /** @class */ (function () {
         var eo1 = u;
         var tem5 = 9999.9;
         var ktr = 1;
-        //    sgp4fix for kepler iteration
-        //    the following iteration needs better limits on corrections
-        var sineo1, coseo1;
+        /*
+         *    Sgp4fix for kepler iteration
+         *    the following iteration needs better limits on corrections
+         */
+        var coseo1, sineo1;
         while (Math.abs(tem5) >= 1.0e-12 && ktr <= 10) {
             sineo1 = Math.sin(eo1);
             coseo1 = Math.cos(eo1);
@@ -2027,9 +2110,9 @@ var Sgp4 = /** @class */ (function () {
         var el2 = axnl * axnl + aynl * aynl;
         var pl = am * (1.0 - el2);
         if (pl < 0.0) {
-            //  printf("// error pl %f\n", pl);
+            //  Printf("// error pl %f\n", pl);
             satrec.error = 4;
-            //  sgp4fix add return
+            //  Sgp4fix add return
             return { position: false, velocity: false };
         }
         var rl = am * (1.0 - ecose);
@@ -2054,9 +2137,9 @@ var Sgp4 = /** @class */ (function () {
         }
         var mrt = rl * (1.0 - 1.5 * temp2 * betal * satrec.con41) + 0.5 * temp1 * satrec.x1mth2 * cos2u;
         /** Moved this up to reduce unnecessary computation if you are going to return false anyway */
-        // sgp4fix for decaying satellites
+        // Sgp4fix for decaying satellites
         if (mrt < 1.0) {
-            // printf("// decay condition %11.6f \n",mrt);
+            // Printf("// decay condition %11.6f \n",mrt);
             satrec.error = 6;
             return {
                 position: false,
@@ -2099,43 +2182,45 @@ var Sgp4 = /** @class */ (function () {
             velocity: v,
         };
     };
-    /* -----------------------------------------------------------------------------
-    *
-    *                           function getgravconst
-    *
-    *  this function gets constants for the propagator. note that mu is identified to
-    *    facilitiate comparisons with newer models. the common useage is wgs72.
-    *
-    *  author        : david vallado                  719-573-2600   21 jul 2006
-    *
-    *  inputs        :
-    *    whichconst  - which set of constants to use  wgs72old, wgs72, wgs84
-    *
-    *  outputs       :
-    *    tumin       - minutes in one time unit
-    *    mu          - earth gravitational parameter
-    *    radiusearthkm - radius of the earth in km
-    *    xke         - reciprocal of tumin
-    *    j2, j3, j4  - un-normalized zonal harmonic values
-    *    j3oj2       - j3 divided by j2
-    *
-    *  locals        :
-    *
-    *  coupling      :
-    *    none
-    *
-    *  references    :
-    *    norad spacetrack report #3
-    *    vallado, crawford, hujsak, kelso  2006
-    --------------------------------------------------------------------------- */
+    /*
+     * -----------------------------------------------------------------------------
+     *
+     *                           function getgravconst
+     *
+     *  this function gets constants for the propagator. note that mu is identified to
+     *    facilitiate comparisons with newer models. the common useage is wgs72.
+     *
+     *  author        : david vallado                  719-573-2600   21 jul 2006
+     *
+     *  inputs        :
+     *    whichconst  - which set of constants to use  wgs72old, wgs72, wgs84
+     *
+     *  outputs       :
+     *    tumin       - minutes in one time unit
+     *    mu          - earth gravitational parameter
+     *    radiusearthkm - radius of the earth in km
+     *    xke         - reciprocal of tumin
+     *    j2, j3, j4  - un-normalized zonal harmonic values
+     *    j3oj2       - j3 divided by j2
+     *
+     *  locals        :
+     *
+     *  coupling      :
+     *    none
+     *
+     *  references    :
+     *    norad spacetrack report #3
+     *    vallado, crawford, hujsak, kelso  2006
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.getgravconst = function (whichconst) {
-        var tumin, mus, radiusearthkm, xke, j2, j3, j4, j3oj2;
+        var j2, j3, j3oj2, j4, mus, radiusearthkm, tumin, xke;
         switch (whichconst) {
             // -- wgs-72 low precision str#3 constants --
             case 'wgs72old':
-                mus = 398600.79964; // in km3 / s2
-                radiusearthkm = 6378.135; // km
-                xke = 0.0743669161; // reciprocal of tumin
+                mus = 398600.79964; // In km3 / s2
+                radiusearthkm = 6378.135; // Km
+                xke = 0.0743669161; // Reciprocal of tumin
                 tumin = 1.0 / xke;
                 j2 = 0.001082616;
                 j3 = -0.00000253881;
@@ -2144,8 +2229,8 @@ var Sgp4 = /** @class */ (function () {
                 break;
             // ------------ wgs-72 constants ------------
             case 'wgs72':
-                mus = 398600.8; // in km3 / s2
-                radiusearthkm = 6378.135; // km
+                mus = 398600.8; // In km3 / s2
+                radiusearthkm = 6378.135; // Km
                 xke = 60.0 / Math.sqrt((radiusearthkm * radiusearthkm * radiusearthkm) / mus);
                 tumin = 1.0 / xke;
                 j2 = 0.001082616;
@@ -2155,8 +2240,8 @@ var Sgp4 = /** @class */ (function () {
                 break;
             case 'wgs84':
                 // ------------ wgs-84 constants ------------
-                mus = 398600.5; // in km3 / s2
-                radiusearthkm = 6378.137; // km
+                mus = 398600.5; // In km3 / s2
+                radiusearthkm = 6378.137; // Km
                 xke = 60.0 / Math.sqrt((radiusearthkm * radiusearthkm * radiusearthkm) / mus);
                 tumin = 1.0 / xke;
                 j2 = 0.00108262998905;
@@ -2177,46 +2262,48 @@ var Sgp4 = /** @class */ (function () {
             j4: j4,
             j3oj2: j3oj2,
         };
-    }; // getgravconst
-    /* -----------------------------------------------------------------------------
-   *
-   *                           function twoline2rv
-   *
-   *  this function converts the two line element set character string data to
-   *    variables and initializes the sgp4 variables. several intermediate varaibles
-   *    and quantities are determined. note that the result is a structure so multiple
-   *    satellites can be processed simultaneously without having to reinitialize. the
-   *    verification mode is an important option that permits quick checks of any
-   *    changes to the underlying technical theory. this option works using a
-   *    modified tle file in which the start, stop, and delta time values are
-   *    included at the end of the second line of data. this only works with the
-   *    verification mode. the catalog mode simply propagates from -1440 to 1440 min
-   *    from epoch and is useful when performing entire catalog runs.
-   *
-   *  author        : david vallado                  719-573-2600    1 mar 2001
-   *
-   *  inputs        :
-   *    longstr1    - first line of the tle
-   *    longstr2    - second line of the tle
-   *    typerun     - type of run                    verification 'v', catalog 'c',
-   *                                                 manual 'm'
-   *    typeinput   - type of manual input           mfe 'm', epoch 'e', dayofyr 'd'
-   *    opsmode     - mode of operation afspc or improved 'a', 'i'
-   *    whichconst  - which set of constants to use  72, 84
-   *
-   *  outputs       :
-   *    satrec      - structure containing all the sgp4 satellite information
-   *
-   *  coupling      :
-   *    getgravconst-
-   *    days2mdhms  - conversion of days to month, day, hour, minute, second
-   *    jday        - convert day month year hour minute second into julian date
-   *    sgp4init    - initialize the sgp4 variables
-   *
-   *  references    :
-   *    norad spacetrack report #3
-   *    vallado, crawford, hujsak, kelso  2006
-   --------------------------------------------------------------------------- */
+    }; // Getgravconst
+    /*
+     * -----------------------------------------------------------------------------
+     *
+     *                           function twoline2rv
+     *
+     *  this function converts the two line element set character string data to
+     *    variables and initializes the sgp4 variables. several intermediate varaibles
+     *    and quantities are determined. note that the result is a structure so multiple
+     *    satellites can be processed simultaneously without having to reinitialize. the
+     *    verification mode is an important option that permits quick checks of any
+     *    changes to the underlying technical theory. this option works using a
+     *    modified tle file in which the start, stop, and delta time values are
+     *    included at the end of the second line of data. this only works with the
+     *    verification mode. the catalog mode simply propagates from -1440 to 1440 min
+     *    from epoch and is useful when performing entire catalog runs.
+     *
+     *  author        : david vallado                  719-573-2600    1 mar 2001
+     *
+     *  inputs        :
+     *    longstr1    - first line of the tle
+     *    longstr2    - second line of the tle
+     *    typerun     - type of run                    verification 'v', catalog 'c',
+     *                                                 manual 'm'
+     *    typeinput   - type of manual input           mfe 'm', epoch 'e', dayofyr 'd'
+     *    opsmode     - mode of operation afspc or improved 'a', 'i'
+     *    whichconst  - which set of constants to use  72, 84
+     *
+     *  outputs       :
+     *    satrec      - structure containing all the sgp4 satellite information
+     *
+     *  coupling      :
+     *    getgravconst-
+     *    days2mdhms  - conversion of days to month, day, hour, minute, second
+     *    jday        - convert day month year hour minute second into julian date
+     *    sgp4init    - initialize the sgp4 variables
+     *
+     *  references    :
+     *    norad spacetrack report #3
+     *    vallado, crawford, hujsak, kelso  2006
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.createSatrec = function (tleLine1, tleLine2, whichconst, opsmode) {
         if (whichconst === void 0) { whichconst = 'wgs72'; }
         if (opsmode === void 0) { opsmode = 'i'; }
@@ -2340,8 +2427,10 @@ var Sgp4 = /** @class */ (function () {
             d2211: null,
             nm: null,
         };
-        // sgp4fix no longer needed
-        // getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+        /*
+         * Sgp4fix no longer needed
+         * getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
+         */
         var xpdotp = 1440.0 / (2.0 * PI); // 229.1831180523293;
         satrec.error = 0;
         satrec.satnum = tleLine1.substring(2, 7);
@@ -2350,7 +2439,7 @@ var Sgp4 = /** @class */ (function () {
         satrec.ndot = parseFloat(tleLine1.substring(33, 43));
         satrec.nddot = parseFloat("".concat(tleLine1.substring(44, 45), ".").concat(tleLine1.substring(45, 50), "E").concat(tleLine1.substring(50, 52)));
         satrec.bstar = parseFloat("".concat(tleLine1.substring(53, 54), ".").concat(tleLine1.substring(54, 59), "E").concat(tleLine1.substring(59, 61)));
-        // satrec.satnum = tleLine2.substring(2, 7);
+        // Satrec.satnum = tleLine2.substring(2, 7);
         satrec.inclo = parseFloat(tleLine2.substring(8, 16));
         satrec.nodeo = parseFloat(tleLine2.substring(17, 25));
         satrec.ecco = parseFloat(".".concat(tleLine2.substring(26, 33)));
@@ -2358,38 +2447,50 @@ var Sgp4 = /** @class */ (function () {
         satrec.mo = parseFloat(tleLine2.substring(43, 51));
         satrec.no = parseFloat(tleLine2.substring(52, 63));
         // ---- find no, ndot, nddot ----
-        satrec.no /= xpdotp; //   rad/min
+        satrec.no /= xpdotp; //   Rad/min
         /** Ootk -- nexp and ibexp are calculated above using template literals */
-        // satrec.nddot = satrec.nddot * Math.pow(10.0, nexp);
-        // satrec.bstar = satrec.bstar * Math.pow(10.0, ibexp);
-        // ---- convert to sgp4 units ----
-        // satrec.a = (satrec.no * tumin) ** (-2.0 / 3.0);
+        /*
+         * Satrec.nddot = satrec.nddot * Math.pow(10.0, nexp);
+         * satrec.bstar = satrec.bstar * Math.pow(10.0, ibexp);
+         */
+        /*
+         * ---- convert to sgp4 units ----
+         * satrec.a = (satrec.no * tumin) ** (-2.0 / 3.0);
+         */
         /** Ootk -- Not sure why the following two lines are added. 1st and 2nd derivatives aren't even used anymore */
-        // satrec.ndot /= xpdotp * 1440.0; // ? * minperday
-        // satrec.nddot /= xpdotp * 1440.0 * 1440;
+        /*
+         * Satrec.ndot /= xpdotp * 1440.0; // ? * minperday
+         * satrec.nddot /= xpdotp * 1440.0 * 1440;
+         */
         // ---- find standard orbital elements ----
         satrec.inclo *= DEG2RAD;
         satrec.nodeo *= DEG2RAD;
         satrec.argpo *= DEG2RAD;
         satrec.mo *= DEG2RAD;
-        // sgp4fix not needed here
-        // satrec.alta = satrec.a * (1.0 + satrec.ecco) - 1.0;
-        // satrec.altp = satrec.a * (1.0 - satrec.ecco) - 1.0;
-        // ----------------------------------------------------------------
-        // find sgp4epoch time of element set
-        // remember that sgp4 uses units of days from 0 jan 1950 (sgp4epoch)
-        // and minutes from the epoch (time)
-        // ----------------------------------------------------------------
-        // ---------------- temp fix for years from 1957-2056 -------------------
-        // --------- correct fix will occur when year is 4-digit in tle ---------
+        /*
+         * Sgp4fix not needed here
+         * satrec.alta = satrec.a * (1.0 + satrec.ecco) - 1.0;
+         * satrec.altp = satrec.a * (1.0 - satrec.ecco) - 1.0;
+         */
+        /*
+         * ----------------------------------------------------------------
+         * find sgp4epoch time of element set
+         * remember that sgp4 uses units of days from 0 jan 1950 (sgp4epoch)
+         * and minutes from the epoch (time)
+         * ----------------------------------------------------------------
+         */
+        /*
+         * ---------------- temp fix for years from 1957-2056 -------------------
+         * --------- correct fix will occur when year is 4-digit in tle ---------
+         */
         if (satrec.epochyr < 57) {
             year = satrec.epochyr + 2000;
         }
         else {
             year = satrec.epochyr + 1900;
         }
-        var _a = Sgp4.days2mdhms(year, satrec.epochdays), mon = _a.mon, day = _a.day, hr = _a.hr, minute = _a.minute, sec = _a.sec;
-        var jdayRes = Sgp4.jday(year, mon, day, hr, minute, sec);
+        var _a = Sgp4.days2mdhms(year, satrec.epochdays), mon = _a.mon, day = _a.day, hr = _a.hr, min = _a.min, sec = _a.sec;
+        var jdayRes = Sgp4.jday(year, mon, day, hr, min, sec);
         satrec.jdsatepoch = jdayRes.jd + jdayRes.jdFrac;
         //  ---------------- initialize the orbit at sgp4epoch -------------------
         Sgp4.sgp4init(satrec, {
@@ -2408,8 +2509,9 @@ var Sgp4 = /** @class */ (function () {
             xnodeo: satrec.nodeo,
         });
         return satrec;
-    }; // twoline2rv
-    /* -----------------------------------------------------------------------------
+    }; // Twoline2rv
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           function gstime
      *
@@ -2433,10 +2535,11 @@ var Sgp4 = /** @class */ (function () {
      *
      *  references    :
      *    vallado       2004, 191, eq 3-45
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.gstime = function (jdut1) {
         var tut1 = (jdut1 - 2451545.0) / 36525.0;
-        var temp = -6.2e-6 * tut1 * tut1 * tut1 + 0.093104 * tut1 * tut1 + (876600.0 * 3600 + 8640184.812866) * tut1 + 67310.54841; // sec
+        var temp = -6.2e-6 * tut1 * tut1 * tut1 + 0.093104 * tut1 * tut1 + (876600.0 * 3600 + 8640184.812866) * tut1 + 67310.54841; // Sec
         temp = ((temp * DEG2RAD) / 240.0) % TAU; // 360/86400 = 1/240, to deg, to rad
         //  ------------------------ check quadrants ---------------------
         if (temp < 0.0) {
@@ -2448,11 +2551,10 @@ var Sgp4 = /** @class */ (function () {
         if (x < 0.0) {
             return -1.0;
         }
-        else {
-            return 1.0;
-        }
-    }; // sgn
-    /* -----------------------------------------------------------------------------
+        return 1.0;
+    }; // Sgn
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           function mag
      *
@@ -2471,39 +2573,43 @@ var Sgp4 = /** @class */ (function () {
      *
      *  coupling      :
      *    none.
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.mag = function (x) {
         return Math.sqrt(x[0] * x[0] + x[1] * x[1] + x[2] * x[2]);
-    }; // mag
-    /* -----------------------------------------------------------------------------
-    *
-    *                           procedure cross_SGP4
-    *
-    *  this procedure crosses two vectors.
-    *
-    *  author        : david vallado                  719-573-2600    1 mar 2001
-    *
-    *  inputs          description                    range / units
-    *    vec1        - vector number 1
-    *    vec2        - vector number 2
-    *
-    *  outputs       :
-    *    outvec      - vector result of a x b
-    *
-    *  locals        :
-    *    none.
-    *
-    *  coupling      :
-    *    mag           magnitude of a vector
-    ---------------------------------------------------------------------------- */
+    }; // Mag
+    /*
+     * -----------------------------------------------------------------------------
+     *
+     *                           procedure cross_SGP4
+     *
+     *  this procedure crosses two vectors.
+     *
+     *  author        : david vallado                  719-573-2600    1 mar 2001
+     *
+     *  inputs          description                    range / units
+     *    vec1        - vector number 1
+     *    vec2        - vector number 2
+     *
+     *  outputs       :
+     *    outvec      - vector result of a x b
+     *
+     *  locals        :
+     *    none.
+     *
+     *  coupling      :
+     *    mag           magnitude of a vector
+     * ----------------------------------------------------------------------------
+     */
     Sgp4.cross = function (vec1, vec2) {
         var outvec = [0, 0, 0];
         outvec[0] = vec1[1] * vec2[2] - vec1[2] * vec2[1];
         outvec[1] = vec1[2] * vec2[0] - vec1[0] * vec2[2];
         outvec[2] = vec1[0] * vec2[1] - vec1[1] * vec2[0];
         return outvec;
-    }; // end cross
-    /* -----------------------------------------------------------------------------
+    }; // End cross
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           function dot_SGP4
      *
@@ -2523,11 +2629,13 @@ var Sgp4 = /** @class */ (function () {
      *
      *  coupling      :
      *    none.
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.dot = function (x, y) {
         return x[0] * y[0] + x[1] * y[1] + x[2] * y[2];
-    }; // dot
-    /* -----------------------------------------------------------------------------
+    }; // Dot
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           procedure angle_SGP4
      *
@@ -2549,7 +2657,8 @@ var Sgp4 = /** @class */ (function () {
      *
      *  coupling      :
      *    dot           dot product of two vectors
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.angle = function (vec1, vec2) {
         var small = 0.00000001;
         var unknown = 999999.1; /** Ootk -- original 'undefined' is protected in JS */
@@ -2557,14 +2666,15 @@ var Sgp4 = /** @class */ (function () {
         var magv2 = Sgp4.mag(vec2);
         if (magv1 * magv2 > small * small) {
             var temp = Sgp4.dot(vec1, vec2) / (magv1 * magv2);
-            if (Math.abs(temp) > 1.0)
-                temp = Sgp4.sgn(temp) * 1.0;
+            if (Math.abs(temp) > 1.0) {
+                temp = Number(Sgp4.sgn(temp));
+            }
             return Math.acos(temp);
         }
-        else
-            return unknown;
-    }; // angle
-    /* -----------------------------------------------------------------------------
+        return unknown;
+    }; // Angle
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           function asinh_SGP4
      *
@@ -2583,11 +2693,13 @@ var Sgp4 = /** @class */ (function () {
      *
      *  coupling      :
      *    none.
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.asinh = function (xval) {
         return Math.log(xval + Math.sqrt(xval * xval + 1.0));
-    }; // asinh
-    /* -----------------------------------------------------------------------------
+    }; // Asinh
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           function newtonnu_SGP4
      *
@@ -2620,52 +2732,55 @@ var Sgp4 = /** @class */ (function () {
      *
      *  references    :
      *    vallado       2013, 77, alg 5
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.newtonnu = function (ecc, nu) {
         // ---------------------  implementation   ---------------------
         var e0 = 999999.9;
         var m = 999999.9;
         var small = 0.00000001;
-        // --------------------------- circular ------------------------
         if (Math.abs(ecc) < small) {
+            // --------------------------- circular ------------------------
             m = nu;
             e0 = nu;
         }
-        // ---------------------- elliptical -----------------------
         else if (ecc < 1.0 - small) {
+            // ---------------------- elliptical -----------------------
             var sine = (Math.sqrt(1.0 - ecc * ecc) * Math.sin(nu)) / (1.0 + ecc * Math.cos(nu));
             var cose = (ecc + Math.cos(nu)) / (1.0 + ecc * Math.cos(nu));
             e0 = Math.atan2(sine, cose);
             m = e0 - ecc * Math.sin(e0);
         }
-        // -------------------- hyperbolic  --------------------
         else if (ecc > 1.0 + small) {
+            // -------------------- hyperbolic  --------------------
             if (ecc > 1.0 && Math.abs(nu) + 0.00001 < PI - Math.acos(1.0 / ecc)) {
                 var sine = (Math.sqrt(ecc * ecc - 1.0) * Math.sin(nu)) / (1.0 + ecc * Math.cos(nu));
                 e0 = Sgp4.asinh(sine);
                 m = ecc * Sgp4.sinh(e0) - e0;
             }
         }
-        // ----------------- parabolic ---------------------
         else if (Math.abs(nu) < (168.0 * PI) / 180.0) {
+            // ----------------- parabolic ---------------------
             e0 = Math.tan(nu * 0.5);
             m = e0 + (e0 * e0 * e0) / 3.0;
         }
         if (ecc < 1.0) {
-            m = m % (2.0 * PI);
-            if (m < 0.0)
-                m = m + 2.0 * PI;
-            e0 = e0 % (2.0 * PI);
+            m %= 2.0 * PI;
+            if (m < 0.0) {
+                m += 2.0 * PI;
+            }
+            e0 %= 2.0 * PI;
         }
         return {
             e0: e0,
             m: m,
         };
-    }; // newtonnu
+    }; // Newtonnu
     Sgp4.sinh = function (x) {
         return (Math.exp(x) - Math.exp(-x)) / 2;
     };
-    /* -----------------------------------------------------------------------------
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           function rv2coe_SGP4
      *
@@ -2720,7 +2835,8 @@ var Sgp4 = /** @class */ (function () {
      *
      *  references    :
      *    vallado       2013, 113, alg 9, ex 2-5
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.rv2coe = function (r, v, mus) {
         var nbar = [0, 0, 0];
         var ebar = [0, 0, 0];
@@ -2740,14 +2856,18 @@ var Sgp4 = /** @class */ (function () {
         var hk;
         var sme;
         var i;
-        // switch this to an integer msvs seems to have probelms with this and strncpy_s
-        //char typeorbit[2];
+        /*
+         *  Switch this to an integer msvs seems to have probelms with this and strncpy_s
+         * char typeorbit[2];
+         */
         var typeorbit;
-        // here
-        // typeorbit = 1 = 'ei'
-        // typeorbit = 2 = 'ce'
-        // typeorbit = 3 = 'ci'
-        // typeorbit = 4 = 'ee'
+        /*
+         * Here
+         * typeorbit = 1 = 'ei'
+         * typeorbit = 2 = 'ce'
+         * typeorbit = 3 = 'ci'
+         * typeorbit = 4 = 'ee'
+         */
         var halfpi = 0.5 * PI;
         var small = 0.00000001;
         var unknown = 999999.1; /** Ootk -- original undefined is illegal in JS */
@@ -2765,21 +2885,26 @@ var Sgp4 = /** @class */ (function () {
             magn = Sgp4.mag(nbar);
             var c1 = magv * magv - mus / magr;
             rdotv = Sgp4.dot(r, v);
-            for (i = 0; i <= 2; i++)
+            for (i = 0; i <= 2; i++) {
                 ebar[i] = (c1 * r[i] - rdotv * v[i]) / mus;
+            }
             ecc = Sgp4.mag(ebar);
             // ------------  find a e and semi-latus rectum   ----------
             sme = magv * magv * 0.5 - mus / magr;
-            if (Math.abs(sme) > small)
+            if (Math.abs(sme) > small) {
                 a = -mus / (2.0 * sme);
-            else
+            }
+            else {
                 a = infinite;
+            }
             p = (magh * magh) / mus;
             // -----------------  find inclination   -------------------
             hk = hbar[2] / magh;
             incl = Math.acos(hk);
-            // --------  determine type of orbit for later use  --------
-            // ------ elliptical, parabolic, hyperbolic inclined -------
+            /*
+             * --------  determine type of orbit for later use  --------
+             * ------ elliptical, parabolic, hyperbolic inclined -------
+             */
             typeorbit = 1;
             if (ecc < small) {
                 // ----------------  circular equatorial ---------------
@@ -2800,68 +2925,86 @@ var Sgp4 = /** @class */ (function () {
             // ----------  find right ascension of the ascending node ------------
             if (magn > small) {
                 var temp = nbar[0] / magn;
-                if (Math.abs(temp) > 1.0)
+                if (Math.abs(temp) > 1.0) {
                     temp = Sgp4.sgn(temp);
+                }
                 omega = Math.acos(temp);
-                if (nbar[1] < 0.0)
+                if (nbar[1] < 0.0) {
                     omega = TAU - omega;
+                }
             }
-            else
+            else {
                 omega = unknown;
+            }
             // ---------------- find argument of perigee ---------------
-            if (typeorbit == 1) {
+            if (typeorbit === 1) {
                 argp = Sgp4.angle(nbar, ebar);
-                if (ebar[2] < 0.0)
+                if (ebar[2] < 0.0) {
                     argp = TAU - argp;
+                }
             }
-            else
+            else {
                 argp = unknown;
+            }
             // ------------  find true anomaly at epoch    -------------
-            if (typeorbit == 1 || typeorbit == 4) {
+            if (typeorbit === 1 || typeorbit === 4) {
                 nu = Sgp4.angle(ebar, r);
-                if (rdotv < 0.0)
+                if (rdotv < 0.0) {
                     nu = TAU - nu;
+                }
             }
-            else
+            else {
                 nu = unknown;
+            }
             // ----  find argument of latitude - circular inclined -----
-            if (typeorbit == 3) {
+            if (typeorbit === 3) {
                 arglat = Sgp4.angle(nbar, r);
-                if (r[2] < 0.0)
+                if (r[2] < 0.0) {
                     arglat = TAU - arglat;
+                }
                 m = arglat;
             }
-            else
+            else {
                 arglat = unknown;
-            if (ecc > small && typeorbit == 4) {
+            }
+            if (ecc > small && typeorbit === 4) {
                 var temp = ebar[0] / ecc;
-                if (Math.abs(temp) > 1.0)
+                if (Math.abs(temp) > 1.0) {
                     temp = Sgp4.sgn(temp);
+                }
                 lonper = Math.acos(temp);
-                if (ebar[1] < 0.0)
+                if (ebar[1] < 0.0) {
                     lonper = TAU - lonper;
-                if (incl > halfpi)
+                }
+                if (incl > halfpi) {
                     lonper = TAU - lonper;
+                }
             }
-            else
+            else {
                 lonper = unknown;
+            }
             // -------- find true longitude - circular equatorial ------
-            if (magr > small && typeorbit == 2) {
+            if (magr > small && typeorbit === 2) {
                 var temp = r[0] / magr;
-                if (Math.abs(temp) > 1.0)
+                if (Math.abs(temp) > 1.0) {
                     temp = Sgp4.sgn(temp);
+                }
                 truelon = Math.acos(temp);
-                if (r[1] < 0.0)
+                if (r[1] < 0.0) {
                     truelon = TAU - truelon;
-                if (incl > halfpi)
+                }
+                if (incl > halfpi) {
                     truelon = TAU - truelon;
+                }
                 m = truelon;
             }
-            else
+            else {
                 truelon = unknown;
+            }
             // ------------ find mean anomaly for all orbits -----------
-            if (typeorbit == 1 || typeorbit == 4)
+            if (typeorbit === 1 || typeorbit === 4) {
                 m = Sgp4.newtonnu(ecc, nu).m;
+            }
         }
         else {
             p = unknown;
@@ -2889,8 +3032,9 @@ var Sgp4 = /** @class */ (function () {
             truelon: truelon,
             lonper: lonper,
         };
-    }; // rv2coe
-    /* -----------------------------------------------------------------------------
+    }; // Rv2coe
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           procedure jday
      *
@@ -2922,24 +3066,41 @@ var Sgp4 = /** @class */ (function () {
      *  references    :
      *    vallado       2013, 183, alg 14, ex 3-4
      *
-     * --------------------------------------------------------------------------- */
-    Sgp4.jday = function (year, mon, day, hr, minute, sec) {
+     * ---------------------------------------------------------------------------
+     */
+    Sgp4.jday = function (year, mon, day, hr, min, sec, ms) {
+        if (mon === void 0) { mon = 0; }
+        if (day === void 0) { day = 0; }
+        if (hr === void 0) { hr = 0; }
+        if (min === void 0) { min = 0; }
+        if (sec === void 0) { sec = 0; }
+        if (ms === void 0) { ms = 0; }
+        if (year instanceof Date) {
+            mon = year.getUTCMonth() + 1;
+            day = year.getUTCDate();
+            hr = year.getUTCHours();
+            min = year.getUTCMinutes();
+            sec = year.getUTCSeconds();
+            ms = year.getUTCMilliseconds();
+            year = year.getUTCFullYear();
+        }
         var jd = 367.0 * year -
             Math.floor(7 * (year + Math.floor((mon + 9) / 12.0)) * 0.25) +
             Math.floor((275 * mon) / 9.0) +
             day +
-            1721013.5; // use - 678987.0 to go to mjd directly
-        var jdFrac = (sec + minute * 60.0 + hr * 3600.0) / 86400.0;
-        // check that the day and fractional day are correct
+            1721013.5; // Use - 678987.0 to go to mjd directly
+        var jdFrac = (ms / 1000 + sec + min * 60.0 + hr * 3600.0) / 86400.0;
+        // Check that the day and fractional day are correct
         if (Math.abs(jdFrac) > 1.0) {
             var dtt = Math.floor(jdFrac);
-            jd = jd + dtt;
-            jdFrac = jdFrac - dtt;
+            jd += dtt;
+            jdFrac -= dtt;
         }
         // - 0.5*sgn(100.0*year + mon - 190002.5) + 0.5;
         return { jd: jd, jdFrac: jdFrac };
-    }; // jday
-    /* -----------------------------------------------------------------------------
+    }; // Jday
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           procedure days2mdhms
      *
@@ -2974,14 +3135,17 @@ var Sgp4 = /** @class */ (function () {
      *
      *  coupling      :
      *    none.
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.days2mdhms = function (year, days) {
         var lmonth = [31, year % 4 === 0 ? 29 : 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31];
         var dayofyr = Math.floor(days);
         //  ----------------- find month and day of month ----------------
         /** Ootk -- Incorporated in the above declaration */
-        // if ((year % 4) == 0)
-        // lmonth[2] = 29;
+        /*
+         * If ((year % 4) == 0)
+         * lmonth[2] = 29;
+         */
         var i = 1;
         var inttemp = 0;
         while (dayofyr > inttemp + lmonth[i - 1] && i < 12) {
@@ -2994,17 +3158,18 @@ var Sgp4 = /** @class */ (function () {
         var temp = (days - dayofyr) * 24.0;
         var hr = Math.floor(temp);
         temp = (temp - hr) * 60.0;
-        var minute = Math.floor(temp);
-        var sec = (temp - minute) * 60.0;
+        var min = Math.floor(temp);
+        var sec = (temp - min) * 60.0;
         return {
             mon: mon,
             day: day,
             hr: hr,
-            minute: minute,
+            min: min,
             sec: sec,
         };
-    }; // days2mdhms
-    /* -----------------------------------------------------------------------------
+    }; // Days2mdhms
+    /*
+     * -----------------------------------------------------------------------------
      *
      *                           procedure invjday
      *
@@ -3043,20 +3208,21 @@ var Sgp4 = /** @class */ (function () {
      *
      *  references    :
      *    vallado       2013, 203, alg 22, ex 3-13
-     * --------------------------------------------------------------------------- */
+     * ---------------------------------------------------------------------------
+     */
     Sgp4.invjday = function (jd, jdfrac) {
         var leapyrs;
         var days;
-        // check jdfrac for multiple days
+        // Check jdfrac for multiple days
         if (Math.abs(jdfrac) >= 1.0) {
-            jd = jd + Math.floor(jdfrac);
-            jdfrac = jdfrac - Math.floor(jdfrac);
+            jd += Math.floor(jdfrac);
+            jdfrac -= Math.floor(jdfrac);
         }
-        // check for fraction of a day included in the jd
+        // Check for fraction of a day included in the jd
         var dt = jd - Math.floor(jd) - 0.5;
         if (Math.abs(dt) > 0.00000001) {
-            jd = jd - dt;
-            jdfrac = jdfrac + dt;
+            jd -= dt;
+            jdfrac += dt;
         }
         /* --------------- find year and days of the year --------------- */
         var temp = jd - 2415019.5;
@@ -3066,21 +3232,21 @@ var Sgp4 = /** @class */ (function () {
         days = Math.floor(temp - ((year - 1900) * 365.0 + leapyrs));
         /* ------------ check for case of beginning of a year ----------- */
         if (days + jdfrac < 1.0) {
-            year = year - 1;
+            year -= 1;
             leapyrs = Math.floor((year - 1901) * 0.25);
             days = Math.floor(temp - ((year - 1900) * 365.0 + leapyrs));
         }
         /* ----------------- find remaining data  ------------------------- */
-        var _a = Sgp4.days2mdhms(year, days + jdfrac), mon = _a.mon, day = _a.day, hr = _a.hr, minute = _a.minute, sec = _a.sec;
+        var _a = Sgp4.days2mdhms(year, days + jdfrac), mon = _a.mon, day = _a.day, hr = _a.hr, min = _a.min, sec = _a.sec;
         return {
             year: year,
             mon: mon,
             day: day,
             hr: hr,
-            minute: minute,
+            min: min,
             sec: sec,
         };
-    }; // invjday
+    }; // Invjday
     return Sgp4;
 }());
 export { Sgp4 };