This allows us to return a GlobalFriction and thus hide some information

* Move coefficient computation into RuinaNonlinearity
* Allow some friction data to vary in space

As a consequence, the nonlinearity is only assembled once

Also fix the order a-mu-V0-stress-b-state-L everywhere

The code for Laursen's model is long gone