• About

Hugh's geophysical blog

~ Occasional musings on meteorology and other bits of Earth science

Monthly Archives: August 2013

Gravity tides

16 Friday Aug 2013

Posted by Hugh in Geophysics

≈ 2 Comments

As a background hobby activity over the summer I have been setting up one of our Lacoste-Romberg gravity meters so that I can show a continuous recording of its output on a web page.

gravplot

The picture shows the first few days of data; if the system is still working then the up-to-date plot is here. In addition to the gravity tides (which sort of agree with the red theoretical curve) the plot shows two earthquakes. Today’s (Friday’s) was in New Zealand; Tuesday’s was in Colombia.

As usual with things involving actual equipment, it all proved to be more complicated to set up than I would have expected. The gravity meter produces a voltage output; the same signal that goes to its beam position galvanometer. I therefore needed a data logging device: The DATAQ DI-145 looked rugged and simple so I got one of those.  Learning to get the data out of it onto a computer running Linux was simple. Understanding why it worked on some computers and not others involved the acquisition of some rather arcane knowledge involving serial port settings and the stty command.

Once I had the data flowing, I then wanted to plot the theoretical tide for comparison. That, too, turned into a bit of a time sink. One of the ways to do the calculation is detailed in Longman (1959); the paper provides all the formulas you need. And there are rather a lot of them. I looked for extant code on the web and found the following:

  • A web form that applies Longman’s formulas
  • A BASIC (yuk!) program that does likewise
  • A C program that purports to do likewise but which seems to get the wrong answer
  • A very sophisticated Windows program, which is not easy to call from a script on a Linux machine

As none of these were what I wanted, I took a deep breath and re-coded the Longman formulas in the very wonderful R, as that is what I used for the rest of the data processing and plotting. Tracking down the last major bug in my code took ages. I was reduced to translating the BASIC program into R as a cross-check. When that worked, I knew that my own  code must have a simple bug; Murphy’s law ensured that it was in almost the last line of the function. It does now seem to be working but it needs a bit of tidying up before I post it on the web.  And I am also itching to annotate the plot with earthquake reports from the USGS.

← Older posts

Subscribe

  • Entries (RSS)
  • Comments (RSS)

Archives

  • December 2019
  • November 2019
  • September 2019
  • June 2019
  • May 2019
  • November 2018
  • October 2018
  • September 2018
  • August 2018
  • July 2018
  • June 2018
  • May 2018
  • March 2018
  • February 2018
  • October 2017
  • September 2017
  • August 2017
  • July 2017
  • June 2017
  • May 2017
  • April 2017
  • March 2017
  • January 2017
  • December 2016
  • November 2016
  • October 2016
  • June 2016
  • May 2016
  • April 2016
  • March 2016
  • February 2016
  • December 2015
  • October 2015
  • September 2015
  • August 2015
  • July 2015
  • June 2015
  • May 2015
  • April 2015
  • March 2015
  • February 2015
  • January 2015
  • December 2014
  • October 2014
  • September 2014
  • August 2014
  • July 2014
  • June 2014
  • May 2014
  • April 2014
  • March 2014
  • February 2014
  • January 2014
  • December 2013
  • November 2013
  • October 2013
  • September 2013
  • August 2013
  • July 2013
  • June 2013
  • May 2013

Categories

  • Classic Lump
  • Garden
  • geology
  • Geophysics
  • Meteorology
  • Uncategorized

Meta

  • Register
  • Log in

Calendar of posts

August 2013
M T W T F S S
« Jul   Sep »
 1234
567891011
12131415161718
19202122232425
262728293031  

Recent posts

  • Above the silvery Tay December 2, 2019
  • Sigh. Again. November 26, 2019
  • Sundial — continued September 7, 2019
  • Sundial June 16, 2019
  • Columbines May 28, 2019

Blog at WordPress.com.

Cancel
Privacy & Cookies: This site uses cookies. By continuing to use this website, you agree to their use.
To find out more, including how to control cookies, see here: Cookie Policy