fn node_geometry(surface: &Surface, i: usize) -> [[f64; 3]; 4]Expand description
The point of node i in the x-z plane, its unnormalized surface normal
r rhat - r' theta-hat (r' is the slope dr/dtheta) and the unit vectors rhat
and -theta-hat, along which the normal changes with the radius and the slope.