Stick to numpy functions when possible
Output runoff rate and instantaneous mass balance
Get rid of division by zero warnings
Add option -l to list available output variables