t50 function

Compute half-emptying time from nlme parameters