@math.gl/sun 3.6.0-alpha.4 → 3.6.0-alpha.7

package/dist/es5/suncalc.js

@@ -32,12 +32,12 @@ function getSunDirection(timestamp, latitude, longitude) {
       azimuth = _getSunPosition.azimuth,
       altitude = _getSunPosition.altitude;
 
-  var azimuthN = azimuth + Math.PI;
-  return [-Math.sin(azimuthN), Math.cos(azimuthN), -Math.sin(altitude)];
+  return [Math.sin(azimuth) * Math.cos(altitude), Math.cos(azimuth) * Math.cos(altitude), -Math.sin(altitude)];
 }
 
 function toJulianDay(timestamp) {
-  return timestamp / DAY_IN_MS - 0.5 + JD1970;
+  var ts = typeof timestamp === 'number' ? timestamp : timestamp.getTime();
+  return ts / DAY_IN_MS - 0.5 + JD1970;
 }
 
 function toDays(timestamp) {

package/dist/es5/suncalc.js.map

@@ -1 +1 @@
-{"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","azimuthN","sin","cos","toJulianDay","getRightAscension","eclipticLongitude","b","lambda","atan2","tan","getDeclination","asin","hourAngle","latitudeInRadians","delta","dates","getSolarMeanAnomaly","days","getEclipticLongitude","meanAnomaly","M","C","P","L"],"mappings":";;;;;;;AAAA,IAAMA,kBAAkB,GAAGC,IAAI,CAACC,EAAL,GAAU,GAArC;AAEA,IAAMC,SAAS,GAAG,OAAO,EAAP,GAAY,EAAZ,GAAiB,EAAnC;AACA,IAAMC,MAAM,GAAG,OAAf;AACA,IAAMC,MAAM,GAAG,OAAf;AAGA,IAAMC,CAAC,GAAGN,kBAAkB,GAAG,OAA/B;AAIA,IAAMO,EAAE,GAAG,QAAX;AACA,IAAMC,EAAE,GAAG,UAAX;AAEA,IAAMC,MAAM,GAAG,OAAf;AACA,IAAMC,MAAM,GAAG,WAAf;;AAkBO,SAASC,cAAT,CACLC,SADK,EAELC,QAFK,EAGLC,SAHK,EAIc;AACnB,MAAMC,sBAAsB,GAAGf,kBAAkB,GAAG,CAACc,SAArD;AACA,MAAME,GAAG,GAAGhB,kBAAkB,GAAGa,QAAjC;AACA,MAAMI,CAAC,GAAGC,MAAM,CAACN,SAAD,CAAhB;AAEA,MAAMO,CAAC,GAAGC,YAAY,CAACH,CAAD,CAAtB;AAEA,MAAMI,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;;AAEM,SAASG,eAAT,CACLjB,SADK,EAELC,QAFK,EAGLC,SAHK,EAIK;AACV,wBAA4BH,cAAc,CAACC,SAAD,EAAYC,QAAZ,EAAsBC,SAAtB,CAA1C;AAAA,MAAOU,OAAP,mBAAOA,OAAP;AAAA,MAAgBG,QAAhB,mBAAgBA,QAAhB;;AAEA,MAAMG,QAAQ,GAAGN,OAAO,GAAGvB,IAAI,CAACC,EAAhC;AAGA,SAAO,CAAC,CAACD,IAAI,CAAC8B,GAAL,CAASD,QAAT,CAAF,EAAsB7B,IAAI,CAAC+B,GAAL,CAASF,QAAT,CAAtB,EAA0C,CAAC7B,IAAI,CAAC8B,GAAL,CAASJ,QAAT,CAA3C,CAAP;AACD;;AAED,SAASM,WAAT,CAAqBrB,SAArB,EAAgC;AAC9B,SAAOA,SAAS,GAAGT,SAAZ,GAAwB,GAAxB,GAA8BC,MAArC;AACD;;AAED,SAASc,MAAT,CAAgBN,SAAhB,EAA2B;AACzB,SAAOqB,WAAW,CAACrB,SAAD,CAAX,GAAyBP,MAAhC;AACD;;AAED,SAAS6B,iBAAT,CAA2BC,iBAA3B,EAA8CC,CAA9C,EAAiD;AAC/C,MAAMC,MAAM,GAAGF,iBAAf;AACA,SAAOlC,IAAI,CAACqC,KAAL,CAAWrC,IAAI,CAAC8B,GAAL,CAASM,MAAT,IAAmBpC,IAAI,CAAC+B,GAAL,CAAS1B,CAAT,CAAnB,GAAiCL,IAAI,CAACsC,GAAL,CAASH,CAAT,IAAcnC,IAAI,CAAC8B,GAAL,CAASzB,CAAT,CAA1D,EAAuEL,IAAI,CAAC+B,GAAL,CAASK,MAAT,CAAvE,CAAP;AACD;;AAED,SAASG,cAAT,CAAwBL,iBAAxB,EAA2CC,CAA3C,EAA8C;AAC5C,MAAMC,MAAM,GAAGF,iBAAf;AACA,SAAOlC,IAAI,CAACwC,IAAL,CAAUxC,IAAI,CAAC8B,GAAL,CAASK,CAAT,IAAcnC,IAAI,CAAC+B,GAAL,CAAS1B,CAAT,CAAd,GAA4BL,IAAI,CAAC+B,GAAL,CAASI,CAAT,IAAcnC,IAAI,CAAC8B,GAAL,CAASzB,CAAT,CAAd,GAA4BL,IAAI,CAAC8B,GAAL,CAASM,MAAT,CAAlE,CAAP;AACD;;AAED,SAASZ,UAAT,CAAoBiB,SAApB,EAA+BC,iBAA/B,EAAkDjB,WAAlD,EAA+D;AAC7D,MAAML,CAAC,GAAGqB,SAAV;AACA,MAAM1B,GAAG,GAAG2B,iBAAZ;AACA,MAAMC,KAAK,GAAGlB,WAAd;AACA,SAAOzB,IAAI,CAACqC,KAAL,CAAWrC,IAAI,CAAC8B,GAAL,CAASV,CAAT,CAAX,EAAwBpB,IAAI,CAAC+B,GAAL,CAASX,CAAT,IAAcpB,IAAI,CAAC8B,GAAL,CAASf,GAAT,CAAd,GAA8Bf,IAAI,CAACsC,GAAL,CAASK,KAAT,IAAkB3C,IAAI,CAAC+B,GAAL,CAAShB,GAAT,CAAxE,CAAP;AACD;;AAED,SAASY,WAAT,CAAqBc,SAArB,EAAgCC,iBAAhC,EAAmDjB,WAAnD,EAAgE;AAC9D,MAAML,CAAC,GAAGqB,SAAV;AACA,MAAM1B,GAAG,GAAG2B,iBAAZ;AACA,MAAMC,KAAK,GAAGlB,WAAd;AACA,SAAOzB,IAAI,CAACwC,IAAL,CAAUxC,IAAI,CAAC8B,GAAL,CAASf,GAAT,IAAgBf,IAAI,CAAC8B,GAAL,CAASa,KAAT,CAAhB,GAAkC3C,IAAI,CAAC+B,GAAL,CAAShB,GAAT,IAAgBf,IAAI,CAAC+B,GAAL,CAASY,KAAT,CAAhB,GAAkC3C,IAAI,CAAC+B,GAAL,CAASX,CAAT,CAA9E,CAAP;AACD;;AAID,SAASC,eAAT,CAAyBuB,KAAzB,EAAgC9B,sBAAhC,EAAwD;AACtD,SAAOf,kBAAkB,IAAIS,MAAM,GAAGC,MAAM,GAAGmC,KAAtB,CAAlB,GAAiD9B,sBAAxD;AACD;;AAED,SAAS+B,mBAAT,CAA6BC,IAA7B,EAAmC;AACjC,SAAO/C,kBAAkB,IAAIO,EAAE,GAAGC,EAAE,GAAGuC,IAAd,CAAzB;AACD;;AAED,SAASC,oBAAT,CAA8BC,WAA9B,EAA2C;AACzC,MAAMC,CAAC,GAAGD,WAAV;AAEA,MAAME,CAAC,GACLnD,kBAAkB,IAAI,SAASC,IAAI,CAAC8B,GAAL,CAASmB,CAAT,CAAT,GAAuB,OAAOjD,IAAI,CAAC8B,GAAL,CAAS,IAAImB,CAAb,CAA9B,GAAgD,SAASjD,IAAI,CAAC8B,GAAL,CAAS,IAAImB,CAAb,CAA7D,CADpB;AAGA,MAAME,CAAC,GAAGpD,kBAAkB,GAAG,QAA/B;AAEA,SAAOkD,CAAC,GAAGC,CAAJ,GAAQC,CAAR,GAAYnD,IAAI,CAACC,EAAxB;AACD;;AAED,SAASkB,YAAT,CAAsByB,KAAtB,EAA6B;AAC3B,MAAMK,CAAC,GAAGJ,mBAAmB,CAACD,KAAD,CAA7B;AACA,MAAMQ,CAAC,GAAGL,oBAAoB,CAACE,CAAD,CAA9B;AAEA,SAAO;AACLxB,IAAAA,WAAW,EAAEc,cAAc,CAACa,CAAD,EAAI,CAAJ,CADtB;AAEL9B,IAAAA,cAAc,EAAEW,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  // convert azimuth from 0 at south to be 0 at north\n  const azimuthN = azimuth + Math.PI;\n\n  // solar position to light direction\n  return [-Math.sin(azimuthN), Math.cos(azimuthN), -Math.sin(altitude)];\n}\n\nfunction toJulianDay(timestamp) {\n  return timestamp / DAY_IN_MS - 0.5 + JD1970;\n}\n\nfunction toDays(timestamp) {\n  return toJulianDay(timestamp) - JD2000;\n}\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}\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}\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}\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\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}\n\nfunction getSolarMeanAnomaly(days) {\n  return DEGREES_TO_RADIANS * (M0 + M1 * days);\n}\n\nfunction getEclipticLongitude(meanAnomaly) {\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) {\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"}
+{"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,IAAMA,kBAAkB,GAAGC,IAAI,CAACC,EAAL,GAAU,GAArC;AAEA,IAAMC,SAAS,GAAG,OAAO,EAAP,GAAY,EAAZ,GAAiB,EAAnC;AACA,IAAMC,MAAM,GAAG,OAAf;AACA,IAAMC,MAAM,GAAG,OAAf;AAGA,IAAMC,CAAC,GAAGN,kBAAkB,GAAG,OAA/B;AAIA,IAAMO,EAAE,GAAG,QAAX;AACA,IAAMC,EAAE,GAAG,UAAX;AAEA,IAAMC,MAAM,GAAG,OAAf;AACA,IAAMC,MAAM,GAAG,WAAf;;AAkBO,SAASC,cAAT,CACLC,SADK,EAELC,QAFK,EAGLC,SAHK,EAIc;AACnB,MAAMC,sBAAsB,GAAGf,kBAAkB,GAAG,CAACc,SAArD;AACA,MAAME,GAAG,GAAGhB,kBAAkB,GAAGa,QAAjC;AACA,MAAMI,CAAC,GAAGC,MAAM,CAACN,SAAD,CAAhB;AAEA,MAAMO,CAAC,GAAGC,YAAY,CAACH,CAAD,CAAtB;AAEA,MAAMI,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;;AAEM,SAASG,eAAT,CACLjB,SADK,EAELC,QAFK,EAGLC,SAHK,EAIK;AACV,wBAA4BH,cAAc,CAACC,SAAD,EAAYC,QAAZ,EAAsBC,SAAtB,CAA1C;AAAA,MAAOU,OAAP,mBAAOA,OAAP;AAAA,MAAgBG,QAAhB,mBAAgBA,QAAhB;;AAGA,SAAO,CACL1B,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,MAAMqB,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,MAAMC,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,MAAMC,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,MAAML,CAAC,GAAGsB,SAAV;AACA,MAAM3B,GAAG,GAAG4B,iBAAZ;AACA,MAAMC,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,MAAML,CAAC,GAAGsB,SAAV;AACA,MAAM3B,GAAG,GAAG4B,iBAAZ;AACA,MAAMC,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,MAAMC,CAAC,GAAGD,WAAV;AAEA,MAAME,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,MAAME,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,MAAMK,CAAC,GAAGJ,mBAAmB,CAACD,KAAD,CAA7B;AACA,MAAMQ,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"}

