@itwin/core-common 4.0.0-dev.46 → 4.0.0-dev.48

package/lib/esm/SolarCalculate.d.ts

@@ -1,28 +1,28 @@
-/** @packageDocumentation
- * @module Utils
- */
-import { Vector3d } from "@itwin/core-geometry";
-import { Cartographic } from "./geometry/Cartographic";
-/** @public
- * calculate solar angles (in radians) based at a given date/time and cartographic location.
- */
-export declare function calculateSolarAngles(date: Date, location: Cartographic): {
-    azimuth: number;
-    elevation: number;
-};
-/** @public
- * calculate solar direction based at a given date/time and cartpgrphic location.
- */
-export declare function calculateSolarDirection(date: Date, location: Cartographic): Vector3d;
-/** @public
- * calculate solar direction corresponding to the given azimuth and elevation (altitude) angles in degrees.
- */
-export declare function calculateSolarDirectionFromAngles(azimuthElevation: {
-    azimuth: number;
-    elevation: number;
-}): Vector3d;
-/** @public
- * calculate solar sunrise or sunset for a given day and cartographic location.
- */
-export declare function calculateSunriseOrSunset(date: Date, location: Cartographic, sunrise: boolean): Date;
+/** @packageDocumentation
+ * @module Utils
+ */
+import { Vector3d } from "@itwin/core-geometry";
+import { Cartographic } from "./geometry/Cartographic";
+/** @public
+ * calculate solar angles (in radians) based at a given date/time and cartographic location.
+ */
+export declare function calculateSolarAngles(date: Date, location: Cartographic): {
+    azimuth: number;
+    elevation: number;
+};
+/** @public
+ * calculate solar direction based at a given date/time and cartpgrphic location.
+ */
+export declare function calculateSolarDirection(date: Date, location: Cartographic): Vector3d;
+/** @public
+ * calculate solar direction corresponding to the given azimuth and elevation (altitude) angles in degrees.
+ */
+export declare function calculateSolarDirectionFromAngles(azimuthElevation: {
+    azimuth: number;
+    elevation: number;
+}): Vector3d;
+/** @public
+ * calculate solar sunrise or sunset for a given day and cartographic location.
+ */
+export declare function calculateSunriseOrSunset(date: Date, location: Cartographic, sunrise: boolean): Date;
 //# sourceMappingURL=SolarCalculate.d.ts.map

package/lib/esm/SolarCalculate.js

