ootk 3.0.0 → 4.0.0-1

package/dist/force/AtmosphericDrag.js

@@ -0,0 +1,73 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { DataHandler, DEG2RAD, Earth, ITRF, J2000, Sun, Vector3D, } from 'ootk-core';
+/**
+ * Harris-Priester atmospheric drag force model.
+ * Atmospheric density model assumes mean solar flux.
+ */
+export class AtmosphericDrag {
+    mass;
+    area;
+    dragCoeff;
+    cosine;
+    constructor(mass, area, dragCoeff, cosine) {
+        this.mass = mass;
+        this.area = area;
+        this.dragCoeff = dragCoeff;
+        this.cosine = cosine;
+    }
+    static _getHPDensity(state, n) {
+        const hpa = DataHandler.getInstance().getHpAtmosphere(state.height);
+        if (hpa === null) {
+            return 0.0;
+        }
+        const sunPos = Sun.positionApparent(state.epoch);
+        const sunVec = new J2000(state.epoch, sunPos, Vector3D.origin).toITRF().position.normalize();
+        const bulVec = sunVec.rotZ(-30.0 * DEG2RAD);
+        const cosPsi = bulVec.normalize().dot(state.position.normalize());
+        const c2Psi2 = 0.5 * (1.0 + cosPsi);
+        const cPsi2 = Math.sqrt(c2Psi2);
+        const cosPow = cPsi2 > 1e-12 ? c2Psi2 * cPsi2 ** (n - 2) : 0.0;
+        const altitude = hpa.height;
+        const [h0, min0, max0] = hpa.hp0;
+        const [h1, min1, max1] = hpa.hp1;
+        const dH = (h0 - altitude) / (h0 - h1);
+        const rhoMin = min0 * (min1 / min0) ** dH;
+        if (cosPow === 0) {
+            return rhoMin;
+        }
+        const rhoMax = max0 * (max1 / max0) ** dH;
+        return rhoMin + (rhoMax - rhoMin) * cosPow;
+    }
+    acceleration(state) {
+        const itrfState = state.toITRF();
+        const density = AtmosphericDrag._getHPDensity(itrfState, this.cosine);
+        if (density === 0) {
+            return Vector3D.origin;
+        }
+        const rotation = new ITRF(state.epoch, 
+        // TODO: #26 This is a bit of a hack to get Radians/Second to work with StateVector class
+        Earth.rotation, Vector3D.origin).toJ2000().position;
+        // TODO: #26 This is a bit of a hack to get Radians/Second to work with StateVector class
+        const vRel = state.velocity.subtract(rotation.cross(state.position))
+            .scale(1000.0);
+        const vm = vRel.magnitude();
+        const fScale = -0.5 * density * ((this.dragCoeff * this.area) / this.mass) * vm;
+        return vRel.scale(fScale / 1000.0);
+    }
+}
+//# sourceMappingURL=AtmosphericDrag.js.map

