Bingo! Many thanks! Bob bg@n4mrv:~/weewx_tonga_browse-main$ sudo python3 ./tonga_barometer.py distance to eruption 12056.6 km arrival at 1642257762 (2022-01-15 09:42:41) opposite pulse arrival at 1642307533 (2022-01-15 23:32:13) second time around pulse arrival at 1642382887 (2022-01-16 20:28:06) select datetime, barometer from archive where datetime > 1642240800 and datetime < 1642395600 and barometer is not null order by dateTime; query returned 171 data points Fitting between 1642240800.0 and 1642395600.0 with 10 knots
[image: Figure_1.png] On Sunday, January 23, 2022 at 11:56:22 AM UTC-5 bgra...@umw.edu wrote: > Hi morrowwm, > I have tried your latest github post and am getting the following error: > ``` > root@n4mrv:/home/bg/weewx_tonga_browse-main# ./tonga_barometer.py > > distance to eruption 12056.6 km > arrival at 1642257762 (2022-01-15 09:42:41) > opposite pulse arrival at 1642307533 (2022-01-15 23:32:13) > second time around pulse arrival at 1642382887 (2022-01-16 20:28:06) > select datetime, barometer from archive where datetime > 1642240800 and > datetime < 1642395600 and barometer is not null order by dateTime; > query returned 171 data points > > Traceback (most recent call last): > File "./tonga_barometer.py", line 96, in <module> > knots = np.linspace(np.min(tdata), np.max(tdata), > (stop_time-start_time)/(3600*smoothing_hours), endpoint=True) # spline > knot every N hours > File "<__array_function__ internals>", line 180, in linspace > File > "/usr/local/lib/python3.8/dist-packages/numpy/core/function_base.py", line > 120, in linspace > num = operator.index(num) > TypeError: 'float' object cannot be interpreted as an integer > ``` > The earlier version worked correctly and I don't think I have changed > anything. I am using the same copy of weewx.sdb and only changed my > lat/lon. > Much appreciated for your work on this software. > Cheers, > Bob > On Saturday, January 22, 2022 at 4:14:51 PM UTC-5 morr...@gmail.com wrote: > >> Also, I found a bug in calculating the smooth (spline fit) curve which is >> removed to highlight the eruption events. Added display of the points used >> to calculate the fitted curve. >> >> On Saturday, January 22, 2022 at 4:32:47 p.m. UTC-4 morr...@gmail.com >> wrote: >> >>> Forgot that, thanks for the reminder. Added now. >>> >>> On Saturday, January 22, 2022 at 4:22:51 p.m. UTC-4 storm...@gmail.com >>> wrote: >>> >>>> Nice work (Morrowwn thanks!). If you have any nulls in the data >>>> during the selected timeframe, change the query to read: (Cameron D >>>> thanks!) >>>> >>>> *query = "select datetime, barometer from archive where datetime > %.0f >>>> and datetime < %.0f and barometer is not null order by dateTime;" % >>>> (start_time, stop_time)* >>>> >>>> >>>> >>>> >>>> [image: hunga_tonga.png] >>>> >>>> If you want the bars closer to the rise and falls of the events, need >>>> to adjust the speed of sound slightly. >>>> On Saturday, January 22, 2022 at 2:37:09 PM UTC-5 pmcg...@gmail.com >>>> wrote: >>>> >>>>> WOW! This is some interesting data, here is my location in Northern >>>>> Utah. Looks like the reverse may be a little off in calculation but this >>>>> made it easy to find instead of hunting through the data. >>>>> >>>>> Great Work! >>>>> [image: hunga_tonga.png] >>>>> On Saturday, January 22, 2022 at 12:16:45 PM UTC-7 morr...@gmail.com >>>>> wrote: >>>>> >>>>>> Did an update to my script at: >>>>>> https://github.com/morrowwm/weewx_tonga_browse >>>>>> >>>>>> Mostly incorporating additions from Cameron (thanks!) >>>>>> >>>>>> All parameters you likely want to set are near the top of the script. >>>>>> >>>>>> 1. support for mysql (untested, I don't have a mysql based weewx >>>>>> installation) >>>>>> 2. calculation and display of all three initial arrival, arrival >>>>>> from opposite direction, and second arrival after one trip around the >>>>>> globe >>>>>> 3. highlighting of the three events with grey bars >>>>>> 4. cubic spline curve fit of raw data and use of that fit to >>>>>> remove the base pressure variations. The new smoothing_hours >>>>>> parameter can >>>>>> be used to control how close the fit is. >>>>>> 5. changed the display of data to highlight the extracted events >>>>>> and subdue the original data, other cosmetic changes >>>>>> >>>>>> I had a data outage about 26 hours after the event, we're quite a >>>>>> distance from Tonga, plus a major storm was clearing out in our area. So >>>>>> not obviously anything unusual happened. >>>>>> >>>>>> Some scientist might be interested in the arrival times of all this >>>>>> data, though. >>>>>> [image: hunga_tonga.png] >>>>>> >>>>>> -- You received this message because you are subscribed to the Google Groups "weewx-user" group. To unsubscribe from this group and stop receiving emails from it, send an email to weewx-user+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/weewx-user/bddac740-5c8c-4617-878a-bf7520018b5dn%40googlegroups.com.