@@ -1,205 +1,205 @@
-/*---------------------------------------------------------------------------------------------
-* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
-* See LICENSE.md in the project root for license terms and full copyright notice.
-*--------------------------------------------------------------------------------------------*/
-/** @packageDocumentation
- * @module Utils
- */
-import { Angle, Vector3d } from "@itwin/core-geometry";
-// cspell:ignore mrad sinm sint aarg
-// Code below loosely translated from  https://www.esrl.noaa.gov/gmd/grad/solcalc/
-function calcTimeJulianCent(jd) {
-    const T = (jd - 2451545.0) / 36525.0;
-    return T;
-}
-function radToDeg(angleRad) {
-    return (180.0 * angleRad / Math.PI);
-}
-function degToRad(angleDeg) {
-    return (Math.PI * angleDeg / 180.0);
-}
-function calcGeomMeanLongSun(t) {
-    let L0 = 280.46646 + t * (36000.76983 + t * (0.0003032));
-    while (L0 > 360.0) {
-        L0 -= 360.0;
-    }
-    while (L0 < 0.0) {
-        L0 += 360.0;
-    }
-    return L0; // in degrees
-}
-function calcGeomMeanAnomalySun(t) {
-    const M = 357.52911 + t * (35999.05029 - 0.0001537 * t);
-    return M; // in degrees
-}
-function calcEccentricityEarthOrbit(t) {
-    const e = 0.016708634 - t * (0.000042037 + 0.0000001267 * t);
-    return e; // unitless
-}
-function calcSunEqOfCenter(t) {
-    const m = calcGeomMeanAnomalySun(t);
-    const mrad = degToRad(m);
-    const sinm = Math.sin(mrad);
-    const sin2m = Math.sin(mrad + mrad);
-    const sin3m = Math.sin(mrad + mrad + mrad);
-    const C = sinm * (1.914602 - t * (0.004817 + 0.000014 * t)) + sin2m * (0.019993 - 0.000101 * t) + sin3m * 0.000289;
-    return C; // in degrees
-}
-function calcSunTrueLong(t) {
-    const l0 = calcGeomMeanLongSun(t);
-    const c = calcSunEqOfCenter(t);
-    const O = l0 + c;
-    return O; // in degrees
-}
-function calcSunApparentLong(t) {
-    const o = calcSunTrueLong(t);
-    const omega = 125.04 - 1934.136 * t;
-    const lambda = o - 0.00569 - 0.00478 * Math.sin(degToRad(omega));
-    return lambda; // in degrees
-}
-function calcMeanObliquityOfEcliptic(t) {
-    const seconds = 21.448 - t * (46.8150 + t * (0.00059 - t * (0.001813)));
-    const e0 = 23.0 + (26.0 + (seconds / 60.0)) / 60.0;
-    return e0; // in degrees
-}
-function calcObliquityCorrection(t) {
-    const e0 = calcMeanObliquityOfEcliptic(t);
-    const omega = 125.04 - 1934.136 * t;
-    const e = e0 + 0.00256 * Math.cos(degToRad(omega));
-    return e; // in degrees
-}
-function calcSunDeclination(t) {
-    const e = calcObliquityCorrection(t);
-    const lambda = calcSunApparentLong(t);
-    const sint = Math.sin(degToRad(e)) * Math.sin(degToRad(lambda));
-    const theta = radToDeg(Math.asin(sint));
-    return theta; // in degrees
-}
-function calcEquationOfTime(t) {
-    const epsilon = calcObliquityCorrection(t);
-    const l0 = calcGeomMeanLongSun(t);
-    const e = calcEccentricityEarthOrbit(t);
-    const m = calcGeomMeanAnomalySun(t);
-    let y = Math.tan(degToRad(epsilon) / 2.0);
-    y *= y;
-    const sin2l0 = Math.sin(2.0 * degToRad(l0));
-    const sinm = Math.sin(degToRad(m));
-    const cos2l0 = Math.cos(2.0 * degToRad(l0));
-    const sin4l0 = Math.sin(4.0 * degToRad(l0));
-    const sin2m = Math.sin(2.0 * degToRad(m));
-    const eTime = y * sin2l0 - 2.0 * e * sinm + 4.0 * e * y * sinm * cos2l0 - 0.5 * y * y * sin4l0 - 1.25 * e * e * sin2m;
-    return radToDeg(eTime) * 4.0; // in minutes of time
-}
-function calcAzEl(t, localTime, latitude, longitude, zone) {
-    const eqTime = calcEquationOfTime(t);
-    const theta = calcSunDeclination(t);
-    const solarTimeFix = eqTime + 4.0 * longitude - 60.0 * zone;
-    let trueSolarTime = localTime + solarTimeFix;
-    while (trueSolarTime > 1440)
-        trueSolarTime -= 1440;
-    let hourAngle = trueSolarTime / 4.0 - 180.0;
-    if (hourAngle < -180) {
-        hourAngle += 360.0;
-    }
-    const haRad = degToRad(hourAngle);
-    let csz = Math.sin(degToRad(latitude)) * Math.sin(degToRad(theta)) + Math.cos(degToRad(latitude)) * Math.cos(degToRad(theta)) * Math.cos(haRad);
-    if (csz > 1.0) {
-        csz = 1.0;
-    }
-    else if (csz < -1.0) {
-        csz = -1.0;
-    }
-    const zenith = radToDeg(Math.acos(csz));
-    const azDenom = (Math.cos(degToRad(latitude)) * Math.sin(degToRad(zenith)));
-    let azimuth;
-    if (Math.abs(azDenom) > 0.001) {
-        let azRad = ((Math.sin(degToRad(latitude)) * Math.cos(degToRad(zenith))) - Math.sin(degToRad(theta))) / azDenom;
-        if (Math.abs(azRad) > 1.0) {
-            if (azRad < 0) {
-                azRad = -1.0;
-            }
-            else {
-                azRad = 1.0;
-            }
-        }
-        azimuth = 180.0 - radToDeg(Math.acos(azRad));
-        if (hourAngle > 0.0) {
-            azimuth = -azimuth;
-        }
-    }
-    else {
-        if (latitude > 0.0) {
-            azimuth = 180.0;
-        }
-        else {
-            azimuth = 0.0;
-        }
-    }
-    if (azimuth < 0.0) {
-        azimuth += 360.0;
-    }
-    return { azimuth, elevation: 90 - zenith };
-}
-function calculateJulianDay(date) {
-    return Math.floor(date.getTime() / 86400000) + 2440587.5; // https://stackoverflow.com/questions/11759992/calculating-jdayjulian-day-in-javascript
-}
-/** @public
- * calculate solar angles (in radians) based at a given date/time and cartographic location.
- */
-export function calculateSolarAngles(date, location) {
-    const jDay = calculateJulianDay(date);
-    const latitude = location.latitudeDegrees;
-    const longitude = location.longitudeDegrees;
-    const utcMinutes = date.getUTCHours() * 60 + date.getUTCMinutes() + date.getUTCSeconds() / 60;
-    const zone = Math.floor(.5 + longitude / 15.0); // date.getTimeZoneOffset mixes in DST.
-    const localMinutes = utcMinutes + zone * 60;
-    const jTotal = jDay + utcMinutes / 1440.0;
-    const T = calcTimeJulianCent(jTotal);
-    return calcAzEl(T, localMinutes, latitude, longitude, zone);
-}
-/** @public
- * calculate solar direction based at a given date/time and cartpgrphic location.
- */
-export function calculateSolarDirection(date, location) {
-    return calculateSolarDirectionFromAngles(calculateSolarAngles(date, location));
-}
-/** @public
- * calculate solar direction corresponding to the given azimuth and elevation (altitude) angles in degrees.
- */
-export function calculateSolarDirectionFromAngles(azimuthElevation) {
-    const azimuth = Angle.degreesToRadians(azimuthElevation.azimuth);
-    const elevation = Angle.degreesToRadians(azimuthElevation.elevation);
-    const cosElevation = Math.cos(elevation);
-    const sinElevation = Math.sin(elevation);
-    return Vector3d.create(-Math.sin(azimuth) * cosElevation, -Math.cos(azimuth) * cosElevation, -sinElevation);
-}
-function dateFromUtcMinutes(date, utcMinutes) {
-    const utcHours = Math.floor(utcMinutes / 60.0);
-    const output = new Date(date);
-    output.setUTCHours(utcHours);
-    output.setUTCMinutes(Math.floor(.5 + utcMinutes - 60.0 * utcHours));
-    output.setUTCSeconds(0);
-    return output;
-}
-function calcSunriseUtcMinutes(rise, lat, longitude, jDay) {
-    const t = calcTimeJulianCent(jDay);
-    const eqTime = calcEquationOfTime(t);
-    const solarDec = calcSunDeclination(t);
-    const latRad = degToRad(lat);
-    const sdRad = degToRad(solarDec);
-    const hAarg = (Math.cos(degToRad(90.833)) / (Math.cos(latRad) * Math.cos(sdRad)) - Math.tan(latRad) * Math.tan(sdRad));
-    const hourAngle = Math.acos(hAarg);
-    const delta = longitude + radToDeg(rise ? hourAngle : -hourAngle);
-    return 720 - (4.0 * delta) - eqTime; // in UTC minutes
-}
-/** @public
- * calculate solar sunrise or sunset for a given day and cartographic location.
- */
-export function calculateSunriseOrSunset(date, location, sunrise) {
-    const jDay = calculateJulianDay(date);
-    const longitude = location.longitudeDegrees;
-    const latitude = location.latitudeDegrees;
-    const utcMinutes = calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay);
-    return sunrise ? dateFromUtcMinutes(date, utcMinutes) : dateFromUtcMinutes(date, calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay + utcMinutes / 1440));
-}
+/*---------------------------------------------------------------------------------------------
+* Copyright (c) Bentley Systems, Incorporated. All rights reserved.
+* See LICENSE.md in the project root for license terms and full copyright notice.
+*--------------------------------------------------------------------------------------------*/
+/** @packageDocumentation
+ * @module Utils
+ */
+import { Angle, Vector3d } from "@itwin/core-geometry";
+// cspell:ignore mrad sinm sint aarg
+// Code below loosely translated from  https://www.esrl.noaa.gov/gmd/grad/solcalc/
+function calcTimeJulianCent(jd) {
+    const T = (jd - 2451545.0) / 36525.0;
+    return T;
+}
+function radToDeg(angleRad) {
+    return (180.0 * angleRad / Math.PI);
+}
+function degToRad(angleDeg) {
+    return (Math.PI * angleDeg / 180.0);
+}
+function calcGeomMeanLongSun(t) {
+    let L0 = 280.46646 + t * (36000.76983 + t * (0.0003032));
+    while (L0 > 360.0) {
+        L0 -= 360.0;
+    }
+    while (L0 < 0.0) {
+        L0 += 360.0;
+    }
+    return L0; // in degrees
+}
+function calcGeomMeanAnomalySun(t) {
+    const M = 357.52911 + t * (35999.05029 - 0.0001537 * t);
+    return M; // in degrees
+}
+function calcEccentricityEarthOrbit(t) {
+    const e = 0.016708634 - t * (0.000042037 + 0.0000001267 * t);
+    return e; // unitless
+}
+function calcSunEqOfCenter(t) {
+    const m = calcGeomMeanAnomalySun(t);
+    const mrad = degToRad(m);
+    const sinm = Math.sin(mrad);
+    const sin2m = Math.sin(mrad + mrad);
+    const sin3m = Math.sin(mrad + mrad + mrad);
+    const C = sinm * (1.914602 - t * (0.004817 + 0.000014 * t)) + sin2m * (0.019993 - 0.000101 * t) + sin3m * 0.000289;
+    return C; // in degrees
+}
+function calcSunTrueLong(t) {
+    const l0 = calcGeomMeanLongSun(t);
+    const c = calcSunEqOfCenter(t);
+    const O = l0 + c;
+    return O; // in degrees
+}
+function calcSunApparentLong(t) {
+    const o = calcSunTrueLong(t);
+    const omega = 125.04 - 1934.136 * t;
+    const lambda = o - 0.00569 - 0.00478 * Math.sin(degToRad(omega));
+    return lambda; // in degrees
+}
+function calcMeanObliquityOfEcliptic(t) {
+    const seconds = 21.448 - t * (46.8150 + t * (0.00059 - t * (0.001813)));
+    const e0 = 23.0 + (26.0 + (seconds / 60.0)) / 60.0;
+    return e0; // in degrees
+}
+function calcObliquityCorrection(t) {
+    const e0 = calcMeanObliquityOfEcliptic(t);
+    const omega = 125.04 - 1934.136 * t;
+    const e = e0 + 0.00256 * Math.cos(degToRad(omega));
+    return e; // in degrees
+}
+function calcSunDeclination(t) {
+    const e = calcObliquityCorrection(t);
+    const lambda = calcSunApparentLong(t);
+    const sint = Math.sin(degToRad(e)) * Math.sin(degToRad(lambda));
+    const theta = radToDeg(Math.asin(sint));
+    return theta; // in degrees
+}
+function calcEquationOfTime(t) {
+    const epsilon = calcObliquityCorrection(t);
+    const l0 = calcGeomMeanLongSun(t);
+    const e = calcEccentricityEarthOrbit(t);
+    const m = calcGeomMeanAnomalySun(t);
+    let y = Math.tan(degToRad(epsilon) / 2.0);
+    y *= y;
+    const sin2l0 = Math.sin(2.0 * degToRad(l0));
+    const sinm = Math.sin(degToRad(m));
+    const cos2l0 = Math.cos(2.0 * degToRad(l0));
+    const sin4l0 = Math.sin(4.0 * degToRad(l0));
+    const sin2m = Math.sin(2.0 * degToRad(m));
+    const eTime = y * sin2l0 - 2.0 * e * sinm + 4.0 * e * y * sinm * cos2l0 - 0.5 * y * y * sin4l0 - 1.25 * e * e * sin2m;
+    return radToDeg(eTime) * 4.0; // in minutes of time
+}
+function calcAzEl(t, localTime, latitude, longitude, zone) {
+    const eqTime = calcEquationOfTime(t);
+    const theta = calcSunDeclination(t);
+    const solarTimeFix = eqTime + 4.0 * longitude - 60.0 * zone;
+    let trueSolarTime = localTime + solarTimeFix;
+    while (trueSolarTime > 1440)
+        trueSolarTime -= 1440;
+    let hourAngle = trueSolarTime / 4.0 - 180.0;
+    if (hourAngle < -180) {
+        hourAngle += 360.0;
+    }
+    const haRad = degToRad(hourAngle);
+    let csz = Math.sin(degToRad(latitude)) * Math.sin(degToRad(theta)) + Math.cos(degToRad(latitude)) * Math.cos(degToRad(theta)) * Math.cos(haRad);
+    if (csz > 1.0) {
+        csz = 1.0;
+    }
+    else if (csz < -1.0) {
+        csz = -1.0;
+    }
+    const zenith = radToDeg(Math.acos(csz));
+    const azDenom = (Math.cos(degToRad(latitude)) * Math.sin(degToRad(zenith)));
+    let azimuth;
+    if (Math.abs(azDenom) > 0.001) {
+        let azRad = ((Math.sin(degToRad(latitude)) * Math.cos(degToRad(zenith))) - Math.sin(degToRad(theta))) / azDenom;
+        if (Math.abs(azRad) > 1.0) {
+            if (azRad < 0) {
+                azRad = -1.0;
+            }
+            else {
+                azRad = 1.0;
+            }
+        }
+        azimuth = 180.0 - radToDeg(Math.acos(azRad));
+        if (hourAngle > 0.0) {
+            azimuth = -azimuth;
+        }
+    }
+    else {
+        if (latitude > 0.0) {
+            azimuth = 180.0;
+        }
+        else {
+            azimuth = 0.0;
+        }
+    }
+    if (azimuth < 0.0) {
+        azimuth += 360.0;
+    }
+    return { azimuth, elevation: 90 - zenith };
+}
+function calculateJulianDay(date) {
+    return Math.floor(date.getTime() / 86400000) + 2440587.5; // https://stackoverflow.com/questions/11759992/calculating-jdayjulian-day-in-javascript
+}
+/** @public
+ * calculate solar angles (in radians) based at a given date/time and cartographic location.
+ */
+export function calculateSolarAngles(date, location) {
+    const jDay = calculateJulianDay(date);
+    const latitude = location.latitudeDegrees;
+    const longitude = location.longitudeDegrees;
+    const utcMinutes = date.getUTCHours() * 60 + date.getUTCMinutes() + date.getUTCSeconds() / 60;
+    const zone = Math.floor(.5 + longitude / 15.0); // date.getTimeZoneOffset mixes in DST.
+    const localMinutes = utcMinutes + zone * 60;
+    const jTotal = jDay + utcMinutes / 1440.0;
+    const T = calcTimeJulianCent(jTotal);
+    return calcAzEl(T, localMinutes, latitude, longitude, zone);
+}
+/** @public
+ * calculate solar direction based at a given date/time and cartpgrphic location.
+ */
+export function calculateSolarDirection(date, location) {
+    return calculateSolarDirectionFromAngles(calculateSolarAngles(date, location));
+}
+/** @public
+ * calculate solar direction corresponding to the given azimuth and elevation (altitude) angles in degrees.
+ */
+export function calculateSolarDirectionFromAngles(azimuthElevation) {
+    const azimuth = Angle.degreesToRadians(azimuthElevation.azimuth);
+    const elevation = Angle.degreesToRadians(azimuthElevation.elevation);
+    const cosElevation = Math.cos(elevation);
+    const sinElevation = Math.sin(elevation);
+    return Vector3d.create(-Math.sin(azimuth) * cosElevation, -Math.cos(azimuth) * cosElevation, -sinElevation);
+}
+function dateFromUtcMinutes(date, utcMinutes) {
+    const utcHours = Math.floor(utcMinutes / 60.0);
+    const output = new Date(date);
+    output.setUTCHours(utcHours);
+    output.setUTCMinutes(Math.floor(.5 + utcMinutes - 60.0 * utcHours));
+    output.setUTCSeconds(0);
+    return output;
+}
+function calcSunriseUtcMinutes(rise, lat, longitude, jDay) {
+    const t = calcTimeJulianCent(jDay);
+    const eqTime = calcEquationOfTime(t);
+    const solarDec = calcSunDeclination(t);
+    const latRad = degToRad(lat);
+    const sdRad = degToRad(solarDec);
+    const hAarg = (Math.cos(degToRad(90.833)) / (Math.cos(latRad) * Math.cos(sdRad)) - Math.tan(latRad) * Math.tan(sdRad));
+    const hourAngle = Math.acos(hAarg);
+    const delta = longitude + radToDeg(rise ? hourAngle : -hourAngle);
+    return 720 - (4.0 * delta) - eqTime; // in UTC minutes
+}
+/** @public
+ * calculate solar sunrise or sunset for a given day and cartographic location.
+ */
+export function calculateSunriseOrSunset(date, location, sunrise) {
+    const jDay = calculateJulianDay(date);
+    const longitude = location.longitudeDegrees;
+    const latitude = location.latitudeDegrees;
+    const utcMinutes = calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay);
+    return sunrise ? dateFromUtcMinutes(date, utcMinutes) : dateFromUtcMinutes(date, calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay + utcMinutes / 1440));
+}
 //# sourceMappingURL=SolarCalculate.js.map