package/dist/force/AtmosphericDrag.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"AtmosphericDrag.js","sourceRoot":"","sources":["../../src/force/AtmosphericDrag.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EACL,WAAW,EACX,OAAO,EACP,KAAK,EACL,IAAI,EACJ,KAAK,EAIL,GAAG,EACH,QAAQ,GACT,MAAM,WAAW,CAAC;AAGnB;;;GAGG;AACH,MAAM,OAAO,eAAe;IAC1B,IAAI,CAAS;IACb,IAAI,CAAS;IACb,SAAS,CAAS;IAClB,MAAM,CAAS;IAEf,YAAY,IAAY,EAAE,IAAY,EAAE,SAAiB,EAAE,MAAc;QACvE,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;IAEO,MAAM,CAAC,aAAa,CAAC,KAAW,EAAE,CAAS;QACjD,MAAM,GAAG,GAAG,WAAW,CAAC,WAAW,EAAE,CAAC,eAAe,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QAEpE,IAAI,GAAG,KAAK,IAAI,EAAE;YAChB,OAAO,GAAG,CAAC;SACZ;QACD,MAAM,MAAM,GAAG,GAAG,CAAC,gBAAgB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QACjD,MAAM,MAAM,GAAG,IAAI,KAAK,CACtB,KAAK,CAAC,KAAK,EACX,MAAM,EACN,QAAQ,CAAC,MAAuC,CACjD,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAChC,MAAM,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,IAAI,GAAG,OAAkB,CAAC,CAAC;QACvD,MAAM,MAAM,GAAG,MAAM,CAAC,SAAS,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC,CAAC;QAClE,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC;QACpC,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAChC,MAAM,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,GAAG,KAAK,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC/D,MAAM,QAAQ,GAAG,GAAG,CAAC,MAAM,CAAC;QAC5B,MAAM,CAAC,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC;QACjC,MAAM,CAAC,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACvC,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;QAE1C,IAAI,MAAM,KAAK,CAAC,EAAE;YAChB,OAAO,MAAM,CAAC;SACf;QACD,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;QAE1C,OAAO,MAAM,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,MAAM,CAAC;IAC7C,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,MAAM,SAAS,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;QACjC,MAAM,OAAO,GAAG,eAAe,CAAC,aAAa,CAAC,SAAS,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAEtE,IAAI,OAAO,KAAK,CAAC,EAAE;YACjB,OAAO,QAAQ,CAAC,MAAM,CAAC;SACxB;QACD,MAAM,QAAQ,GAAG,IAAI,IAAI,CACvB,KAAK,CAAC,KAAK;QACX,yFAAyF;QACzF,KAAK,CAAC,QAA2C,EACjD,QAAQ,CAAC,MAAuC,CACjD,CAAC,OAAO,EAAE,CAAC,QAAQ,CAAC;QACrB,yFAAyF;QACzF,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,KAAK,CAAC,KAAK,CAAC,QAAQ,CAA6C,CAAC;aAC7G,KAAK,CAAC,MAA6B,CAAC,CAAC;QACxC,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;QAC5B,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC;QAEhF,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC;IACrC,CAAC;CACF"}

package/dist/force/EarthGravity.d.ts

@@ -0,0 +1,27 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class EarthGravity implements Force {
+    degree: number;
+    order: number;
+    _asphericalFlag: boolean;
+    constructor(degree: number, order: number);
+    _spherical(state: J2000): Vector3D;
+    _aspherical(state: J2000): Vector3D;
+    acceleration(state: J2000): Vector3D;
+}

package/dist/force/EarthGravity.js

