from dolfin import *

coarse_mesh = UnitSquareMesh(3, 3)
fine_mesh = refine(refine(refine(refine(coarse_mesh))))

P2 = FunctionSpace(coarse_mesh, "CG", 2)
P1 = FunctionSpace(fine_mesh, "CG", 1)

f = interpolate(Expression("exp(-(pow(x[0] - 0.5, 2) + pow(x[1] - 0.5, 2)) / 0.05)"), P2)
g = interpolate(f, P1)

plot(f, interactive=True)
plot(g, title="what one wants", interactive=True)