package/lib/esm/SolarCalculate.js.map

@@ -1 +1 @@
-{"version":3,"file":"SolarCalculate.js","sourceRoot":"","sources":["../../src/SolarCalculate.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,KAAK,EAAE,QAAQ,EAAE,MAAM,sBAAsB,CAAC;AAGvD,oCAAoC;AAEpC,kFAAkF;AAClF,SAAS,kBAAkB,CAAC,EAAU;IACpC,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC,GAAG,OAAO,CAAC;IACrC,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,QAAQ,CAAC,QAAgB;IAChC,OAAO,CAAC,KAAK,GAAG,QAAQ,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;AACtC,CAAC;AAED,SAAS,QAAQ,CAAC,QAAgB;IAChC,OAAO,CAAC,IAAI,CAAC,EAAE,GAAG,QAAQ,GAAG,KAAK,CAAC,CAAC;AACtC,CAAC;AAED,SAAS,mBAAmB,CAAC,CAAS;IACpC,IAAI,EAAE,GAAG,SAAS,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC;IACzD,OAAO,EAAE,GAAG,KAAK,EAAE;QACjB,EAAE,IAAI,KAAK,CAAC;KACb;IACD,OAAO,EAAE,GAAG,GAAG,EAAE;QACf,EAAE,IAAI,KAAK,CAAC;KACb;IACD,OAAO,EAAE,CAAC,CAAE,aAAa;AAC3B,CAAC;AAED,SAAS,sBAAsB,CAAC,CAAS;IACvC,MAAM,CAAC,GAAG,SAAS,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,SAAS,GAAG,CAAC,CAAC,CAAC;IACxD,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,0BAA0B,CAAC,CAAS;IAC3C,MAAM,CAAC,GAAG,WAAW,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,YAAY,GAAG,CAAC,CAAC,CAAC;IAC7D,OAAO,CAAC,CAAC,CAAE,WAAW;AACxB,CAAC;AAED,SAAS,iBAAiB,CAAC,CAAS;IAClC,MAAM,CAAC,GAAG,sBAAsB,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,IAAI,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAC5B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;IAC3C,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,QAAQ,CAAC;IACnH,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,eAAe,CAAC,CAAS;IAChC,MAAM,EAAE,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,CAAC,GAAG,iBAAiB,CAAC,CAAC,CAAC,CAAC;IAC/B,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACjB,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,mBAAmB,CAAC,CAAS;IACpC,MAAM,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,KAAK,GAAG,MAAM,GAAG,QAAQ,GAAG,CAAC,CAAC;IACpC,MAAM,MAAM,GAAG,CAAC,GAAG,OAAO,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC;IACjE,OAAO,MAAM,CAAC,CAAE,aAAa;AAC/B,CAAC;AAED,SAAS,2BAA2B,CAAC,CAAS;IAC5C,MAAM,OAAO,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IACxE,MAAM,EAAE,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC;IACnD,OAAO,EAAE,CAAC,CAAE,aAAa;AAC3B,CAAC;AAED,SAAS,uBAAuB,CAAC,CAAS;IACxC,MAAM,EAAE,GAAG,2BAA2B,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,KAAK,GAAG,MAAM,GAAG,QAAQ,GAAG,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,EAAE,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC;IACnD,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,kBAAkB,CAAC,CAAS;IACnC,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,MAAM,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAEtC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;IAChE,MAAM,KAAK,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;IACxC,OAAO,KAAK,CAAC,CAAE,aAAa;AAC9B,CAAC;AAED,SAAS,kBAAkB,CAAC,CAAS;IACnC,MAAM,OAAO,GAAG,uBAAuB,CAAC,CAAC,CAAC,CAAC;IAC3C,MAAM,EAAE,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,CAAC,GAAG,0BAA0B,CAAC,CAAC,CAAC,CAAC;IACxC,MAAM,CAAC,GAAG,sBAAsB,CAAC,CAAC,CAAC,CAAC;IAEpC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC;IAC1C,CAAC,IAAI,CAAC,CAAC;IAEP,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAE1C,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;IACtH,OAAO,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,qBAAqB;AACrD,CAAC;AAED,SAAS,QAAQ,CAAC,CAAS,EAAE,SAAiB,EAAE,QAAgB,EAAE,SAAiB,EAAE,IAAY;IAC/F,MAAM,MAAM,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,YAAY,GAAG,MAAM,GAAG,GAAG,GAAG,SAAS,GAAG,IAAI,GAAG,IAAI,CAAC;IAC5D,IAAI,aAAa,GAAG,SAAS,GAAG,YAAY,CAAC;IAC7C,OAAO,aAAa,GAAG,IAAI;QACzB,aAAa,IAAI,IAAI,CAAC;IAExB,IAAI,SAAS,GAAG,aAAa,GAAG,GAAG,GAAG,KAAK,CAAC;IAC5C,IAAI,SAAS,GAAG,CAAC,GAAG,EAAE;QACpB,SAAS,IAAI,KAAK,CAAC;KACpB;IACD,MAAM,KAAK,GAAG,QAAQ,CAAC,SAAS,CAAC,CAAC;IAClC,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAChJ,IAAI,GAAG,GAAG,GAAG,EAAE;QACb,GAAG,GAAG,GAAG,CAAC;KACX;SAAM,IAAI,GAAG,GAAG,CAAC,GAAG,EAAE;QACrB,GAAG,GAAG,CAAC,GAAG,CAAC;KACZ;IACD,MAAM,MAAM,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAC5E,IAAI,OAAO,CAAC;IACZ,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,KAAK,EAAE;QAC7B,IAAI,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QAChH,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,EAAE;YACzB,IAAI,KAAK,GAAG,CAAC,EAAE;gBACb,KAAK,GAAG,CAAC,GAAG,CAAC;aACd;iBAAM;gBACL,KAAK,GAAG,GAAG,CAAC;aACb;SACF;QACD,OAAO,GAAG,KAAK,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QAC7C,IAAI,SAAS,GAAG,GAAG,EAAE;YACnB,OAAO,GAAG,CAAC,OAAO,CAAC;SACpB;KACF;SAAM;QACL,IAAI,QAAQ,GAAG,GAAG,EAAE;YAClB,OAAO,GAAG,KAAK,CAAC;SACjB;aAAM;YACL,OAAO,GAAG,GAAG,CAAC;SACf;KACF;IACD,IAAI,OAAO,GAAG,GAAG,EAAE;QACjB,OAAO,IAAI,KAAK,CAAC;KAClB;IACD,OAAO,EAAE,OAAO,EAAE,SAAS,EAAE,EAAE,GAAG,MAAM,EAAE,CAAC;AAC7C,CAAC;AAED,SAAS,kBAAkB,CAAC,IAAU;IACpC,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,QAAQ,CAAC,GAAG,SAAS,CAAC,CAAI,wFAAwF;AACvJ,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,oBAAoB,CAAC,IAAU,EAAE,QAAsB;IACrE,MAAM,IAAI,GAAG,kBAAkB,CAAC,IAAI,CAAC,CAAC;IACtC,MAAM,QAAQ,GAAG,QAAQ,CAAC,eAAe,CAAC;IAC1C,MAAM,SAAS,GAAG,QAAQ,CAAC,gBAAgB,CAAC;IAC5C,MAAM,UAAU,GAAG,IAAI,CAAC,WAAW,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,aAAa,EAAE,GAAG,IAAI,CAAC,aAAa,EAAE,GAAG,EAAE,CAAC;IAC9F,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,uCAAuC;IACvF,MAAM,YAAY,GAAG,UAAU,GAAG,IAAI,GAAG,EAAE,CAAC;IAC5C,MAAM,MAAM,GAAG,IAAI,GAAG,UAAU,GAAG,MAAM,CAAC;IAC1C,MAAM,CAAC,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;IACrC,OAAO,QAAQ,CAAC,CAAC,EAAE,YAAY,EAAE,QAAQ,EAAE,SAAS,EAAE,IAAI,CAAC,CAAC;AAC9D,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,uBAAuB,CAAC,IAAU,EAAE,QAAsB;IACxE,OAAO,iCAAiC,CAAC,oBAAoB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC,CAAC;AACjF,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,iCAAiC,CAAC,gBAAwD;IACxG,MAAM,OAAO,GAAG,KAAK,CAAC,gBAAgB,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;IACjE,MAAM,SAAS,GAAG,KAAK,CAAC,gBAAgB,CAAC,gBAAgB,CAAC,SAAS,CAAC,CAAC;IACrE,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACzC,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACzC,OAAO,QAAQ,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,YAAY,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,YAAY,EAAE,CAAC,YAAY,CAAC,CAAC;AAC9G,CAAC;AAED,SAAS,kBAAkB,CAAC,IAAU,EAAE,UAAkB;IACxD,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,CAAC;IAC/C,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC;IAC9B,MAAM,CAAC,WAAW,CAAC,QAAQ,CAAC,CAAC;IAC7B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,GAAG,QAAQ,CAAC,CAAC,CAAC;IACpE,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;IACxB,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,qBAAqB,CAAC,IAAa,EAAE,GAAW,EAAE,SAAiB,EAAE,IAAY;IACxF,MAAM,CAAC,GAAG,kBAAkB,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,MAAM,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,QAAQ,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACvC,MAAM,MAAM,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,KAAK,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;IACvH,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,SAAS,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAE,SAAS,CAAC,CAAC;IACnE,OAAO,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC,iBAAiB;AACxD,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,wBAAwB,CAAC,IAAU,EAAE,QAAsB,EAAE,OAAgB;IAC3F,MAAM,IAAI,GAAG,kBAAkB,CAAC,IAAI,CAAC,CAAC;IACtC,MAAM,SAAS,GAAG,QAAQ,CAAC,gBAAgB,CAAC;IAC5C,MAAM,QAAQ,GAAG,QAAQ,CAAC,eAAe,CAAC;IAC1C,MAAM,UAAU,GAAG,qBAAqB,CAAC,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,IAAI,CAAC,CAAC;IAC7E,OAAO,OAAO,CAAC,CAAC,CAAC,kBAAkB,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,kBAAkB,CAAC,IAAI,EAAE,qBAAqB,CAAC,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,IAAI,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC,CAAC;AAClK,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\r\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\r\n* See LICENSE.md in the project root for license terms and full copyright notice.\r\n*--------------------------------------------------------------------------------------------*/\r\n/** @packageDocumentation\r\n * @module Utils\r\n */\r\n\r\nimport { Angle, Vector3d } from \"@itwin/core-geometry\";\r\nimport { Cartographic } from \"./geometry/Cartographic\";\r\n\r\n// cspell:ignore mrad sinm sint aarg\r\n\r\n// Code below loosely translated from  https://www.esrl.noaa.gov/gmd/grad/solcalc/\r\nfunction calcTimeJulianCent(jd: number) {\r\n  const T = (jd - 2451545.0) / 36525.0;\r\n  return T;\r\n}\r\n\r\nfunction radToDeg(angleRad: number) {\r\n  return (180.0 * angleRad / Math.PI);\r\n}\r\n\r\nfunction degToRad(angleDeg: number) {\r\n  return (Math.PI * angleDeg / 180.0);\r\n}\r\n\r\nfunction calcGeomMeanLongSun(t: number) {\r\n  let L0 = 280.46646 + t * (36000.76983 + t * (0.0003032));\r\n  while (L0 > 360.0) {\r\n    L0 -= 360.0;\r\n  }\r\n  while (L0 < 0.0) {\r\n    L0 += 360.0;\r\n  }\r\n  return L0;\t\t// in degrees\r\n}\r\n\r\nfunction calcGeomMeanAnomalySun(t: number) {\r\n  const M = 357.52911 + t * (35999.05029 - 0.0001537 * t);\r\n  return M;\t\t// in degrees\r\n}\r\n\r\nfunction calcEccentricityEarthOrbit(t: number) {\r\n  const e = 0.016708634 - t * (0.000042037 + 0.0000001267 * t);\r\n  return e;\t\t// unitless\r\n}\r\n\r\nfunction calcSunEqOfCenter(t: number) {\r\n  const m = calcGeomMeanAnomalySun(t);\r\n  const mrad = degToRad(m);\r\n  const sinm = Math.sin(mrad);\r\n  const sin2m = Math.sin(mrad + mrad);\r\n  const sin3m = Math.sin(mrad + mrad + mrad);\r\n  const C = sinm * (1.914602 - t * (0.004817 + 0.000014 * t)) + sin2m * (0.019993 - 0.000101 * t) + sin3m * 0.000289;\r\n  return C;\t\t// in degrees\r\n}\r\n\r\nfunction calcSunTrueLong(t: number) {\r\n  const l0 = calcGeomMeanLongSun(t);\r\n  const c = calcSunEqOfCenter(t);\r\n  const O = l0 + c;\r\n  return O;\t\t// in degrees\r\n}\r\n\r\nfunction calcSunApparentLong(t: number) {\r\n  const o = calcSunTrueLong(t);\r\n  const omega = 125.04 - 1934.136 * t;\r\n  const lambda = o - 0.00569 - 0.00478 * Math.sin(degToRad(omega));\r\n  return lambda;\t\t// in degrees\r\n}\r\n\r\nfunction calcMeanObliquityOfEcliptic(t: number) {\r\n  const seconds = 21.448 - t * (46.8150 + t * (0.00059 - t * (0.001813)));\r\n  const e0 = 23.0 + (26.0 + (seconds / 60.0)) / 60.0;\r\n  return e0;\t\t// in degrees\r\n}\r\n\r\nfunction calcObliquityCorrection(t: number) {\r\n  const e0 = calcMeanObliquityOfEcliptic(t);\r\n  const omega = 125.04 - 1934.136 * t;\r\n  const e = e0 + 0.00256 * Math.cos(degToRad(omega));\r\n  return e;\t\t// in degrees\r\n}\r\n\r\nfunction calcSunDeclination(t: number) {\r\n  const e = calcObliquityCorrection(t);\r\n  const lambda = calcSunApparentLong(t);\r\n\r\n  const sint = Math.sin(degToRad(e)) * Math.sin(degToRad(lambda));\r\n  const theta = radToDeg(Math.asin(sint));\r\n  return theta;\t\t// in degrees\r\n}\r\n\r\nfunction calcEquationOfTime(t: number) {\r\n  const epsilon = calcObliquityCorrection(t);\r\n  const l0 = calcGeomMeanLongSun(t);\r\n  const e = calcEccentricityEarthOrbit(t);\r\n  const m = calcGeomMeanAnomalySun(t);\r\n\r\n  let y = Math.tan(degToRad(epsilon) / 2.0);\r\n  y *= y;\r\n\r\n  const sin2l0 = Math.sin(2.0 * degToRad(l0));\r\n  const sinm = Math.sin(degToRad(m));\r\n  const cos2l0 = Math.cos(2.0 * degToRad(l0));\r\n  const sin4l0 = Math.sin(4.0 * degToRad(l0));\r\n  const sin2m = Math.sin(2.0 * degToRad(m));\r\n\r\n  const eTime = y * sin2l0 - 2.0 * e * sinm + 4.0 * e * y * sinm * cos2l0 - 0.5 * y * y * sin4l0 - 1.25 * e * e * sin2m;\r\n  return radToDeg(eTime) * 4.0;\t// in minutes of time\r\n}\r\n\r\nfunction calcAzEl(t: number, localTime: number, latitude: number, longitude: number, zone: number): { azimuth: number, elevation: number } {\r\n  const eqTime = calcEquationOfTime(t);\r\n  const theta = calcSunDeclination(t);\r\n  const solarTimeFix = eqTime + 4.0 * longitude - 60.0 * zone;\r\n  let trueSolarTime = localTime + solarTimeFix;\r\n  while (trueSolarTime > 1440)\r\n    trueSolarTime -= 1440;\r\n\r\n  let hourAngle = trueSolarTime / 4.0 - 180.0;\r\n  if (hourAngle < -180) {\r\n    hourAngle += 360.0;\r\n  }\r\n  const haRad = degToRad(hourAngle);\r\n  let csz = Math.sin(degToRad(latitude)) * Math.sin(degToRad(theta)) + Math.cos(degToRad(latitude)) * Math.cos(degToRad(theta)) * Math.cos(haRad);\r\n  if (csz > 1.0) {\r\n    csz = 1.0;\r\n  } else if (csz < -1.0) {\r\n    csz = -1.0;\r\n  }\r\n  const zenith = radToDeg(Math.acos(csz));\r\n  const azDenom = (Math.cos(degToRad(latitude)) * Math.sin(degToRad(zenith)));\r\n  let azimuth;\r\n  if (Math.abs(azDenom) > 0.001) {\r\n    let azRad = ((Math.sin(degToRad(latitude)) * Math.cos(degToRad(zenith))) - Math.sin(degToRad(theta))) / azDenom;\r\n    if (Math.abs(azRad) > 1.0) {\r\n      if (azRad < 0) {\r\n        azRad = -1.0;\r\n      } else {\r\n        azRad = 1.0;\r\n      }\r\n    }\r\n    azimuth = 180.0 - radToDeg(Math.acos(azRad));\r\n    if (hourAngle > 0.0) {\r\n      azimuth = -azimuth;\r\n    }\r\n  } else {\r\n    if (latitude > 0.0) {\r\n      azimuth = 180.0;\r\n    } else {\r\n      azimuth = 0.0;\r\n    }\r\n  }\r\n  if (azimuth < 0.0) {\r\n    azimuth += 360.0;\r\n  }\r\n  return { azimuth, elevation: 90 - zenith };\r\n}\r\n\r\nfunction calculateJulianDay(date: Date) {\r\n  return Math.floor(date.getTime() / 86400000) + 2440587.5;    // https://stackoverflow.com/questions/11759992/calculating-jdayjulian-day-in-javascript\r\n}\r\n\r\n/** @public\r\n * calculate solar angles (in radians) based at a given date/time and cartographic location.\r\n */\r\nexport function calculateSolarAngles(date: Date, location: Cartographic): { azimuth: number, elevation: number } {\r\n  const jDay = calculateJulianDay(date);\r\n  const latitude = location.latitudeDegrees;\r\n  const longitude = location.longitudeDegrees;\r\n  const utcMinutes = date.getUTCHours() * 60 + date.getUTCMinutes() + date.getUTCSeconds() / 60;\r\n  const zone = Math.floor(.5 + longitude / 15.0); // date.getTimeZoneOffset mixes in DST.\r\n  const localMinutes = utcMinutes + zone * 60;\r\n  const jTotal = jDay + utcMinutes / 1440.0;\r\n  const T = calcTimeJulianCent(jTotal);\r\n  return calcAzEl(T, localMinutes, latitude, longitude, zone);\r\n}\r\n\r\n/** @public\r\n * calculate solar direction based at a given date/time and cartpgrphic location.\r\n */\r\nexport function calculateSolarDirection(date: Date, location: Cartographic): Vector3d {\r\n  return calculateSolarDirectionFromAngles(calculateSolarAngles(date, location));\r\n}\r\n\r\n/** @public\r\n * calculate solar direction corresponding to the given azimuth and elevation (altitude) angles in degrees.\r\n */\r\nexport function calculateSolarDirectionFromAngles(azimuthElevation: { azimuth: number, elevation: number }): Vector3d {\r\n  const azimuth = Angle.degreesToRadians(azimuthElevation.azimuth);\r\n  const elevation = Angle.degreesToRadians(azimuthElevation.elevation);\r\n  const cosElevation = Math.cos(elevation);\r\n  const sinElevation = Math.sin(elevation);\r\n  return Vector3d.create(-Math.sin(azimuth) * cosElevation, -Math.cos(azimuth) * cosElevation, -sinElevation);\r\n}\r\n\r\nfunction dateFromUtcMinutes(date: Date, utcMinutes: number) {\r\n  const utcHours = Math.floor(utcMinutes / 60.0);\r\n  const output = new Date(date);\r\n  output.setUTCHours(utcHours);\r\n  output.setUTCMinutes(Math.floor(.5 + utcMinutes - 60.0 * utcHours));\r\n  output.setUTCSeconds(0);\r\n  return output;\r\n}\r\n\r\nfunction calcSunriseUtcMinutes(rise: boolean, lat: number, longitude: number, jDay: number) {\r\n  const t = calcTimeJulianCent(jDay);\r\n  const eqTime = calcEquationOfTime(t);\r\n  const solarDec = calcSunDeclination(t);\r\n  const latRad = degToRad(lat);\r\n  const sdRad = degToRad(solarDec);\r\n  const hAarg = (Math.cos(degToRad(90.833)) / (Math.cos(latRad) * Math.cos(sdRad)) - Math.tan(latRad) * Math.tan(sdRad));\r\n  const hourAngle = Math.acos(hAarg);\r\n  const delta = longitude + radToDeg(rise ? hourAngle : - hourAngle);\r\n  return 720 - (4.0 * delta) - eqTime;\t// in UTC minutes\r\n}\r\n\r\n/** @public\r\n * calculate solar sunrise or sunset for a given day and cartographic location.\r\n */\r\nexport function calculateSunriseOrSunset(date: Date, location: Cartographic, sunrise: boolean): Date {\r\n  const jDay = calculateJulianDay(date);\r\n  const longitude = location.longitudeDegrees;\r\n  const latitude = location.latitudeDegrees;\r\n  const utcMinutes = calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay);\r\n  return sunrise ? dateFromUtcMinutes(date, utcMinutes) : dateFromUtcMinutes(date, calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay + utcMinutes / 1440));\r\n}\r\n"]}
+{"version":3,"file":"SolarCalculate.js","sourceRoot":"","sources":["../../src/SolarCalculate.ts"],"names":[],"mappings":"AAAA;;;+FAG+F;AAC/F;;GAEG;AAEH,OAAO,EAAE,KAAK,EAAE,QAAQ,EAAE,MAAM,sBAAsB,CAAC;AAGvD,oCAAoC;AAEpC,kFAAkF;AAClF,SAAS,kBAAkB,CAAC,EAAU;IACpC,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,SAAS,CAAC,GAAG,OAAO,CAAC;IACrC,OAAO,CAAC,CAAC;AACX,CAAC;AAED,SAAS,QAAQ,CAAC,QAAgB;IAChC,OAAO,CAAC,KAAK,GAAG,QAAQ,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;AACtC,CAAC;AAED,SAAS,QAAQ,CAAC,QAAgB;IAChC,OAAO,CAAC,IAAI,CAAC,EAAE,GAAG,QAAQ,GAAG,KAAK,CAAC,CAAC;AACtC,CAAC;AAED,SAAS,mBAAmB,CAAC,CAAS;IACpC,IAAI,EAAE,GAAG,SAAS,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC;IACzD,OAAO,EAAE,GAAG,KAAK,EAAE;QACjB,EAAE,IAAI,KAAK,CAAC;KACb;IACD,OAAO,EAAE,GAAG,GAAG,EAAE;QACf,EAAE,IAAI,KAAK,CAAC;KACb;IACD,OAAO,EAAE,CAAC,CAAE,aAAa;AAC3B,CAAC;AAED,SAAS,sBAAsB,CAAC,CAAS;IACvC,MAAM,CAAC,GAAG,SAAS,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,SAAS,GAAG,CAAC,CAAC,CAAC;IACxD,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,0BAA0B,CAAC,CAAS;IAC3C,MAAM,CAAC,GAAG,WAAW,GAAG,CAAC,GAAG,CAAC,WAAW,GAAG,YAAY,GAAG,CAAC,CAAC,CAAC;IAC7D,OAAO,CAAC,CAAC,CAAE,WAAW;AACxB,CAAC;AAED,SAAS,iBAAiB,CAAC,CAAS;IAClC,MAAM,CAAC,GAAG,sBAAsB,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,IAAI,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;IACzB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IAC5B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC;IAC3C,MAAM,CAAC,GAAG,IAAI,GAAG,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,QAAQ,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,KAAK,GAAG,QAAQ,CAAC;IACnH,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,eAAe,CAAC,CAAS;IAChC,MAAM,EAAE,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,CAAC,GAAG,iBAAiB,CAAC,CAAC,CAAC,CAAC;IAC/B,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACjB,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,mBAAmB,CAAC,CAAS;IACpC,MAAM,CAAC,GAAG,eAAe,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,KAAK,GAAG,MAAM,GAAG,QAAQ,GAAG,CAAC,CAAC;IACpC,MAAM,MAAM,GAAG,CAAC,GAAG,OAAO,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC;IACjE,OAAO,MAAM,CAAC,CAAE,aAAa;AAC/B,CAAC;AAED,SAAS,2BAA2B,CAAC,CAAS;IAC5C,MAAM,OAAO,GAAG,MAAM,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IACxE,MAAM,EAAE,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC;IACnD,OAAO,EAAE,CAAC,CAAE,aAAa;AAC3B,CAAC;AAED,SAAS,uBAAuB,CAAC,CAAS;IACxC,MAAM,EAAE,GAAG,2BAA2B,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,KAAK,GAAG,MAAM,GAAG,QAAQ,GAAG,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,EAAE,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC;IACnD,OAAO,CAAC,CAAC,CAAE,aAAa;AAC1B,CAAC;AAED,SAAS,kBAAkB,CAAC,CAAS;IACnC,MAAM,CAAC,GAAG,uBAAuB,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,MAAM,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAEtC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;IAChE,MAAM,KAAK,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;IACxC,OAAO,KAAK,CAAC,CAAE,aAAa;AAC9B,CAAC;AAED,SAAS,kBAAkB,CAAC,CAAS;IACnC,MAAM,OAAO,GAAG,uBAAuB,CAAC,CAAC,CAAC,CAAC;IAC3C,MAAM,EAAE,GAAG,mBAAmB,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,CAAC,GAAG,0BAA0B,CAAC,CAAC,CAAC,CAAC;IACxC,MAAM,CAAC,GAAG,sBAAsB,CAAC,CAAC,CAAC,CAAC;IAEpC,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC;IAC1C,CAAC,IAAI,CAAC,CAAC;IAEP,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC;IAC5C,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAE1C,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC;IACtH,OAAO,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,qBAAqB;AACrD,CAAC;AAED,SAAS,QAAQ,CAAC,CAAS,EAAE,SAAiB,EAAE,QAAgB,EAAE,SAAiB,EAAE,IAAY;IAC/F,MAAM,MAAM,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,YAAY,GAAG,MAAM,GAAG,GAAG,GAAG,SAAS,GAAG,IAAI,GAAG,IAAI,CAAC;IAC5D,IAAI,aAAa,GAAG,SAAS,GAAG,YAAY,CAAC;IAC7C,OAAO,aAAa,GAAG,IAAI;QACzB,aAAa,IAAI,IAAI,CAAC;IAExB,IAAI,SAAS,GAAG,aAAa,GAAG,GAAG,GAAG,KAAK,CAAC;IAC5C,IAAI,SAAS,GAAG,CAAC,GAAG,EAAE;QACpB,SAAS,IAAI,KAAK,CAAC;KACpB;IACD,MAAM,KAAK,GAAG,QAAQ,CAAC,SAAS,CAAC,CAAC;IAClC,IAAI,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAChJ,IAAI,GAAG,GAAG,GAAG,EAAE;QACb,GAAG,GAAG,GAAG,CAAC;KACX;SAAM,IAAI,GAAG,GAAG,CAAC,GAAG,EAAE;QACrB,GAAG,GAAG,CAAC,GAAG,CAAC;KACZ;IACD,MAAM,MAAM,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACxC,MAAM,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAC5E,IAAI,OAAO,CAAC;IACZ,IAAI,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,KAAK,EAAE;QAC7B,IAAI,KAAK,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC;QAChH,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,EAAE;YACzB,IAAI,KAAK,GAAG,CAAC,EAAE;gBACb,KAAK,GAAG,CAAC,GAAG,CAAC;aACd;iBAAM;gBACL,KAAK,GAAG,GAAG,CAAC;aACb;SACF;QACD,OAAO,GAAG,KAAK,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;QAC7C,IAAI,SAAS,GAAG,GAAG,EAAE;YACnB,OAAO,GAAG,CAAC,OAAO,CAAC;SACpB;KACF;SAAM;QACL,IAAI,QAAQ,GAAG,GAAG,EAAE;YAClB,OAAO,GAAG,KAAK,CAAC;SACjB;aAAM;YACL,OAAO,GAAG,GAAG,CAAC;SACf;KACF;IACD,IAAI,OAAO,GAAG,GAAG,EAAE;QACjB,OAAO,IAAI,KAAK,CAAC;KAClB;IACD,OAAO,EAAE,OAAO,EAAE,SAAS,EAAE,EAAE,GAAG,MAAM,EAAE,CAAC;AAC7C,CAAC;AAED,SAAS,kBAAkB,CAAC,IAAU;IACpC,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,QAAQ,CAAC,GAAG,SAAS,CAAC,CAAI,wFAAwF;AACvJ,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,oBAAoB,CAAC,IAAU,EAAE,QAAsB;IACrE,MAAM,IAAI,GAAG,kBAAkB,CAAC,IAAI,CAAC,CAAC;IACtC,MAAM,QAAQ,GAAG,QAAQ,CAAC,eAAe,CAAC;IAC1C,MAAM,SAAS,GAAG,QAAQ,CAAC,gBAAgB,CAAC;IAC5C,MAAM,UAAU,GAAG,IAAI,CAAC,WAAW,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,aAAa,EAAE,GAAG,IAAI,CAAC,aAAa,EAAE,GAAG,EAAE,CAAC;IAC9F,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,SAAS,GAAG,IAAI,CAAC,CAAC,CAAC,uCAAuC;IACvF,MAAM,YAAY,GAAG,UAAU,GAAG,IAAI,GAAG,EAAE,CAAC;IAC5C,MAAM,MAAM,GAAG,IAAI,GAAG,UAAU,GAAG,MAAM,CAAC;IAC1C,MAAM,CAAC,GAAG,kBAAkB,CAAC,MAAM,CAAC,CAAC;IACrC,OAAO,QAAQ,CAAC,CAAC,EAAE,YAAY,EAAE,QAAQ,EAAE,SAAS,EAAE,IAAI,CAAC,CAAC;AAC9D,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,uBAAuB,CAAC,IAAU,EAAE,QAAsB;IACxE,OAAO,iCAAiC,CAAC,oBAAoB,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC,CAAC;AACjF,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,iCAAiC,CAAC,gBAAwD;IACxG,MAAM,OAAO,GAAG,KAAK,CAAC,gBAAgB,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;IACjE,MAAM,SAAS,GAAG,KAAK,CAAC,gBAAgB,CAAC,gBAAgB,CAAC,SAAS,CAAC,CAAC;IACrE,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACzC,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACzC,OAAO,QAAQ,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,YAAY,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,GAAG,YAAY,EAAE,CAAC,YAAY,CAAC,CAAC;AAC9G,CAAC;AAED,SAAS,kBAAkB,CAAC,IAAU,EAAE,UAAkB;IACxD,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,UAAU,GAAG,IAAI,CAAC,CAAC;IAC/C,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,CAAC;IAC9B,MAAM,CAAC,WAAW,CAAC,QAAQ,CAAC,CAAC;IAC7B,MAAM,CAAC,aAAa,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,UAAU,GAAG,IAAI,GAAG,QAAQ,CAAC,CAAC,CAAC;IACpE,MAAM,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC;IACxB,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,SAAS,qBAAqB,CAAC,IAAa,EAAE,GAAW,EAAE,SAAiB,EAAE,IAAY;IACxF,MAAM,CAAC,GAAG,kBAAkB,CAAC,IAAI,CAAC,CAAC;IACnC,MAAM,MAAM,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,QAAQ,GAAG,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACvC,MAAM,MAAM,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;IAC7B,MAAM,KAAK,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;IACjC,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;IACvH,MAAM,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IACnC,MAAM,KAAK,GAAG,SAAS,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,CAAE,SAAS,CAAC,CAAC;IACnE,OAAO,GAAG,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,CAAC,iBAAiB;AACxD,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,wBAAwB,CAAC,IAAU,EAAE,QAAsB,EAAE,OAAgB;IAC3F,MAAM,IAAI,GAAG,kBAAkB,CAAC,IAAI,CAAC,CAAC;IACtC,MAAM,SAAS,GAAG,QAAQ,CAAC,gBAAgB,CAAC;IAC5C,MAAM,QAAQ,GAAG,QAAQ,CAAC,eAAe,CAAC;IAC1C,MAAM,UAAU,GAAG,qBAAqB,CAAC,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,IAAI,CAAC,CAAC;IAC7E,OAAO,OAAO,CAAC,CAAC,CAAC,kBAAkB,CAAC,IAAI,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,kBAAkB,CAAC,IAAI,EAAE,qBAAqB,CAAC,OAAO,EAAE,QAAQ,EAAE,SAAS,EAAE,IAAI,GAAG,UAAU,GAAG,IAAI,CAAC,CAAC,CAAC;AAClK,CAAC","sourcesContent":["/*---------------------------------------------------------------------------------------------\n* Copyright (c) Bentley Systems, Incorporated. All rights reserved.\n* See LICENSE.md in the project root for license terms and full copyright notice.\n*--------------------------------------------------------------------------------------------*/\n/** @packageDocumentation\n * @module Utils\n */\n\nimport { Angle, Vector3d } from \"@itwin/core-geometry\";\nimport { Cartographic } from \"./geometry/Cartographic\";\n\n// cspell:ignore mrad sinm sint aarg\n\n// Code below loosely translated from  https://www.esrl.noaa.gov/gmd/grad/solcalc/\nfunction calcTimeJulianCent(jd: number) {\n  const T = (jd - 2451545.0) / 36525.0;\n  return T;\n}\n\nfunction radToDeg(angleRad: number) {\n  return (180.0 * angleRad / Math.PI);\n}\n\nfunction degToRad(angleDeg: number) {\n  return (Math.PI * angleDeg / 180.0);\n}\n\nfunction calcGeomMeanLongSun(t: number) {\n  let L0 = 280.46646 + t * (36000.76983 + t * (0.0003032));\n  while (L0 > 360.0) {\n    L0 -= 360.0;\n  }\n  while (L0 < 0.0) {\n    L0 += 360.0;\n  }\n  return L0;\t\t// in degrees\n}\n\nfunction calcGeomMeanAnomalySun(t: number) {\n  const M = 357.52911 + t * (35999.05029 - 0.0001537 * t);\n  return M;\t\t// in degrees\n}\n\nfunction calcEccentricityEarthOrbit(t: number) {\n  const e = 0.016708634 - t * (0.000042037 + 0.0000001267 * t);\n  return e;\t\t// unitless\n}\n\nfunction calcSunEqOfCenter(t: number) {\n  const m = calcGeomMeanAnomalySun(t);\n  const mrad = degToRad(m);\n  const sinm = Math.sin(mrad);\n  const sin2m = Math.sin(mrad + mrad);\n  const sin3m = Math.sin(mrad + mrad + mrad);\n  const C = sinm * (1.914602 - t * (0.004817 + 0.000014 * t)) + sin2m * (0.019993 - 0.000101 * t) + sin3m * 0.000289;\n  return C;\t\t// in degrees\n}\n\nfunction calcSunTrueLong(t: number) {\n  const l0 = calcGeomMeanLongSun(t);\n  const c = calcSunEqOfCenter(t);\n  const O = l0 + c;\n  return O;\t\t// in degrees\n}\n\nfunction calcSunApparentLong(t: number) {\n  const o = calcSunTrueLong(t);\n  const omega = 125.04 - 1934.136 * t;\n  const lambda = o - 0.00569 - 0.00478 * Math.sin(degToRad(omega));\n  return lambda;\t\t// in degrees\n}\n\nfunction calcMeanObliquityOfEcliptic(t: number) {\n  const seconds = 21.448 - t * (46.8150 + t * (0.00059 - t * (0.001813)));\n  const e0 = 23.0 + (26.0 + (seconds / 60.0)) / 60.0;\n  return e0;\t\t// in degrees\n}\n\nfunction calcObliquityCorrection(t: number) {\n  const e0 = calcMeanObliquityOfEcliptic(t);\n  const omega = 125.04 - 1934.136 * t;\n  const e = e0 + 0.00256 * Math.cos(degToRad(omega));\n  return e;\t\t// in degrees\n}\n\nfunction calcSunDeclination(t: number) {\n  const e = calcObliquityCorrection(t);\n  const lambda = calcSunApparentLong(t);\n\n  const sint = Math.sin(degToRad(e)) * Math.sin(degToRad(lambda));\n  const theta = radToDeg(Math.asin(sint));\n  return theta;\t\t// in degrees\n}\n\nfunction calcEquationOfTime(t: number) {\n  const epsilon = calcObliquityCorrection(t);\n  const l0 = calcGeomMeanLongSun(t);\n  const e = calcEccentricityEarthOrbit(t);\n  const m = calcGeomMeanAnomalySun(t);\n\n  let y = Math.tan(degToRad(epsilon) / 2.0);\n  y *= y;\n\n  const sin2l0 = Math.sin(2.0 * degToRad(l0));\n  const sinm = Math.sin(degToRad(m));\n  const cos2l0 = Math.cos(2.0 * degToRad(l0));\n  const sin4l0 = Math.sin(4.0 * degToRad(l0));\n  const sin2m = Math.sin(2.0 * degToRad(m));\n\n  const eTime = y * sin2l0 - 2.0 * e * sinm + 4.0 * e * y * sinm * cos2l0 - 0.5 * y * y * sin4l0 - 1.25 * e * e * sin2m;\n  return radToDeg(eTime) * 4.0;\t// in minutes of time\n}\n\nfunction calcAzEl(t: number, localTime: number, latitude: number, longitude: number, zone: number): { azimuth: number, elevation: number } {\n  const eqTime = calcEquationOfTime(t);\n  const theta = calcSunDeclination(t);\n  const solarTimeFix = eqTime + 4.0 * longitude - 60.0 * zone;\n  let trueSolarTime = localTime + solarTimeFix;\n  while (trueSolarTime > 1440)\n    trueSolarTime -= 1440;\n\n  let hourAngle = trueSolarTime / 4.0 - 180.0;\n  if (hourAngle < -180) {\n    hourAngle += 360.0;\n  }\n  const haRad = degToRad(hourAngle);\n  let csz = Math.sin(degToRad(latitude)) * Math.sin(degToRad(theta)) + Math.cos(degToRad(latitude)) * Math.cos(degToRad(theta)) * Math.cos(haRad);\n  if (csz > 1.0) {\n    csz = 1.0;\n  } else if (csz < -1.0) {\n    csz = -1.0;\n  }\n  const zenith = radToDeg(Math.acos(csz));\n  const azDenom = (Math.cos(degToRad(latitude)) * Math.sin(degToRad(zenith)));\n  let azimuth;\n  if (Math.abs(azDenom) > 0.001) {\n    let azRad = ((Math.sin(degToRad(latitude)) * Math.cos(degToRad(zenith))) - Math.sin(degToRad(theta))) / azDenom;\n    if (Math.abs(azRad) > 1.0) {\n      if (azRad < 0) {\n        azRad = -1.0;\n      } else {\n        azRad = 1.0;\n      }\n    }\n    azimuth = 180.0 - radToDeg(Math.acos(azRad));\n    if (hourAngle > 0.0) {\n      azimuth = -azimuth;\n    }\n  } else {\n    if (latitude > 0.0) {\n      azimuth = 180.0;\n    } else {\n      azimuth = 0.0;\n    }\n  }\n  if (azimuth < 0.0) {\n    azimuth += 360.0;\n  }\n  return { azimuth, elevation: 90 - zenith };\n}\n\nfunction calculateJulianDay(date: Date) {\n  return Math.floor(date.getTime() / 86400000) + 2440587.5;    // https://stackoverflow.com/questions/11759992/calculating-jdayjulian-day-in-javascript\n}\n\n/** @public\n * calculate solar angles (in radians) based at a given date/time and cartographic location.\n */\nexport function calculateSolarAngles(date: Date, location: Cartographic): { azimuth: number, elevation: number } {\n  const jDay = calculateJulianDay(date);\n  const latitude = location.latitudeDegrees;\n  const longitude = location.longitudeDegrees;\n  const utcMinutes = date.getUTCHours() * 60 + date.getUTCMinutes() + date.getUTCSeconds() / 60;\n  const zone = Math.floor(.5 + longitude / 15.0); // date.getTimeZoneOffset mixes in DST.\n  const localMinutes = utcMinutes + zone * 60;\n  const jTotal = jDay + utcMinutes / 1440.0;\n  const T = calcTimeJulianCent(jTotal);\n  return calcAzEl(T, localMinutes, latitude, longitude, zone);\n}\n\n/** @public\n * calculate solar direction based at a given date/time and cartpgrphic location.\n */\nexport function calculateSolarDirection(date: Date, location: Cartographic): Vector3d {\n  return calculateSolarDirectionFromAngles(calculateSolarAngles(date, location));\n}\n\n/** @public\n * calculate solar direction corresponding to the given azimuth and elevation (altitude) angles in degrees.\n */\nexport function calculateSolarDirectionFromAngles(azimuthElevation: { azimuth: number, elevation: number }): Vector3d {\n  const azimuth = Angle.degreesToRadians(azimuthElevation.azimuth);\n  const elevation = Angle.degreesToRadians(azimuthElevation.elevation);\n  const cosElevation = Math.cos(elevation);\n  const sinElevation = Math.sin(elevation);\n  return Vector3d.create(-Math.sin(azimuth) * cosElevation, -Math.cos(azimuth) * cosElevation, -sinElevation);\n}\n\nfunction dateFromUtcMinutes(date: Date, utcMinutes: number) {\n  const utcHours = Math.floor(utcMinutes / 60.0);\n  const output = new Date(date);\n  output.setUTCHours(utcHours);\n  output.setUTCMinutes(Math.floor(.5 + utcMinutes - 60.0 * utcHours));\n  output.setUTCSeconds(0);\n  return output;\n}\n\nfunction calcSunriseUtcMinutes(rise: boolean, lat: number, longitude: number, jDay: number) {\n  const t = calcTimeJulianCent(jDay);\n  const eqTime = calcEquationOfTime(t);\n  const solarDec = calcSunDeclination(t);\n  const latRad = degToRad(lat);\n  const sdRad = degToRad(solarDec);\n  const hAarg = (Math.cos(degToRad(90.833)) / (Math.cos(latRad) * Math.cos(sdRad)) - Math.tan(latRad) * Math.tan(sdRad));\n  const hourAngle = Math.acos(hAarg);\n  const delta = longitude + radToDeg(rise ? hourAngle : - hourAngle);\n  return 720 - (4.0 * delta) - eqTime;\t// in UTC minutes\n}\n\n/** @public\n * calculate solar sunrise or sunset for a given day and cartographic location.\n */\nexport function calculateSunriseOrSunset(date: Date, location: Cartographic, sunrise: boolean): Date {\n  const jDay = calculateJulianDay(date);\n  const longitude = location.longitudeDegrees;\n  const latitude = location.latitudeDegrees;\n  const utcMinutes = calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay);\n  return sunrise ? dateFromUtcMinutes(date, utcMinutes) : dateFromUtcMinutes(date, calcSunriseUtcMinutes(sunrise, latitude, longitude, jDay + utcMinutes / 1440));\n}\n"]}

