[CIG-SHORT] embedded fault plane

Rebecca Bendick bendick at mso.umt.edu
Fri Jan 15 13:12:50 PST 2010


Hi, everyone.  This may be a very rookie question; I'm just trying to  
learn how to use PyLith.  Are there any examples using a dislocation  
on a fault plane embedded in a larger volume?  I've had no trouble   
reproducing results for a throughgoing fault (cutting though the whole  
model volume) using the examples provided in the distribution and in  
the latest manual.  However, now I've used the demos from LaGRiT, in  
particular variations on the example at http://meshing.lanl.gov/proj/crustal_dynamics_CFEM/insert_fault/gallery.html 
  to create a fully embedded dislocation, and I can't get through the  
PyLith run.  Here's the screen dump, but (given my newbie-ness) maybe  
an already-worked example would be more helpful than debugging because  
I might have made many nested mistakes.

  -- meshiolagrit(info)
  -- Reading finite-element mesh
  >> meshio/PsetFileBinary.cc:70:<unknown>
  -- psetfile(info)
  -- Reading binary Pset file 'tet4_1000m_binary.pset'.
  >> meshio/PsetFileBinary.cc:86:<unknown>
  -- psetfile(info)
  -- Reading 7 point sets from file.
  >> meshio/PsetFileBinary.cc:203:<unknown>
  -- psetfile(info)
  -- Reading point set 'fault' with 0 points.
  >> meshio/PsetFileBinary.cc:215:<unknown>
  -- psetfile(info)
  -- Done.
  >> meshio/PsetFileBinary.cc:203:<unknown>
  -- psetfile(info)
  -- Reading point set 'boundary_xm' with 287 points.
  >> meshio/PsetFileBinary.cc:215:<unknown>
  -- psetfile(info)
  -- Done.
  >> meshio/PsetFileBinary.cc:203:<unknown>
  -- psetfile(info)
  -- Reading point set 'boundary_xp' with 287 points.
  >> meshio/PsetFileBinary.cc:215:<unknown>
  -- psetfile(info)
  -- Done.
  >> meshio/PsetFileBinary.cc:203:<unknown>
  -- psetfile(info)
  -- Reading point set 'boundary_ym' with 287 points.
  >> meshio/PsetFileBinary.cc:215:<unknown>
  -- psetfile(info)
  -- Done.
  >> meshio/PsetFileBinary.cc:203:<unknown>
  -- psetfile(info)
  -- Reading point set 'boundary_yp' with 287 points.
  >> meshio/PsetFileBinary.cc:215:<unknown>
  -- psetfile(info)
  -- Done.
  >> meshio/PsetFileBinary.cc:203:<unknown>
  -- psetfile(info)
  -- Reading point set 'boundary_zm' with 49 points.
  >> meshio/PsetFileBinary.cc:215:<unknown>
  -- psetfile(info)
  -- Done.
  >> meshio/PsetFileBinary.cc:203:<unknown>
  -- psetfile(info)
  -- Reading point set 'boundary_zp' with 49 points.
  >> meshio/PsetFileBinary.cc:215:<unknown>
  -- psetfile(info)
  -- Done.
faults/TopologyOps.cc:330: failed assertion `numFaultCorners ==  
faceSize'
mpirun noticed that job rank 0 with PID 7750 on node kham.gs.umt.edu  
exited on signal 6 (Abort trap).
/Users/rbk/pylith/pylith-1.4.2-darwin-powerpc/bin/nemesis: mpirun:  
exit 134
/Users/rbk/pylith/pylith-1.4.2-darwin-powerpc/bin/pylith: /Users/rbk/ 
pylith/pylith-1.4.2-darwin-powerpc/bin/nemesis: exit 1



Thanks in advance,
Becks
---------------------------------------------------------------
Rebecca Bendick
SC 331
Department of Geosciences
University of Montana
Missoula, MT 59812
(406) 243-5774

bendick at mso.umt.edu





More information about the CIG-SHORT mailing list