2D/1D Flory Huggins works in principle. Currently output only via paraview, need to figure out plotting from fenics.