package/dist/esm/suncalc.js

@@ -23,12 +23,12 @@ export function getSunDirection(timestamp, latitude, longitude) {
     azimuth,
     altitude
   } = getSunPosition(timestamp, latitude, longitude);
-  const azimuthN = azimuth + Math.PI;
-  return [-Math.sin(azimuthN), Math.cos(azimuthN), -Math.sin(altitude)];
+  return [Math.sin(azimuth) * Math.cos(altitude), Math.cos(azimuth) * Math.cos(altitude), -Math.sin(altitude)];
 }
 
 function toJulianDay(timestamp) {
-  return timestamp / 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) {

package/dist/esm/suncalc.js.map

@@ -1 +1 @@
-{"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","azimuthN","sin","cos","toJulianDay","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;AAEA,QAAMgB,QAAQ,GAAGN,OAAO,GAAGvB,IAAI,CAACC,EAAhC;AAGA,SAAO,CAAC,CAACD,IAAI,CAAC8B,GAAL,CAASD,QAAT,CAAF,EAAsB7B,IAAI,CAAC+B,GAAL,CAASF,QAAT,CAAtB,EAA0C,CAAC7B,IAAI,CAAC8B,GAAL,CAASJ,QAAT,CAA3C,CAAP;AACD;;AAED,SAASM,WAAT,CAAqBrB,SAArB,EAAgC;AAC9B,SAAOA,SAAS,GAAGT,SAAZ,GAAwB,GAAxB,GAA8BC,MAArC;AACD;;AAED,SAASc,MAAT,CAAgBN,SAAhB,EAA2B;AACzB,SAAOqB,WAAW,CAACrB,SAAD,CAAX,GAAyBP,MAAhC;AACD;;AAED,SAAS6B,iBAAT,CAA2BC,iBAA3B,EAA8CC,CAA9C,EAAiD;AAC/C,QAAMC,MAAM,GAAGF,iBAAf;AACA,SAAOlC,IAAI,CAACqC,KAAL,CAAWrC,IAAI,CAAC8B,GAAL,CAASM,MAAT,IAAmBpC,IAAI,CAAC+B,GAAL,CAAS1B,CAAT,CAAnB,GAAiCL,IAAI,CAACsC,GAAL,CAASH,CAAT,IAAcnC,IAAI,CAAC8B,GAAL,CAASzB,CAAT,CAA1D,EAAuEL,IAAI,CAAC+B,GAAL,CAASK,MAAT,CAAvE,CAAP;AACD;;AAED,SAASG,cAAT,CAAwBL,iBAAxB,EAA2CC,CAA3C,EAA8C;AAC5C,QAAMC,MAAM,GAAGF,iBAAf;AACA,SAAOlC,IAAI,CAACwC,IAAL,CAAUxC,IAAI,CAAC8B,GAAL,CAASK,CAAT,IAAcnC,IAAI,CAAC+B,GAAL,CAAS1B,CAAT,CAAd,GAA4BL,IAAI,CAAC+B,GAAL,CAASI,CAAT,IAAcnC,IAAI,CAAC8B,GAAL,CAASzB,CAAT,CAAd,GAA4BL,IAAI,CAAC8B,GAAL,CAASM,MAAT,CAAlE,CAAP;AACD;;AAED,SAASZ,UAAT,CAAoBiB,SAApB,EAA+BC,iBAA/B,EAAkDjB,WAAlD,EAA+D;AAC7D,QAAML,CAAC,GAAGqB,SAAV;AACA,QAAM1B,GAAG,GAAG2B,iBAAZ;AACA,QAAMC,KAAK,GAAGlB,WAAd;AACA,SAAOzB,IAAI,CAACqC,KAAL,CAAWrC,IAAI,CAAC8B,GAAL,CAASV,CAAT,CAAX,EAAwBpB,IAAI,CAAC+B,GAAL,CAASX,CAAT,IAAcpB,IAAI,CAAC8B,GAAL,CAASf,GAAT,CAAd,GAA8Bf,IAAI,CAACsC,GAAL,CAASK,KAAT,IAAkB3C,IAAI,CAAC+B,GAAL,CAAShB,GAAT,CAAxE,CAAP;AACD;;AAED,SAASY,WAAT,CAAqBc,SAArB,EAAgCC,iBAAhC,EAAmDjB,WAAnD,EAAgE;AAC9D,QAAML,CAAC,GAAGqB,SAAV;AACA,QAAM1B,GAAG,GAAG2B,iBAAZ;AACA,QAAMC,KAAK,GAAGlB,WAAd;AACA,SAAOzB,IAAI,CAACwC,IAAL,CAAUxC,IAAI,CAAC8B,GAAL,CAASf,GAAT,IAAgBf,IAAI,CAAC8B,GAAL,CAASa,KAAT,CAAhB,GAAkC3C,IAAI,CAAC+B,GAAL,CAAShB,GAAT,IAAgBf,IAAI,CAAC+B,GAAL,CAASY,KAAT,CAAhB,GAAkC3C,IAAI,CAAC+B,GAAL,CAASX,CAAT,CAA9E,CAAP;AACD;;AAID,SAASC,eAAT,CAAyBuB,KAAzB,EAAgC9B,sBAAhC,EAAwD;AACtD,SAAOf,kBAAkB,IAAIS,MAAM,GAAGC,MAAM,GAAGmC,KAAtB,CAAlB,GAAiD9B,sBAAxD;AACD;;AAED,SAAS+B,mBAAT,CAA6BC,IAA7B,EAAmC;AACjC,SAAO/C,kBAAkB,IAAIO,EAAE,GAAGC,EAAE,GAAGuC,IAAd,CAAzB;AACD;;AAED,SAASC,oBAAT,CAA8BC,WAA9B,EAA2C;AACzC,QAAMC,CAAC,GAAGD,WAAV;AAEA,QAAME,CAAC,GACLnD,kBAAkB,IAAI,SAASC,IAAI,CAAC8B,GAAL,CAASmB,CAAT,CAAT,GAAuB,OAAOjD,IAAI,CAAC8B,GAAL,CAAS,IAAImB,CAAb,CAA9B,GAAgD,SAASjD,IAAI,CAAC8B,GAAL,CAAS,IAAImB,CAAb,CAA7D,CADpB;AAGA,QAAME,CAAC,GAAGpD,kBAAkB,GAAG,QAA/B;AAEA,SAAOkD,CAAC,GAAGC,CAAJ,GAAQC,CAAR,GAAYnD,IAAI,CAACC,EAAxB;AACD;;AAED,SAASkB,YAAT,CAAsByB,KAAtB,EAA6B;AAC3B,QAAMK,CAAC,GAAGJ,mBAAmB,CAACD,KAAD,CAA7B;AACA,QAAMQ,CAAC,GAAGL,oBAAoB,CAACE,CAAD,CAA9B;AAEA,SAAO;AACLxB,IAAAA,WAAW,EAAEc,cAAc,CAACa,CAAD,EAAI,CAAJ,CADtB;AAEL9B,IAAAA,cAAc,EAAEW,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  // convert azimuth from 0 at south to be 0 at north\n  const azimuthN = azimuth + Math.PI;\n\n  // solar position to light direction\n  return [-Math.sin(azimuthN), Math.cos(azimuthN), -Math.sin(altitude)];\n}\n\nfunction toJulianDay(timestamp) {\n  return timestamp / DAY_IN_MS - 0.5 + JD1970;\n}\n\nfunction toDays(timestamp) {\n  return toJulianDay(timestamp) - JD2000;\n}\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}\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}\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}\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\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}\n\nfunction getSolarMeanAnomaly(days) {\n  return DEGREES_TO_RADIANS * (M0 + M1 * days);\n}\n\nfunction getEclipticLongitude(meanAnomaly) {\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) {\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"}
+{"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"}

package/dist/suncalc.d.ts.map

@@ -1 +1 @@
-{"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,CAOV"}
+{"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

@@ -29,13 +29,16 @@ export function getSunPosition(timestamp, latitude, longitude) {
 }
 export function getSunDirection(timestamp, latitude, longitude) {
     const { azimuth, altitude } = getSunPosition(timestamp, latitude, longitude);
-    // convert azimuth from 0 at south to be 0 at north
-    const azimuthN = azimuth + Math.PI;
     // solar position to light direction
-    return [-Math.sin(azimuthN), Math.cos(azimuthN), -Math.sin(altitude)];
+    return [
+        Math.sin(azimuth) * Math.cos(altitude),
+        Math.cos(azimuth) * Math.cos(altitude),
+        -Math.sin(altitude)
+    ];
 }
 function toJulianDay(timestamp) {
-    return timestamp / 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;

package/package.json

@@ -5,7 +5,7 @@
   "publishConfig": {
     "access": "public"
   },
-  "version": "3.6.0-alpha.4",
+  "version": "3.6.0-alpha.7",
   "keywords": [
     "javascript",
     "math",
@@ -28,5 +28,5 @@
   "dependencies": {
     "@babel/runtime": "^7.12.0"
   },
-  "gitHead": "d2ce86da8b6f07f0af0236818138442f68e1d8e4"
+  "gitHead": "166f7411d07d545e6a3f247b55b56f5cb8efccde"
 }

package/src/suncalc.ts

@@ -56,39 +56,42 @@ export function getSunDirection(
   longitude: number
 ): number[] {
   const {azimuth, altitude} = getSunPosition(timestamp, latitude, longitude);
-  // convert azimuth from 0 at south to be 0 at north
-  const azimuthN = azimuth + Math.PI;
 
   // solar position to light direction
-  return [-Math.sin(azimuthN), Math.cos(azimuthN), -Math.sin(altitude)];
+  return [
+    Math.sin(azimuth) * Math.cos(altitude),
+    Math.cos(azimuth) * Math.cos(altitude),
+    -Math.sin(altitude)
+  ];
 }
 
-function toJulianDay(timestamp) {
-  return timestamp / DAY_IN_MS - 0.5 + JD1970;
+function toJulianDay(timestamp: number | Date): number {
+  const ts = typeof timestamp === 'number' ? timestamp : timestamp.getTime();
+  return ts / DAY_IN_MS - 0.5 + JD1970;
 }
 
-function toDays(timestamp) {
+function toDays(timestamp: number | Date): number {
   return toJulianDay(timestamp) - JD2000;
 }
 
-function getRightAscension(eclipticLongitude, b) {
+function getRightAscension(eclipticLongitude: number, b: number): number {
   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) {
+function getDeclination(eclipticLongitude: number, b: number): number {
   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) {
+function getAzimuth(hourAngle: number, latitudeInRadians: number, declination: number): number {
   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) {
+function getAltitude(hourAngle: number, latitudeInRadians: number, declination: number): number {
   const H = hourAngle;
   const phi = latitudeInRadians;
   const delta = declination;
@@ -97,15 +100,15 @@ function getAltitude(hourAngle, latitudeInRadians, declination) {
 
 // https://www.aa.quae.nl/en/reken/zonpositie.html
 // "The Observer section"
-function getSiderealTime(dates, longitudeWestInRadians) {
+function getSiderealTime(dates: number, longitudeWestInRadians: number): number {
   return DEGREES_TO_RADIANS * (THETA0 + THETA1 * dates) - longitudeWestInRadians;
 }
 
-function getSolarMeanAnomaly(days) {
+function getSolarMeanAnomaly(days: number): number {
   return DEGREES_TO_RADIANS * (M0 + M1 * days);
 }
 
-function getEclipticLongitude(meanAnomaly) {
+function getEclipticLongitude(meanAnomaly: number): number {
   const M = meanAnomaly;
   // equation of center
   const C =
@@ -116,7 +119,10 @@ function getEclipticLongitude(meanAnomaly) {
   return M + C + P + Math.PI;
 }
 
-function getSunCoords(dates) {
+function getSunCoords(dates: number): {
+  declination: number;
+  rightAscension: number;
+} {
   const M = getSolarMeanAnomaly(dates);
   const L = getEclipticLongitude(M);