Implement substrate decay rate

The code currently does not work when KDEGS is set to a non-zero value.