I am struggling with converting ASPECT mesh to OpenFOAM, I know the ASPECT vertices can be accessed using the following codes,

```
for (size_t i = 0; i < triangulation.n_vertices(); i++)
{
for (size_t j = 0; j < dim; j++)
{
std::cout<<triangulation.get_vertices()[i][j]<<" ";
}
std::cout<<"\n";
}
```

I noticed that this code can only access the initial mesh. For a deformed mesh (e.g. free surface case), I want to access the vertices coordinate (x,y, or z) but I still don’t manage it. I tried to add “mesh_displacements” to vertices(see code below), but the result seems not correct. **I would appreciate it very much if anyone can tell me how to access vertices coordinate of a deformed mesh!**

```
for (size_t i = 0; i < triangulation.n_vertices(); i++)
{
for (size_t j = 0; j < dim; j++)
{
std::cout<<triangulation.get_vertices()[i][j] + mesh_deformation->mesh_displacements[i*dim+j]<<" ";
}
std::cout<<"\n";
}
```

Variables `triangulation`

and `mesh_deformation`

could be accessed in `aspect/source/simulator/core.cc`