@neaps/tide-predictor 0.0.5 → 0.1.1

package/dist/{tide-predictor.js → web/tide-predictor.js}

@@ -1,8 +1,8 @@
 (function (global, factory) {
   typeof exports === 'object' && typeof module !== 'undefined' ? module.exports = factory() :
   typeof define === 'function' && define.amd ? define(factory) :
-  (global = global || self, global.tidePredictor = factory());
-}(this, function () { 'use strict';
+  (global = typeof globalThis !== 'undefined' ? globalThis : global || self, global.tidePredictor = factory());
+})(this, (function () { 'use strict';
 
   const d2r = Math.PI / 180.0;
   const r2d = 180.0 / Math.PI;
@@ -36,7 +36,7 @@
       sexagesimalToDecimal(0, 0, 7.12),
       sexagesimalToDecimal(0, 0, 27.87),
       sexagesimalToDecimal(0, 0, 5.79),
-      sexagesimalToDecimal(0, 0, 2.45)
+      sexagesimalToDecimal(0, 0, 2.45),
     ].map((number, index) => {
       return number * Math.pow(1e-2, index)
     }),
@@ -45,7 +45,7 @@
       280.46645 - 357.5291,
       36000.76932 - 35999.0503,
       0.0003032 + 0.0001559,
-      0.00000048
+      0.00000048,
     ],
 
     solarLongitude: [280.46645, 36000.76983, 0.0003032],
@@ -56,7 +56,7 @@
       218.3164591,
       481267.88134236,
       -0.0013268,
-      1 / 538841.0 - 1 / 65194000.0
+      1 / 538841.0 - 1 / 65194000.0,
     ],
 
     lunarNode: [
@@ -64,7 +64,7 @@
       -1934.1361849,
       0.0020762,
       1 / 467410.0,
-      -1 / 60616000.0
+      -1 / 60616000.0,
     ],
 
     lunarPerigee: [
@@ -72,8 +72,8 @@
       4069.0137111,
       -0.0103238,
       -1 / 80053.0,
-      1 / 18999000.0
-    ]
+      1 / 18999000.0,
+    ],
   };
 
   // Evaluates a polynomial at argument
@@ -99,12 +99,12 @@
   };
 
   // Meeus formula 11.1
