[CIG-SHORT] Output Material Subset

[Matthew Knepley]

On Sun, Dec 1, 2013 at 4:37 PM, Scott Henderson <sth54 at cornell.edu> wrote:

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

I assume you can do this with OutputSolnPoints, just like its done in the
GreenFunc examples.

   Matt


> 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
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/cig-short/attachments/20131201/ceac32c8/attachment.html>