[CIG-SHORT] (no subject)

[Brad Aagaard]

Xiao,

The stresses are computed at each of the quadrature points within a 
cell. This is true for all cell related data. By default the output 
includes the values at all of the quadrature points. You can use the 
CellFilterAvgMesh filter (see the tutorials in the manual and the 
examples) to have PyLith average the values over the quadrature points 
so you end up with the values (e.g., stress tensor) at just point per cell.

Regards,
Brad


On 12/04/2013 12:45 PM, Ma, Xiao wrote:
> Hi Brad,
> I am able to use the h5py , put the stress into a variable called stress, however the shape of this stress variable is [500,10000,12], I know the first tow 500 is the total times steps, and 10000 is the total element numbers, however why there is 12 components for the stress  tensor?
> Best,
> Xiao
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>