@@ -0,0 +1,117 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { DataHandler, Earth, ITRF, Vector3D } from 'ootk-core';
+// / tesseral geopotential perturbations.
+export class EarthGravity {
+    degree;
+    order;
+    _asphericalFlag;
+    constructor(degree, order) {
+        this.degree = Math.min(Math.max(degree, 0), 36);
+        this.order = Math.min(Math.max(order, 0), 36);
+        this._asphericalFlag = degree >= 2;
+    }
+    _spherical(state) {
+        const rMag = state.position.magnitude();
+        return state.position.scale(-Earth.mu / (rMag * rMag * rMag));
+    }
+    // eslint-disable-next-line max-statements
+    _aspherical(state) {
+        const posEcef = state.toITRF().position;
+        const ri = 1.0 / posEcef.magnitude();
+        const xor = posEcef.x * ri;
+        const yor = posEcef.y * ri;
+        const zor = posEcef.z * ri;
+        const ep = zor;
+        const reor = Earth.radiusEquator * ri;
+        let reorn = reor;
+        const muor2 = Earth.mu * ri * ri;
+        let sumH = 0.0;
+        let sumGm = 0.0;
+        let sumJ = 0.0;
+        let sumK = 0.0;
+        const cTil = new Float64Array(this.order + 4);
+        const sTil = new Float64Array(this.order + 4);
+        const pN = new Float64Array(this.order + 4);
+        const pNm1 = new Float64Array(this.order + 4);
+        const pNm2 = new Float64Array(this.order + 4);
+        pNm2[0] = 1.0;
+        pNm1[0] = ep;
+        pNm1[1] = 1.0;
+        cTil[0] = 1.0;
+        cTil[1] = xor;
+        sTil[1] = yor;
+        const dh = DataHandler.getInstance();
+        for (let n = 2, nm1 = 1, nm2 = 0, np1 = 3; n <= this.degree; nm2++, nm1++, n++, np1++) {
+            const twonm1 = 2.0 * n - 1.0;
+            reorn *= reor;
+            const cN0 = dh.getEgm96Coeffs(n, 0)[2];
+            pN[0] = (twonm1 * ep * pNm1[0] - nm1 * pNm2[0]) / n;
+            pN[1] = pNm2[1] + twonm1 * pNm1[0];
+            pN[2] = pNm2[2] + twonm1 * pNm1[1];
+            let sumHn = pN[1] * cN0;
+            let sumGmn = pN[0] * cN0 * np1;
+            if (this.order > 0) {
+                let sumJn = 0.0;
+                let sumKn = 0.0;
+                cTil[n] = cTil[1] * cTil[nm1] - sTil[1] * sTil[nm1];
+                sTil[n] = sTil[1] * cTil[nm1] + cTil[1] * sTil[nm1];
+                const lim = n < this.order ? n : this.order;
+                for (let m = 1, mm1 = 0, mm2 = -1, mp1 = 2, mp2 = 3; m <= lim; mm2++, mm1++, m++, mp1++, mp2++) {
+                    pN[mp1] = pNm2[mp1] + twonm1 * pNm1[m];
+                    const dm = m;
+                    const npmp1 = n + mp1;
+                    const pNm = pN[m];
+                    const pNmp1 = pN[mp1];
+                    const coefs = dh.getEgm96Coeffs(n, m);
+                    const cNm = coefs[2];
+                    const sNm = coefs[3];
+                    const mxPnm = dm * pNm;
+                    const bNmtil = cNm * cTil[m] + sNm * sTil[m];
+                    const pNmBnm = pNm * bNmtil;
+                    const bNmtm1 = cNm * cTil[mm1] + sNm * sTil[mm1];
+                    const aNmtm1 = cNm * sTil[mm1] - sNm * cTil[mm1];
+                    sumHn += pNmp1 * bNmtil;
+                    sumGmn += npmp1 * pNmBnm;
+                    sumJn += mxPnm * bNmtm1;
+                    sumKn -= mxPnm * aNmtm1;
+                }
+                sumJ += reorn * sumJn;
+                sumK += reorn * sumKn;
+            }
+            sumH += reorn * sumHn;
+            sumGm += reorn * sumGmn;
+            if (n < this.degree) {
+                for (let i = 0; i <= n; i++) {
+                    pNm2[i] = pNm1[i];
+                    pNm1[i] = pN[i];
+                }
+            }
+        }
+        const lambda = sumGm + ep * sumH;
+        const g = new Vector3D(-muor2 * (lambda * xor - sumJ), -muor2 * (lambda * yor - sumK), -muor2 * (lambda * zor - sumH));
+        return new ITRF(state.epoch, g, Vector3D.origin).toJ2000().position;
+    }
+    acceleration(state) {
+        let accVec = this._spherical(state);
+        if (this._asphericalFlag) {
+            accVec = accVec.add(this._aspherical(state));
+        }
+        return accVec;
+    }
+}
+//# sourceMappingURL=EarthGravity.js.map

