On 11/29/2017 06:13 AM, firstname.lastname@example.org wrote:
I think I found the key, I can call cmesh.get_incident to get the conn of a specific region of facets and I assume that the quad nodes are constructed based on the order of this conn, right?
The surface connectivity is directly available in the SurfaceData of the field
sd = virtual.field.surface_data[self.region.name]
Note that the connectivity from CMesh is not the same as the field connectivity
It should be also possible to construct a surface field for the region, define a corresponding variable, set its DOFs to the nodal pressure values, and then set the material parameter of dw_surface_ltr by a function, where you just evaluate the variable in the required coordinates. Let me know if you need more advice (send the example script).
Otherwise for your purpose, maybe you can use dw_point_load term?