Add keyword to allow specific r to be calculated
Add keyword to allow specific r to be calculated

Before the r was always given by r_sing + offset.
This is a lazy addition to my function it would be better to rewrite the
function for the general case and then a specific function calling the
general fro the old case.

jensv/fluxtubestability,jensv/fluxtubestability