package/dist/force/EarthGravity.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"EarthGravity.js","sourceRoot":"","sources":["../../src/force/EarthGravity.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,WAAW,EAAE,KAAK,EAAE,IAAI,EAA0C,QAAQ,EAAE,MAAM,WAAW,CAAC;AAMvG,yCAAyC;AACzC,MAAM,OAAO,YAAY;IACvB,MAAM,CAAS;IACf,KAAK,CAAS;IACd,eAAe,CAAU;IAEzB,YAAY,MAAc,EAAE,KAAa;QACvC,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAChD,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC9C,IAAI,CAAC,eAAe,GAAG,MAAM,IAAI,CAAC,CAAC;IACrC,CAAC;IAED,UAAU,CAAC,KAAY;QACrB,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAExC,OAAO,KAAK,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC;IAChE,CAAC;IAED,0CAA0C;IAC1C,WAAW,CAAC,KAAY;QACtB,MAAM,OAAO,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC;QACxC,MAAM,EAAE,GAAG,GAAG,GAAG,OAAO,CAAC,SAAS,EAAE,CAAC;QACrC,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC;QAE3B,MAAM,EAAE,GAAG,GAAG,CAAC;QACf,MAAM,IAAI,GAAG,KAAK,CAAC,aAAa,GAAG,EAAE,CAAC;QACtC,IAAI,KAAK,GAAG,IAAI,CAAC;QACjB,MAAM,KAAK,GAAG,KAAK,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAEjC,IAAI,IAAI,GAAG,GAAG,CAAC;QACf,IAAI,KAAK,GAAG,GAAG,CAAC;QAChB,IAAI,IAAI,GAAG,GAAG,CAAC;QACf,IAAI,IAAI,GAAG,GAAG,CAAC;QAEf,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC9C,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAE9C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC5C,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC9C,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAE9C,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAEd,MAAM,EAAE,GAAG,WAAW,CAAC,WAAW,EAAE,CAAC;QAErC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,MAAM,EAAE,GAAG,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE;YACrF,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;YAE7B,KAAK,IAAI,IAAI,CAAC;YACd,MAAM,GAAG,GAAG,EAAE,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAEvC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACpD,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACnC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YAEnC,IAAI,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YACxB,IAAI,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;YAE/B,IAAI,IAAI,CAAC,KAAK,GAAG,CAAC,EAAE;gBAClB,IAAI,KAAK,GAAG,GAAG,CAAC;gBAChB,IAAI,KAAK,GAAG,GAAG,CAAC;gBAEhB,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;gBACpD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;gBAEpD,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;gBAE5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,GAAG,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,GAAG,EAAE,EAAE;oBAC9F,EAAE,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBAEvC,MAAM,EAAE,GAAG,CAAC,CAAC;oBACb,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,CAAC;oBAEtB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBAClB,MAAM,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;oBAEtB,MAAM,KAAK,GAAG,EAAE,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;oBACtC,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBACrB,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAErB,MAAM,KAAK,GAAG,EAAE,GAAG,GAAG,CAAC;oBACvB,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBAC7C,MAAM,MAAM,GAAG,GAAG,GAAG,MAAM,CAAC;oBAC5B,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;oBACjD,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;oBAEjD,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;oBACxB,MAAM,IAAI,KAAK,GAAG,MAAM,CAAC;oBACzB,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;oBACxB,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;iBACzB;gBAED,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC;gBACtB,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC;aACvB;YAED,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC;YACtB,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;YAExB,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE;gBACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;oBAC3B,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBAClB,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;iBACjB;aACF;SACF;QAED,MAAM,MAAM,GAAG,KAAK,GAAG,EAAE,GAAG,IAAI,CAAC;QACjC,MAAM,CAAC,GAAG,IAAI,QAAQ,CACpB,CAAC,KAAK,GAAG,CAAC,MAAM,GAAG,GAAG,GAAG,IAAI,CAAe,EAC5C,CAAC,KAAK,GAAG,CAAC,MAAM,GAAG,GAAG,GAAG,IAAI,CAAe,EAC5C,CAAC,KAAK,GAAG,CAAC,MAAM,GAAG,GAAG,GAAG,IAAI,CAAe,CAC7C,CAAC;QAEF,OAAO,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,EAAE,QAAQ,CAAC,MAAuC,CAAC,CAAC,OAAO,EAAE,CAAC,QAAQ,CAAC;IACvG,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,IAAI,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;QAEpC,IAAI,IAAI,CAAC,eAAe,EAAE;YACxB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC;SAC9C;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF"}

