[CIG-SHORT] Output Material Subset

[Scott Henderson]

In Pylith is it possible to output stresses resolved onto an arbitrary surface (not a defined fault) for a 3D model? 

I’m imagining something like OutputSolnSubset for materials where you specify a Cubit nodeset based on a material surface. For example, this would be useful for plotting stress contours in a profile through a 3D domain:

#---------------
[pylithapp.problem.materials.elastic]
output =  pylith.meshio.OutputSolnSubset

[pylithapp.problem.materials.elastic.output]
label = face_y
cell_filter = pylith.meshio.CellFilterAvgSubMesh
cell_data_fields = [total-strain, stress] 
writer = pylith.meshio.DataWriterHDF5SubMesh
writer.filename = output/step01/NorthSouthMaterialCrossSection.h5
#———————

This is similar to the question of getting stresses at specified points (http://www.geodynamics.org/pipermail/cig-short/2012-October/001165.html). So I know it’s possible to do this as post-processing in Python or Matlab. Maybe there is also an easy way to export ’slices' of 3D data from Paraview too? Any suggestions would be appreciated,

Thanks,
Scott