[CIG-SHORT] [EXTERNAL] Re: defining fault_edge nodes

[Brad Aagaard]

Niloufar,

The cells (numCells, numCorners) are the indices (zero based) of the 
vertices.

So cells[iCell, iCorner) is the index (zero based) into the vertices 
array for the iCorner vertex of iCell cell. So the coordinates of that 
vertex is vertices[cells[iCell, iCorner, :] (using Python syntax). For 
one-based indexing, use vertices[1+cells[iCell,iCorner],:].

Regards,
Brad


On 6/29/18 3:29 PM, Niloufar Abolfathian wrote:
> Hi,
> 
> I had another question regarding reading the output files. I am using 
> matlab for reading the h5 files and further plottings.
> For reading stress and strain related to each material I do have the 
> stress/strain regarding each gird cell plus the topology information 
> (only number of vertices not the geometrical values). I could not figure 
> out how to relate the topology cells with the geometrical vertices.
> 
> Thanks,
> Niloufar
> 
> 
> On Thu, Jun 21, 2018 at 7:14 PM, Brad Aagaard <baagaard at usgs.gov 
> <mailto:baagaard at usgs.gov>> wrote:
> 
>     On 6/21/18 5:50 PM, Niloufar Abolfathian wrote:
> 
>         Thanks, Brad. By my second question, I mean in cases where the
>         fault is extended to the external boundaries, and we want to set
>         boundary conditions on these surfaces(external boundaries),
>         fault nodes should be all excluded. Is this statement right?
> 
> 
>     Yes, that is correct.
> 
> 
> 
> 
>         On Thu, Jun 21, 2018 at 5:40 PM, Brad Aagaard <baagaard at usgs.gov
>         <mailto:baagaard at usgs.gov> <mailto:baagaard at usgs.gov
>         <mailto:baagaard at usgs.gov>>> wrote:
> 
>              On 6/21/18 12:38 PM, Niloufar Abolfathian wrote:
> 
>                  Hi,
> 
>                  I am modeling a 3D vertical fault in an upper crust,
>         lying on a
>                  lower crust. It is as your 3d/hex example instead of
>         that the
>                  fault is only located in the upper crust.  I set up the
>         nodes
>                  related to the fault and fault_edge. From my meshing
>         and node
>                  selection the fault should only be located in the upper
>         crust,
>                  but after running the model the fault is extended below the
>                  upper crust. I guess it may be due to the selection of
>         fault and
>                  fault_edge nodes. Hope you can help me with this issue.
> 
>                  Question: Should the fault_edge nodes be also included
>         in fault
>                  nodes? (I suppose the answer is yes, due to the figure
>         6.4 in
>                  the manual.)
>                  I tried excluding fault_edge nodes from the fault but
>         then the
>                  model will not be converged.
> 
> 
>              Yes. The fault edge nodes are included in the fault nodes.
> 
>              The edge nodes are any buried edges of the fault. That is,
>         any nodes
>              on edges that are not on the external domain boundary should be
>              included in the fault edges.
> 
>                  I have another question. In the example 3d/hex model,
>         if we want
>                  to define a boundary condition on y_neg/y_pos , the
>         fault nodes
>                  be excluded az z_neg, right?
> 
> 
>              I don't understand this question.
> 
>              Regards,
>              Brad
>              _______________________________________________
>              CIG-SHORT mailing list
>         CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org>
>         <mailto:CIG-SHORT at geodynamics.org
>         <mailto:CIG-SHORT at geodynamics.org>>
>         http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>         <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>
>             
>         <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>         <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>>
> 
> 
>     _______________________________________________
>     CIG-SHORT mailing list
>     CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org>
>     http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>     <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>
> 
>