Losing subdomain after dolfin.refine()

You need to call the adapt-command as for instance done for facet markers here Adapt for mesh and MeshFunctions - #2 by dokken