question
The calculated geoid magnitude is incorrect and differs significantly from the actual result. The document says that no additional scaling is needed, and I’m not sure what’s wrong with that.
Strangely, the deflection of vertical based on geoid is correct. geoid’s magnitude is not correct but its difference (deflection of vertical) is correct, I guess it could be removing the trend term, but using N+l is still not correct.
Hi,
I don’t think you can do this comparison. The EGM is a global grid and you cannot expect to capture all the wavelengths when computing the the geoid from a small region as your faa grid. And extending that grids area is limited to the extent that a FFT is valid when applied to a geographical grid. I also don’t think there is an error in grdfft. As the main page says:
This is equivalent to divide by kr in the frequency domain (kr is radial wave number).
and that has been tested during years (ofc nothing is bug-free but …)