WBL with speed relaxation

`Speed relaxation' is the simplest form of an acceleration relation:
dD/dt = [Dasym(K) - D] / tau(K)
type name
wbl_relax
parameters
As for WBL, with in addition
tau(K)
- another object of type Realfunction1d