package/{lib/ootk.d.ts → dist/force/Force.d.ts}

@@ -1,10 +1,7 @@
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
- * with satellites and other orbital objects.
- *
+ * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @Copyright (c) 2020-2023 Theodore Kruczek
+ * @copyright (c) 2020-2024 Theodore Kruczek
  *
  * 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
@@ -17,9 +14,13 @@
  * 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/>.
  */
-export { Sgp4 } from './sgp4/sgp4';
-export { Transforms } from './transforms/transforms';
-export { Utils } from './utils/utils';
-export { Tle } from './tle/tle';
-export * from './objects/index';
-export * from './types/types';
+import { J2000, Vector3D } from 'ootk-core';
+export declare abstract class Force {
+    /**
+     * Calculate the acceleration due to the perturbing force on a given
+     * state vector.
+     * @param state The state vector.
+     * @throws If the force cannot be calculated.
+     */
+    acceleration(state: J2000): Vector3D;
+}

package/dist/force/Force.js

@@ -0,0 +1,31 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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/>.
+ */
+/* eslint-disable @typescript-eslint/no-unused-vars */
+/* eslint-disable class-methods-use-this */
+// / Base class for perturbation forces.
+export class Force {
+    /**
+     * Calculate the acceleration due to the perturbing force on a given
+     * state vector.
+     * @param state The state vector.
+     * @throws If the force cannot be calculated.
+     */
+    acceleration(state) {
+        throw Error('Not implemented');
+    }
+}
+//# sourceMappingURL=Force.js.map

package/dist/force/Force.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Force.js","sourceRoot":"","sources":["../../src/force/Force.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAGH,sDAAsD;AACtD,2CAA2C;AAE3C,wCAAwC;AACxC,MAAM,OAAgB,KAAK;IACzB;;;;;OAKG;IACH,YAAY,CAAC,KAAY;QACvB,MAAM,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACjC,CAAC;CACF"}

package/dist/force/ForceModel.d.ts

@@ -0,0 +1,37 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { J2000, Vector, Vector3D } from 'ootk-core';
+import { Thrust } from './Thrust';
+export declare class ForceModel {
+    private _centralGravity?;
+    private _thirdBodyGravity?;
+    private _solarRadiationPressure?;
+    private _atmosphericDrag?;
+    private _maneuverThrust;
+    setGravity(mu?: number): this;
+    setEarthGravity(degree: number, order: number): void;
+    setThirdBodyGravity({ moon, sun }: {
+        moon?: boolean | undefined;
+        sun?: boolean | undefined;
+    }): void;
+    setSolarRadiationPressure(mass: number, area: number, coeff?: number): void;
+    setAtmosphericDrag(mass: number, area: number, coeff?: number, cosine?: number): void;
+    loadManeuver(maneuver: Thrust): void;
+    clearManeuver(): void;
+    acceleration(state: J2000): Vector3D;
+    derivative(state: J2000): Vector;
+}

package/dist/force/ForceModel.js

