Straight field line angle not yet implemented

We should add the map S:(s, th, ze) -> (s, th*, ze) with th* = th + lambda to straight field line coordinates (s, th*, ze).

For this, one has to compute the inverse S^-1 by solving the equation th + lambda(s_i, th, ze_k) - th^*_j = 0 for th.