Dear vuong,
If the plane z=z0 is part of the mesh as a boundary you can integrate on the boundary using a proper region flag:
Code:
defined=mesh.Boundaries("myplane")
.
Otherwise there is no core functionality in NGSolve. However, you can use ngs-xfem, which is a small add-on module for things like that ( see
ngsxfemhttps://github.com/ngsxfem/ngsxfem
). There you have to provide a level set function where the zero level corresponds to your plane, i.e. this would look like
Code:
Integrate(levelset_domain = { "levelset" : z, "domain_type" : IF }, cf = your_integrand)
.
Best,
Christoph