@@ -0,0 +1,75 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { Earth, Vector3D } from 'ootk-core';
+import { AtmosphericDrag } from './AtmosphericDrag';
+import { EarthGravity } from './EarthGravity';
+import { Gravity } from './Gravity';
+import { SolarRadiationPressure } from './SolarRadiationPressure';
+import { ThirdBodyGravity } from './ThirdBodyGravity';
+// / Force model for spacecraft propagation.
+export class ForceModel {
+    _centralGravity;
+    _thirdBodyGravity;
+    _solarRadiationPressure;
+    _atmosphericDrag;
+    _maneuverThrust = null;
+    setGravity(mu = Earth.mu) {
+        this._centralGravity = new Gravity(mu);
+        return this;
+    }
+    setEarthGravity(degree, order) {
+        this._centralGravity = new EarthGravity(degree, order);
+    }
+    setThirdBodyGravity({ moon = false, sun = false }) {
+        this._thirdBodyGravity = new ThirdBodyGravity(moon, sun);
+    }
+    setSolarRadiationPressure(mass, area, coeff = 1.2) {
+        this._solarRadiationPressure = new SolarRadiationPressure(mass, area, coeff);
+    }
+    setAtmosphericDrag(mass, area, coeff = 2.2, cosine = 4) {
+        this._atmosphericDrag = new AtmosphericDrag(mass, area, coeff, cosine);
+    }
+    loadManeuver(maneuver) {
+        this._maneuverThrust = maneuver;
+    }
+    clearManeuver() {
+        this._maneuverThrust = null;
+    }
+    acceleration(state) {
+        let accVec = Vector3D.origin;
+        if (this._centralGravity) {
+            accVec = accVec.add(this._centralGravity.acceleration(state));
+        }
+        if (this._thirdBodyGravity) {
+            accVec = accVec.add(this._thirdBodyGravity.acceleration(state));
+        }
+        if (this._solarRadiationPressure) {
+            accVec = accVec.add(this._solarRadiationPressure.acceleration(state));
+        }
+        if (this._atmosphericDrag) {
+            accVec = accVec.add(this._atmosphericDrag.acceleration(state));
+        }
+        if (this._maneuverThrust) {
+            accVec = accVec.add(this._maneuverThrust.acceleration(state));
+        }
+        return accVec;
+    }
+    derivative(state) {
+        return state.velocity.join(this.acceleration(state));
+    }
+}
+//# sourceMappingURL=ForceModel.js.map

package/dist/force/ForceModel.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ForceModel.js","sourceRoot":"","sources":["../../src/force/ForceModel.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,KAAK,EAAiB,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AACpD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,OAAO,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AACpC,OAAO,EAAE,sBAAsB,EAAE,MAAM,0BAA0B,CAAC;AAClE,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAGtD,4CAA4C;AAE5C,MAAM,OAAO,UAAU;IACb,eAAe,CAAS;IACxB,iBAAiB,CAAS;IAC1B,uBAAuB,CAAS;IAChC,gBAAgB,CAAS;IACzB,eAAe,GAAiB,IAAI,CAAC;IAE7C,UAAU,CAAC,KAAa,KAAK,CAAC,EAAE;QAC9B,IAAI,CAAC,eAAe,GAAG,IAAI,OAAO,CAAC,EAAE,CAAC,CAAC;QAEvC,OAAO,IAAI,CAAC;IACd,CAAC;IAED,eAAe,CAAC,MAAc,EAAE,KAAa;QAC3C,IAAI,CAAC,eAAe,GAAG,IAAI,YAAY,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;IACzD,CAAC;IAED,mBAAmB,CAAC,EAAE,IAAI,GAAG,KAAK,EAAE,GAAG,GAAG,KAAK,EAAE;QAC/C,IAAI,CAAC,iBAAiB,GAAG,IAAI,gBAAgB,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IAC3D,CAAC;IAED,yBAAyB,CAAC,IAAY,EAAE,IAAY,EAAE,KAAK,GAAG,GAAG;QAC/D,IAAI,CAAC,uBAAuB,GAAG,IAAI,sBAAsB,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;IAC/E,CAAC;IAED,kBAAkB,CAAC,IAAY,EAAE,IAAY,EAAE,KAAK,GAAG,GAAG,EAAE,MAAM,GAAG,CAAC;QACpE,IAAI,CAAC,gBAAgB,GAAG,IAAI,eAAe,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IACzE,CAAC;IAED,YAAY,CAAC,QAAgB;QAC3B,IAAI,CAAC,eAAe,GAAG,QAAQ,CAAC;IAClC,CAAC;IAED,aAAa;QACX,IAAI,CAAC,eAAe,GAAG,IAAI,CAAC;IAC9B,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,IAAI,MAAM,GAAG,QAAQ,CAAC,MAAM,CAAC;QAE7B,IAAI,IAAI,CAAC,eAAe,EAAE;YACxB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,eAAe,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SAC/D;QACD,IAAI,IAAI,CAAC,iBAAiB,EAAE;YAC1B,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,iBAAiB,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SACjE;QACD,IAAI,IAAI,CAAC,uBAAuB,EAAE;YAChC,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,uBAAuB,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SACvE;QACD,IAAI,IAAI,CAAC,gBAAgB,EAAE;YACzB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,gBAAgB,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SAChE;QACD,IAAI,IAAI,CAAC,eAAe,EAAE;YACxB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,eAAe,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SAC/D;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,UAAU,CAAC,KAAY;QACrB,OAAO,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;IACvD,CAAC;CACF"}

