Movemesh23 and move solution for visualisation

[fb77 (François Bouchut)]

There is a method that works here

Transfer scalar field to deformed mesh with movemesh2d

Dear Angel Daniel, Indeed the following lines make it work for P2 or any finite element space: int[int] renum(Ch.ndof); for (int k=0;k<Th2.nt;++k){ for (int j=0;j<Ch.ndofK;++j){ renum(Ch3(k,j))=Ch(k,j); } } temp3d[]=Temp[](renum); full code: load "msh3" real r=1.; real R=2.; real x0 = pi/2.; real x1 = 3.*pi/2.; real y0 = pi; real y1 = 2.*pi; int n = 5; int m = 5; mesh Th2 = square(n, m, [x0+(x1-x0)*x, y0+(y1-y0)*y]); fespace Ch(Th2,P2); Ch Temp; //calculations of Temp follow Temp=x; plot…

but it is related to the property that the mesh nodes are purely transported by th transformation.
In your case you have “removeduplicate=1” that breaks this property (the two meshes do not have the same number of nodes).
Maybe you can do it if you set “removeduplicate=0”