GET_RADIUS_NORMAL normal along the prime vertical section ellipsoidal radius of curvature

Contents

Inputs

Outputs

function N = get_radius_normal(lat, E)
arguments
  lat {mustBeReal}
  E (1,1) matmap3d.referenceEllipsoid
end

if isempty(E)
  E = matmap3d.wgs84Ellipsoid();
end

N = E.SemimajorAxis^2 ./ sqrt( E.SemimajorAxis^2 .* cos(lat).^2 + E.SemiminorAxis^2 .* sin(lat).^2 );
end