Debugging Numerical ice sheet modelling using FEniCS

To me it seems like your StokesSolver does not return anything.

It should return the solution w.

Please follow the guidelines about posting a minimal code example: