(49g 50g) Theoretical Earth gravity g = g(latitude, height), WGS84, GRS80/67

10172021, 04:22 PM
(This post was last modified: 10192021 07:42 PM by Gil.)
Post: #12




RE: HP4950G:Theorical Earth gravity g = g(latitude, height), WGS84, GRS80/67
New version 12b
Minor changes ("nice to have"), above all relative to gE (g at Equator), gE.STR, gP (g at Pole), gP.STR, k and k.STR values in GRS80 and WGS84 directory. I'm particularly indebted to Charles Karney for his patience and insight regarding some of the formulae and results in relation to the above and to the topic in general. Version 11f With latitude and height h, calculates automatically g2C67 g8C67 gS67 g2C80 g8C80 gS80 gS84 gSFR g±SFR gClFR. g 2C67: g2C67: gravity result g2C67: 2nd order appro g2C67: Chebyshev g2C67: GRS67 g8C80: 8th order appro g8C80: GRS80 ±S Smodif gS84 : Somigliana gS84 : WGS84 gµ:gMean gFR(EE):gyour FREE ellips Cl : closed form for g calculation, without having to introduce gE / GP (self calculated) As previously, you can calculate your own FREE ellipsoid, changing the values of GM w f a in FREE Dir (no need to look for g at Equator or at Pole). Two methods automatically calculated : closed form or modified Somigliana's. J2.CAL: is to get a close value of J2 in GRS67 WGS84 and FREE. It is not used, however, to further calculations (seemingly to many rounding intermediate errors). F.CAL: is to get a close value of f in WGS80. It is not used, however, to further calculations (seemingly to many rounding intermediate errors). Only f and f.STR values (more accurate) are used instead for further calculations. Many selfexplaining NOTES at the end of each directory. 

« Next Oldest  Next Newest »

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