package/lib/esm/SolarShadows.d.ts

@@ -1,36 +1,36 @@
-/** @packageDocumentation
- * @module DisplayStyles
- */
-import { ColorDefProps } from "./ColorDef";
-import { RgbColor } from "./RgbColor";
-/** JSON representation of [[SolarShadowSettings]].
- * @public
- * @extensions
- */
-export interface SolarShadowSettingsProps {
-    /** Shadow color. Default: [[ColorByName.grey]]. */
-    color?: ColorDefProps;
-    /** @internal */
-    bias?: number;
-}
-/** Settings controlling display of solar shadows for a [[DisplayStyle3dSettings]].
- * Solar shadows are imposed as a color scaling on geometry occluded from solar lighting.
- * @public
- */
-export declare class SolarShadowSettings {
-    /** Shadow color. */
-    readonly color: RgbColor;
-    /** @internal */
-    readonly bias: number;
-    private constructor();
-    static defaults: SolarShadowSettings;
-    static fromJSON(props?: SolarShadowSettingsProps): SolarShadowSettings;
-    toJSON(): SolarShadowSettingsProps | undefined;
-    equals(rhs: SolarShadowSettings): boolean;
-    /** Create a copy of these settings.
-     * @param changedProps Any property explicitly defined will be overridden in the copy.
-     * @returns A settings object equivalent to this one except for any properties explicitly overridden by `changedProps`.
-     */
-    clone(changedProps?: SolarShadowSettingsProps): SolarShadowSettings;
-}
+/** @packageDocumentation
+ * @module DisplayStyles
+ */
+import { ColorDefProps } from "./ColorDef";
+import { RgbColor } from "./RgbColor";
+/** JSON representation of [[SolarShadowSettings]].
+ * @public
+ * @extensions
+ */
+export interface SolarShadowSettingsProps {
+    /** Shadow color. Default: [[ColorByName.grey]]. */
+    color?: ColorDefProps;
+    /** @internal */
+    bias?: number;
+}
+/** Settings controlling display of solar shadows for a [[DisplayStyle3dSettings]].
+ * Solar shadows are imposed as a color scaling on geometry occluded from solar lighting.
+ * @public
+ */
+export declare class SolarShadowSettings {
+    /** Shadow color. */
+    readonly color: RgbColor;
+    /** @internal */
+    readonly bias: number;
+    private constructor();
+    static defaults: SolarShadowSettings;
+    static fromJSON(props?: SolarShadowSettingsProps): SolarShadowSettings;
+    toJSON(): SolarShadowSettingsProps | undefined;
+    equals(rhs: SolarShadowSettings): boolean;
+    /** Create a copy of these settings.
+     * @param changedProps Any property explicitly defined will be overridden in the copy.
+     * @returns A settings object equivalent to this one except for any properties explicitly overridden by `changedProps`.
+     */
+    clone(changedProps?: SolarShadowSettingsProps): SolarShadowSettings;
+}
 //# sourceMappingURL=SolarShadows.d.ts.map