Ocean loading for terrestrial gravity corrections

[craigmillernz]

Hi, I've just come across pyTMD and thanks for making a really comprehensive package! I've downloaded the FES2022 model and i'd like to calculate a timeseries of ocean loading (in mGal or nm/s^2 units) for a list of long/lat/elevation. This is to add to the solid earth tide correction that i already calculate (using the simple Longman formula) to correct terrestrial gravity data for solid earth tides. The examples show calculating displacements, currents, elevations etc, but is there a gravity ouput as well for ocean loading? Any assistance appreciated!

[tsutterley]

Hey @craigmillernz, unfortunately outputting in terms of mGals is not in the current functionality. I'd have to think about the best approach to getting it implemented (e.g. spherical harmonics versus Green's functions). Would gladly take some help getting it put together (it might be a more long-term task).

[craigmillernz]

Hi, Thanks for the prompt reply, agree more of a long-term task for me a the moment. Thanks again for an interesting package.