ootk 1.0.4 → 1.0.8

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

@@ -30,6 +30,7 @@
  * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
  * SOFTWARE.
  */
+import { DEG2RAD, PI, TAU, temp4, x2o3 } from './utils/constants';
 /*     ----------------------------------------------------------------
  *
  *                               sgp4unit.cpp
@@ -86,47 +87,6 @@
  *                     80  norad
  *                           original baseline
  *       ----------------------------------------------------------------      */
-var PI = Math.PI;
-var TAU = PI * 2;
-var DEG2RAD = PI / 180.0;
-var x2o3 = 2.0 / 3.0;
-var xpdotp = 1440.0 / (2.0 * PI); // 229.1831180523293;
-var zes = 0.01675;
-var zel = 0.0549;
-var c1ss = 2.9864797e-6;
-var c1l = 4.7968065e-7;
-var zsinis = 0.39785416;
-var zcosis = 0.91744867;
-var zcosgs = 0.1945905;
-var zsings = -0.98088458;
-var q22 = 1.7891679e-6;
-var q31 = 2.1460748e-6;
-var q33 = 2.2123015e-7;
-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
-var root32 = 3.7393792e-7;
-var root52 = 1.1428639e-7;
-var znl = 1.5835218e-4;
-var zns = 1.19459e-5;
-var temp4 = 1.5e-12;
-var fasx2 = 0.13130908;
-var fasx4 = 2.8843198;
-var fasx6 = 0.37448087;
-var g22 = 5.7686396;
-var g32 = 0.95240898;
-var g44 = 1.8014998;
-var g52 = 1.050833;
-var g54 = 4.4108898;
-var stepp = 720.0;
-var stepn = -720.0;
-var step2 = 259200.0;
-/**
- * Ootk code is run synchronously and often on super tight loops, so lets reuse the variables as much as possible.
- */
-// prettier-ignore
-var axnl, aynl, xl, u, ktr, vkmpersec, ecose, esine, el2, pl, rl, rdotl, rvdotl, betal, sinu, cosu, sin2u, coseo1, sineo1, cosip, sinip, cosisq, delm, delomg, eo1, argpm, argpp, su, t3, t4, tc, tem5, temp, tempa, tempe, templ, inclm, mm, nm, nodem, xincp, xlm, mp, nodep, xmdf, argpdf, nodedf, t2, delmtemp, em, dspaceOptions, dspaceResult, am, ep, cosim, sinim, dpperResult, dpperParameters, method, ts70, ds70, c1, tfrac, thgr70, fk5r, c1p2p, dndt, ft, delt, xndt, xldot, xnddt, x2omi, xomi, x2li, theta, tut1, snodm, cnodm, sinomm, cosomm, betasq, rtemsq, peo, a1, a2, a3, a4, a5, a6, a7, a8, a9, a10, cc, x1, x2, x3, x4, x5, x6, x7, x8, zcosg, zsing, zcosh, zsinh, zcosi, zsini, ss1, ss2, ss3, ss4, ss5, ss6, ss7, sz1, sz2, sz3, sz11, sz12, sz13, sz21, sz22, sz23, sz31, sz32, sz33, s1, s2, s3, s4, s5, s6, s7, z1, z2, z3, z11, z12, z13, z21, z22, z23, z31, z32, z33, 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, temp1, xno2, ainv2, aonv, cos2u, mrt, temp2, xnode, xinc, mvt, rvdot, sinsu, cossu, snod, cnod, cosi, sini, xmx, xmy, ux, uy, uz, vx, vy, vz, r, v, year, satrec, dsinitResult, dsinitOptions, ao, mon, day, hr, minute, sec, ss, qzms24temp, qzms2ttemp, initlOptions, initlResult, cosio, delmotemp, dscomOptions, dscomResult, dpperOptions, sinio, qzms2t, con42, cosio2, lmonth, dayofyr, inttemp, i, eccsq, omeosq, rteosq, ak, d1, adel, po, ses, sis, sls, sghs, shs, sel, sil, sll, con41, ainv, posq, rp, sghl, shll, PInco, plo, pho, ctem, zcosil, zsinhl, stem, zsingl, zsinil, gam, zy, pgho, zcoshl, xnoi, zcosgl, zmol, zmos, se2, se3, si2, si3, sl2, sl3, sl4, sgh2, xnodce, sh3, sgh4, sh2, ee2, e3, xi2, xi3, xl2, xl3, xl4, sgh3, xgh2, xgh3, xgh4, xh2, xh3, sgs, emo, emsqo, eoc, alfdp, betdp, cosop, sinop, dalf, dbet, dls, f2, f3, pe, pgh, ph, PInc, sinzf, xls, xnoh, zf, zm, cc1sq, cc2, cc3, coef, coef1, cosio4, emsq, eeta, etasq, perige, PInvsq, psisq, qzms24, sfour, temp3, tsi, xPIdot, xhdot1, leadingChar, gsto, lsflg, j2, j3, j3oj2, j4, xke, mus, radiusearthkm, tumin;
 var Sgp4 = /** @class */ (function () {
     function Sgp4() {
     }
@@ -199,28 +159,28 @@ var Sgp4 = /** @class */ (function () {
     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;
-        var _a = options.opsmode, opsmode = _a === void 0 ? i : _a, init = options.init;
+        var _a = options.opsmode, opsmode = _a === void 0 ? 'i' : _a, init = options.init;
         //  ---------------------- constants -----------------------------
-        /** Ootk -- Declared outside the class at the top */
-        // zns = 1.19459e-5;
-        // zes = 0.01675;
-        // znl = 1.5835218e-4;
-        // zel = 0.0549;
+        /** Ootk -- Some variables imported from outside the class at the top */
+        var zns = 1.19459e-5;
+        var zes = 0.01675;
+        var znl = 1.5835218e-4;
+        var zel = 0.0549;
         //  --------------- calculate time varying periodics -----------
-        zm = zmos + zns * t;
+        var zm = zmos + zns * t;
         // be sure that the initial call has time set to zero
         if (init === 'y') {
             zm = zmos;
         }
-        zf = zm + 2.0 * zes * Math.sin(zm);
-        sinzf = Math.sin(zf);
-        f2 = 0.5 * sinzf * sinzf - 0.25;
-        f3 = -0.5 * sinzf * Math.cos(zf);
-        ses = se2 * f2 + se3 * f3;
-        sis = si2 * f2 + si3 * f3;
-        sls = sl2 * f2 + sl3 * f3 + sl4 * sinzf;
-        sghs = sgh2 * f2 + sgh3 * f3 + sgh4 * sinzf;
-        shs = sh2 * f2 + sh3 * f3;
+        var zf = zm + 2.0 * zes * Math.sin(zm);
+        var sinzf = Math.sin(zf);
+        var f2 = 0.5 * sinzf * sinzf - 0.25;
+        var f3 = -0.5 * sinzf * Math.cos(zf);
+        var ses = se2 * f2 + se3 * f3;
+        var sis = si2 * f2 + si3 * f3;
+        var sls = sl2 * f2 + sl3 * f3 + sl4 * sinzf;
+        var sghs = sgh2 * f2 + sgh3 * f3 + sgh4 * sinzf;
+        var shs = sh2 * f2 + sh3 * f3;
         zm = zmol + znl * t;
         if (init === 'y') {
             zm = zmol;
@@ -229,16 +189,16 @@ var Sgp4 = /** @class */ (function () {
         sinzf = Math.sin(zf);
         f2 = 0.5 * sinzf * sinzf - 0.25;
         f3 = -0.5 * sinzf * Math.cos(zf);
-        sel = ee2 * f2 + e3 * f3;
-        sil = xi2 * f2 + xi3 * f3;
-        sll = xl2 * f2 + xl3 * f3 + xl4 * sinzf;
-        sghl = xgh2 * f2 + xgh3 * f3 + xgh4 * sinzf;
-        shll = xh2 * f2 + xh3 * f3;
-        pe = ses + sel;
-        PInc = sis + sil;
-        pl = sls + sll;
-        pgh = sghs + sghl;
-        ph = shs + shll;
+        var sel = ee2 * f2 + e3 * f3;
+        var sil = xi2 * f2 + xi3 * f3;
+        var sll = xl2 * f2 + xl3 * f3 + xl4 * sinzf;
+        var sghl = xgh2 * f2 + xgh3 * f3 + xgh4 * sinzf;
+        var shll = xh2 * f2 + xh3 * f3;
+        var pe = ses + sel;
+        var PInc = sis + sil;
+        var pl = sls + sll;
+        var pgh = sghs + sghl;
+        var ph = shs + shll;
         if (init === 'n') {
             pe -= peo;
             PInc -= PInco;
@@ -247,8 +207,8 @@ var Sgp4 = /** @class */ (function () {
             ph -= pho;
             inclp += PInc;
             ep += pe;
-            sinip = Math.sin(inclp);
-            cosip = Math.cos(inclp);
+            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
@@ -267,12 +227,12 @@ var Sgp4 = /** @class */ (function () {
             }
             else {
                 //  ---- apply periodics with lyddane modification ----
-                sinop = Math.sin(nodep);
-                cosop = Math.cos(nodep);
-                alfdp = sinip * sinop;
-                betdp = sinip * cosop;
-                dalf = ph * cosop + PInc * cosip * sinop;
-                dbet = -ph * sinop + PInc * cosip * cosop;
+                var sinop = Math.sin(nodep);
+                var cosop = Math.cos(nodep);
+                var alfdp = sinip * sinop;
+                var betdp = sinip * cosop;
+                var dalf = ph * cosop + PInc * cosip * sinop;
+                var dbet = -ph * sinop + PInc * cosip * cosop;
                 alfdp += dalf;
                 betdp += dbet;
                 nodep %= TAU;
@@ -282,10 +242,10 @@ var Sgp4 = /** @class */ (function () {
                 if (nodep < 0.0 && opsmode === 'a') {
                     nodep += TAU;
                 }
-                xls = mp + argpp + cosip * nodep;
-                dls = pl + pgh - PInc * nodep * sinip;
+                var xls = mp + argpp + cosip * nodep;
+                var dls = pl + pgh - PInc * nodep * sinip;
                 xls += dls;
-                xnoh = nodep;
+                var xnoh = nodep;
                 nodep = Math.atan2(alfdp, betdp);
                 //  sgp4fix for afspc written intrinsic functions
                 //  nodep used without a trigonometric function ahead
@@ -384,77 +344,78 @@ var Sgp4 = /** @class */ (function () {
     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 -------------------------
-        /** Ootk -- These are declared outside the class at the top of the file */
-        // zes = 0.01675;
-        // zel = 0.0549;
-        // c1ss = 2.9864797e-6;
-        // c1l = 4.7968065e-7;
-        // zsinis = 0.39785416;
-        // zcosis = 0.91744867;
-        // zcosgs = 0.1945905;
-        // zsings = -0.98088458;
-        // TAU = 2.0 * Math.PI;
+        /** Ootk -- Some variables imported from outside the class at the top */
+        var zes = 0.01675;
+        var zel = 0.0549;
+        var c1ss = 2.9864797e-6;
+        var c1l = 4.7968065e-7;
+        var zsinis = 0.39785416;
+        var zcosis = 0.91744867;
+        var zcosgs = 0.1945905;
+        var zsings = -0.98088458;
+        var TAU = 2.0 * Math.PI;
         //  --------------------- local variables ------------------------
-        nm = np;
-        em = ep;
-        snodm = Math.sin(nodep);
-        cnodm = Math.cos(nodep);
-        sinomm = Math.sin(argpp);
-        cosomm = Math.cos(argpp);
-        sinim = Math.sin(inclp);
-        cosim = Math.cos(inclp);
-        emsq = em * em;
-        betasq = 1.0 - emsq;
-        rtemsq = Math.sqrt(betasq);
+        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 nm = np;
+        var em = ep;
+        var snodm = Math.sin(nodep);
+        var cnodm = Math.cos(nodep);
+        var sinomm = Math.sin(argpp);
+        var cosomm = Math.cos(argpp);
+        var sinim = Math.sin(inclp);
+        var cosim = Math.cos(inclp);
+        var emsq = em * em;
+        var betasq = 1.0 - emsq;
+        var rtemsq = Math.sqrt(betasq);
         //  ----------------- initialize lunar solar terms ---------------
-        peo = 0.0;
-        PInco = 0.0;
-        plo = 0.0;
-        pgho = 0.0;
-        pho = 0.0;
-        day = epoch + 18261.5 + tc / 1440.0;
-        xnodce = (4.523602 - 9.2422029e-4 * day) % TAU;
-        stem = Math.sin(xnodce);
-        ctem = Math.cos(xnodce);
-        zcosil = 0.91375164 - 0.03568096 * ctem;
-        zsinil = Math.sqrt(1.0 - zcosil * zcosil);
-        zsinhl = (0.089683511 * stem) / zsinil;
-        zcoshl = Math.sqrt(1.0 - zsinhl * zsinhl);
-        gam = 5.8351514 + 0.001944368 * day;
+        var peo = 0.0;
+        var PInco = 0.0;
+        var plo = 0.0;
+        var pgho = 0.0;
+        var pho = 0.0;
+        var day = epoch + 18261.5 + tc / 1440.0;
+        var xnodce = (4.523602 - 9.2422029e-4 * day) % TAU;
+        var stem = Math.sin(xnodce);
+        var ctem = Math.cos(xnodce);
+        var zcosil = 0.91375164 - 0.03568096 * ctem;
+        var zsinil = Math.sqrt(1.0 - zcosil * zcosil);
+        var zsinhl = (0.089683511 * stem) / zsinil;
+        var zcoshl = Math.sqrt(1.0 - zsinhl * zsinhl);
+        var gam = 5.8351514 + 0.001944368 * day;
         var zx = (0.39785416 * stem) / zsinil;
-        zy = zcoshl * ctem + 0.91744867 * zsinhl * stem;
+        var zy = zcoshl * ctem + 0.91744867 * zsinhl * stem;
         zx = Math.atan2(zx, zy);
         zx += gam - xnodce;
-        zcosgl = Math.cos(zx);
-        zsingl = Math.sin(zx);
+        var zcosgl = Math.cos(zx);
+        var zsingl = Math.sin(zx);
         //  ------------------------- do solar terms ---------------------
-        zcosg = zcosgs;
-        zsing = zsings;
-        zcosi = zcosis;
-        zsini = zsinis;
-        zcosh = cnodm;
-        zsinh = snodm;
-        cc = c1ss;
-        xnoi = 1.0 / nm;
-        for (lsflg = 1; lsflg <= 2; lsflg++) {
-            a1 = zcosg * zcosh + zsing * zcosi * zsinh;
-            a3 = -zsing * zcosh + zcosg * zcosi * zsinh;
-            a7 = -zcosg * zsinh + zsing * zcosi * zcosh;
-            a8 = zsing * zsini;
-            a9 = zsing * zsinh + zcosg * zcosi * zcosh;
-            a10 = zcosg * zsini;
-            a2 = cosim * a7 + sinim * a8;
-            a4 = cosim * a9 + sinim * a10;
-            a5 = -sinim * a7 + cosim * a8;
-            a6 = -sinim * a9 + cosim * a10;
-            x1 = a1 * cosomm + a2 * sinomm;
-            x2 = a3 * cosomm + a4 * sinomm;
-            x3 = -a1 * sinomm + a2 * cosomm;
-            x4 = -a3 * sinomm + a4 * cosomm;
-            x5 = a5 * sinomm;
-            x6 = a6 * sinomm;
-            x7 = a5 * cosomm;
-            x8 = a6 * cosomm;
+        var zcosg = zcosgs;
+        var zsing = zsings;
+        var zcosi = zcosis;
+        var zsini = zsinis;
+        var zcosh = cnodm;
+        var zsinh = snodm;
+        var cc = c1ss;
+        var xnoi = 1.0 / nm;
+        for (var lsflg = 1; lsflg <= 2; lsflg++) {
+            var a1 = zcosg * zcosh + zsing * zcosi * zsinh;
+            var a3 = -zsing * zcosh + zcosg * zcosi * zsinh;
+            var a7 = -zcosg * zsinh + zsing * zcosi * zcosh;
+            var a8 = zsing * zsini;
+            var a9 = zsing * zsinh + zcosg * zcosi * zcosh;
+            var a10 = zcosg * zsini;
+            var a2 = cosim * a7 + sinim * a8;
+            var a4 = cosim * a9 + sinim * a10;
+            var a5 = -sinim * a7 + cosim * a8;
+            var a6 = -sinim * a9 + cosim * a10;
+            var x1 = a1 * cosomm + a2 * sinomm;
+            var x2 = a3 * cosomm + a4 * sinomm;
+            var x3 = -a1 * sinomm + a2 * cosomm;
+            var x4 = -a3 * sinomm + a4 * cosomm;
+            var x5 = a5 * sinomm;
+            var x6 = a6 * sinomm;
+            var x7 = a5 * cosomm;
+            var x8 = a6 * cosomm;
             z31 = 12.0 * x1 * x1 - 3.0 * x3 * x3;
             z32 = 24.0 * x1 * x2 - 6.0 * x3 * x4;
             z33 = 12.0 * x2 * x2 - 3.0 * x4 * x4;
@@ -462,13 +423,10 @@ var Sgp4 = /** @class */ (function () {
             z2 = 6.0 * (a1 * a3 + a2 * a4) + z32 * emsq;
             z3 = 3.0 * (a3 * a3 + a4 * a4) + z33 * emsq;
             z11 = -6.0 * a1 * a5 + emsq * (-24.0 * x1 * x7 - 6.0 * x3 * x5);
-            z12 =
-                -6.0 * (a1 * a6 + a3 * a5) +
-                    emsq * (-24.0 * (x2 * x7 + x1 * x8) + -6.0 * (x3 * x6 + x4 * x5));
+            z12 = -6.0 * (a1 * a6 + a3 * a5) + emsq * (-24.0 * (x2 * x7 + x1 * x8) + -6.0 * (x3 * x6 + x4 * x5));
             z13 = -6.0 * a3 * a6 + emsq * (-24.0 * x2 * x8 - 6.0 * x4 * x6);
             z21 = 6.0 * a2 * a5 + emsq * (24.0 * x1 * x5 - 6.0 * x3 * x7);
-            z22 =
-                6.0 * (a4 * a5 + a2 * a6) + emsq * (24.0 * (x2 * x5 + x1 * x6) - 6.0 * (x4 * x7 + x3 * x8));
+            z22 = 6.0 * (a4 * a5 + a2 * a6) + emsq * (24.0 * (x2 * x5 + x1 * x6) - 6.0 * (x4 * x7 + x3 * x8));
             z23 = 6.0 * a4 * a6 + emsq * (24.0 * x2 * x6 - 6.0 * x4 * x8);
             z1 = z1 + z1 + betasq * z31;
             z2 = z2 + z2 + betasq * z32;
@@ -510,34 +468,34 @@ var Sgp4 = /** @class */ (function () {
                 cc = c1l;
             }
         }
-        zmol = (4.7199672 + (0.2299715 * day - gam)) % TAU;
-        zmos = (6.2565837 + 0.017201977 * day) % TAU;
+        var zmol = (4.7199672 + (0.2299715 * day - gam)) % TAU;
+        var zmos = (6.2565837 + 0.017201977 * day) % TAU;
         //  ------------------------ do solar terms ----------------------
-        se2 = 2.0 * ss1 * ss6;
-        se3 = 2.0 * ss1 * ss7;
-        si2 = 2.0 * ss2 * sz12;
-        si3 = 2.0 * ss2 * (sz13 - sz11);
-        sl2 = -2.0 * ss3 * sz2;
-        sl3 = -2.0 * ss3 * (sz3 - sz1);
-        sl4 = -2.0 * ss3 * (-21.0 - 9.0 * emsq) * zes;
-        sgh2 = 2.0 * ss4 * sz32;
-        sgh3 = 2.0 * ss4 * (sz33 - sz31);
-        sgh4 = -18.0 * ss4 * zes;
-        sh2 = -2.0 * ss2 * sz22;
-        sh3 = -2.0 * ss2 * (sz23 - sz21);
+        var se2 = 2.0 * ss1 * ss6;
+        var se3 = 2.0 * ss1 * ss7;
+        var si2 = 2.0 * ss2 * sz12;
+        var si3 = 2.0 * ss2 * (sz13 - sz11);
+        var sl2 = -2.0 * ss3 * sz2;
+        var sl3 = -2.0 * ss3 * (sz3 - sz1);
+        var sl4 = -2.0 * ss3 * (-21.0 - 9.0 * emsq) * zes;
+        var sgh2 = 2.0 * ss4 * sz32;
+        var sgh3 = 2.0 * ss4 * (sz33 - sz31);
+        var sgh4 = -18.0 * ss4 * zes;
+        var sh2 = -2.0 * ss2 * sz22;
+        var sh3 = -2.0 * ss2 * (sz23 - sz21);
         //  ------------------------ do lunar terms ----------------------
-        ee2 = 2.0 * s1 * s6;
-        e3 = 2.0 * s1 * s7;
-        xi2 = 2.0 * s2 * z12;
-        xi3 = 2.0 * s2 * (z13 - z11);
-        xl2 = -2.0 * s3 * z2;
-        xl3 = -2.0 * s3 * (z3 - z1);
-        xl4 = -2.0 * s3 * (-21.0 - 9.0 * emsq) * zel;
-        xgh2 = 2.0 * s4 * z32;
-        xgh3 = 2.0 * s4 * (z33 - z31);
-        xgh4 = -18.0 * s4 * zel;
-        xh2 = -2.0 * s2 * z22;
-        xh3 = -2.0 * s2 * (z23 - z21);
+        var ee2 = 2.0 * s1 * s6;
+        var e3 = 2.0 * s1 * s7;
+        var xi2 = 2.0 * s2 * z12;
+        var xi3 = 2.0 * s2 * (z13 - z11);
+        var xl2 = -2.0 * s3 * z2;
+        var xl3 = -2.0 * s3 * (z3 - z1);
+        var xl4 = -2.0 * s3 * (-21.0 - 9.0 * emsq) * zel;
+        var xgh2 = 2.0 * s4 * z32;
+        var xgh3 = 2.0 * s4 * (z33 - z31);
+        var xgh4 = -18.0 * s4 * zel;
+        var xh2 = -2.0 * s2 * z22;
+        var xh3 = -2.0 * s2 * (z23 - z21);
         return {
             snodm: snodm,
             cnodm: cnodm,
@@ -708,19 +666,19 @@ var Sgp4 = /** @class */ (function () {
         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 -- these are declared at the top instead */
-        // q22 = 1.7891679e-6;
-        // q31 = 2.1460748e-6;
-        // q33 = 2.2123015e-7;
-        // root22 = 1.7891679e-6;
-        // root44 = 7.3636953e-9;
-        // root54 = 2.1765803e-9;
-        // rptim = 4.37526908801129966e-3; // equates to 7.29211514668855e-5 rad/sec
-        // root32 = 3.7393792e-7;
-        // root52 = 1.1428639e-7;
-        // x2o3 = 2.0 / 3.0;
-        // znl = 1.5835218e-4;
-        // zns = 1.19459e-5;
+        /** Ootk -- Some variables imported from outside the class at the top */
+        var q22 = 1.7891679e-6;
+        var q31 = 2.1460748e-6;
+        var q33 = 2.2123015e-7;
+        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
+        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 );
@@ -733,10 +691,10 @@ var Sgp4 = /** @class */ (function () {
             irez = 2;
         }
         // ------------------------ do solar terms -------------------
-        ses = ss1 * zns * ss5;
-        sis = ss2 * zns * (sz11 + sz13);
-        sls = -zns * ss3 * (sz1 + sz3 - 14.0 - 6.0 * emsq);
-        sghs = ss4 * zns * (sz31 + sz33 - 6.0);
+        var ses = ss1 * zns * ss5;
+        var sis = ss2 * zns * (sz11 + sz13);
+        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
         if (inclm < 5.2359877e-2 || inclm > PI - 5.2359877e-2) {
@@ -745,12 +703,12 @@ var Sgp4 = /** @class */ (function () {
         if (sinim !== 0.0) {
             shs /= sinim;
         }
-        sgs = sghs - cosim * shs;
+        var sgs = sghs - cosim * shs;
         // ------------------------- do lunar terms ------------------
         dedt = ses + s1 * znl * s5;
         didt = sis + s2 * znl * (z11 + z13);
         dmdt = sls - znl * s3 * (z1 + z3 - 14.0 - 6.0 * emsq);
-        sghl = s4 * znl * (z31 + z33 - 6.0);
+        var sghl = s4 * znl * (z31 + z33 - 6.0);
         var shll = -znl * s2 * (z21 + z23);
         // sgp4fix for 180 deg incl
         if (inclm < 5.2359877e-2 || inclm > PI - 5.2359877e-2) {
@@ -763,8 +721,8 @@ var Sgp4 = /** @class */ (function () {
             dnodt += shll / sinim;
         }
         // ----------- calculate deep space resonance effects --------
-        dndt = 0.0;
-        theta = (gsto + tc * rptim) % TAU;
+        var dndt = 0.0;
+        var theta = (gsto + tc * rptim) % TAU;
         em += dedt * t;
         inclm += didt * t;
         argpm += domdt * t;
@@ -780,16 +738,17 @@ var Sgp4 = /** @class */ (function () {
         // }
         // -------------- initialize the resonance terms -------------
         if (irez !== 0) {
-            aonv = Math.pow((nm / xke), x2o3);
+            var aonv = Math.pow((nm / xke), x2o3);
             // ---------- geopotential resonance for 12 hour orbits ------
             if (irez === 2) {
-                cosisq = cosim * cosim;
-                emo = em;
+                var cosisq = cosim * cosim;
+                var emo = em;
                 em = ecco;
-                emsqo = emsq;
+                var emsqo = emsq;
                 emsq = eccsq;
-                eoc = em * emsq;
-                g201 = -0.306 - (em - 0.64) * 0.44;
+                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;
                 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;
@@ -821,30 +780,24 @@ var Sgp4 = /** @class */ (function () {
                     g521 = -51752.104 + 218913.95 * em - 309468.16 * emsq + 146349.42 * eoc;
                     g532 = -40023.88 + 170470.89 * em - 242699.48 * emsq + 115605.82 * eoc;
                 }
-                sini2 = sinim * sinim;
-                f220 = 0.75 * (1.0 + 2.0 * cosim + cosisq);
-                f221 = 1.5 * sini2;
-                f321 = 1.875 * sinim * (1.0 - 2.0 * cosim - 3.0 * cosisq);
-                f322 = -1.875 * sinim * (1.0 + 2.0 * cosim - 3.0 * cosisq);
-                f441 = 35.0 * sini2 * f220;
-                f442 = 39.375 * sini2 * sini2;
-                f522 =
-                    9.84375 *
-                        sinim *
-                        (sini2 * (1.0 - 2.0 * cosim - 5.0 * cosisq) +
-                            0.33333333 * (-2.0 + 4.0 * cosim + 6.0 * cosisq));
-                f523 =
+                var sini2 = sinim * sinim;
+                var f220 = 0.75 * (1.0 + 2.0 * cosim + cosisq);
+                var f221 = 1.5 * sini2;
+                var f321 = 1.875 * sinim * (1.0 - 2.0 * cosim - 3.0 * cosisq);
+                var f322 = -1.875 * sinim * (1.0 + 2.0 * cosim - 3.0 * cosisq);
+                var f441 = 35.0 * sini2 * f220;
+                var f442 = 39.375 * sini2 * sini2;
+                var f522 = 9.84375 *
                     sinim *
-                        (4.92187512 * sini2 * (-2.0 - 4.0 * cosim + 10.0 * cosisq) +
-                            6.56250012 * (1.0 + 2.0 * cosim - 3.0 * cosisq));
-                f542 =
-                    29.53125 * sinim * (2.0 - 8.0 * cosim + cosisq * (-12.0 + 8.0 * cosim + 10.0 * cosisq));
-                f543 =
-                    29.53125 * sinim * (-2.0 - 8.0 * cosim + cosisq * (12.0 + 8.0 * cosim - 10.0 * cosisq));
-                xno2 = nm * nm;
-                ainv2 = aonv * aonv;
-                temp1 = 3.0 * xno2 * ainv2;
-                temp = temp1 * root22;
+                    (sini2 * (1.0 - 2.0 * cosim - 5.0 * cosisq) + 0.33333333 * (-2.0 + 4.0 * cosim + 6.0 * cosisq));
+                var f523 = sinim *
+                    (4.92187512 * sini2 * (-2.0 - 4.0 * cosim + 10.0 * cosisq) + 6.56250012 * (1.0 + 2.0 * cosim - 3.0 * cosisq));
+                var f542 = 29.53125 * sinim * (2.0 - 8.0 * cosim + cosisq * (-12.0 + 8.0 * cosim + 10.0 * cosisq));
+                var f543 = 29.53125 * sinim * (-2.0 - 8.0 * cosim + cosisq * (12.0 + 8.0 * cosim - 10.0 * cosisq));
+                var xno2 = nm * nm;
+                var ainv2 = aonv * aonv;
+                var temp1 = 3.0 * xno2 * ainv2;
+                var temp = temp1 * root22;
                 d2201 = temp * f220 * g201;
                 d2211 = temp * f221 * g211;
                 temp1 *= aonv;
@@ -869,12 +822,12 @@ var Sgp4 = /** @class */ (function () {
             }
             //  ---------------- synchronous resonance terms --------------
             if (irez === 1) {
-                g200 = 1.0 + emsq * (-2.5 + 0.8125 * emsq);
-                g310 = 1.0 + 2.0 * emsq;
-                g300 = 1.0 + emsq * (-6.0 + 6.60937 * emsq);
-                f220 = 0.75 * (1.0 + cosim) * (1.0 + cosim);
-                f311 = 0.9375 * sinim * sinim * (1.0 + 3.0 * cosim) - 0.75 * (1.0 + cosim);
-                f330 = 1.0 + cosim;
+                var g200 = 1.0 + emsq * (-2.5 + 0.8125 * emsq);
+                var g310 = 1.0 + 2.0 * emsq;
+                var g300 = 1.0 + emsq * (-6.0 + 6.60937 * emsq);
+                var f220 = 0.75 * (1.0 + cosim) * (1.0 + cosim);
+                var f311 = 0.9375 * sinim * sinim * (1.0 + 3.0 * cosim) - 0.75 * (1.0 + cosim);
+                var f330 = 1.0 + cosim;
                 f330 *= 1.875 * f330 * f330;
                 del1 = 3.0 * nm * nm * aonv * aonv;
                 del2 = 2.0 * del1 * f220 * g200 * q22;
@@ -995,25 +948,23 @@ var Sgp4 = /** @class */ (function () {
    *    hoots, schumacher and glover 2004
    *    vallado, crawford, hujsak, kelso  2006
    ----------------------------------------------------------------------------*/
-    Sgp4.dspace = function (options) {
-        var irez = options.irez, 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, del1 = options.del1, del2 = options.del2, del3 = options.del3, didt = options.didt, dmdt = options.dmdt, dnodt = options.dnodt, domdt = options.domdt, argpo = options.argpo, argpdot = options.argpdot, t = options.t, tc = options.tc, gsto = options.gsto, xfact = options.xfact, xlamo = options.xlamo, no = options.no;
-        var atime = options.atime, em = options.em, argpm = options.argpm, inclm = options.inclm, xli = options.xli, mm = options.mm, xni = options.xni, nodem = options.nodem, nm = options.nm;
-        /** Ootk -- Declared at the top of the file instead */
-        // fasx2 = 0.13130908;
-        // fasx4 = 2.8843198;
-        // fasx6 = 0.37448087;
-        // g22 = 5.7686396;
-        // g32 = 0.95240898;
-        // g44 = 1.8014998;
-        // g52 = 1.050833;
-        // g54 = 4.4108898;
-        // rptim = 4.37526908801129966e-3; // equates to 7.29211514668855e-5 rad/sec
-        // stepp = 720.0;
-        // stepn = -720.0;
-        // step2 = 259200.0;
+    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;
+        var fasx4 = 2.8843198;
+        var fasx6 = 0.37448087;
+        var g22 = 5.7686396;
+        var g32 = 0.95240898;
+        var g44 = 1.8014998;
+        var g52 = 1.050833;
+        var g54 = 4.4108898;
+        var rptim = 4.37526908801129966e-3; // equates to 7.29211514668855e-5 rad/sec
+        var stepp = 720.0;
+        var stepn = -720.0;
+        var step2 = 259200.0;
         //  ----------- calculate deep space resonance effects -----------
-        dndt = 0.0;
-        theta = (gsto + tc * rptim) % TAU;
+        var dndt = 0.0;
+        var theta = (gsto + tc * rptim) % TAU;
         em += dedt * t;
         inclm += didt * t;
         argpm += domdt * t;
@@ -1042,21 +993,21 @@ var Sgp4 = /** @class */ (function () {
                 xli = xlamo;
             }
             // sgp4fix move check outside loop
+            var delt = void 0;
             if (t > 0.0) {
                 delt = stepp;
             }
             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
             while (iretn === 381) {
                 //  ------------------- 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));
+                        del1 * Math.sin(xli - fasx2) + del2 * Math.sin(2.0 * (xli - fasx4)) + del3 * Math.sin(3.0 * (xli - fasx6));
                     xldot = xni + xfact;
                     xnddt =
                         del1 * Math.cos(xli - fasx2) +
@@ -1111,7 +1062,7 @@ var Sgp4 = /** @class */ (function () {
                 }
             } // while iretn = 381
             nm = xni + xndt * ft + xnddt * ft * ft * 0.5;
-            xl = xli + xldot * ft + xndt * ft * ft * 0.5;
+            var xl = xli + xldot * ft + xndt * ft * ft * 0.5;
             if (irez !== 1) {
                 mm = xl - 2.0 * nodem + 2.0 * theta;
                 dndt = nm - no;
@@ -1122,18 +1073,7 @@ var Sgp4 = /** @class */ (function () {
             }
             nm = no + dndt;
         }
-        return {
-            atime: atime,
-            em: em,
-            argpm: argpm,
-            inclm: inclm,
-            xli: xli,
-            mm: mm,
-            xni: xni,
-            nodem: nodem,
-            dndt: dndt,
-            nm: nm,
-        };
+        return [em, argpm, inclm, mm, nodem, nm];
     }; // dspace
     /*-----------------------------------------------------------------------------
    *
@@ -1195,50 +1135,51 @@ var Sgp4 = /** @class */ (function () {
         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 ------------------------ */
-        /** Ootk -- Declared at the top of the file */
+        var PI = Math.PI;
+        var TAU = PI * 2;
+        var x2o3 = 2.0 / 3.0;
         // sgp4fix use old way of finding gst
-        /** Ootk -- Declared at the top of the file */
+        /** 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 )
         // ------------- calculate auxillary epoch quantities ----------
-        eccsq = ecco * ecco;
-        omeosq = 1.0 - eccsq;
-        rteosq = Math.sqrt(omeosq);
-        cosio = Math.cos(inclo);
-        cosio2 = cosio * cosio;
+        var eccsq = ecco * ecco;
+        var omeosq = 1.0 - eccsq;
+        var rteosq = Math.sqrt(omeosq);
+        var cosio = Math.cos(inclo);
+        var cosio2 = cosio * cosio;
         // ------------------ un-kozai the mean motion -----------------
-        ak = Math.pow((xke / no), x2o3);
-        d1 = (0.75 * j2 * (3.0 * cosio2 - 1.0)) / (rteosq * omeosq);
+        var ak = Math.pow((xke / no), x2o3);
+        var d1 = (0.75 * j2 * (3.0 * cosio2 - 1.0)) / (rteosq * omeosq);
         var delPrime = d1 / (ak * ak);
-        adel =
-            ak *
-                (1.0 - delPrime * delPrime - delPrime * (1.0 / 3.0 + (134.0 * delPrime * delPrime) / 81.0));
+        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);
-        ao = Math.pow((xke / no), x2o3);
-        sinio = Math.sin(inclo);
-        po = ao * omeosq;
-        con42 = 1.0 - 5.0 * cosio2;
-        con41 = -con42 - cosio2 - cosio2;
-        ainv = 1.0 / ao;
-        posq = po * po;
-        rp = ao * (1.0 - ecco);
-        method = 'n';
+        var ao = Math.pow((xke / no), x2o3);
+        var sinio = Math.sin(inclo);
+        var po = ao * omeosq;
+        var con42 = 1.0 - 5.0 * cosio2;
+        var con41 = -con42 - cosio2 - cosio2;
+        var ainv = 1.0 / ao;
+        var posq = po * po;
+        var rp = ao * (1.0 - ecco);
+        var method = 'n';
         //  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
-            ts70 = epoch - 7305.0;
-            ds70 = Math.floor(ts70 + 1.0e-8);
-            tfrac = ts70 - ds70;
+            var ts70 = epoch - 7305.0;
+            var ds70 = Math.floor(ts70 + 1.0e-8);
+            var tfrac = ts70 - ds70;
             //  find greenwich location at epoch
-            c1 = 1.72027916940703639e-2;
-            thgr70 = 1.7321343856509374;
-            fk5r = 5.07551419432269442e-15;
-            c1p2p = c1 + TAU;
+            var c1 = 1.72027916940703639e-2;
+            var thgr70 = 1.7321343856509374;
+            var fk5r = 5.07551419432269442e-15;
+            var c1p2p = c1 + TAU;
             gsto = (thgr70 + c1 * ds70 + c1p2p * tfrac + ts70 * ts70 * fk5r) % TAU;
             /* istanbul ignore next | This is no longer possible*/
             if (gsto < 0.0) {
@@ -1349,14 +1290,11 @@ var Sgp4 = /** @class */ (function () {
    ----------------------------------------------------------------------------*/
     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;
-        /* --------------------- local variables ------------------------ */
-        /** Ootk -- Declared at the top of the file */
+        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
-        // temp4 = 1.5e-12;
         // ----------- set all near earth variables to zero ------------
         satrec.isimp = 0;
         satrec.method = 'n';
@@ -1454,6 +1392,10 @@ var Sgp4 = /** @class */ (function () {
         satrec.j3 = gravResults.j3;
         satrec.j4 = gravResults.j4;
         satrec.j3oj2 = gravResults.j3oj2;
+        var j2 = gravResults.j2;
+        var j4 = gravResults.j4;
+        var xke = gravResults.xke;
+        var j3oj2 = gravResults.j3oj2;
         //-------------------------------------------------------------------------
         satrec.error = 0;
         satrec.operationmode = opsmode;
@@ -1462,8 +1404,8 @@ var Sgp4 = /** @class */ (function () {
          * Ex. A2525 = 102525
          * Ex. Z1234 = 351234
          */
-        leadingChar = satn.split('')[0].toLowerCase(); // Using uppercase will break the -96 math.
-        if (isNaN(leadingChar)) {
+        var leadingChar = satn.split('')[0].toLowerCase(); // Using uppercase will break the -96 math.
+        if (isNaN(parseInt(leadingChar))) {
             satrec.satnum = parseInt(leadingChar.charCodeAt(0) - 96 + 9 + satrec.satnum.slice(1, 5)).toString();
         }
         else {
@@ -1489,14 +1431,14 @@ var Sgp4 = /** @class */ (function () {
         // ------------------------ earth constants -----------------------
         // sgp4fix identify constants and allow alternate values
         // getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
-        ss = 78.0 / satrec.radiusearthkm + 1.0;
+        var ss = 78.0 / satrec.radiusearthkm + 1.0;
         // sgp4fix use multiply for speed instead of pow
-        qzms2ttemp = (120.0 - 78.0) / satrec.radiusearthkm;
-        qzms2t = qzms2ttemp * qzms2ttemp * qzms2ttemp * qzms2ttemp;
+        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
-        initlOptions = {
+        var initlOptions = {
             ecco: satrec.ecco,
             epoch: epoch,
             inclo: satrec.inclo,
@@ -1506,7 +1448,7 @@ var Sgp4 = /** @class */ (function () {
             xke: satrec.xke,
             j2: satrec.j2,
         };
-        initlResult = Sgp4.initl(initlOptions);
+        var initlResult = Sgp4.initl(initlOptions);
         var ao = initlResult.ao, con42 = initlResult.con42, cosio = initlResult.cosio, cosio2 = initlResult.cosio2, eccsq = initlResult.eccsq, omeosq = initlResult.omeosq, posq = initlResult.posq, rp = initlResult.rp, rteosq = initlResult.rteosq, sinio = initlResult.sinio;
         satrec.no = initlResult.no;
         satrec.con41 = initlResult.con41;
@@ -1528,9 +1470,9 @@ var Sgp4 = /** @class */ (function () {
             if (rp < 220.0 / satrec.radiusearthkm + 1.0) {
                 satrec.isimp = 1;
             }
-            sfour = ss;
-            qzms24 = qzms2t;
-            perige = (rp - 1.0) * satrec.radiusearthkm;
+            var sfour = ss;
+            var qzms24 = qzms2t;
+            var perige = (rp - 1.0) * satrec.radiusearthkm;
             // - for perigees below 156 km, s and qoms2t are altered -
             if (perige < 156.0) {
                 sfour = perige - 78.0;
@@ -1538,25 +1480,24 @@ var Sgp4 = /** @class */ (function () {
                     sfour = 20.0;
                 }
                 // sgp4fix use multiply for speed instead of pow
-                qzms24temp = (120.0 - sfour) / satrec.radiusearthkm;
+                var qzms24temp = (120.0 - sfour) / satrec.radiusearthkm;
                 qzms24 = qzms24temp * qzms24temp * qzms24temp * qzms24temp;
                 sfour = sfour / satrec.radiusearthkm + 1.0;
             }
-            PInvsq = 1.0 / posq;
-            tsi = 1.0 / (ao - sfour);
+            var PInvsq = 1.0 / posq;
+            var tsi = 1.0 / (ao - sfour);
             satrec.eta = ao * satrec.ecco * tsi;
-            etasq = satrec.eta * satrec.eta;
-            eeta = satrec.ecco * satrec.eta;
-            psisq = Math.abs(1.0 - etasq);
-            coef = qzms24 * Math.pow(tsi, 4.0);
-            coef1 = coef / Math.pow(psisq, 3.5);
-            cc2 =
-                coef1 *
-                    satrec.no *
-                    (ao * (1.0 + 1.5 * etasq + eeta * (4.0 + etasq)) +
-                        ((0.375 * j2 * tsi) / psisq) * satrec.con41 * (8.0 + 3.0 * etasq * (8.0 + etasq)));
+            var etasq = satrec.eta * satrec.eta;
+            var eeta = satrec.ecco * satrec.eta;
+            var psisq = Math.abs(1.0 - etasq);
+            var coef = qzms24 * Math.pow(tsi, 4.0);
+            var coef1 = coef / Math.pow(psisq, 3.5);
+            var cc2 = coef1 *
+                satrec.no *
+                (ao * (1.0 + 1.5 * etasq + eeta * (4.0 + etasq)) +
+                    ((0.375 * j2 * tsi) / psisq) * satrec.con41 * (8.0 + 3.0 * etasq * (8.0 + etasq)));
             satrec.cc1 = satrec.bstar * cc2;
-            cc3 = 0.0;
+            var cc3 = 0.0;
             if (satrec.ecco > 1.0e-4) {
                 cc3 = (-2.0 * coef * tsi * j3oj2 * satrec.no * sinio) / satrec.ecco;
             }
@@ -1571,15 +1512,12 @@ var Sgp4 = /** @class */ (function () {
                         satrec.ecco * (0.5 + 2.0 * etasq) -
                         ((j2 * tsi) / (ao * psisq)) *
                             (-3.0 * satrec.con41 * (1.0 - 2.0 * eeta + etasq * (1.5 - 0.5 * eeta)) +
-                                0.75 *
-                                    satrec.x1mth2 *
-                                    (2.0 * etasq - eeta * (1.0 + etasq)) *
-                                    Math.cos(2.0 * satrec.argpo)));
+                                0.75 * satrec.x1mth2 * (2.0 * etasq - eeta * (1.0 + etasq)) * Math.cos(2.0 * satrec.argpo)));
             satrec.cc5 = 2.0 * coef1 * ao * omeosq * (1.0 + 2.75 * (etasq + eeta) + eeta * etasq);
-            cosio4 = cosio2 * cosio2;
-            temp1 = 1.5 * j2 * PInvsq * satrec.no;
-            temp2 = 0.5 * temp1 * j2 * PInvsq;
-            temp3 = -0.46875 * j4 * PInvsq * PInvsq * satrec.no;
+            var cosio4 = cosio2 * cosio2;
+            var temp1 = 1.5 * j2 * PInvsq * satrec.no;
+            var temp2 = 0.5 * temp1 * j2 * PInvsq;
+            var temp3 = -0.46875 * j4 * PInvsq * PInvsq * satrec.no;
             satrec.mdot =
                 satrec.no +
                     0.5 * temp1 * rteosq * satrec.con41 +
@@ -1588,10 +1526,9 @@ var Sgp4 = /** @class */ (function () {
                 -0.5 * temp1 * con42 +
                     0.0625 * temp2 * (7.0 - 114.0 * cosio2 + 395.0 * cosio4) +
                     temp3 * (3.0 - 36.0 * cosio2 + 49.0 * cosio4);
-            xhdot1 = -temp1 * cosio;
-            satrec.nodedot =
-                xhdot1 + (0.5 * temp2 * (4.0 - 19.0 * cosio2) + 2.0 * temp3 * (3.0 - 7.0 * cosio2)) * cosio;
-            xPIdot = satrec.argpdot + satrec.nodedot;
+            var xhdot1 = -temp1 * cosio;
+            satrec.nodedot = xhdot1 + (0.5 * temp2 * (4.0 - 19.0 * cosio2) + 2.0 * temp3 * (3.0 - 7.0 * cosio2)) * cosio;
+            var xPIdot = satrec.argpdot + satrec.nodedot;
             satrec.omgcof = satrec.bstar * cc3 * Math.cos(satrec.argpo);
             satrec.xmcof = 0.0;
             if (satrec.ecco > 1.0e-4) {
@@ -1608,7 +1545,7 @@ var Sgp4 = /** @class */ (function () {
             }
             satrec.aycof = -0.5 * j3oj2 * sinio;
             // sgp4fix use multiply for speed instead of pow
-            delmotemp = 1.0 + satrec.eta * Math.cos(satrec.mo);
+            var delmotemp = 1.0 + satrec.eta * Math.cos(satrec.mo);
             satrec.delmo = delmotemp * delmotemp * delmotemp;
             satrec.sinmao = Math.sin(satrec.mo);
             satrec.x7thm1 = 7.0 * cosio2 - 1.0;
@@ -1616,9 +1553,9 @@ var Sgp4 = /** @class */ (function () {
             if (TAU / satrec.no >= 225.0) {
                 satrec.method = 'd';
                 satrec.isimp = 1;
-                tc = 0.0;
-                inclm = satrec.inclo;
-                dscomOptions = {
+                var tc = 0.0;
+                var inclm = satrec.inclo;
+                var dscomOptions = {
                     epoch: epoch,
                     ep: satrec.ecco,
                     argpp: satrec.argpo,
@@ -1658,7 +1595,7 @@ var Sgp4 = /** @class */ (function () {
                     zmol: satrec.zmol,
                     zmos: satrec.zmos,
                 };
-                dscomResult = Sgp4.dscom(dscomOptions);
+                var dscomResult = Sgp4.dscom(dscomOptions);
                 satrec.e3 = dscomResult.e3;
                 satrec.ee2 = dscomResult.ee2;
                 satrec.peo = dscomResult.peo;
@@ -1678,7 +1615,7 @@ var Sgp4 = /** @class */ (function () {
                 satrec.sl2 = dscomResult.sl2;
                 satrec.sl3 = dscomResult.sl3;
                 satrec.sl4 = dscomResult.sl4;
-                var sinim_1 = dscomResult.sinim, cosim_1 = dscomResult.cosim, em_1 = dscomResult.em, emsq_1 = dscomResult.emsq, s1_1 = dscomResult.s1, s2_1 = dscomResult.s2, s3_1 = dscomResult.s3, s4_1 = dscomResult.s4, s5_1 = dscomResult.s5, ss1_1 = dscomResult.ss1, ss2_1 = dscomResult.ss2, ss3_1 = dscomResult.ss3, ss4_1 = dscomResult.ss4, ss5_1 = dscomResult.ss5, sz1_1 = dscomResult.sz1, sz3_1 = dscomResult.sz3, sz11_1 = dscomResult.sz11, sz13_1 = dscomResult.sz13, sz21_1 = dscomResult.sz21, sz23_1 = dscomResult.sz23, sz31_1 = dscomResult.sz31, sz33_1 = dscomResult.sz33;
+                var sinim = dscomResult.sinim, cosim = dscomResult.cosim, em = dscomResult.em, emsq = dscomResult.emsq, s1 = dscomResult.s1, s2 = dscomResult.s2, s3 = dscomResult.s3, s4 = dscomResult.s4, s5 = dscomResult.s5, ss1 = dscomResult.ss1, ss2 = dscomResult.ss2, ss3 = dscomResult.ss3, ss4 = dscomResult.ss4, ss5 = dscomResult.ss5, sz1 = dscomResult.sz1, sz3 = dscomResult.sz3, sz11 = dscomResult.sz11, sz13 = dscomResult.sz13, sz21 = dscomResult.sz21, sz23 = dscomResult.sz23, sz31 = dscomResult.sz31, sz33 = dscomResult.sz33;
                 satrec.xgh2 = dscomResult.xgh2;
                 satrec.xgh3 = dscomResult.xgh3;
                 satrec.xgh4 = dscomResult.xgh4;
@@ -1691,8 +1628,8 @@ var Sgp4 = /** @class */ (function () {
                 satrec.xl4 = dscomResult.xl4;
                 satrec.zmol = dscomResult.zmol;
                 satrec.zmos = dscomResult.zmos;
-                (nm = dscomResult.nm, z1 = dscomResult.z1, z3 = dscomResult.z3, z11 = dscomResult.z11, z13 = dscomResult.z13, z21 = dscomResult.z21, z23 = dscomResult.z23, z31 = dscomResult.z31, z33 = dscomResult.z33);
-                dpperOptions = {
+                var nm = dscomResult.nm, z1 = dscomResult.z1, z3 = dscomResult.z3, z11 = dscomResult.z11, z13 = dscomResult.z13, z21 = dscomResult.z21, z23 = dscomResult.z23, z31 = dscomResult.z31, z33 = dscomResult.z33;
+                var dpperOptions = {
                     inclo: inclm,
                     init: satrec.init,
                     ep: satrec.ecco,
@@ -1702,39 +1639,39 @@ var Sgp4 = /** @class */ (function () {
                     mp: satrec.mo,
                     opsmode: satrec.operationmode,
                 };
-                dpperResult = Sgp4.dpper(satrec, dpperOptions);
+                var dpperResult = Sgp4.dpper(satrec, dpperOptions);
                 satrec.ecco = dpperResult.ep;
                 satrec.inclo = dpperResult.inclp;
                 satrec.nodeo = dpperResult.nodep;
                 satrec.argpo = dpperResult.argpp;
                 satrec.mo = dpperResult.mp;
-                argpm = 0.0;
-                nodem = 0.0;
-                mm = 0.0;
-                dsinitOptions = {
+                var argpm = 0.0;
+                var nodem = 0.0;
+                var mm = 0.0;
+                var dsinitOptions = {
                     xke: xke,
-                    cosim: cosim_1,
-                    emsq: emsq_1,
+                    cosim: cosim,
+                    emsq: emsq,
                     argpo: satrec.argpo,
-                    s1: s1_1,
-                    s2: s2_1,
-                    s3: s3_1,
-                    s4: s4_1,
-                    s5: s5_1,
-                    sinim: sinim_1,
-                    ss1: ss1_1,
-                    ss2: ss2_1,
-                    ss3: ss3_1,
-                    ss4: ss4_1,
-                    ss5: ss5_1,
-                    sz1: sz1_1,
-                    sz3: sz3_1,
-                    sz11: sz11_1,
-                    sz13: sz13_1,
-                    sz21: sz21_1,
-                    sz23: sz23_1,
-                    sz31: sz31_1,
-                    sz33: sz33_1,
+                    s1: s1,
+                    s2: s2,
+                    s3: s3,
+                    s4: s4,
+                    s5: s5,
+                    sinim: sinim,
+                    ss1: ss1,
+                    ss2: ss2,
+                    ss3: ss3,
+                    ss4: ss4,
+                    ss5: ss5,
+                    sz1: sz1,
+                    sz3: sz3,
+                    sz11: sz11,
+                    sz13: sz13,
+                    sz21: sz21,
+                    sz23: sz23,
+                    sz31: sz31,
+                    sz33: sz33,
                     t: satrec.t,
                     tc: tc,
                     gsto: satrec.gsto,
@@ -1754,7 +1691,7 @@ var Sgp4 = /** @class */ (function () {
                     z33: z33,
                     ecco: satrec.ecco,
                     eccsq: eccsq,
-                    em: em_1,
+                    em: em,
                     argpm: argpm,
                     inclm: inclm,
                     mm: mm,
@@ -1785,7 +1722,7 @@ var Sgp4 = /** @class */ (function () {
                     xli: satrec.xli,
                     xni: satrec.xni,
                 };
-                dsinitResult = Sgp4.dsinit(dsinitOptions);
+                var dsinitResult = Sgp4.dsinit(dsinitOptions);
                 satrec.irez = dsinitResult.irez;
                 satrec.atime = dsinitResult.atime;
                 satrec.d2201 = dsinitResult.d2201;
@@ -1813,9 +1750,9 @@ var Sgp4 = /** @class */ (function () {
             }
             // ----------- set variables if not deep space -----------
             if (satrec.isimp !== 1) {
-                cc1sq = satrec.cc1 * satrec.cc1;
+                var cc1sq = satrec.cc1 * satrec.cc1;
                 satrec.d2 = 4.0 * ao * tsi * cc1sq;
-                temp = (satrec.d2 * tsi * satrec.cc1) / 3.0;
+                var temp = (satrec.d2 * tsi * satrec.cc1) / 3.0;
                 satrec.d3 = (17.0 * ao + sfour) * temp;
                 satrec.d4 = 0.5 * temp * ao * tsi * (221.0 * ao + 31.0 * sfour) * satrec.cc1;
                 satrec.t3cof = satrec.d2 + 2.0 * cc1sq;
@@ -1927,84 +1864,45 @@ var Sgp4 = /** @class */ (function () {
         // 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
-        /** Ootk -- Declared at the top of the page */
-        // temp4 = 1.5e-12;
         // sgp4fix identify constants and allow alternate values
         // getgravconst( whichconst, tumin, mu, radiusearthkm, xke, j2, j3, j4, j3oj2 );
-        vkmpersec = (satrec.radiusearthkm * satrec.xke) / 60.0;
+        var xke = satrec.xke, j2 = satrec.j2, j3oj2 = satrec.j3oj2;
+        var vkmpersec = (satrec.radiusearthkm * satrec.xke) / 60.0;
         // --------------------- clear sgp4 error flag -----------------
         satrec.t = tsince;
         satrec.error = 0;
         //  ------- update for secular gravity and atmospheric drag -----
-        xmdf = satrec.mo + satrec.mdot * satrec.t;
-        argpdf = satrec.argpo + satrec.argpdot * satrec.t;
-        nodedf = satrec.nodeo + satrec.nodedot * satrec.t;
-        argpm = argpdf;
-        mm = xmdf;
-        t2 = satrec.t * satrec.t;
-        nodem = nodedf + satrec.nodecf * t2;
-        tempa = 1.0 - satrec.cc1 * satrec.t;
-        tempe = satrec.bstar * satrec.cc4 * satrec.t;
-        templ = satrec.t2cof * t2;
+        var xmdf = satrec.mo + satrec.mdot * satrec.t;
+        var argpdf = satrec.argpo + satrec.argpdot * satrec.t;
+        var nodedf = satrec.nodeo + satrec.nodedot * satrec.t;
+        var argpm = argpdf;
+        var mm = xmdf;
+        var t2 = satrec.t * satrec.t;
+        var nodem = nodedf + satrec.nodecf * t2;
+        var tempa = 1.0 - satrec.cc1 * satrec.t;
+        var tempe = satrec.bstar * satrec.cc4 * satrec.t;
+        var templ = satrec.t2cof * t2;
         if (satrec.isimp !== 1) {
-            delomg = satrec.omgcof * satrec.t;
+            var delomg = satrec.omgcof * satrec.t;
             //  sgp4fix use mutliply for speed instead of pow
-            delmtemp = 1.0 + satrec.eta * Math.cos(xmdf);
-            delm = satrec.xmcof * (delmtemp * delmtemp * delmtemp - satrec.delmo);
-            temp = delomg + delm;
-            mm = xmdf + temp;
-            argpm = argpdf - temp;
-            t3 = t2 * satrec.t;
-            t4 = t3 * satrec.t;
+            var delmtemp = 1.0 + satrec.eta * Math.cos(xmdf);
+            var delm = satrec.xmcof * (delmtemp * delmtemp * delmtemp - satrec.delmo);
+            var temp_1 = delomg + delm;
+            mm = xmdf + temp_1;
+            argpm = argpdf - temp_1;
+            var t3 = t2 * satrec.t;
+            var t4 = t3 * satrec.t;
             tempa = tempa - satrec.d2 * t2 - satrec.d3 * t3 - satrec.d4 * t4;
             tempe += satrec.bstar * satrec.cc5 * (Math.sin(mm) - satrec.sinmao);
             templ = templ + satrec.t3cof * t3 + t4 * (satrec.t4cof + satrec.t * satrec.t5cof);
         }
-        nm = satrec.no;
-        em = satrec.ecco;
-        inclm = satrec.inclo;
+        var nm = satrec.no;
+        var em = satrec.ecco;
+        var inclm = satrec.inclo;
         if (satrec.method === 'd') {
-            tc = satrec.t;
-            dspaceOptions = {
-                irez: satrec.irez,
-                d2201: satrec.d2201,
-                d2211: satrec.d2211,
-                d3210: satrec.d3210,
-                d3222: satrec.d3222,
-                d4410: satrec.d4410,
-                d4422: satrec.d4422,
-                d5220: satrec.d5220,
-                d5232: satrec.d5232,
-                d5421: satrec.d5421,
-                d5433: satrec.d5433,
-                dedt: satrec.dedt,
-                del1: satrec.del1,
-                del2: satrec.del2,
-                del3: satrec.del3,
-                didt: satrec.didt,
-                dmdt: satrec.dmdt,
-                dnodt: satrec.dnodt,
-                domdt: satrec.domdt,
-                argpo: satrec.argpo,
-                argpdot: satrec.argpdot,
-                t: satrec.t,
-                tc: tc,
-                gsto: satrec.gsto,
-                xfact: satrec.xfact,
-                xlamo: satrec.xlamo,
-                no: satrec.no,
-                atime: satrec.atime,
-                em: em,
-                argpm: argpm,
-                inclm: inclm,
-                xli: satrec.xli,
-                mm: mm,
-                xni: satrec.xni,
-                nodem: nodem,
-                nm: nm,
-            };
-            dspaceResult = Sgp4.dspace(dspaceOptions);
-            (em = dspaceResult.em, argpm = dspaceResult.argpm, inclm = dspaceResult.inclm, mm = dspaceResult.mm, nodem = dspaceResult.nodem, nm = dspaceResult.nm);
+            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 (nm <= 0.0) {
             // printf("// error nm %f\n", nm);
@@ -2012,7 +1910,7 @@ var Sgp4 = /** @class */ (function () {
             // sgp4fix add return
             return { position: false, velocity: false };
         }
-        am = Math.pow((xke / nm), x2o3) * tempa * tempa;
+        var am = Math.pow((xke / nm), x2o3) * tempa * tempa;
         nm = xke / Math.pow(am, 1.5);
         em -= tempe;
         // fix tolerance for error recognition
@@ -2030,32 +1928,33 @@ var Sgp4 = /** @class */ (function () {
             em = 1.0e-6;
         }
         mm += satrec.no * templ;
-        xlm = mm + argpm + nodem;
+        var xlm = mm + argpm + nodem;
         nodem %= TAU;
         argpm %= TAU;
         xlm %= TAU;
         mm = (xlm - argpm - nodem) % TAU;
         // sgp4fix recover singly averaged mean elements
-        satrec.am = am;
-        satrec.em = em;
-        satrec.im = inclm;
-        satrec.Om = nodem;
-        satrec.om = argpm;
-        satrec.mm = mm;
-        satrec.nm = nm;
+        // 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 -------------
-        sinim = Math.sin(inclm);
-        cosim = Math.cos(inclm);
+        var sinim = Math.sin(inclm);
+        var cosim = Math.cos(inclm);
         // -------------------- add lunar-solar periodics --------------
-        ep = em;
-        xincp = inclm;
-        argpp = argpm;
-        nodep = nodem;
-        mp = mm;
-        sinip = sinim;
-        cosip = cosim;
+        var ep = em;
+        var xincp = inclm;
+        var argpp = argpm;
+        var nodep = nodem;
+        var mp = mm;
+        var sinip = sinim;
+        var cosip = cosim;
         if (satrec.method === 'd') {
-            dpperParameters = {
+            var dpperParameters = {
                 inclo: satrec.inclo,
                 init: 'n',
                 ep: ep,
@@ -2065,7 +1964,7 @@ var Sgp4 = /** @class */ (function () {
                 mp: mp,
                 opsmode: satrec.operationmode,
             };
-            dpperResult = Sgp4.dpper(satrec, dpperParameters);
+            var dpperResult = Sgp4.dpper(satrec, dpperParameters);
             (ep = dpperResult.ep, nodep = dpperResult.nodep, argpp = dpperResult.argpp, mp = dpperResult.mp);
             xincp = dpperResult.inclp;
             if (xincp < 0.0) {
@@ -2094,17 +1993,18 @@ var Sgp4 = /** @class */ (function () {
                 satrec.xlcof = (-0.25 * j3oj2 * sinip * (3.0 + 5.0 * cosip)) / temp4;
             }
         }
-        axnl = ep * Math.cos(argpp);
-        temp = 1.0 / (am * (1.0 - ep * ep));
-        aynl = ep * Math.sin(argpp) + temp * satrec.aycof;
-        xl = mp + argpp + nodep + temp * satrec.xlcof * axnl;
+        var axnl = ep * Math.cos(argpp);
+        var temp = 1.0 / (am * (1.0 - ep * ep));
+        var aynl = ep * Math.sin(argpp) + temp * satrec.aycof;
+        var xl = mp + argpp + nodep + temp * satrec.xlcof * axnl;
         // --------------------- solve kepler's equation ---------------
-        u = (xl - nodep) % TAU;
-        eo1 = u;
-        tem5 = 9999.9;
-        ktr = 1;
+        var u = (xl - nodep) % TAU;
+        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;
         while (Math.abs(tem5) >= 1.0e-12 && ktr <= 10) {
             sineo1 = Math.sin(eo1);
             coseo1 = Math.cos(eo1);
@@ -2122,37 +2022,37 @@ var Sgp4 = /** @class */ (function () {
             ktr += 1;
         }
         //  ------------- short period preliminary quantities -----------
-        ecose = axnl * coseo1 + aynl * sineo1;
-        esine = axnl * sineo1 - aynl * coseo1;
-        el2 = axnl * axnl + aynl * aynl;
-        pl = am * (1.0 - el2);
+        var ecose = axnl * coseo1 + aynl * sineo1;
+        var esine = axnl * sineo1 - aynl * coseo1;
+        var el2 = axnl * axnl + aynl * aynl;
+        var pl = am * (1.0 - el2);
         if (pl < 0.0) {
             //  printf("// error pl %f\n", pl);
             satrec.error = 4;
             //  sgp4fix add return
             return { position: false, velocity: false };
         }
-        rl = am * (1.0 - ecose);
-        rdotl = (Math.sqrt(am) * esine) / rl;
-        rvdotl = Math.sqrt(pl) / rl;
-        betal = Math.sqrt(1.0 - el2);
+        var rl = am * (1.0 - ecose);
+        var rdotl = (Math.sqrt(am) * esine) / rl;
+        var rvdotl = Math.sqrt(pl) / rl;
+        var betal = Math.sqrt(1.0 - el2);
         temp = esine / (1.0 + betal);
-        sinu = (am / rl) * (sineo1 - aynl - axnl * temp);
-        cosu = (am / rl) * (coseo1 - axnl + aynl * temp);
-        su = Math.atan2(sinu, cosu);
-        sin2u = (cosu + cosu) * sinu;
-        cos2u = 1.0 - 2.0 * sinu * sinu;
+        var sinu = (am / rl) * (sineo1 - aynl - axnl * temp);
+        var cosu = (am / rl) * (coseo1 - axnl + aynl * temp);
+        var su = Math.atan2(sinu, cosu);
+        var sin2u = (cosu + cosu) * sinu;
+        var cos2u = 1.0 - 2.0 * sinu * sinu;
         temp = 1.0 / pl;
-        temp1 = 0.5 * j2 * temp;
-        temp2 = temp1 * temp;
+        var temp1 = 0.5 * j2 * temp;
+        var temp2 = temp1 * temp;
         // -------------- update for short period periodics ------------
         if (satrec.method === 'd') {
-            cosisq = cosip * cosip;
+            var cosisq = cosip * cosip;
             satrec.con41 = 3.0 * cosisq - 1.0;
             satrec.x1mth2 = 1.0 - cosisq;
             satrec.x7thm1 = 7.0 * cosisq - 1.0;
         }
-        mrt = rl * (1.0 - 1.5 * temp2 * betal * satrec.con41) + 0.5 * temp1 * satrec.x1mth2 * cos2u;
+        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
         if (mrt < 1.0) {
@@ -2164,32 +2064,32 @@ var Sgp4 = /** @class */ (function () {
             };
         }
         su -= 0.25 * temp2 * satrec.x7thm1 * sin2u;
-        xnode = nodep + 1.5 * temp2 * cosip * sin2u;
-        xinc = xincp + 1.5 * temp2 * cosip * sinip * cos2u;
-        mvt = rdotl - (nm * temp1 * satrec.x1mth2 * sin2u) / xke;
-        rvdot = rvdotl + (nm * temp1 * (satrec.x1mth2 * cos2u + 1.5 * satrec.con41)) / xke;
+        var xnode = nodep + 1.5 * temp2 * cosip * sin2u;
+        var xinc = xincp + 1.5 * temp2 * cosip * sinip * cos2u;
+        var mvt = rdotl - (nm * temp1 * satrec.x1mth2 * sin2u) / xke;
+        var rvdot = rvdotl + (nm * temp1 * (satrec.x1mth2 * cos2u + 1.5 * satrec.con41)) / xke;
         // --------------------- orientation vectors -------------------
-        sinsu = Math.sin(su);
-        cossu = Math.cos(su);
-        snod = Math.sin(xnode);
-        cnod = Math.cos(xnode);
-        sini = Math.sin(xinc);
-        cosi = Math.cos(xinc);
-        xmx = -snod * cosi;
-        xmy = cnod * cosi;
-        ux = xmx * sinsu + cnod * cossu;
-        uy = xmy * sinsu + snod * cossu;
-        uz = sini * sinsu;
-        vx = xmx * cossu - cnod * sinsu;
-        vy = xmy * cossu - snod * sinsu;
-        vz = sini * cossu;
+        var sinsu = Math.sin(su);
+        var cossu = Math.cos(su);
+        var snod = Math.sin(xnode);
+        var cnod = Math.cos(xnode);
+        var sini = Math.sin(xinc);
+        var cosi = Math.cos(xinc);
+        var xmx = -snod * cosi;
+        var xmy = cnod * cosi;
+        var ux = xmx * sinsu + cnod * cossu;
+        var uy = xmy * sinsu + snod * cossu;
+        var uz = sini * sinsu;
+        var vx = xmx * cossu - cnod * sinsu;
+        var vy = xmy * cossu - snod * sinsu;
+        var vz = sini * cossu;
         // --------- position and velocity (in km and km/sec) ----------
-        r = {
+        var r = {
             x: mrt * ux * satrec.radiusearthkm,
             y: mrt * uy * satrec.radiusearthkm,
             z: mrt * uz * satrec.radiusearthkm,
         };
-        v = {
+        var v = {
             x: (mvt * ux + rvdot * vx) * vkmpersec,
             y: (mvt * uy + rvdot * vy) * vkmpersec,
             z: (mvt * uz + rvdot * vz) * vkmpersec,
@@ -2229,6 +2129,7 @@ var Sgp4 = /** @class */ (function () {
     *    vallado, crawford, hujsak, kelso  2006
     --------------------------------------------------------------------------- */
     Sgp4.getgravconst = function (whichconst) {
+        var tumin, mus, radiusearthkm, xke, j2, j3, j4, j3oj2;
         switch (whichconst) {
             // -- wgs-72 low precision str#3 constants --
             case 'wgs72old':
@@ -2264,7 +2165,7 @@ var Sgp4 = /** @class */ (function () {
                 j3oj2 = j3 / j2;
                 break;
             default:
-                throw new Error("unknown gravity option " + whichconst);
+                throw new Error("unknown gravity option ".concat(whichconst));
         }
         return {
             tumin: tumin,
@@ -2319,22 +2220,140 @@ var Sgp4 = /** @class */ (function () {
     Sgp4.createSatrec = function (tleLine1, tleLine2, whichconst, opsmode) {
         if (whichconst === void 0) { whichconst = 'wgs72'; }
         if (opsmode === void 0) { opsmode = 'i'; }
-        // xpdotp = 1440.0 / (2.0 * PI); // 229.1831180523293;
-        year = 0;
-        satrec = {};
+        var year = 0;
+        var satrec = {
+            a: null,
+            am: null,
+            alta: null,
+            altp: null,
+            argpdot: null,
+            argpo: null,
+            aycof: null,
+            bstar: null,
+            cc1: null,
+            cc4: null,
+            cc5: null,
+            con41: null,
+            d2: null,
+            d3: null,
+            d4: null,
+            d5232: null,
+            d5421: null,
+            d5433: null,
+            dedt: null,
+            delmo: null,
+            del1: null,
+            ecco: null,
+            em: null,
+            epochdays: null,
+            epochyr: null,
+            error: null,
+            eta: null,
+            gsto: null,
+            im: null,
+            inclo: null,
+            init: null,
+            isimp: null,
+            jdsatepoch: null,
+            mdot: null,
+            method: null,
+            mo: null,
+            mm: null,
+            nddot: null,
+            ndot: null,
+            no: null,
+            nodecf: null,
+            nodedot: null,
+            nodeo: null,
+            om: null,
+            Om: null,
+            omgcof: null,
+            operationmode: null,
+            satnum: null,
+            sinmao: null,
+            t: null,
+            t2cof: null,
+            t3cof: null,
+            t4cof: null,
+            t5cof: null,
+            x1mth2: null,
+            x7thm1: null,
+            xlcof: null,
+            xmcof: null,
+            xfact: null,
+            xlamo: null,
+            xli: null,
+            xgh4: null,
+            xgh3: null,
+            xh2: null,
+            xi2: null,
+            xi3: null,
+            xl2: null,
+            xl3: null,
+            xl4: null,
+            zmol: null,
+            zmos: null,
+            dmdt: null,
+            dnodt: null,
+            domdt: null,
+            e3: null,
+            ee2: null,
+            peo: null,
+            pgho: null,
+            pho: null,
+            PInco: null,
+            plo: null,
+            se2: null,
+            se3: null,
+            sgh2: null,
+            sgh3: null,
+            sgh4: null,
+            sh2: null,
+            sh3: null,
+            si2: null,
+            si3: null,
+            sl2: null,
+            sl3: null,
+            sl4: null,
+            xgh2: null,
+            xh3: null,
+            tumin: null,
+            radiusearthkm: null,
+            irez: null,
+            d3210: null,
+            d3222: null,
+            d4410: null,
+            d4422: null,
+            d5220: null,
+            del2: null,
+            del3: null,
+            didt: null,
+            atime: null,
+            j2: null,
+            j3: null,
+            j4: null,
+            mus: null,
+            xke: null,
+            j3oj2: null,
+            xni: null,
+            d2201: null,
+            d2211: null,
+            nm: null,
+        };
         // 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);
         satrec.epochyr = parseInt(tleLine1.substring(18, 20));
         satrec.epochdays = parseFloat(tleLine1.substring(20, 32));
         satrec.ndot = parseFloat(tleLine1.substring(33, 43));
-        satrec.nddot = parseFloat(tleLine1.substring(44, 45) + "." + tleLine1.substring(45, 50) + "E" + tleLine1.substring(50, 52));
-        satrec.bstar = parseFloat(tleLine1.substring(53, 54) + "." + tleLine1.substring(54, 59) + "E" + tleLine1.substring(59, 61));
+        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.inclo = parseFloat(tleLine2.substring(8, 16));
         satrec.nodeo = parseFloat(tleLine2.substring(17, 25));
-        satrec.ecco = parseFloat("." + tleLine2.substring(26, 33));
+        satrec.ecco = parseFloat(".".concat(tleLine2.substring(26, 33)));
         satrec.argpo = parseFloat(tleLine2.substring(34, 42));
         satrec.mo = parseFloat(tleLine2.substring(43, 51));
         satrec.no = parseFloat(tleLine2.substring(52, 63));
@@ -2416,11 +2435,8 @@ var Sgp4 = /** @class */ (function () {
      *    vallado       2004, 191, eq 3-45
      * --------------------------------------------------------------------------- */
     Sgp4.gstime = function (jdut1) {
-        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 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
         temp = ((temp * DEG2RAD) / 240.0) % TAU; // 360/86400 = 1/240, to deg, to rad
         //  ------------------------ check quadrants ---------------------
         if (temp < 0.0) {
@@ -2540,7 +2556,7 @@ var Sgp4 = /** @class */ (function () {
         var magv1 = Sgp4.mag(vec1);
         var magv2 = Sgp4.mag(vec2);
         if (magv1 * magv2 > small * small) {
-            temp = Sgp4.dot(vec1, vec2) / (magv1 * magv2);
+            var temp = Sgp4.dot(vec1, vec2) / (magv1 * magv2);
             if (Math.abs(temp) > 1.0)
                 temp = Sgp4.sgn(temp) * 1.0;
             return Math.acos(temp);
@@ -2747,7 +2763,7 @@ var Sgp4 = /** @class */ (function () {
             nbar[1] = hbar[0];
             nbar[2] = 0.0;
             magn = Sgp4.mag(nbar);
-            c1 = magv * magv - mus / magr;
+            var c1 = magv * magv - mus / magr;
             rdotv = Sgp4.dot(r, v);
             for (i = 0; i <= 2; i++)
                 ebar[i] = (c1 * r[i] - rdotv * v[i]) / mus;
@@ -2783,7 +2799,7 @@ var Sgp4 = /** @class */ (function () {
             }
             // ----------  find right ascension of the ascending node ------------
             if (magn > small) {
-                temp = nbar[0] / magn;
+                var temp = nbar[0] / magn;
                 if (Math.abs(temp) > 1.0)
                     temp = Sgp4.sgn(temp);
                 omega = Math.acos(temp);
@@ -2818,7 +2834,7 @@ var Sgp4 = /** @class */ (function () {
             else
                 arglat = unknown;
             if (ecc > small && typeorbit == 4) {
-                temp = ebar[0] / ecc;
+                var temp = ebar[0] / ecc;
                 if (Math.abs(temp) > 1.0)
                     temp = Sgp4.sgn(temp);
                 lonper = Math.acos(temp);
@@ -2831,7 +2847,7 @@ var Sgp4 = /** @class */ (function () {
                 lonper = unknown;
             // -------- find true longitude - circular equatorial ------
             if (magr > small && typeorbit == 2) {
-                temp = r[0] / magr;
+                var temp = r[0] / magr;
                 if (Math.abs(temp) > 1.0)
                     temp = Sgp4.sgn(temp);
                 truelon = Math.acos(temp);
@@ -2960,26 +2976,26 @@ var Sgp4 = /** @class */ (function () {
      *    none.
      * --------------------------------------------------------------------------- */
     Sgp4.days2mdhms = function (year, days) {
-        lmonth = [31, year % 4 === 0 ? 29 : 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31];
-        dayofyr = Math.floor(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;
-        i = 1;
-        inttemp = 0;
+        var i = 1;
+        var inttemp = 0;
         while (dayofyr > inttemp + lmonth[i - 1] && i < 12) {
             inttemp += lmonth[i - 1];
             i += 1;
         }
-        mon = i;
-        day = dayofyr - inttemp;
+        var mon = i;
+        var day = dayofyr - inttemp;
         //  ----------------- find hours minutes and seconds -------------
-        temp = (days - dayofyr) * 24.0;
-        hr = Math.floor(temp);
+        var temp = (days - dayofyr) * 24.0;
+        var hr = Math.floor(temp);
         temp = (temp - hr) * 60.0;
-        minute = Math.floor(temp);
-        sec = (temp - minute) * 60.0;
+        var minute = Math.floor(temp);
+        var sec = (temp - minute) * 60.0;
         return {
             mon: mon,
             day: day,
@@ -3045,7 +3061,7 @@ var Sgp4 = /** @class */ (function () {
         /* --------------- find year and days of the year --------------- */
         var temp = jd - 2415019.5;
         var tu = temp / 365.25;
-        year = 1900 + Math.floor(tu);
+        var year = 1900 + Math.floor(tu);
         leapyrs = Math.floor((year - 1901) * 0.25);
         days = Math.floor(temp - ((year - 1900) * 365.0 + leapyrs));
         /* ------------ check for case of beginning of a year ----------- */
@@ -3068,4 +3084,4 @@ var Sgp4 = /** @class */ (function () {
     return Sgp4;
 }());
 export { Sgp4 };
-//# sourceMappingURL=ootk-sgp4.js.map
+//# sourceMappingURL=sgp4.js.map