import numpy as np
# Constantes
a = 6378.137
j2 = 0.0010826265
c1 = 0.0034134
f = 1/298.257222
phi = -0.9784
# Radio del geoide
r_0 = a*(1-(3/2*j2+0.5*c1)*(np.sin(phi)**2))
# Radio del elipsoide
r_c = a*(1 + (((2*f - f**2)/(1-f)**2)*(np.sin(phi))**2))**(-0.5)
# comparacion
r_0, r_c
# Transformar las unidades a metros
com = r_0-r_c
com
metros = com*1000
com, metros