@math.gl/sun 4.0.0-beta.1 → 4.0.1

package/dist/index.cjs

@@ -16,15 +16,15 @@ var __copyProps = (to, from, except, desc) => {
 };
 var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
 
-// src/index.ts
-var src_exports = {};
-__export(src_exports, {
+// dist/index.js
+var dist_exports = {};
+__export(dist_exports, {
   getSunDirection: () => getSunDirection,
   getSunPosition: () => getSunPosition
 });
-module.exports = __toCommonJS(src_exports);
+module.exports = __toCommonJS(dist_exports);
 
-// src/suncalc.ts
+// dist/suncalc.js
 var DEGREES_TO_RADIANS = Math.PI / 180;
 var DAY_IN_MS = 1e3 * 60 * 60 * 24;
 var JD1970 = 2440588;
@@ -100,3 +100,4 @@ function getSunCoords(dates) {
     rightAscension: getRightAscension(L, 0)
   };
 }
+//# sourceMappingURL=index.cjs.map

package/dist/index.cjs.map

@@ -0,0 +1,7 @@
+{
+  "version": 3,
+  "sources": ["index.js", "suncalc.js"],
+  "sourcesContent": ["export { getSunPosition, getSunDirection } from \"./suncalc.js\";\n", "const DEGREES_TO_RADIANS = Math.PI / 180;\nconst DAY_IN_MS = 1000 * 60 * 60 * 24;\nconst JD1970 = 2440588; // Julian Day year 1970\nconst JD2000 = 2451545; // Julian Day year 2000\n// This angle \u03B5 [epsilon] is called the obliquity of the ecliptic and its value at the beginning of 2000 was 23.4397\u00B0\nconst e = DEGREES_TO_RADIANS * 23.4397; // obliquity of the Earth\n// Refer https://www.aa.quae.nl/en/reken/zonpositie.html\n// \"The Mean Anomaly\" section for explanation\nconst M0 = 357.5291; // Earth mean anomaly on start day\nconst M1 = 0.98560028; // Earth angle traverses on average per day seen from the sun\nconst THETA0 = 280.147; // The sidereal time (in degrees) at longitude 0\u00B0 at the instant defined by JD2000\nconst THETA1 = 360.9856235; // The rate of change of the sidereal time, in degrees per day.\n/**\n * Calculate sun position\n * based on https://www.aa.quae.nl/en/reken/zonpositie.html\n * inspired by https://github.com/mourner/suncalc/blob/master/suncalc.js\n */\nexport function getSunPosition(timestamp, latitude, longitude) {\n    const longitudeWestInRadians = DEGREES_TO_RADIANS * -longitude;\n    const phi = DEGREES_TO_RADIANS * latitude;\n    const d = toDays(timestamp);\n    const c = getSunCoords(d);\n    // hour angle\n    const H = getSiderealTime(d, longitudeWestInRadians) - c.rightAscension;\n    return {\n        azimuth: getAzimuth(H, phi, c.declination),\n        altitude: getAltitude(H, phi, c.declination)\n    };\n}\nexport function getSunDirection(timestamp, latitude, longitude) {\n    const { azimuth, altitude } = getSunPosition(timestamp, latitude, longitude);\n    // solar position to light direction\n    return [\n        Math.sin(azimuth) * Math.cos(altitude),\n        Math.cos(azimuth) * Math.cos(altitude),\n        -Math.sin(altitude)\n    ];\n}\nfunction toJulianDay(timestamp) {\n    const ts = typeof timestamp === 'number' ? timestamp : timestamp.getTime();\n    return ts / DAY_IN_MS - 0.5 + JD1970;\n}\nfunction toDays(timestamp) {\n    return toJulianDay(timestamp) - JD2000;\n}\nfunction getRightAscension(eclipticLongitude, b) {\n    const lambda = eclipticLongitude;\n    return Math.atan2(Math.sin(lambda) * Math.cos(e) - Math.tan(b) * Math.sin(e), Math.cos(lambda));\n}\nfunction getDeclination(eclipticLongitude, b) {\n    const lambda = eclipticLongitude;\n    return Math.asin(Math.sin(b) * Math.cos(e) + Math.cos(b) * Math.sin(e) * Math.sin(lambda));\n}\nfunction getAzimuth(hourAngle, latitudeInRadians, declination) {\n    const H = hourAngle;\n    const phi = latitudeInRadians;\n    const delta = declination;\n    return Math.atan2(Math.sin(H), Math.cos(H) * Math.sin(phi) - Math.tan(delta) * Math.cos(phi));\n}\nfunction getAltitude(hourAngle, latitudeInRadians, declination) {\n    const H = hourAngle;\n    const phi = latitudeInRadians;\n    const delta = declination;\n    return Math.asin(Math.sin(phi) * Math.sin(delta) + Math.cos(phi) * Math.cos(delta) * Math.cos(H));\n}\n// https://www.aa.quae.nl/en/reken/zonpositie.html\n// \"The Observer section\"\nfunction getSiderealTime(dates, longitudeWestInRadians) {\n    return DEGREES_TO_RADIANS * (THETA0 + THETA1 * dates) - longitudeWestInRadians;\n}\nfunction getSolarMeanAnomaly(days) {\n    return DEGREES_TO_RADIANS * (M0 + M1 * days);\n}\nfunction getEclipticLongitude(meanAnomaly) {\n    const M = meanAnomaly;\n    // equation of center\n    const C = DEGREES_TO_RADIANS * (1.9148 * Math.sin(M) + 0.02 * Math.sin(2 * M) + 0.0003 * Math.sin(3 * M));\n    // perihelion of the Earth\n    const P = DEGREES_TO_RADIANS * 102.9372;\n    return M + C + P + Math.PI;\n}\nfunction getSunCoords(dates) {\n    const M = getSolarMeanAnomaly(dates);\n    const L = getEclipticLongitude(M);\n    return {\n        declination: getDeclination(L, 0),\n        rightAscension: getRightAscension(L, 0)\n    };\n}\n"],
+  "mappings": ";;;;;;;;;;;;;;;;;;;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;;;ACAA,IAAM,qBAAqB,KAAK,KAAK;AACrC,IAAM,YAAY,MAAO,KAAK,KAAK;AACnC,IAAM,SAAS;AACf,IAAM,SAAS;AAEf,IAAM,IAAI,qBAAqB;AAG/B,IAAM,KAAK;AACX,IAAM,KAAK;AACX,IAAM,SAAS;AACf,IAAM,SAAS;AAMR,SAAS,eAAe,WAAW,UAAU,WAAW;AAC3D,QAAM,yBAAyB,qBAAqB,CAAC;AACrD,QAAM,MAAM,qBAAqB;AACjC,QAAM,IAAI,OAAO,SAAS;AAC1B,QAAM,IAAI,aAAa,CAAC;AAExB,QAAM,IAAI,gBAAgB,GAAG,sBAAsB,IAAI,EAAE;AACzD,SAAO;AAAA,IACH,SAAS,WAAW,GAAG,KAAK,EAAE,WAAW;AAAA,IACzC,UAAU,YAAY,GAAG,KAAK,EAAE,WAAW;AAAA,EAC/C;AACJ;AACO,SAAS,gBAAgB,WAAW,UAAU,WAAW;AAC5D,QAAM,EAAE,SAAS,SAAS,IAAI,eAAe,WAAW,UAAU,SAAS;AAE3E,SAAO;AAAA,IACH,KAAK,IAAI,OAAO,IAAI,KAAK,IAAI,QAAQ;AAAA,IACrC,KAAK,IAAI,OAAO,IAAI,KAAK,IAAI,QAAQ;AAAA,IACrC,CAAC,KAAK,IAAI,QAAQ;AAAA,EACtB;AACJ;AACA,SAAS,YAAY,WAAW;AAC5B,QAAM,KAAK,OAAO,cAAc,WAAW,YAAY,UAAU,QAAQ;AACzE,SAAO,KAAK,YAAY,MAAM;AAClC;AACA,SAAS,OAAO,WAAW;AACvB,SAAO,YAAY,SAAS,IAAI;AACpC;AACA,SAAS,kBAAkB,mBAAmB,GAAG;AAC7C,QAAM,SAAS;AACf,SAAO,KAAK,MAAM,KAAK,IAAI,MAAM,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,GAAG,KAAK,IAAI,MAAM,CAAC;AAClG;AACA,SAAS,eAAe,mBAAmB,GAAG;AAC1C,QAAM,SAAS;AACf,SAAO,KAAK,KAAK,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,MAAM,CAAC;AAC7F;AACA,SAAS,WAAW,WAAW,mBAAmB,aAAa;AAC3D,QAAM,IAAI;AACV,QAAM,MAAM;AACZ,QAAM,QAAQ;AACd,SAAO,KAAK,MAAM,KAAK,IAAI,CAAC,GAAG,KAAK,IAAI,CAAC,IAAI,KAAK,IAAI,GAAG,IAAI,KAAK,IAAI,KAAK,IAAI,KAAK,IAAI,GAAG,CAAC;AAChG;AACA,SAAS,YAAY,WAAW,mBAAmB,aAAa;AAC5D,QAAM,IAAI;AACV,QAAM,MAAM;AACZ,QAAM,QAAQ;AACd,SAAO,KAAK,KAAK,KAAK,IAAI,GAAG,IAAI,KAAK,IAAI,KAAK,IAAI,KAAK,IAAI,GAAG,IAAI,KAAK,IAAI,KAAK,IAAI,KAAK,IAAI,CAAC,CAAC;AACpG;AAGA,SAAS,gBAAgB,OAAO,wBAAwB;AACpD,SAAO,sBAAsB,SAAS,SAAS,SAAS;AAC5D;AACA,SAAS,oBAAoB,MAAM;AAC/B,SAAO,sBAAsB,KAAK,KAAK;AAC3C;AACA,SAAS,qBAAqB,aAAa;AACvC,QAAM,IAAI;AAEV,QAAM,IAAI,sBAAsB,SAAS,KAAK,IAAI,CAAC,IAAI,OAAO,KAAK,IAAI,IAAI,CAAC,IAAI,OAAS,KAAK,IAAI,IAAI,CAAC;AAEvG,QAAM,IAAI,qBAAqB;AAC/B,SAAO,IAAI,IAAI,IAAI,KAAK;AAC5B;AACA,SAAS,aAAa,OAAO;AACzB,QAAM,IAAI,oBAAoB,KAAK;AACnC,QAAM,IAAI,qBAAqB,CAAC;AAChC,SAAO;AAAA,IACH,aAAa,eAAe,GAAG,CAAC;AAAA,IAChC,gBAAgB,kBAAkB,GAAG,CAAC;AAAA,EAC1C;AACJ;",
+  "names": []
+}

package/dist/index.d.ts

@@ -1,2 +1 @@
-export { getSunPosition, getSunDirection } from './suncalc';
-//# sourceMappingURL=index.d.ts.map
+export { getSunPosition, getSunDirection } from "./suncalc.js";

package/dist/index.js

@@ -1,2 +1 @@
 export { getSunPosition, getSunDirection } from "./suncalc.js";
-//# sourceMappingURL=index.js.map

package/dist/suncalc.d.ts

@@ -4,7 +4,7 @@
  * The azimuth is the direction along the horizon, which we measure from south to west.
  * South has azimuth 0°, west +90°, north +180°, and east +270° (or −90°, that's the same thing).
  */
-export declare type CelestialPosition = {
+export type CelestialPosition = {
     azimuth: number;
     altitude: number;
 };
@@ -15,4 +15,3 @@ export declare type CelestialPosition = {
  */
 export declare function getSunPosition(timestamp: number | Date, latitude: number, longitude: number): CelestialPosition;
 export declare function getSunDirection(timestamp: number | Date, latitude: number, longitude: number): number[];
-//# sourceMappingURL=suncalc.d.ts.map

package/dist/suncalc.js

@@ -1,85 +1,89 @@
 const DEGREES_TO_RADIANS = Math.PI / 180;
 const DAY_IN_MS = 1000 * 60 * 60 * 24;
-const JD1970 = 2440588;
-const JD2000 = 2451545;
-const e = DEGREES_TO_RADIANS * 23.4397;
-const M0 = 357.5291;
-const M1 = 0.98560028;
-const THETA0 = 280.147;
-const THETA1 = 360.9856235;
+const JD1970 = 2440588; // Julian Day year 1970
+const JD2000 = 2451545; // Julian Day year 2000
+// This angle ε [epsilon] is called the obliquity of the ecliptic and its value at the beginning of 2000 was 23.4397°
+const e = DEGREES_TO_RADIANS * 23.4397; // obliquity of the Earth
+// Refer https://www.aa.quae.nl/en/reken/zonpositie.html
+// "The Mean Anomaly" section for explanation
+const M0 = 357.5291; // Earth mean anomaly on start day
+const M1 = 0.98560028; // Earth angle traverses on average per day seen from the sun
+const THETA0 = 280.147; // The sidereal time (in degrees) at longitude 0° at the instant defined by JD2000
+const THETA1 = 360.9856235; // The rate of change of the sidereal time, in degrees per day.
+/**
+ * Calculate sun position
+ * based on https://www.aa.quae.nl/en/reken/zonpositie.html
+ * inspired by https://github.com/mourner/suncalc/blob/master/suncalc.js
+ */
 export function getSunPosition(timestamp, latitude, longitude) {
-  const longitudeWestInRadians = DEGREES_TO_RADIANS * -longitude;
-  const phi = DEGREES_TO_RADIANS * latitude;
-  const d = toDays(timestamp);
-  const c = getSunCoords(d);
-  const H = getSiderealTime(d, longitudeWestInRadians) - c.rightAscension;
-  return {
-    azimuth: getAzimuth(H, phi, c.declination),
-    altitude: getAltitude(H, phi, c.declination)
-  };
+    const longitudeWestInRadians = DEGREES_TO_RADIANS * -longitude;
+    const phi = DEGREES_TO_RADIANS * latitude;
+    const d = toDays(timestamp);
+    const c = getSunCoords(d);
+    // hour angle
+    const H = getSiderealTime(d, longitudeWestInRadians) - c.rightAscension;
+    return {
+        azimuth: getAzimuth(H, phi, c.declination),
+        altitude: getAltitude(H, phi, c.declination)
+    };
 }
 export function getSunDirection(timestamp, latitude, longitude) {
-  const {
-    azimuth,
-    altitude
-  } = getSunPosition(timestamp, latitude, longitude);
-  return [Math.sin(azimuth) * Math.cos(altitude), Math.cos(azimuth) * Math.cos(altitude), -Math.sin(altitude)];
+    const { azimuth, altitude } = getSunPosition(timestamp, latitude, longitude);
+    // solar position to light direction
+    return [
+        Math.sin(azimuth) * Math.cos(altitude),
+        Math.cos(azimuth) * Math.cos(altitude),
+        -Math.sin(altitude)
+    ];
 }
-
 function toJulianDay(timestamp) {
-  const ts = typeof timestamp === 'number' ? timestamp : timestamp.getTime();
-  return ts / DAY_IN_MS - 0.5 + JD1970;
+    const ts = typeof timestamp === 'number' ? timestamp : timestamp.getTime();
+    return ts / DAY_IN_MS - 0.5 + JD1970;
 }
-
 function toDays(timestamp) {
-  return toJulianDay(timestamp) - JD2000;
+    return toJulianDay(timestamp) - JD2000;
 }
-
 function getRightAscension(eclipticLongitude, b) {
-  const lambda = eclipticLongitude;
-  return Math.atan2(Math.sin(lambda) * Math.cos(e) - Math.tan(b) * Math.sin(e), Math.cos(lambda));
+    const lambda = eclipticLongitude;
+    return Math.atan2(Math.sin(lambda) * Math.cos(e) - Math.tan(b) * Math.sin(e), Math.cos(lambda));
 }
-
 function getDeclination(eclipticLongitude, b) {
-  const lambda = eclipticLongitude;
-  return Math.asin(Math.sin(b) * Math.cos(e) + Math.cos(b) * Math.sin(e) * Math.sin(lambda));
+    const lambda = eclipticLongitude;
+    return Math.asin(Math.sin(b) * Math.cos(e) + Math.cos(b) * Math.sin(e) * Math.sin(lambda));
 }
-
 function getAzimuth(hourAngle, latitudeInRadians, declination) {
-  const H = hourAngle;
-  const phi = latitudeInRadians;
-  const delta = declination;
-  return Math.atan2(Math.sin(H), Math.cos(H) * Math.sin(phi) - Math.tan(delta) * Math.cos(phi));
+    const H = hourAngle;
+    const phi = latitudeInRadians;
+    const delta = declination;
+    return Math.atan2(Math.sin(H), Math.cos(H) * Math.sin(phi) - Math.tan(delta) * Math.cos(phi));
 }
-
 function getAltitude(hourAngle, latitudeInRadians, declination) {
-  const H = hourAngle;
-  const phi = latitudeInRadians;
-  const delta = declination;
-  return Math.asin(Math.sin(phi) * Math.sin(delta) + Math.cos(phi) * Math.cos(delta) * Math.cos(H));
+    const H = hourAngle;
+    const phi = latitudeInRadians;
+    const delta = declination;
+    return Math.asin(Math.sin(phi) * Math.sin(delta) + Math.cos(phi) * Math.cos(delta) * Math.cos(H));
 }
-
+// https://www.aa.quae.nl/en/reken/zonpositie.html
+// "The Observer section"
 function getSiderealTime(dates, longitudeWestInRadians) {
-  return DEGREES_TO_RADIANS * (THETA0 + THETA1 * dates) - longitudeWestInRadians;
+    return DEGREES_TO_RADIANS * (THETA0 + THETA1 * dates) - longitudeWestInRadians;
 }
-
 function getSolarMeanAnomaly(days) {
-  return DEGREES_TO_RADIANS * (M0 + M1 * days);
+    return DEGREES_TO_RADIANS * (M0 + M1 * days);
 }
-
 function getEclipticLongitude(meanAnomaly) {
-  const M = meanAnomaly;
-  const C = DEGREES_TO_RADIANS * (1.9148 * Math.sin(M) + 0.02 * Math.sin(2 * M) + 0.0003 * Math.sin(3 * M));
-  const P = DEGREES_TO_RADIANS * 102.9372;
-  return M + C + P + Math.PI;
+    const M = meanAnomaly;
+    // equation of center
+    const C = DEGREES_TO_RADIANS * (1.9148 * Math.sin(M) + 0.02 * Math.sin(2 * M) + 0.0003 * Math.sin(3 * M));
+    // perihelion of the Earth
+    const P = DEGREES_TO_RADIANS * 102.9372;
+    return M + C + P + Math.PI;
 }
-
 function getSunCoords(dates) {
-  const M = getSolarMeanAnomaly(dates);
-  const L = getEclipticLongitude(M);
-  return {
-    declination: getDeclination(L, 0),
-    rightAscension: getRightAscension(L, 0)
-  };
+    const M = getSolarMeanAnomaly(dates);
+    const L = getEclipticLongitude(M);
+    return {
+        declination: getDeclination(L, 0),
+        rightAscension: getRightAscension(L, 0)
+    };
 }
-//# sourceMappingURL=suncalc.js.map

package/package.json

@@ -6,7 +6,7 @@
   "publishConfig": {
     "access": "public"
   },
-  "version": "4.0.0-beta.1",
+  "version": "4.0.1",
   "keywords": [
     "javascript",
     "math",
@@ -24,17 +24,14 @@
   "module": "dist/index.js",
   "exports": {
     ".": {
+      "types": "./dist/index.d.ts",
       "import": "./dist/index.js",
-      "require": "./dist/index.cjs",
-      "types": "./dist/index.d.ts"
+      "require": "./dist/index.cjs"
     }
   },
   "files": [
     "dist",
     "src"
   ],
-  "dependencies": {
-    "@babel/runtime": "^7.12.0"
-  },
-  "gitHead": "5a2364216b6aadbece690e05b9b5b71753aee472"
+  "gitHead": "33f369ba3a259f79acc3fa8181190c9da8841648"
 }

package/dist/index.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,cAAc,EAAE,eAAe,EAAC,MAAM,WAAW,CAAC"}

package/dist/index.js.map

@@ -1 +0,0 @@
-{"version":3,"sources":["../src/index.ts"],"names":["getSunPosition","getSunDirection"],"mappings":"SAAQA,c,EAAgBC,e","sourcesContent":["export {getSunPosition, getSunDirection} from './suncalc';\n"],"file":"index.js"}

package/dist/suncalc.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"suncalc.d.ts","sourceRoot":"","sources":["../src/suncalc.ts"],"names":[],"mappings":"AAiBA;;;;;GAKG;AACH,oBAAY,iBAAiB,GAAG;IAC9B,OAAO,EAAE,MAAM,CAAC;IAChB,QAAQ,EAAE,MAAM,CAAC;CAClB,CAAC;AAEF;;;;GAIG;AACH,wBAAgB,cAAc,CAC5B,SAAS,EAAE,MAAM,GAAG,IAAI,EACxB,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,GAChB,iBAAiB,CAanB;AAED,wBAAgB,eAAe,CAC7B,SAAS,EAAE,MAAM,GAAG,IAAI,EACxB,QAAQ,EAAE,MAAM,EAChB,SAAS,EAAE,MAAM,GAChB,MAAM,EAAE,CASV"}

package/dist/suncalc.js.map

@@ -1 +0,0 @@
-{"version":3,"sources":["../src/suncalc.ts"],"names":["DEGREES_TO_RADIANS","Math","PI","DAY_IN_MS","JD1970","JD2000","e","M0","M1","THETA0","THETA1","getSunPosition","timestamp","latitude","longitude","longitudeWestInRadians","phi","d","toDays","c","getSunCoords","H","getSiderealTime","rightAscension","azimuth","getAzimuth","declination","altitude","getAltitude","getSunDirection","sin","cos","toJulianDay","ts","getTime","getRightAscension","eclipticLongitude","b","lambda","atan2","tan","getDeclination","asin","hourAngle","latitudeInRadians","delta","dates","getSolarMeanAnomaly","days","getEclipticLongitude","meanAnomaly","M","C","P","L"],"mappings":"AAAA,MAAMA,kBAAkB,GAAGC,IAAI,CAACC,EAAL,GAAU,GAArC;AAEA,MAAMC,SAAS,GAAG,OAAO,EAAP,GAAY,EAAZ,GAAiB,EAAnC;AACA,MAAMC,MAAM,GAAG,OAAf;AACA,MAAMC,MAAM,GAAG,OAAf;AAGA,MAAMC,CAAC,GAAGN,kBAAkB,GAAG,OAA/B;AAIA,MAAMO,EAAE,GAAG,QAAX;AACA,MAAMC,EAAE,GAAG,UAAX;AAEA,MAAMC,MAAM,GAAG,OAAf;AACA,MAAMC,MAAM,GAAG,WAAf;AAkBA,OAAO,SAASC,cAAT,CACLC,SADK,EAELC,QAFK,EAGLC,SAHK,EAIc;AACnB,QAAMC,sBAAsB,GAAGf,kBAAkB,GAAG,CAACc,SAArD;AACA,QAAME,GAAG,GAAGhB,kBAAkB,GAAGa,QAAjC;AACA,QAAMI,CAAC,GAAGC,MAAM,CAACN,SAAD,CAAhB;AAEA,QAAMO,CAAC,GAAGC,YAAY,CAACH,CAAD,CAAtB;AAEA,QAAMI,CAAC,GAAGC,eAAe,CAACL,CAAD,EAAIF,sBAAJ,CAAf,GAA6CI,CAAC,CAACI,cAAzD;AAEA,SAAO;AACLC,IAAAA,OAAO,EAAEC,UAAU,CAACJ,CAAD,EAAIL,GAAJ,EAASG,CAAC,CAACO,WAAX,CADd;AAELC,IAAAA,QAAQ,EAAEC,WAAW,CAACP,CAAD,EAAIL,GAAJ,EAASG,CAAC,CAACO,WAAX;AAFhB,GAAP;AAID;AAED,OAAO,SAASG,eAAT,CACLjB,SADK,EAELC,QAFK,EAGLC,SAHK,EAIK;AACV,QAAM;AAACU,IAAAA,OAAD;AAAUG,IAAAA;AAAV,MAAsBhB,cAAc,CAACC,SAAD,EAAYC,QAAZ,EAAsBC,SAAtB,CAA1C;AAGA,SAAO,CACLb,IAAI,CAAC6B,GAAL,CAASN,OAAT,IAAoBvB,IAAI,CAAC8B,GAAL,CAASJ,QAAT,CADf,EAEL1B,IAAI,CAAC8B,GAAL,CAASP,OAAT,IAAoBvB,IAAI,CAAC8B,GAAL,CAASJ,QAAT,CAFf,EAGL,CAAC1B,IAAI,CAAC6B,GAAL,CAASH,QAAT,CAHI,CAAP;AAKD;;AAED,SAASK,WAAT,CAAqBpB,SAArB,EAAuD;AACrD,QAAMqB,EAAE,GAAG,OAAOrB,SAAP,KAAqB,QAArB,GAAgCA,SAAhC,GAA4CA,SAAS,CAACsB,OAAV,EAAvD;AACA,SAAOD,EAAE,GAAG9B,SAAL,GAAiB,GAAjB,GAAuBC,MAA9B;AACD;;AAED,SAASc,MAAT,CAAgBN,SAAhB,EAAkD;AAChD,SAAOoB,WAAW,CAACpB,SAAD,CAAX,GAAyBP,MAAhC;AACD;;AAED,SAAS8B,iBAAT,CAA2BC,iBAA3B,EAAsDC,CAAtD,EAAyE;AACvE,QAAMC,MAAM,GAAGF,iBAAf;AACA,SAAOnC,IAAI,CAACsC,KAAL,CAAWtC,IAAI,CAAC6B,GAAL,CAASQ,MAAT,IAAmBrC,IAAI,CAAC8B,GAAL,CAASzB,CAAT,CAAnB,GAAiCL,IAAI,CAACuC,GAAL,CAASH,CAAT,IAAcpC,IAAI,CAAC6B,GAAL,CAASxB,CAAT,CAA1D,EAAuEL,IAAI,CAAC8B,GAAL,CAASO,MAAT,CAAvE,CAAP;AACD;;AAED,SAASG,cAAT,CAAwBL,iBAAxB,EAAmDC,CAAnD,EAAsE;AACpE,QAAMC,MAAM,GAAGF,iBAAf;AACA,SAAOnC,IAAI,CAACyC,IAAL,CAAUzC,IAAI,CAAC6B,GAAL,CAASO,CAAT,IAAcpC,IAAI,CAAC8B,GAAL,CAASzB,CAAT,CAAd,GAA4BL,IAAI,CAAC8B,GAAL,CAASM,CAAT,IAAcpC,IAAI,CAAC6B,GAAL,CAASxB,CAAT,CAAd,GAA4BL,IAAI,CAAC6B,GAAL,CAASQ,MAAT,CAAlE,CAAP;AACD;;AAED,SAASb,UAAT,CAAoBkB,SAApB,EAAuCC,iBAAvC,EAAkElB,WAAlE,EAA+F;AAC7F,QAAML,CAAC,GAAGsB,SAAV;AACA,QAAM3B,GAAG,GAAG4B,iBAAZ;AACA,QAAMC,KAAK,GAAGnB,WAAd;AACA,SAAOzB,IAAI,CAACsC,KAAL,CAAWtC,IAAI,CAAC6B,GAAL,CAAST,CAAT,CAAX,EAAwBpB,IAAI,CAAC8B,GAAL,CAASV,CAAT,IAAcpB,IAAI,CAAC6B,GAAL,CAASd,GAAT,CAAd,GAA8Bf,IAAI,CAACuC,GAAL,CAASK,KAAT,IAAkB5C,IAAI,CAAC8B,GAAL,CAASf,GAAT,CAAxE,CAAP;AACD;;AAED,SAASY,WAAT,CAAqBe,SAArB,EAAwCC,iBAAxC,EAAmElB,WAAnE,EAAgG;AAC9F,QAAML,CAAC,GAAGsB,SAAV;AACA,QAAM3B,GAAG,GAAG4B,iBAAZ;AACA,QAAMC,KAAK,GAAGnB,WAAd;AACA,SAAOzB,IAAI,CAACyC,IAAL,CAAUzC,IAAI,CAAC6B,GAAL,CAASd,GAAT,IAAgBf,IAAI,CAAC6B,GAAL,CAASe,KAAT,CAAhB,GAAkC5C,IAAI,CAAC8B,GAAL,CAASf,GAAT,IAAgBf,IAAI,CAAC8B,GAAL,CAASc,KAAT,CAAhB,GAAkC5C,IAAI,CAAC8B,GAAL,CAASV,CAAT,CAA9E,CAAP;AACD;;AAID,SAASC,eAAT,CAAyBwB,KAAzB,EAAwC/B,sBAAxC,EAAgF;AAC9E,SAAOf,kBAAkB,IAAIS,MAAM,GAAGC,MAAM,GAAGoC,KAAtB,CAAlB,GAAiD/B,sBAAxD;AACD;;AAED,SAASgC,mBAAT,CAA6BC,IAA7B,EAAmD;AACjD,SAAOhD,kBAAkB,IAAIO,EAAE,GAAGC,EAAE,GAAGwC,IAAd,CAAzB;AACD;;AAED,SAASC,oBAAT,CAA8BC,WAA9B,EAA2D;AACzD,QAAMC,CAAC,GAAGD,WAAV;AAEA,QAAME,CAAC,GACLpD,kBAAkB,IAAI,SAASC,IAAI,CAAC6B,GAAL,CAASqB,CAAT,CAAT,GAAuB,OAAOlD,IAAI,CAAC6B,GAAL,CAAS,IAAIqB,CAAb,CAA9B,GAAgD,SAASlD,IAAI,CAAC6B,GAAL,CAAS,IAAIqB,CAAb,CAA7D,CADpB;AAGA,QAAME,CAAC,GAAGrD,kBAAkB,GAAG,QAA/B;AAEA,SAAOmD,CAAC,GAAGC,CAAJ,GAAQC,CAAR,GAAYpD,IAAI,CAACC,EAAxB;AACD;;AAED,SAASkB,YAAT,CAAsB0B,KAAtB,EAGE;AACA,QAAMK,CAAC,GAAGJ,mBAAmB,CAACD,KAAD,CAA7B;AACA,QAAMQ,CAAC,GAAGL,oBAAoB,CAACE,CAAD,CAA9B;AAEA,SAAO;AACLzB,IAAAA,WAAW,EAAEe,cAAc,CAACa,CAAD,EAAI,CAAJ,CADtB;AAEL/B,IAAAA,cAAc,EAAEY,iBAAiB,CAACmB,CAAD,EAAI,CAAJ;AAF5B,GAAP;AAID","sourcesContent":["const DEGREES_TO_RADIANS = Math.PI / 180;\n\nconst DAY_IN_MS = 1000 * 60 * 60 * 24;\nconst JD1970 = 2440588; // Julian Day year 1970\nconst JD2000 = 2451545; // Julian Day year 2000\n\n// This angle ε [epsilon] is called the obliquity of the ecliptic and its value at the beginning of 2000 was 23.4397°\nconst e = DEGREES_TO_RADIANS * 23.4397; // obliquity of the Earth\n\n// Refer https://www.aa.quae.nl/en/reken/zonpositie.html\n// \"The Mean Anomaly\" section for explanation\nconst M0 = 357.5291; // Earth mean anomaly on start day\nconst M1 = 0.98560028; // Earth angle traverses on average per day seen from the sun\n\nconst THETA0 = 280.147; // The sidereal time (in degrees) at longitude 0° at the instant defined by JD2000\nconst THETA1 = 360.9856235; // The rate of change of the sidereal time, in degrees per day.\n\n/**\n * A position in the sky defined by two angles\n * The altitude is 0° at the horizon, +90° in the zenith (straight over your head), and −90° in the nadir (straight down).\n * The azimuth is the direction along the horizon, which we measure from south to west.\n * South has azimuth 0°, west +90°, north +180°, and east +270° (or −90°, that's the same thing).\n */\nexport type CelestialPosition = {\n  azimuth: number;\n  altitude: number;\n};\n\n/**\n * Calculate sun position\n * based on https://www.aa.quae.nl/en/reken/zonpositie.html\n * inspired by https://github.com/mourner/suncalc/blob/master/suncalc.js\n */\nexport function getSunPosition(\n  timestamp: number | Date,\n  latitude: number,\n  longitude: number\n): CelestialPosition {\n  const longitudeWestInRadians = DEGREES_TO_RADIANS * -longitude;\n  const phi = DEGREES_TO_RADIANS * latitude;\n  const d = toDays(timestamp);\n\n  const c = getSunCoords(d);\n  // hour angle\n  const H = getSiderealTime(d, longitudeWestInRadians) - c.rightAscension;\n\n  return {\n    azimuth: getAzimuth(H, phi, c.declination),\n    altitude: getAltitude(H, phi, c.declination)\n  };\n}\n\nexport function getSunDirection(\n  timestamp: number | Date,\n  latitude: number,\n  longitude: number\n): number[] {\n  const {azimuth, altitude} = getSunPosition(timestamp, latitude, longitude);\n\n  // solar position to light direction\n  return [\n    Math.sin(azimuth) * Math.cos(altitude),\n    Math.cos(azimuth) * Math.cos(altitude),\n    -Math.sin(altitude)\n  ];\n}\n\nfunction toJulianDay(timestamp: number | Date): number {\n  const ts = typeof timestamp === 'number' ? timestamp : timestamp.getTime();\n  return ts / DAY_IN_MS - 0.5 + JD1970;\n}\n\nfunction toDays(timestamp: number | Date): number {\n  return toJulianDay(timestamp) - JD2000;\n}\n\nfunction getRightAscension(eclipticLongitude: number, b: number): number {\n  const lambda = eclipticLongitude;\n  return Math.atan2(Math.sin(lambda) * Math.cos(e) - Math.tan(b) * Math.sin(e), Math.cos(lambda));\n}\n\nfunction getDeclination(eclipticLongitude: number, b: number): number {\n  const lambda = eclipticLongitude;\n  return Math.asin(Math.sin(b) * Math.cos(e) + Math.cos(b) * Math.sin(e) * Math.sin(lambda));\n}\n\nfunction getAzimuth(hourAngle: number, latitudeInRadians: number, declination: number): number {\n  const H = hourAngle;\n  const phi = latitudeInRadians;\n  const delta = declination;\n  return Math.atan2(Math.sin(H), Math.cos(H) * Math.sin(phi) - Math.tan(delta) * Math.cos(phi));\n}\n\nfunction getAltitude(hourAngle: number, latitudeInRadians: number, declination: number): number {\n  const H = hourAngle;\n  const phi = latitudeInRadians;\n  const delta = declination;\n  return Math.asin(Math.sin(phi) * Math.sin(delta) + Math.cos(phi) * Math.cos(delta) * Math.cos(H));\n}\n\n// https://www.aa.quae.nl/en/reken/zonpositie.html\n// \"The Observer section\"\nfunction getSiderealTime(dates: number, longitudeWestInRadians: number): number {\n  return DEGREES_TO_RADIANS * (THETA0 + THETA1 * dates) - longitudeWestInRadians;\n}\n\nfunction getSolarMeanAnomaly(days: number): number {\n  return DEGREES_TO_RADIANS * (M0 + M1 * days);\n}\n\nfunction getEclipticLongitude(meanAnomaly: number): number {\n  const M = meanAnomaly;\n  // equation of center\n  const C =\n    DEGREES_TO_RADIANS * (1.9148 * Math.sin(M) + 0.02 * Math.sin(2 * M) + 0.0003 * Math.sin(3 * M));\n  // perihelion of the Earth\n  const P = DEGREES_TO_RADIANS * 102.9372;\n\n  return M + C + P + Math.PI;\n}\n\nfunction getSunCoords(dates: number): {\n  declination: number;\n  rightAscension: number;\n} {\n  const M = getSolarMeanAnomaly(dates);\n  const L = getEclipticLongitude(M);\n\n  return {\n    declination: getDeclination(L, 0),\n    rightAscension: getRightAscension(L, 0)\n  };\n}\n"],"file":"suncalc.js"}