A semi-Lagrangian semi-implicit numerical method for models of the urine concentrating mechanism