[multidomain,Beavers-Joseph] Move evalBeaversJoseph to LOP which prevents a reimplementation in the non-isothermal coupling local residual.

Move beaversJospehCoeff to SpatialParamters, where it belongs.
Changes in reference solution was triggered, comments see below.

2cstokes2p2c:
Checked storage.out (e.g. the evaporation rate):
  - changes were lass than 1 per mile.
Checked vertical velocity profiles in the middle of the domain:
  - v_0: no visible difference
  - v_1: main deviation (5%) in middle, where the profile is oscillating anyway and NOT at the interface
Checked horizontal velocity profiles:
  - v_0: small differences towards the outflow boundary
  - v_1: difference during the profile development

2cnizeroeq2p2cni (freeflow):
Checked vertical velocity profiles in the middle of the domain:
  - v_0: no visible difference
  - v_1: no visible difference
Checked horizontal velocity profiles:
  - v_0: no visible difference
  - v_1: small differences towards the outflow boundary

2cnizeroeq2p2cni (porousmedium):
Checked vertical velocity profiles in the middle of the domain:
  - v_0: no visible difference
  - v_1: very small deviation towards the interface (3%)