package/dist/force/Gravity.d.ts

@@ -0,0 +1,38 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class Gravity implements Force {
+    mu: number;
+    /**
+     * Create a new [Gravity] object with optional gravitational
+     * @param mu Gravitational parameter _(km²/s³)_.
+     */
+    constructor(mu?: number);
+    /**
+     * Calculates the gravitational force in spherical coordinates.
+     * @param state The J2000 state containing the position and velocity vectors.
+     * @returns The gravitational force vector in spherical coordinates.
+     */
+    private _spherical;
+    /**
+     * Calculates the acceleration due to gravity at a given state.
+     * @param state The J2000 state at which to calculate the acceleration.
+     * @returns The acceleration vector.
+     */
+    acceleration(state: J2000): Vector3D;
+}

package/dist/force/Gravity.js

@@ -0,0 +1,47 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { Earth } from 'ootk-core';
+// / Simple central-body gravity model.
+export class Gravity {
+    // / Gravitational parameter _(km²/s³)_.
+    mu;
+    /**
+     * Create a new [Gravity] object with optional gravitational
+     * @param mu Gravitational parameter _(km²/s³)_.
+     */
+    constructor(mu = Earth.mu) {
+        this.mu = mu;
+    }
+    /**
+     * Calculates the gravitational force in spherical coordinates.
+     * @param state The J2000 state containing the position and velocity vectors.
+     * @returns The gravitational force vector in spherical coordinates.
+     */
+    _spherical(state) {
+        const rMag = state.position.magnitude();
+        return state.position.scale(-this.mu / (rMag * rMag * rMag));
+    }
+    /**
+     * Calculates the acceleration due to gravity at a given state.
+     * @param state The J2000 state at which to calculate the acceleration.
+     * @returns The acceleration vector.
+     */
+    acceleration(state) {
+        return this._spherical(state);
+    }
+}
+//# sourceMappingURL=Gravity.js.map

package/dist/force/Gravity.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Gravity.js","sourceRoot":"","sources":["../../src/force/Gravity.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,KAAK,EAAmB,MAAM,WAAW,CAAC;AAGnD,uCAAuC;AACvC,MAAM,OAAO,OAAO;IAClB,wCAAwC;IACxC,EAAE,CAAS;IAEX;;;OAGG;IACH,YAAY,KAAa,KAAK,CAAC,EAAE;QAC/B,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACf,CAAC;IAED;;;;OAIG;IACK,UAAU,CAAC,KAAY;QAC7B,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAExC,OAAO,KAAK,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC;IAC/D,CAAC;IAED;;;;OAIG;IACH,YAAY,CAAC,KAAY;QACvB,OAAO,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;IAChC,CAAC;CACF"}

package/dist/{ootk.d.ts → force/SolarRadiationPressure.d.ts}

