diff options
author | Mattias Andrée <m@maandree.se> | 2024-10-22 22:22:41 +0200 |
---|---|---|
committer | Mattias Andrée <m@maandree.se> | 2024-10-22 22:22:41 +0200 |
commit | bb5de3aa2ee118df78f0347cffd4e58f846dc1fb (patch) | |
tree | b313a20b79d67f35a5e4e05bc9106b22dc17f451 /libtellurian_transverse_radius_radians.c | |
parent | ... (diff) | |
download | libtellurian-bb5de3aa2ee118df78f0347cffd4e58f846dc1fb.tar.gz libtellurian-bb5de3aa2ee118df78f0347cffd4e58f846dc1fb.tar.bz2 libtellurian-bb5de3aa2ee118df78f0347cffd4e58f846dc1fb.tar.xz |
Signed-off-by: Mattias Andrée <m@maandree.se>
Diffstat (limited to 'libtellurian_transverse_radius_radians.c')
-rw-r--r-- | libtellurian_transverse_radius_radians.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/libtellurian_transverse_radius_radians.c b/libtellurian_transverse_radius_radians.c index c4f80aa..3e4950e 100644 --- a/libtellurian_transverse_radius_radians.c +++ b/libtellurian_transverse_radius_radians.c @@ -9,7 +9,7 @@ libtellurian_transverse_radius_radians(double latitude) double a = LIBTELLURIAN_EQUATORIAL_RADIUS; double b = LIBTELLURIAN_POLAR_RADIUS; double neg_e2 = fma(b / a, b / a, -1.0); - double s = sin(latitude); + double s = sin(geodetic(latitude)); return a * pow(fma(neg_e2, s * s, 1.0), -0.5); } |