1function [r, z] = cymbal_point(basis, sfrac)
2% CYMBAL_POINT Meridian coordinates (r, z) at an arclength fraction.
3nn = length(basis.snode);
4s = min(0.999999, max(0, sfrac)) * basis.S;
5idx = sum(basis.snode <= s);
6idx = min(max(idx, 1), nn - 1);
7xi = (s - basis.snode(idx)) / (basis.snode(idx + 1) - basis.snode(idx));
8r = basis.rnode(idx) + xi * (basis.rnode(idx + 1) - basis.rnode(idx));
9z = basis.znode(idx) + xi * (basis.znode(idx + 1) - basis.znode(idx));
10end