@@ -1,10 +1,7 @@
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
- * with satellites and other orbital objects.
- *
+ * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @Copyright (c) 2020-2023 Theodore Kruczek
+ * @copyright (c) 2020-2024 Theodore Kruczek
  *
  * 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
@@ -17,9 +14,13 @@
  * 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/>.
  */
-export { Sgp4 } from './sgp4/sgp4';
-export { Transforms } from './transforms/transforms';
-export { Utils } from './utils/utils';
-export { Tle } from './tle/tle';
-export * from './objects/index';
-export * from './types/types';
+import { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class SolarRadiationPressure extends Force {
+    mass: number;
+    area: number;
+    reflectCoeff: number;
+    constructor(mass: number, area: number, reflectCoeff: number);
+    private static readonly _kRef;
+    acceleration(state: J2000): Vector3D;
+}

package/dist/force/SolarRadiationPressure.js

@@ -0,0 +1,44 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { astronomicalUnit, Sun } from 'ootk-core';
+import { Force } from './Force';
+// / Solar radiation pressure model.
+export class SolarRadiationPressure extends Force {
+    mass;
+    area;
+    reflectCoeff;
+    // / Create a new [SolarRadiationPressure] object.
+    constructor(mass, area, reflectCoeff) {
+        super();
+        this.mass = mass;
+        this.area = area;
+        this.reflectCoeff = reflectCoeff;
+    }
+    // / Solar pressure _(N/m²)_;
+    static _kRef = 4.56e-6 * astronomicalUnit ** 2;
+    acceleration(state) {
+        const rSun = Sun.positionApparent(state.epoch);
+        const r = state.position.subtract(rSun);
+        const rMag = r.magnitude();
+        const r2 = rMag * rMag;
+        const ratio = Sun.lightingRatio(state.position, rSun);
+        const p = (ratio * SolarRadiationPressure._kRef) / r2;
+        const flux = r.scale(p / rMag);
+        return flux.scale(((this.area * this.reflectCoeff) / this.mass) * 1e-3);
+    }
+}
+//# sourceMappingURL=SolarRadiationPressure.js.map

package/dist/force/SolarRadiationPressure.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SolarRadiationPressure.js","sourceRoot":"","sources":["../../src/force/SolarRadiationPressure.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,gBAAgB,EAAS,GAAG,EAAY,MAAM,WAAW,CAAC;AACnE,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAEhC,oCAAoC;AACpC,MAAM,OAAO,sBAAuB,SAAQ,KAAK;IAE5B;IAAqB;IAAqB;IAD7D,kDAAkD;IAClD,YAAmB,IAAY,EAAS,IAAY,EAAS,YAAoB;QAC/E,KAAK,EAAE,CAAC;QADS,SAAI,GAAJ,IAAI,CAAQ;QAAS,SAAI,GAAJ,IAAI,CAAQ;QAAS,iBAAY,GAAZ,YAAY,CAAQ;IAEjF,CAAC;IAED,6BAA6B;IACrB,MAAM,CAAU,KAAK,GAAW,OAAO,GAAG,gBAAgB,IAAI,CAAC,CAAC;IAExE,YAAY,CAAC,KAAY;QACvB,MAAM,IAAI,GAAG,GAAG,CAAC,gBAAgB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QAC/C,MAAM,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC;QACvB,MAAM,KAAK,GAAG,GAAG,CAAC,aAAa,CAAC,KAAK,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;QACtD,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,sBAAsB,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;QACtD,MAAM,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QAE/B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,YAAY,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC;IAC1E,CAAC"}

package/dist/force/ThirdBodyGravity.d.ts

@@ -0,0 +1,26 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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.
+ *
+ * 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.
+ *
+ * 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 { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class ThirdBodyGravity implements Force {
+    moon: boolean;
+    sun: boolean;
+    constructor(moon?: boolean, sun?: boolean);
+    private static _moonGravity;
+    private static _sunGravity;
+    acceleration(state: J2000): Vector3D;
+}