Boundary condition zero

This branch changes the behaviour of the nonlinear solver when applying boundary conditions to the RHS. Instead of setting the boundary condition nodes to bc-value, they are set to zero. This is consistent with the way we treat the LHS and with the documentation of boundary conditions.