Add new function TWI to compute Beven and Kirkby 1979 topographic wetness index as ln(sca/slp).