It looks like when ex12 clones the DM, it reuses the mesh data structure but copies the coordinate vector. Is there a reason it can't reuse that storage? Given that the PetscFEM routines accept the same position array for dependent and auxiliary fields, does it even make sense to have distinct coordinate arrays? Geoffrey