On 05/07/18 18:03, Amireh Mousavi wrote:

Dear all, in Firedrake how can I plot the solution as a 3D surface which has volume? thank you in advance.

For anything more than very simple plotting (and 3D falls into this category), we recommend that you use Paraview for visualisation. You can save data to a pvd file, readable by Paraview with: output = File("output.pvd") output.write(function) Cheers, Lawrence