-  const T = t => {
+  const T = (t) => {
     return (JD(t) - 2451545.0) / 36525
   };
 
   // Meeus formula 7.1
-  const JD = t => {
+  const JD = (t) => {
     let Y = t.getFullYear();
     let M = t.getMonth() + 1;
     const D =
@@ -200,11 +200,11 @@
     return r2d * 0.5 * Math.atan(tan2nupp)
   };
 
-  const modulus = (a, b) => {
+  const modulus$1 = (a, b) => {
     return ((a % b) + b) % b
   };
 
-  const astro = time => {
+  const astro = (time) => {
     const result = {};
     const polynomials = {
       s: coefficients.lunarLongitude,
@@ -220,9 +220,9 @@
     // Polynomials are in T, that is Julian Centuries; we want our speeds to be
     // in the more convenient unit of degrees per hour.
     const dTdHour = 1 / (24 * 365.25 * 100);
-    Object.keys(polynomials).forEach(name => {
+    Object.keys(polynomials).forEach((name) => {
       result[name] = {
-        value: modulus(polynomial(polynomials[name], T(time)), 360.0),
+        value: modulus$1(polynomial(polynomials[name], T(time)), 360.0),
         speed: derivativePolynomial(polynomials[name], T(time)) * dTdHour
       };
     });
@@ -237,10 +237,10 @@
       nup: _nup,
       nupp: _nupp
     };
-    Object.keys(functions).forEach(name => {
+    Object.keys(functions).forEach((name) => {
       const functionCall = functions[name];
       result[name] = {
-        value: modulus(
+        value: modulus$1(
           functionCall(result.N.value, result.i.value, result.omega.value),
           360.0
         ),
@@ -272,7 +272,7 @@
     return result
   };
 
-  const modulus$1 = (a, b) => {
+  const modulus = (a, b) => {
     return ((a % b) + b) % b
   };
 
@@ -318,7 +318,7 @@
       const amplitudes = [];
       let result = 0;
 
-      constituents.forEach(constituent => {
+      constituents.forEach((constituent) => {
         const amplitude = constituent.amplitude;
         const phase = constituent._phase;
         const f = modelF[constituent.name];
@@ -328,7 +328,7 @@
         amplitudes.push(amplitude * f * Math.cos(speed * hour + (V0 + u) - phase));
       });
       // sum up each row
-      amplitudes.forEach(item => {
+      amplitudes.forEach((item) => {
         result += item;
       });
       return result
@@ -336,7 +336,7 @@
 
     const prediction = {};
 
-    prediction.getExtremesPrediction = options => {
+    prediction.getExtremesPrediction = (options) => {
       const { labels, offsets } = typeof options !== 'undefined' ? options : {};
       const results = [];
       const { baseSpeed, u, f, baseValue } = prepare();
@@ -395,8 +395,8 @@
       timeline.items.forEach((time, index) => {
         const hour = timeline.hours[index];
         const prediction = {
-          time: time,
-          hour: hour,
+          time,
+          hour,
           level: getLevel(hour, baseSpeed, u[index], f[index], baseValue)
         };
 
@@ -412,31 +412,31 @@
       const baseSpeed = {};
       const u = [];
       const f = [];
-      constituents.forEach(constituent => {
+      constituents.forEach((constituent) => {
         const value = constituent._model.value(baseAstro);
         const speed = constituent._model.speed(baseAstro);
         baseValue[constituent.name] = d2r * value;
         baseSpeed[constituent.name] = d2r * speed;
       });
-      timeline.items.forEach(time => {
+      timeline.items.forEach((time) => {
         const uItem = {};
         const fItem = {};
         const itemAstro = astro(time);
-        constituents.forEach(constituent => {
-          const constituentU = modulus$1(constituent._model.u(itemAstro), 360);
+        constituents.forEach((constituent) => {
+          const constituentU = modulus(constituent._model.u(itemAstro), 360);
 
           uItem[constituent.name] = d2r * constituentU;
-          fItem[constituent.name] = modulus$1(constituent._model.f(itemAstro), 360);
+          fItem[constituent.name] = modulus(constituent._model.f(itemAstro), 360);
         });
         u.push(uItem);
         f.push(fItem);
       });
 
       return {
-        baseValue: baseValue,
-        baseSpeed: baseSpeed,
-        u: u,
-        f: f
+        baseValue,
+        baseSpeed,
+        u,
+        f
       }
     };
 
@@ -588,7 +588,7 @@
 
     // Node factors u, see Table 2 of Schureman.
 
-    uZero(a) {
+    uZero() {
       return 0.0
     },
 
@@ -663,7 +663,7 @@
     })
   };
 
-  const astronimicDoodsonNumber = astro => {
+  const astronimicDoodsonNumber = (astro) => {
     return [
       astro['T+h-s'],
       astro.s,
@@ -675,17 +675,17 @@
     ]
   };
 
-  const astronomicSpeed = astro => {
+  const astronomicSpeed = (astro) => {
     const results = [];
-    astronimicDoodsonNumber(astro).forEach(number => {
+    astronimicDoodsonNumber(astro).forEach((number) => {
       results.push(number.speed);
     });
     return results
   };
 
-  const astronomicValues = astro => {
+  const astronomicValues = (astro) => {
     const results = [];
-    astronimicDoodsonNumber(astro).forEach(number => {
+    astronimicDoodsonNumber(astro).forEach((number) => {
       results.push(number.value);
     });
     return results
@@ -697,11 +697,10 @@
     }
 
     const constituent = {
-      name: name,
+      name,
+      coefficients,
 
-      coefficients: coefficients,
-
-      value: astro => {
+      value: (astro) => {
         return dotArray(coefficients, astronomicValues(astro))
       },
 
@@ -729,11 +728,10 @@
     });
 
     const compoundConstituent = {
-      name: name,
-
-      coefficients: coefficients,
+      name,
+      coefficients,
 
-      speed: astro => {
+      speed: (astro) => {
         let speed = 0;
         members.forEach(({ constituent, factor }) => {
           speed += constituent.speed(astro) * factor;
@@ -741,7 +739,7 @@
         return speed
       },
 
-      value: astro => {
+      value: (astro) => {
         let value = 0;
         members.forEach(({ constituent, factor }) => {
           value += constituent.value(astro) * factor;
@@ -749,7 +747,7 @@
         return value
       },
 
-      u: astro => {
+      u: (astro) => {
         let u = 0;
         members.forEach(({ constituent, factor }) => {
           u += constituent.u(astro) * factor;
@@ -757,7 +755,7 @@
         return u
       },
 
-      f: astro => {
+      f: (astro) => {
         const f = [];
         members.forEach(({ constituent, factor }) => {
           f.push(Math.pow(constituent.f(astro), Math.abs(factor)));
@@ -812,10 +810,10 @@
   constituents.M3 = constituentFactory(
     'M3',
     [3, 0, 0, 0, 0, 0, 0],
-    a => {
+    (a) => {
       return corrections.uModd(a, 3)
     },
-    a => {
+    (a) => {
       return corrections.fModd(a, 3)
     }
   );
@@ -885,7 +883,7 @@
     { constituent: constituents.M2, factor: 4 }
   ]);
 
-  const getDate = time => {
+  const getDate = (time) => {
     if (time instanceof Date) {
       return time
     }
@@ -897,20 +895,20 @@
 
   const getTimeline = (start, end, seconds) => {
     seconds = typeof seconds !== 'undefined' ? seconds : 10 * 60;
-    const timeline = [];
+    const items = [];
     const endTime = end.getTime() / 1000;
     let lastTime = start.getTime() / 1000;
     const startTime = lastTime;
     const hours = [];
     while (lastTime <= endTime) {
-      timeline.push(new Date(lastTime * 1000));
+      items.push(new Date(lastTime * 1000));
       hours.push((lastTime - startTime) / (60 * 60));
       lastTime += seconds;
     }
 
     return {
-      items: timeline,
-      hours: hours
+      items,
+      hours
     }
   };
 
@@ -919,7 +917,7 @@
       throw new Error('Harmonic constituents are not an array')
     }
     const constituents$1 = [];
-    harmonicConstituents.forEach((constituent, index) => {
+    harmonicConstituents.forEach((constituent) => {
       if (typeof constituent.name === 'undefined') {
         throw new Error('Harmonic constituents must have a name property')
       }
@@ -953,13 +951,13 @@
       return harmonics
     };
 
-    harmonics.prediction = options => {
+    harmonics.prediction = (options) => {
       options =
         typeof options !== 'undefined' ? options : { timeFidelity: 10 * 60 };
       return predictionFactory({
         timeline: getTimeline(start, end, options.timeFidelity),
         constituents: constituents$1,
-        start: start
+        start
       })
     };
 
@@ -974,7 +972,7 @@
     };
 
     if (typeof options !== 'undefined') {
-      Object.keys(harmonicsOptions).forEach(key => {
+      Object.keys(harmonicsOptions).forEach((key) => {
         if (typeof options[key] !== 'undefined') {
           harmonicsOptions[key] = options[key];
         }
@@ -992,7 +990,7 @@
       getExtremesPrediction: ({ start, end, labels, offsets, timeFidelity }) => {
         return harmonicsFactory(harmonicsOptions)
           .setTimeSpan(start, end)
-          .prediction({ timeFidelity: timeFidelity })
+          .prediction({ timeFidelity })
           .getExtremesPrediction(labels, offsets)
       },
 

package/dist/web/tide-predictor.min.js

@@ -0,0 +1 @@
+(function(global,factory){typeof exports==="object"&&typeof module!=="undefined"?module.exports=factory():typeof define==="function"&&define.amd?define(factory):(global=typeof globalThis!=="undefined"?globalThis:global||self,global.tidePredictor=factory())})(this,function(){"use strict";const d2r=Math.PI/180;const r2d=180/Math.PI;const sexagesimalToDecimal=(degrees,arcmins,arcsecs,mas,muas)=>{arcmins=typeof arcmins!=="undefined"?arcmins:0;arcsecs=typeof arcsecs!=="undefined"?arcsecs:0;mas=typeof mas!=="undefined"?mas:0;muas=typeof muas!=="undefined"?muas:0;return degrees+arcmins/60+arcsecs/(60*60)+mas/(60*60*1e3)+muas/(60*60*1e6)};const coefficients={terrestrialObliquity:[sexagesimalToDecimal(23,26,21.448),-sexagesimalToDecimal(0,0,4680.93),-sexagesimalToDecimal(0,0,1.55),sexagesimalToDecimal(0,0,1999.25),-sexagesimalToDecimal(0,0,51.38),-sexagesimalToDecimal(0,0,249.67),-sexagesimalToDecimal(0,0,39.05),sexagesimalToDecimal(0,0,7.12),sexagesimalToDecimal(0,0,27.87),sexagesimalToDecimal(0,0,5.79),sexagesimalToDecimal(0,0,2.45)].map((number,index)=>{return number*Math.pow(.01,index)}),solarPerigee:[280.46645-357.5291,36000.76932-35999.0503,3032e-7+1559e-7,48e-8],solarLongitude:[280.46645,36000.76983,3032e-7],lunarInclination:[5.145],lunarLongitude:[218.3164591,481267.88134236,-.0013268,1/538841-1/65194e3],lunarNode:[125.044555,-1934.1361849,.0020762,1/467410,-1/60616e3],lunarPerigee:[83.353243,4069.0137111,-.0103238,-1/80053,1/18999e3]};const polynomial=(coefficients,argument)=>{const result=[];coefficients.forEach((coefficient,index)=>{result.push(coefficient*Math.pow(argument,index))});return result.reduce((a,b)=>{return a+b})};const derivativePolynomial=(coefficients,argument)=>{const result=[];coefficients.forEach((coefficient,index)=>{result.push(coefficient*index*Math.pow(argument,index-1))});return result.reduce((a,b)=>{return a+b})};const T=t=>{return(JD(t)-2451545)/36525};const JD=t=>{let Y=t.getFullYear();let M=t.getMonth()+1;const D=t.getDate()+t.getHours()/24+t.getMinutes()/(24*60)+t.getSeconds()/(24*60*60)+t.getMilliseconds()/(24*60*60*1e6);if(M<=2){Y=Y-1;M=M+12}const A=Math.floor(Y/100);const B=2-A+Math.floor(A/4);return Math.floor(365.25*(Y+4716))+Math.floor(30.6001*(M+1))+D+B-1524.5};const _I=(N,i,omega)=>{N=d2r*N;i=d2r*i;omega=d2r*omega;const cosI=Math.cos(i)*Math.cos(omega)-Math.sin(i)*Math.sin(omega)*Math.cos(N);return r2d*Math.acos(cosI)};const _xi=(N,i,omega)=>{N=d2r*N;i=d2r*i;omega=d2r*omega;let e1=Math.cos(.5*(omega-i))/Math.cos(.5*(omega+i))*Math.tan(.5*N);let e2=Math.sin(.5*(omega-i))/Math.sin(.5*(omega+i))*Math.tan(.5*N);e1=Math.atan(e1);e2=Math.atan(e2);e1=e1-.5*N;e2=e2-.5*N;return-(e1+e2)*r2d};const _nu=(N,i,omega)=>{N=d2r*N;i=d2r*i;omega=d2r*omega;let e1=Math.cos(.5*(omega-i))/Math.cos(.5*(omega+i))*Math.tan(.5*N);let e2=Math.sin(.5*(omega-i))/Math.sin(.5*(omega+i))*Math.tan(.5*N);e1=Math.atan(e1);e2=Math.atan(e2);e1=e1-.5*N;e2=e2-.5*N;return(e1-e2)*r2d};const _nup=(N,i,omega)=>{const I=d2r*_I(N,i,omega);const nu=d2r*_nu(N,i,omega);return r2d*Math.atan(Math.sin(2*I)*Math.sin(nu)/(Math.sin(2*I)*Math.cos(nu)+.3347))};const _nupp=(N,i,omega)=>{const I=d2r*_I(N,i,omega);const nu=d2r*_nu(N,i,omega);const tan2nupp=Math.sin(I)**2*Math.sin(2*nu)/(Math.sin(I)**2*Math.cos(2*nu)+.0727);return r2d*.5*Math.atan(tan2nupp)};const modulus$1=(a,b)=>{return(a%b+b)%b};const astro=time=>{const result={};const polynomials={s:coefficients.lunarLongitude,h:coefficients.solarLongitude,p:coefficients.lunarPerigee,N:coefficients.lunarNode,pp:coefficients.solarPerigee,90:[90],omega:coefficients.terrestrialObliquity,i:coefficients.lunarInclination};const dTdHour=1/(24*365.25*100);Object.keys(polynomials).forEach(name=>{result[name]={value:modulus$1(polynomial(polynomials[name],T(time)),360),speed:derivativePolynomial(polynomials[name],T(time))*dTdHour}});const functions={I:_I,xi:_xi,nu:_nu,nup:_nup,nupp:_nupp};Object.keys(functions).forEach(name=>{const functionCall=functions[name];result[name]={value:modulus$1(functionCall(result.N.value,result.i.value,result.omega.value),360),speed:null}});const hour={value:(JD(time)-Math.floor(JD(time)))*360,speed:15};result["T+h-s"]={value:hour.value+result.h.value-result.s.value,speed:hour.speed+result.h.speed-result.s.speed};result.P={value:result.p.value-result.xi.value%360,speed:null};return result};const modulus=(a,b)=>{return(a%b+b)%b};const addExtremesOffsets=(extreme,offsets)=>{if(typeof offsets==="undefined"||!offsets){return extreme}if(extreme.high&&offsets.height_offset&&offsets.height_offset.high){extreme.level*=offsets.height_offset.high}if(extreme.low&&offsets.height_offset&&offsets.height_offset.low){extreme.level*=offsets.height_offset.low}if(extreme.high&&offsets.time_offset&&offsets.time_offset.high){extreme.time=new Date(extreme.time.getTime()+offsets.time_offset.high*60*1e3)}if(extreme.low&&offsets.time_offset&&offsets.time_offset.low){extreme.time=new Date(extreme.time.getTime()+offsets.time_offset.low*60*1e3)}return extreme};const getExtremeLabel=(label,highLowLabels)=>{if(typeof highLowLabels!=="undefined"&&typeof highLowLabels[label]!=="undefined"){return highLowLabels[label]}const labels={high:"High",low:"Low"};return labels[label]};const predictionFactory=({timeline,constituents,start})=>{const getLevel=(hour,modelBaseSpeed,modelU,modelF,modelBaseValue)=>{const amplitudes=[];let result=0;constituents.forEach(constituent=>{const amplitude=constituent.amplitude;const phase=constituent._phase;const f=modelF[constituent.name];const speed=modelBaseSpeed[constituent.name];const u=modelU[constituent.name];const V0=modelBaseValue[constituent.name];amplitudes.push(amplitude*f*Math.cos(speed*hour+(V0+u)-phase))});amplitudes.forEach(item=>{result+=item});return result};const prediction={};prediction.getExtremesPrediction=options=>{const{labels,offsets}=typeof options!=="undefined"?options:{};const results=[];const{baseSpeed,u,f,baseValue}=prepare();let goingUp=false;let goingDown=false;let lastLevel=getLevel(0,baseSpeed,u[0],f[0],baseValue);timeline.items.forEach((time,index)=>{const hour=timeline.hours[index];const level=getLevel(hour,baseSpeed,u[index],f[index],baseValue);if(level>lastLevel&&goingDown){results.push(addExtremesOffsets({time:timeline.items[index-1],level:lastLevel,high:false,low:true,label:getExtremeLabel("low",labels)},offsets))}if(level<lastLevel&&goingUp){results.push(addExtremesOffsets({time:timeline.items[index-1],level:lastLevel,high:true,low:false,label:getExtremeLabel("high",labels)},offsets))}if(level>lastLevel){goingUp=true;goingDown=false}if(level<lastLevel){goingUp=false;goingDown=true}lastLevel=level});return results};prediction.getTimelinePrediction=()=>{const results=[];const{baseSpeed,u,f,baseValue}=prepare();timeline.items.forEach((time,index)=>{const hour=timeline.hours[index];const prediction={time:time,hour:hour,level:getLevel(hour,baseSpeed,u[index],f[index],baseValue)};results.push(prediction)});return results};const prepare=()=>{const baseAstro=astro(start);const baseValue={};const baseSpeed={};const u=[];const f=[];constituents.forEach(constituent=>{const value=constituent._model.value(baseAstro);const speed=constituent._model.speed(baseAstro);baseValue[constituent.name]=d2r*value;baseSpeed[constituent.name]=d2r*speed});timeline.items.forEach(time=>{const uItem={};const fItem={};const itemAstro=astro(time);constituents.forEach(constituent=>{const constituentU=modulus(constituent._model.u(itemAstro),360);uItem[constituent.name]=d2r*constituentU;fItem[constituent.name]=modulus(constituent._model.f(itemAstro),360)});u.push(uItem);f.push(fItem)});return{baseValue:baseValue,baseSpeed:baseSpeed,u:u,f:f}};return Object.freeze(prediction)};const corrections={fUnity(){return 1},fMm(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const mean=(2/3-Math.pow(Math.sin(omega),2))*(1-3/2*Math.pow(Math.sin(i),2));return(2/3-Math.pow(Math.sin(I),2))/mean},fMf(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const mean=Math.pow(Math.sin(omega),2)*Math.pow(Math.cos(.5*i),4);return Math.pow(Math.sin(I),2)/mean},fO1(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const mean=Math.sin(omega)*Math.pow(Math.cos(.5*omega),2)*Math.pow(Math.cos(.5*i),4);return Math.sin(I)*Math.pow(Math.cos(.5*I),2)/mean},fJ1(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const mean=Math.sin(2*omega)*(1-3/2*Math.pow(Math.sin(i),2));return Math.sin(2*I)/mean},fOO1(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const mean=Math.sin(omega)*Math.pow(Math.sin(.5*omega),2)*Math.pow(Math.cos(.5*i),4);return Math.sin(I)*Math.pow(Math.sin(.5*I),2)/mean},fM2(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const mean=Math.pow(Math.cos(.5*omega),4)*Math.pow(Math.cos(.5*i),4);return Math.pow(Math.cos(.5*I),4)/mean},fK1(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const nu=d2r*a.nu.value;const sin2IcosnuMean=Math.sin(2*omega)*(1-3/2*Math.pow(Math.sin(i),2));const mean=.5023*sin2IcosnuMean+.1681;return Math.pow(.2523*Math.pow(Math.sin(2*I),2)+.1689*Math.sin(2*I)*Math.cos(nu)+.0283,.5)/mean},fL2(a){const P=d2r*a.P.value;const I=d2r*a.I.value;const rAInv=Math.pow(1-12*Math.pow(Math.tan(.5*I),2)*Math.cos(2*P)+36*Math.pow(Math.tan(.5*I),4),.5);return corrections.fM2(a)*rAInv},fK2(a){const omega=d2r*a.omega.value;const i=d2r*a.i.value;const I=d2r*a.I.value;const nu=d2r*a.nu.value;const sinsqIcos2nuMean=Math.sin(omega)**2*(1-3/2*Math.sin(i)**2);const mean=.5023*sinsqIcos2nuMean+.0365;return Math.pow(.2523*Math.pow(Math.sin(I),4)+.0367*Math.pow(Math.sin(I),2)*Math.cos(2*nu)+.0013,.5)/mean},fM1(a){const P=d2r*a.P.value;const I=d2r*a.I.value;const qAInv=Math.pow(.25+1.5*Math.cos(I)*Math.cos(2*P)*Math.pow(Math.cos(.5*I),-.5)+2.25*Math.pow(Math.cos(I),2)*Math.pow(Math.cos(.5*I),-4),.5);return corrections.fO1(a)*qAInv},fModd(a,n){return Math.pow(corrections.fM2(a),n/2)},uZero(){return 0},uMf(a){return-2*a.xi.value},uO1(a){return 2*a.xi.value-a.nu.value},uJ1(a){return-a.nu.value},uOO1(a){return-2*a.xi.value-a.nu.value},uM2(a){return 2*a.xi.value-2*a.nu.value},uK1(a){return-a.nup.value},uL2(a){const I=d2r*a.I.value;const P=d2r*a.P.value;const R=r2d*Math.atan(Math.sin(2*P)/(1/6*Math.pow(Math.tan(.5*I),-2)-Math.cos(2*P)));return 2*a.xi.value-2*a.nu.value-R},uK2(a){return-2*a.nupp.value},uM1(a){const I=d2r*a.I.value;const P=d2r*a.P.value;const Q=r2d*Math.atan((5*Math.cos(I)-1)/(7*Math.cos(I)+1)*Math.tan(P));return a.xi.value-a.nu.value+Q},uModd(a,n){return n/2*corrections.uM2(a)}};const dotArray=(a,b)=>{const results=[];a.forEach((value,index)=>{results.push(value*b[index])});return results.reduce((total,value)=>{return total+value})};const astronimicDoodsonNumber=astro=>{return[astro["T+h-s"],astro.s,astro.h,astro.p,astro.N,astro.pp,astro["90"]]};const astronomicSpeed=astro=>{const results=[];astronimicDoodsonNumber(astro).forEach(number=>{results.push(number.speed)});return results};const astronomicValues=astro=>{const results=[];astronimicDoodsonNumber(astro).forEach(number=>{results.push(number.value)});return results};const constituentFactory=(name,coefficients,u,f)=>{if(!coefficients){throw new Error("Coefficient must be defined for a constituent")}const constituent={name:name,coefficients:coefficients,value:astro=>{return dotArray(coefficients,astronomicValues(astro))},speed(astro){return dotArray(coefficients,astronomicSpeed(astro))},u:typeof u!=="undefined"?u:corrections.uZero,f:typeof f!=="undefined"?f:corrections.fUnity};return Object.freeze(constituent)};const compoundConstituentFactory=(name,members)=>{const coefficients=[];members.forEach(({constituent,factor})=>{constituent.coefficients.forEach((coefficient,index)=>{if(typeof coefficients[index]==="undefined"){coefficients[index]=0}coefficients[index]+=coefficient*factor})});const compoundConstituent={name:name,coefficients:coefficients,speed:astro=>{let speed=0;members.forEach(({constituent,factor})=>{speed+=constituent.speed(astro)*factor});return speed},value:astro=>{let value=0;members.forEach(({constituent,factor})=>{value+=constituent.value(astro)*factor});return value},u:astro=>{let u=0;members.forEach(({constituent,factor})=>{u+=constituent.u(astro)*factor});return u},f:astro=>{const f=[];members.forEach(({constituent,factor})=>{f.push(Math.pow(constituent.f(astro),Math.abs(factor)))});return f.reduce((previous,value)=>{return previous*value})}};return Object.freeze(compoundConstituent)};const constituents={};constituents.Z0=constituentFactory("Z0",[0,0,0,0,0,0,0],corrections.uZero,corrections.fUnity);constituents.SA=constituentFactory("Sa",[0,0,1,0,0,0,0],corrections.uZero,corrections.fUnity);constituents.SSA=constituentFactory("Ssa",[0,0,2,0,0,0,0],corrections.uZero,corrections.fUnity);constituents.MM=constituentFactory("MM",[0,1,0,-1,0,0,0],corrections.uZero,corrections.fMm);constituents.MF=constituentFactory("MF",[0,2,0,0,0,0,0],corrections.uMf,corrections.fMf);constituents.Q1=constituentFactory("Q1",[1,-2,0,1,0,0,1],corrections.uO1,corrections.fO1);constituents.O1=constituentFactory("O1",[1,-1,0,0,0,0,1],corrections.uO1,corrections.fO1);constituents.K1=constituentFactory("K1",[1,1,0,0,0,0,-1],corrections.uK1,corrections.fK1);constituents.J1=constituentFactory("J1",[1,2,0,-1,0,0,-1],corrections.uJ1,corrections.fJ1);constituents.M1=constituentFactory("M1",[1,0,0,0,0,0,1],corrections.uM1,corrections.fM1);constituents.P1=constituentFactory("P1",[1,1,-2,0,0,0,1],corrections.uZero,corrections.fUnity);constituents.S1=constituentFactory("S1",[1,1,-1,0,0,0,0],corrections.uZero,corrections.fUnity);constituents.OO1=constituentFactory("OO1",[1,3,0,0,0,0,-1],corrections.uOO1,corrections.fOO1);constituents["2N2"]=constituentFactory("2N2",[2,-2,0,2,0,0,0],corrections.uM2,corrections.fM2);constituents.N2=constituentFactory("N2",[2,-1,0,1,0,0,0],corrections.uM2,corrections.fM2);constituents.NU2=constituentFactory("NU2",[2,-1,2,-1,0,0,0],corrections.uM2,corrections.fM2);constituents.M2=constituentFactory("M2",[2,0,0,0,0,0,0],corrections.uM2,corrections.fM2);constituents.LAM2=constituentFactory("LAM2",[2,1,-2,1,0,0,2],corrections.uM2,corrections.fM2);constituents.L2=constituentFactory("L2",[2,1,0,-1,0,0,2],corrections.uL2,corrections.fL2);constituents.T2=constituentFactory("T2",[2,2,-3,0,0,1,0],corrections.uZero,corrections.fUnity);constituents.S2=constituentFactory("S2",[2,2,-2,0,0,0,0],corrections.uZero,corrections.fUnity);constituents.R2=constituentFactory("R2",[2,2,-1,0,0,-1,2],corrections.uZero,corrections.fUnity);constituents.K2=constituentFactory("K2",[2,2,0,0,0,0,0],corrections.uK2,corrections.fK2);constituents.M3=constituentFactory("M3",[3,0,0,0,0,0,0],a=>{return corrections.uModd(a,3)},a=>{return corrections.fModd(a,3)});constituents.MSF=compoundConstituentFactory("MSF",[{constituent:constituents.S2,factor:1},{constituent:constituents.M2,factor:-1}]);constituents["2Q1"]=compoundConstituentFactory("2Q1",[{constituent:constituents.N2,factor:1},{constituent:constituents.J1,factor:-1}]);constituents.RHO=compoundConstituentFactory("RHO",[{constituent:constituents.NU2,factor:1},{constituent:constituents.K1,factor:-1}]);constituents.MU2=compoundConstituentFactory("MU2",[{constituent:constituents.M2,factor:2},{constituent:constituents.S2,factor:-1}]);constituents["2SM2"]=compoundConstituentFactory("2SM2",[{constituent:constituents.S2,factor:2},{constituent:constituents.M2,factor:-1}]);constituents["2MK3"]=compoundConstituentFactory("2MK3",[{constituent:constituents.M2,factor:1},{constituent:constituents.O1,factor:1}]);constituents.MK3=compoundConstituentFactory("MK3",[{constituent:constituents.M2,factor:1},{constituent:constituents.K1,factor:1}]);constituents.MN4=compoundConstituentFactory("MN4",[{constituent:constituents.M2,factor:1},{constituent:constituents.N2,factor:1}]);constituents.M4=compoundConstituentFactory("M4",[{constituent:constituents.M2,factor:2}]);constituents.MS4=compoundConstituentFactory("MS4",[{constituent:constituents.M2,factor:1},{constituent:constituents.S2,factor:1}]);constituents.S4=compoundConstituentFactory("S4",[{constituent:constituents.S2,factor:2}]);constituents.M6=compoundConstituentFactory("M6",[{constituent:constituents.M2,factor:3}]);constituents.S6=compoundConstituentFactory("S6",[{constituent:constituents.S2,factor:3}]);constituents.M8=compoundConstituentFactory("M8",[{constituent:constituents.M2,factor:4}]);const getDate=time=>{if(time instanceof Date){return time}if(typeof time==="number"){return new Date(time*1e3)}throw new Error("Invalid date format, should be a Date object, or timestamp")};const getTimeline=(start,end,seconds)=>{seconds=typeof seconds!=="undefined"?seconds:10*60;const items=[];const endTime=end.getTime()/1e3;let lastTime=start.getTime()/1e3;const startTime=lastTime;const hours=[];while(lastTime<=endTime){items.push(new Date(lastTime*1e3));hours.push((lastTime-startTime)/(60*60));lastTime+=seconds}return{items:items,hours:hours}};const harmonicsFactory=({harmonicConstituents,phaseKey,offset})=>{if(!Array.isArray(harmonicConstituents)){throw new Error("Harmonic constituents are not an array")}const constituents$1=[];harmonicConstituents.forEach(constituent=>{if(typeof constituent.name==="undefined"){throw new Error("Harmonic constituents must have a name property")}if(typeof constituents[constituent.name]!=="undefined"){constituent._model=constituents[constituent.name];constituent._phase=d2r*constituent[phaseKey];constituents$1.push(constituent)}});if(offset!==false){constituents$1.push({name:"Z0",_model:constituents.Z0,_phase:0,amplitude:offset})}let start=new Date;let end=new Date;const harmonics={};harmonics.setTimeSpan=(startTime,endTime)=>{start=getDate(startTime);end=getDate(endTime);if(start.getTime()>=end.getTime()){throw new Error("Start time must be before end time")}return harmonics};harmonics.prediction=options=>{options=typeof options!=="undefined"?options:{timeFidelity:10*60};return predictionFactory({timeline:getTimeline(start,end,options.timeFidelity),constituents:constituents$1,start:start})};return Object.freeze(harmonics)};const tidePredictionFactory=(constituents,options)=>{const harmonicsOptions={harmonicConstituents:constituents,phaseKey:"phase_GMT",offset:false};if(typeof options!=="undefined"){Object.keys(harmonicsOptions).forEach(key=>{if(typeof options[key]!=="undefined"){harmonicsOptions[key]=options[key]}})}const tidePrediction={getTimelinePrediction:({start,end})=>{return harmonicsFactory(harmonicsOptions).setTimeSpan(start,end).prediction().getTimelinePrediction()},getExtremesPrediction:({start,end,labels,offsets,timeFidelity})=>{return harmonicsFactory(harmonicsOptions).setTimeSpan(start,end).prediction({timeFidelity:timeFidelity}).getExtremesPrediction(labels,offsets)},getWaterLevelAtTime:({time})=>{const endDate=new Date(time.getTime()+10*60*1e3);return harmonicsFactory(harmonicsOptions).setTimeSpan(time,endDate).prediction().getTimelinePrediction()[0]}};return tidePrediction};return tidePredictionFactory});

package/package.json

@@ -1,45 +1,36 @@
 {
   "name": "@neaps/tide-predictor",
-  "version": "0.0.5",
+  "version": "0.1.1",
   "description": "Tide predictor",
-  "main": "",
   "repository": "https://github.com/neaps/tide-predictor",
   "author": "Kevin Miller <keveemiller@gmail.com>",
   "license": "MIT",
-  "module": "./lib/index.es6.js",
-  "browser": "./dist/tide-predictor.js",
-  "dependencies": {},
+  "type": "module",
+  "module": "./src/index.js",
+  "browser": "./dist/web/tide-predictor.js",
+  "commonjs": "./dist/commonjs/index.js",
+  "exports": {
+    ".": {
+      "import": "./src/index.js",
+      "require": "./dist/commonjs/index.js"
+    }
+  },
   "devDependencies": {
-    "@babel/core": "^7.6.2",
-    "@babel/preset-env": "^7.6.2",
-    "babel-jest": "^27.2.0",
-    "babel-preset-env": "^1.7.0",
-    "eslint": "^7.10.0",
-    "eslint-config-standard": "^14.1.0",
-    "eslint-plugin-import": "^2.18.2",
-    "eslint-plugin-jest": "^24.0.2",
-    "eslint-plugin-node": "^11.1.0",
-    "eslint-plugin-promise": "^4.2.1",
-    "eslint-plugin-standard": "^4.0.1",
-    "grunt": "^1.4.1",
-    "grunt-babel": "^8.0.0",
-    "grunt-contrib-clean": "^2.0.0",
-    "grunt-contrib-uglify": "^5.0.0",
-    "grunt-contrib-watch": "^1.1.0",
-    "grunt-coveralls": "^2.0.0",
-    "grunt-eslint": "^23.0.0",
-    "grunt-exec": "^3.0.0",
-    "jest": "^27.2.0",
-    "jest-junit": "^11.0.1",
-    "kind-of": "^6.0.3",
-    "node-fetch": "^2.6.0",
-    "rollup-plugin-babel": "^4.3.3",
-    "rollup-plugin-node-resolve": "^5.2.0"
+    "c8": "^7.11.3",
+    "eslint": "^8.18.0",
+    "mocha": "^10.0.0",
+    "npm-run-all": "^4.1.5",
+    "rollup-plugin-node-resolve": "^5.2.0",
+    "uglify-js": "^3.16.1"
   },
   "scripts": {
-    "test": "BABEL_ENV=test jest",
-    "coverage": "jest --coverage; grunt coveralls;",
+    "test": "mocha --recursive",
     "lint": "eslint ./src",
-    "ci": "yarn run lint; jest --ci --runInBand --reporters=default --reporters=jest-junit"
+    "coverage": "c8 --reporter=lcov mocha --recursive && codecov",
+    "ci": "run-s lint test coverage",
+    "build:rollup": "rollup --config rollup.config.js",
+    "build:uglify": "uglifyjs ./dist/web/tide-predictor.js -o ./dist/web/tide-predictor.min.js",
+    "build:commonDist": "echo \"{\\\"type\\\":\\\"commonjs\\\"}\" > ./dist/commonjs/package.json",
+    "build": "run-s build:rollup build:commonDist build:uglify"
   }
 }

package/rollup.config.js

@@ -9,13 +9,13 @@ export default {
       file: pkg.browser,
       format: 'umd'
     },
-    { file: pkg.module, format: 'es' }
+    { file: pkg.commonjs, format: 'commonjs', exports: 'default' }
   ],
   plugins: [
     resolve({
       mainFields: ['module', 'main'],
 
-      jail: '/src' // Default: '/'
+      jail: '/src'
     })
   ]
 }

package/src/astronomy/index.js

@@ -1,5 +1,5 @@
-import { d2r, r2d } from './constants'
-import coefficients from './coefficients'
+import { d2r, r2d } from './constants.js'
+import coefficients from './coefficients.js'
 
 // Evaluates a polynomial at argument
 const polynomial = (coefficients, argument) => {
@@ -139,7 +139,7 @@ const astro = (time) => {
     pp: coefficients.solarPerigee,
     90: [90.0],
     omega: coefficients.terrestrialObliquity,
-    i: coefficients.lunarInclination,
+    i: coefficients.lunarInclination
   }
 
   // Polynomials are in T, that is Julian Centuries; we want our speeds to be
@@ -148,7 +148,7 @@ const astro = (time) => {
   Object.keys(polynomials).forEach((name) => {
     result[name] = {
       value: modulus(polynomial(polynomials[name], T(time)), 360.0),
-      speed: derivativePolynomial(polynomials[name], T(time)) * dTdHour,
+      speed: derivativePolynomial(polynomials[name], T(time)) * dTdHour
     }
   })
 
@@ -160,7 +160,7 @@ const astro = (time) => {
     xi: _xi,
     nu: _nu,
     nup: _nup,
-    nupp: _nupp,
+    nupp: _nupp
   }
   Object.keys(functions).forEach((name) => {
     const functionCall = functions[name]
@@ -169,7 +169,7 @@ const astro = (time) => {
         functionCall(result.N.value, result.i.value, result.omega.value),
         360.0
       ),
-      speed: null,
+      speed: null
     }
   })
 
@@ -178,12 +178,12 @@ const astro = (time) => {
   // This is in line with convention.
   const hour = {
     value: (JD(time) - Math.floor(JD(time))) * 360.0,
-    speed: 15.0,
+    speed: 15.0
   }
 
   result['T+h-s'] = {
     value: hour.value + result.h.value - result.s.value,
-    speed: hour.speed + result.h.speed - result.s.speed,
+    speed: hour.speed + result.h.speed - result.s.speed
   }
 
   // It is convenient to calculate Schureman's P here since several node
@@ -191,7 +191,7 @@ const astro = (time) => {
   // (along with I, xi, nu etc) belong somewhere else.
   result.P = {
     value: result.p.value - (result.xi.value % 360.0),
-    speed: null,
+    speed: null
   }
 
   return result

package/src/constituents/compound-constituent.js

@@ -10,9 +10,8 @@ const compoundConstituentFactory = (name, members) => {
   })
 
   const compoundConstituent = {
-    name: name,
-
-    coefficients: coefficients,
+    name,
+    coefficients,
 
     speed: (astro) => {
       let speed = 0
@@ -46,7 +45,7 @@ const compoundConstituentFactory = (name, members) => {
       return f.reduce((previous, value) => {
         return previous * value
       })
-    },
+    }
   }
 
   return Object.freeze(compoundConstituent)

package/src/constituents/constituent.js

@@ -1,4 +1,4 @@
-import nodeCorrections from '../node-corrections/index'
+import nodeCorrections from '../node-corrections/index.js'
 
 /**
  * Computes the dot notation of two arrays
@@ -23,7 +23,7 @@ const astronimicDoodsonNumber = (astro) => {
     astro.p,
     astro.N,
     astro.pp,
-    astro['90'],
+    astro['90']
   ]
 }
 
@@ -49,9 +49,8 @@ const constituentFactory = (name, coefficients, u, f) => {
   }
 
   const constituent = {
-    name: name,
-
-    coefficients: coefficients,
+    name,
+    coefficients,
 
     value: (astro) => {
       return dotArray(coefficients, astronomicValues(astro))
@@ -63,7 +62,7 @@ const constituentFactory = (name, coefficients, u, f) => {
 
     u: typeof u !== 'undefined' ? u : nodeCorrections.uZero,
 
-    f: typeof f !== 'undefined' ? f : nodeCorrections.fUnity,
+    f: typeof f !== 'undefined' ? f : nodeCorrections.fUnity
   }
 
   return Object.freeze(constituent)