Post Reply 
Gravity g WGS84 84 at pole
09-28-2021, 08:53 AM
Post: #1
Gravity g WGS84 84 at pole
Sure that here, in that HP calculator forum, is not the normal place to place my question.

But that interrogation occurred now when developing my gravity program for the HP49-HP50.

And as I know that the level of some HP user is quite high, counting a lot of maths people or ingenieurs, I dare and submit my problem, hoping for your indulgence... and help.

In the Wikipedia "Theoretical Gravity" article,
choose the Chinese language.

There is to be found a formula j.e
(indeed j.p) for the gravity g at the pole.

j.pole =
'GM/(a*a)*(1+m+3/7*e'²*m)' equation 1

In WGS84
a=exactly 6 378 137
with w: = exactly .00007292115
b=a-a*f and f=exactly 298.257223563
so that b=about 6356752. 31424517949756396659963365515679817131108549733884857165128320852208683510093940​57094507132374770633
And GM=exactly 3.986004418E14
so that m=about 0.003449786506840845307027496505078670004841648268835534929013642504025303575739​68704757009899238851678016698284958722982651012885967016569924935793184562396036​96294749064174268559478550984385788003910837612

so that e'²=about 0.006739496742276434954782158956759376656122065858844009981905028258322797320891​4914189860809667349634

So that (1+m+3/7*e'²*m)=

And finally g. Pole (without the dot/comma) should be exactly equal to, according to equation 1:

Or j. Pole, placing back the dot/comma at the right place=
9.832185104404435416623584524111001588952012585446633794280230591536194799286886304178567025​98679431255291172411572286558495901685221296632175331598046360918175599536950480​56826641872205769685950655291862172166302389558892330844759850013860675472423221​62331131594199840054395148802370298815134288969023537376906669191507694470181419​8886293763195146108008175411982135184646619151544169766231068819297264703012453.​ (Full digit calculated result)

But all the formulae given for the gravity at the pole write that it is exactly? =
9.83218 49378 (General given 11 digits result)

My questions
Does a term, or several, miss in the equation 1?
Or, on the the contrary, Full digit calculated result is correct and General given 11 digits result was I then incorrectly rounded (with the 7th digit already false)?

As a curious layman, I would appreciate your precious insights.

Thanks for your help.

Gil Campart
Find all posts by this user
Quote this message in a reply
09-28-2021, 10:48 AM
Post: #2
RE: Gravity g WGS84 84 at pole
(09-28-2021 08:53 AM)Gil Wrote:  e'²='(a*a-b*b)/(b*b)'

e²= 1 - (b/a)^2 = (a*a-b*b)/(a*a)

The typo may still not make up the difference ...
Calculated gravity is derived from assumed standards, thus the need to keep saying WGS84
Find all posts by this user
Quote this message in a reply
09-28-2021, 11:16 AM (This post was last modified: 09-28-2021 11:17 AM by Gil.)
Post: #3
RE: Gravity g WGS84 84 at pole
In fact, using the reference of
H. Moritz Geodesic Reference System 80,

we get quite closer, better results.

The differences are certainly due to rounding errors by tan^-1 and sqrt functions.


Find all posts by this user
Quote this message in a reply
09-28-2021, 03:21 PM
Post: #4
RE: Gravity g WGS84 84 at pole
RE: HP49-50G : —>g gravity calculation = g(latitude, height) with WGS84
To check,
the paper variables of Geodetic Reference System 1980, by H.Moritz, should be taken, instead of the simplifications given in Wikipedia "Theorical Gravity", Chinese page.


j.e =
9.78032533482, from the above formulae
9.7803253359 official
—> almost equal value

j.p =
9.83218494001, from the above formulae
9.8321849378, official
—> almost equal value!

ƒé² = sqrt (é²) =e'


q0´ =
'3*(1+1/é²)*(1-1/ƒé²*ATAN(ƒé²))-1' .00268804118

q0 =

é² =

GM =

m =

w =
a 6378137

a =

b =

The differences are now quite small, due to the roundings of the calculator.

Find all posts by this user
Quote this message in a reply
Post Reply 

User(s) browsing this thread: 1 Guest(s)