[cig-commits] r11654 - in short/3D/PyLith/trunk/examples/3d: . bar_hex8 bar_tet4

brad at geodynamics.org brad at geodynamics.org
Mon Mar 31 00:04:29 PDT 2008


Author: brad
Date: 2008-03-31 00:04:29 -0700 (Mon, 31 Mar 2008)
New Revision: 11654

Added:
   short/3D/PyLith/trunk/examples/3d/bar_tet4/
   short/3D/PyLith/trunk/examples/3d/bar_tet4/README
   short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.gmv
   short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.pset
   short/3D/PyLith/trunk/examples/3d/bar_tet4/matprops.spatialdb
   short/3D/PyLith/trunk/examples/3d/bar_tet4/mesh_tet4_200m.lagrit
   short/3D/PyLith/trunk/examples/3d/bar_tet4/output/
   short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave.cfg
   short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_slip.spatialdb
   short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliprate.spatialdb
   short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliptime.spatialdb
Modified:
   short/3D/PyLith/trunk/examples/3d/bar_hex8/README
   short/3D/PyLith/trunk/examples/3d/bar_hex8/shearwave.cfg
Log:
Added bar_tet4 example (same as bar_hex8 but with tet4 cells).

Modified: short/3D/PyLith/trunk/examples/3d/bar_hex8/README
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_hex8/README	2008-03-31 06:39:25 UTC (rev 11653)
+++ short/3D/PyLith/trunk/examples/3d/bar_hex8/README	2008-03-31 07:04:29 UTC (rev 11654)
@@ -4,7 +4,8 @@
 vertices are fixed in the longitudinal and other shear direction)
 generates a shear wave that propagates down the bar. Absorbing
 boundary conditions on each end perfectly absorb the shear wave.
-The bar is discretized with hexahedral cells.
+The bar is discretized with hexahedral cells (see examples/3d/bar_tet4)
+for the same problem discretized with tetrahedral cells.
 
 Files associated with the example problem:
 
@@ -18,3 +19,7 @@
     peak slip rate
   - shearwave_sliptime.spatialdb: Spatial database with values for the
     slip initiation time
+
+To run the problem, enter
+
+  pylith shearwave.cfg

Modified: short/3D/PyLith/trunk/examples/3d/bar_hex8/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_hex8/shearwave.cfg	2008-03-31 06:39:25 UTC (rev 11653)
+++ short/3D/PyLith/trunk/examples/3d/bar_hex8/shearwave.cfg	2008-03-31 07:04:29 UTC (rev 11654)
@@ -48,7 +48,7 @@
 formulation = pylith.problems.Explicit
 
 # Reuse solution from previous time step as initial guess for next time step
-#formulation.solver.initial_guess_zero = False
+formulation.solver.initial_guess_zero = False
 
 # Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
 bc = [x_pos,x_neg,y_pos]

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/README
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/README	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/README	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,25 @@
+This directory contains a wave propagation example problem. We create
+a shear dislocation in a rectangular, elastic bar 8 km long with a
+square cross-section 400 m on a side. The shear dislocation (the
+vertices are fixed in the longitudinal and other shear direction)
+generates a shear wave that propagates down the bar. Absorbing
+boundary conditions on each end perfectly absorb the shear wave.
+The bar is discretized with tetrahedral cells (see examples/3d/bar_hex8)
+for the same problem discretized with hexahedral cells.
+
+Files associated with the example problem:
+
+  - mesh_tet4_200m.lagrit: LaGriT script file for meshing the problem
+  - shearwave.cfg: Parameter settings for the shear wave dislocation problem
+  - matprops.spatialdb: Material database of physical properties
+  - bar_tet4_200m.gmv: GMV file with mesh generated by LaGriT
+  - bar_tet4_200m.pset: Pset file with point sets for mesh generated by LaGriT
+  - shearwave_slip.spatialdb: Spatial database with values for final slip
+  - shearwave_sliprate.spatialdb: Spatial database with values for the
+    peak slip rate
+  - shearwave_sliptime.spatialdb: Spatial database with values for the
+    slip initiation time
+
+To run the problem, enter
+
+  pylith shearwave.cfg

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.gmv
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.gmv	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.gmv	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,1300 @@
+gmvinput ascii
+codename LaGriT
+simdate ô/14/07
+nodes          369
+  0.00000E+000  0.00000E+000  0.00000E+000  2.00000E+002  2.00000E+002  2.00000E+002  4.00000E+002  4.00000E+002  4.00000E+002  6.00000E+002
+  6.00000E+002  6.00000E+002  8.00000E+002  8.00000E+002  8.00000E+002  1.00000E+003  1.00000E+003  1.00000E+003  1.20000E+003  1.20000E+003
+  1.20000E+003  1.40000E+003  1.40000E+003  1.40000E+003  1.60000E+003  1.60000E+003  1.60000E+003  1.80000E+003  1.80000E+003  1.80000E+003
+  2.00000E+003  2.00000E+003  2.00000E+003  2.20000E+003  2.20000E+003  2.20000E+003  2.40000E+003  2.40000E+003  2.40000E+003  2.60000E+003
+  2.60000E+003  2.60000E+003  2.80000E+003  2.80000E+003  2.80000E+003  3.00000E+003  3.00000E+003  3.00000E+003  3.20000E+003  3.20000E+003
+  3.20000E+003  3.40000E+003  3.40000E+003  3.40000E+003  3.60000E+003  3.60000E+003  3.60000E+003  3.80000E+003  3.80000E+003  3.80000E+003
+  4.00000E+003  4.00000E+003  4.00000E+003  0.00000E+000  0.00000E+000  0.00000E+000  2.00000E+002  2.00000E+002  2.00000E+002  4.00000E+002
+  4.00000E+002  4.00000E+002  6.00000E+002  6.00000E+002  6.00000E+002  8.00000E+002  8.00000E+002  8.00000E+002  1.00000E+003  1.00000E+003
+  1.00000E+003  1.20000E+003  1.20000E+003  1.20000E+003  1.40000E+003  1.40000E+003  1.40000E+003  1.60000E+003  1.60000E+003  1.60000E+003
+  1.80000E+003  1.80000E+003  1.80000E+003  2.00000E+003  2.00000E+003  2.00000E+003  2.20000E+003  2.20000E+003  2.20000E+003  2.40000E+003
+  2.40000E+003  2.40000E+003  2.60000E+003  2.60000E+003  2.60000E+003  2.80000E+003  2.80000E+003  2.80000E+003  3.00000E+003  3.00000E+003
+  3.00000E+003  3.20000E+003  3.20000E+003  3.20000E+003  3.40000E+003  3.40000E+003  3.40000E+003  3.60000E+003  3.60000E+003  3.60000E+003
+  3.80000E+003  3.80000E+003  3.80000E+003  4.00000E+003  4.00000E+003  4.00000E+003  0.00000E+000  0.00000E+000  0.00000E+000  2.00000E+002
+  2.00000E+002  2.00000E+002  4.00000E+002  4.00000E+002  4.00000E+002  6.00000E+002  6.00000E+002  6.00000E+002  8.00000E+002  8.00000E+002
+  8.00000E+002  1.00000E+003  1.00000E+003  1.00000E+003  1.20000E+003  1.20000E+003  1.20000E+003  1.40000E+003  1.40000E+003  1.40000E+003
+  1.60000E+003  1.60000E+003  1.60000E+003  1.80000E+003  1.80000E+003  1.80000E+003  2.00000E+003  2.00000E+003  2.00000E+003  2.20000E+003
+  2.20000E+003  2.20000E+003  2.40000E+003  2.40000E+003  2.40000E+003  2.60000E+003  2.60000E+003  2.60000E+003  2.80000E+003  2.80000E+003
+  2.80000E+003  3.00000E+003  3.00000E+003  3.00000E+003  3.20000E+003  3.20000E+003  3.20000E+003  3.40000E+003  3.40000E+003  3.40000E+003
+  3.60000E+003  3.60000E+003  3.60000E+003  3.80000E+003  3.80000E+003  3.80000E+003  4.00000E+003  4.00000E+003  4.00000E+003 -4.00000E+003
+ -4.00000E+003 -4.00000E+003 -3.80000E+003 -3.80000E+003 -3.80000E+003 -3.60000E+003 -3.60000E+003 -3.60000E+003 -3.40000E+003 -3.40000E+003
+ -3.40000E+003 -3.20000E+003 -3.20000E+003 -3.20000E+003 -3.00000E+003 -3.00000E+003 -3.00000E+003 -2.80000E+003 -2.80000E+003 -2.80000E+003
+ -2.60000E+003 -2.60000E+003 -2.60000E+003 -2.40000E+003 -2.40000E+003 -2.40000E+003 -2.20000E+003 -2.20000E+003 -2.20000E+003 -2.00000E+003
+ -2.00000E+003 -2.00000E+003 -1.80000E+003 -1.80000E+003 -1.80000E+003 -1.60000E+003 -1.60000E+003 -1.60000E+003 -1.40000E+003 -1.40000E+003
+ -1.40000E+003 -1.20000E+003 -1.20000E+003 -1.20000E+003 -1.00000E+003 -1.00000E+003 -1.00000E+003 -8.00000E+002 -8.00000E+002 -8.00000E+002
+ -6.00000E+002 -6.00000E+002 -6.00000E+002 -4.00000E+002 -4.00000E+002 -4.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -4.00000E+003
+ -4.00000E+003 -4.00000E+003 -3.80000E+003 -3.80000E+003 -3.80000E+003 -3.60000E+003 -3.60000E+003 -3.60000E+003 -3.40000E+003 -3.40000E+003
+ -3.40000E+003 -3.20000E+003 -3.20000E+003 -3.20000E+003 -3.00000E+003 -3.00000E+003 -3.00000E+003 -2.80000E+003 -2.80000E+003 -2.80000E+003
+ -2.60000E+003 -2.60000E+003 -2.60000E+003 -2.40000E+003 -2.40000E+003 -2.40000E+003 -2.20000E+003 -2.20000E+003 -2.20000E+003 -2.00000E+003
+ -2.00000E+003 -2.00000E+003 -1.80000E+003 -1.80000E+003 -1.80000E+003 -1.60000E+003 -1.60000E+003 -1.60000E+003 -1.40000E+003 -1.40000E+003
+ -1.40000E+003 -1.20000E+003 -1.20000E+003 -1.20000E+003 -1.00000E+003 -1.00000E+003 -1.00000E+003 -8.00000E+002 -8.00000E+002 -8.00000E+002
+ -6.00000E+002 -6.00000E+002 -6.00000E+002 -4.00000E+002 -4.00000E+002 -4.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -4.00000E+003
+ -4.00000E+003 -4.00000E+003 -3.80000E+003 -3.80000E+003 -3.80000E+003 -3.60000E+003 -3.60000E+003 -3.60000E+003 -3.40000E+003 -3.40000E+003
+ -3.40000E+003 -3.20000E+003 -3.20000E+003 -3.20000E+003 -3.00000E+003 -3.00000E+003 -3.00000E+003 -2.80000E+003 -2.80000E+003 -2.80000E+003
+ -2.60000E+003 -2.60000E+003 -2.60000E+003 -2.40000E+003 -2.40000E+003 -2.40000E+003 -2.20000E+003 -2.20000E+003 -2.20000E+003 -2.00000E+003
+ -2.00000E+003 -2.00000E+003 -1.80000E+003 -1.80000E+003 -1.80000E+003 -1.60000E+003 -1.60000E+003 -1.60000E+003 -1.40000E+003 -1.40000E+003
+ -1.40000E+003 -1.20000E+003 -1.20000E+003 -1.20000E+003 -1.00000E+003 -1.00000E+003 -1.00000E+003 -8.00000E+002 -8.00000E+002 -8.00000E+002
+ -6.00000E+002 -6.00000E+002 -6.00000E+002 -4.00000E+002 -4.00000E+002 -4.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002
+ -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002 -2.00000E+002  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002
+  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000
+  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+ -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002 -2.00000E+002  0.00000E+000  2.00000E+002
+cells          950
+tet         4    354    293    356    353
+tet         4    208    209    269    211
+tet         4    153    155    152     89
+tet         4     19     82     83     80
+tet         4    353    293    356    352
+tet         4    160    158    161     97
+tet         4     80     84     18     81
+tet         4    335    276    333    336
+tet         4    335    276    336    339
+tet         4    303    242    305    302
+tet         4    298    296    356    359
+tet         4    158     95     96     99
+tet         4    302    242    305    301
+tet         4    303    242    302    299
+tet         4     82     86     83    146
+tet         4    299    303    362    302
+tet         4     96     92    155    156
+tet         4     96     92    158    155
+tet         4    113    179    177    176
+tet         4    302    242    301    299
+tet         4     89     91    152    151
+tet         4    264    324    261    320
+tet         4    284    285    345    288
+tet         4    201    261    198    257
+tet         4    172    175    109    113
+tet         4    160    163     97    101
+tet         4    284    225    285    288
+tet         4    208    268    269    265
+tet         4    354    293    353    350
+tet         4    359    298    361    358
+tet         4    165    168    167    104
+tet         4    101    165    164    167
+tet         4    111    113    177    174
+tet         4    268    328    269    265
+tet         4    130    133     67    131
+tet         4    101    102    165    104
+tet         4    252    254    311    251
+tet         4    101     35     99     38
+tet         4    119    123    120    182
+tet         4    162    101     99    165
+tet         4     55    121    119    118
+tet         4    109    170    173    107
+tet         4    177    116    180    179
+tet         4    121    125    122    185
+tet         4     44     48     45    107
+tet         4    179    118    182    181
+tet         4    113    176    177    174
+tet         4    251    254    311    250
+tet         4    233    237    234    294
+tet         4    162    161     99    101
+tet         4    131     70     68    134
+tet         4     79    140     77    139
+tet         4    158     92     94    155
+tet         4    264    323    260    263
+tet         4    168    108    167    104
+tet         4     83     86     84    146
+tet         4    201    260    197    200
+tet         4    252    254    251    191
+tet         4     89     85    148     86
+tet         4    281    342    341    344
+tet         4    139    140     77    136
+tet         4    200    260    197    199
+tet         4    315    257    255    318
+tet         4    203    260    201    200
+tet         4    203    260    200    199
+tet         4     41    107     45    104
+tet         4    288    230    291    290
+tet         4    163    103    167    166
+tet         4    296    300    356    359
+tet         4    271    334    332    331
+tet         4    166    103    167    169
+tet         4    156     92    155    153
+tet         4    146    148     86    149
+tet         4     23     26     87     89
+tet         4    232    230    290    293
+tet         4    133     70     67    131
+tet         4    148    152     86    149
+tet         4    155     91     94    157
+tet         4    305    366    303    363
+tet         4    299    242    301    298
+tet         4    361    305    362    364
+tet         4    232    293    292    295
+tet         4    149    152     86    150
+tet         4    266    269    210    267
+tet         4    351    291    293    290
+tet         4    311    254    310    250
+tet         4    235    233    293    296
+tet         4    278    281    341    337
+tet         4    220    218    278    281
+tet         4     23     84     21     83
+tet         4    178    112    116    115
+tet         4    286    230    289    226
+tet         4    178    115    116    118
+tet         4     95     35     99     98
+tet         4     95     35     98     97
+tet         4    353    293    352    350
+tet         4    131    134     68     72
+tet         4      5     68      7     71
+tet         4     16     77     80     14
+tet         4     97     35     98    101
+tet         4      7     71     11      8
+tet         4    164    161    162    101
+tet         4    114    116    113     50
+tet         4    131     65     69     68
+tet         4     98     35     99    101
+tet         4    158     94     95     97
+tet         4     80     81     78    144
+tet         4     80    146     84    143
+tet         4     80     81     18     78
+tet         4    109    172    170    106
+tet         4     97    100     38    101
+tet         4    173    170    111    107
+tet         4    122    125    123    185
+tet         4    113    116    114    177
+tet         4     80     84     81    144
+tet         4     75     77    141    138
+tet         4    220    281    223    221
+tet         4    247     65      2    308
+tet         4    245    306    303    305
+tet         4    138     74     77    140
+tet         4      1     64      4      2
+tet         4    247    307     64     65
+tet         4     79     82    143    145
+tet         4     44     47     48    110
+tet         4    164    160    161    101
+tet         4     92     91    155     89
+tet         4     44     46     47    110
+tet         4     30     29     33     93
+tet         4    362    305    366    365
+tet         4    362    305    365    364
+tet         4    315    314    254    317
+tet         4     95     92     31     94
+tet         4     62    126     60    123
+tet         4    314    313    254    317
+tet         4    173    174    111    170
+tet         4    220    281    277    280
+tet         4     96     35     33     99
+tet         4    131     65     68     67
+tet         4    280    281    277    340
+tet         4     96    155    159    156
+tet         4    203    263    264    260
+tet         4    203    262    263    260
+tet         4    150     89    152     86
+tet         4    262    266    263    323
+tet         4    263    266    264    323
+tet         4    121    124     62    125
+tet         4    226    230    289    229
+tet         4    229    230    289    232
+tet         4    211    215    212    272
+tet         4    123    189    186    125
+tet         4    119    185    123    182
+tet         4    212    215    213    272
+tet         4    265    328    269    325
+tet         4    278    222    218    279
+tet         4    315    311    254    314
+tet         4    192    255    194    195
+tet         4    195    255    194    198
+tet         4    206    204    264    267
+tet         4    227    230    290    286
+tet         4    178    118    179    181
+tet         4    300    360    357    359
+tet         4    289    293    349    352
+tet         4    223    281    283    284
+tet         4    158     95     99    161
+tet         4    349    350    290    347
+tet         4    158     95    161     97
+tet         4     74    136    140    137
+tet         4    301    305    361    364
+tet         4     23     22     26     85
+tet         4    161     97     98    101
+tet         4    160     97    161    101
+tet         4    335    276    332    333
+tet         4    161     98     99    101
+tet         4    317    320    261    318
+tet         4    334    278    337    274
+tet         4    211    272    269    271
+tet         4    162    158     99    161
+tet         4    277    281    337    340
+tet         4    143     84    144    147
+tet         4      4     64     67     65
+tet         4    266    269    267    327
+tet         4    295    298    356    358
+tet         4     64     65      4      2
+tet         4    205    266    262    265
+tet         4    300    363    360    359
+tet         4     46     50    110    109
+tet         4    107    110     48    111
+tet         4    110     50     48    111
+tet         4    314    311    254    313
+tet         4    107     46    110    109
+tet         4    109    173    110    107
+tet         4    276    335    272    275
+tet         4    274    218    278    277
+tet         4    281    282    342    284
+tet         4    110    173    111    107
+tet         4     27     29     30     93
+tet         4    310    313    250    254
+tet         4     90     29     93     92
+tet         4     93     29     33     92
+tet         4    308    368    369    366
+tet         4     13     76     14     10
+tet         4    223    284    283    227
+tet         4    274    278    335    334
+tet         4     16     76     14     13
+tet         4    350    351    290    347
+tet         4    295    358    356    355
+tet         4     76     73     14     10
+tet         4      5     71      7      8
+tet         4      5     71      8      9
+tet         4      4     65     67     68
+tet         4     82    148     86    145
+tet         4    266    326    325    329
+tet         4    266    327    326    329
+tet         4     80     77     18     14
+tet         4    161     95     99     98
+tet         4     23     87     21     84
+tet         4    108    170    171    167
+tet         4    308    309    249    306
+tet         4     35     34     38     97
+tet         4     19     85     23     82
+tet         4    308    369    309    306
+tet         4    232    293    289    292
+tet         4    292    293    289    352
+tet         4    215    275    276    272
+tet         4     23     86     84     83
+tet         4     80    144     78    141
+tet         4     19     80     17     16
+tet         4    215    274    275    272
+tet         4    308     66    249    309
+tet         4     79    139     77     76
+tet         4    274    278    275    335
+tet         4    275    278    276    335
+tet         4    255    257    198    258
+tet         4    101    104    167    103
+tet         4    255    257    258    318
+tet         4    244    242    241    301
+tet         4     40    106    104    103
+tet         4    180    116    120    179
+tet         4     41     43    107    104
+tet         4    241    242    238    301
+tet         4    300    303    363    359
+tet         4    223    224    284    227
+tet         4    133    136     74    137
+tet         4     80     82    146    143
+tet         4    224    225    284    227
+tet         4    206    207    204    267
+tet         4    206    210    207    267
+tet         4    190    193    191    250
+tet         4     65     68      6     69
+tet         4    131     72     68     69
+tet         4    113    176    174    173
+tet         4    113    176    173    172
+tet         4    109    172    173    170
+tet         4    106    169    167    103
+tet         4    113    174    111    173
+tet         4     91    154    155    152
+tet         4     93     92     33     96
+tet         4    156     92    153     90
+tet         4    175    112    116    178
+tet         4     27     29     93     90
+tet         4     79    142    140    139
+tet         4    139    136     77     73
+tet         4    164    101    162    165
+tet         4    266    327    264    324
+tet         4     56     58    122    119
+tet         4     56     58     59    122
+tet         4     58     62    122    121
+tet         4    119    122     60    123
+tet         4    119    122    123    185
+tet         4    119    121    122    185
+tet         4    262    266    323    322
+tet         4    203    206    264    266
+tet         4    105     41     45    104
+tet         4    266    327    323    326
+tet         4     39     41    105    102
+tet         4    235    239    295    298
+tet         4    295    355    356    352
+tet         4    105    104     45    108
+tet         4     85     88     26     89
+tet         4     88     28     26     89
+tet         4    308     65    369    368
+tet         4     88     91     28     89
+tet         4     22     88     26     85
+tet         4     11     71     74     72
+tet         4      7     70     11     71
+tet         4     68     72     71      9
+tet         4    214    218    277    217
+tet         4    134     70     68     71
+tet         4    223    221    284    224
+tet         4    134     71     68     72
+tet         4    275    335    272    274
+tet         4    278    279    218    276
+tet         4    337    278    338    341
+tet         4    338    278    339    341
+tet         4    318    320    261    321
+tet         4    321    320    261    324
+tet         4    299    302    362    301
+tet         4    286    349    290    346
+tet         4    244    305    301    304
+tet         4    304    305    301    364
+tet         4    347    286    284    287
+tet         4    266    326    323    325
+tet         4    347    287    284    288
+tet         4    227    286    290    287
+tet         4    227    287    290    288
+tet         4    203    259    260    199
+tet         4    269    213    210    273
+tet         4    335    334    332    274
+tet         4    335    278    339    338
+tet         4    269    270    210    267
+tet         4    269    330    270    267
+tet         4    250    253    193    254
+tet         4    250    313    253    254
+tet         4    278    281    339    341
+tet         4    235    236    233    296
+tet         4    128     67    131    130
+tet         4    128     65    369    129
+tet         4      5     68     71      9
+tet         4    368     65    369    128
+tet         4    364    307    305    368
+tet         4    364    307    368    367
+tet         4    236    237    233    296
+tet         4     73    136     77     74
+tet         4    235    239    236    296
+tet         4    134     74     70     71
+tet         4    133    136     70     74
+tet         4     90     92    153     89
+tet         4     23     83     21     20
+tet         4     23     83     20     19
+tet         4     79     80     77    143
+tet         4    236    239    237    296
+tet         4    104     43    107    106
+tet         4     53     52    116     50
+tet         4    205    262    203    202
+tet         4    202    262    203    199
+tet         4    127    307     65     64
+tet         4    182    121    185    184
+tet         4    135     72    131     69
+tet         4     10     73     14     11
+tet         4     56    122     60    119
+tet         4     95     35     96     99
+tet         4    121    125    185    184
+tet         4    182    185    123    186
+tet         4     89    151    152    148
+tet         4    106    170    107    104
+tet         4    102     41    105    104
+tet         4    102    168    165    104
+tet         4    101     40    104    103
+tet         4     31     97     35     94
+tet         4    160    163    101    164
+tet         4    163    103    101    164
+tet         4     46    112     50    109
+tet         4     42     41     45    105
+tet         4     52     55     53    116
+tet         4    168    108    171    167
+tet         4     89    151    148     85
+tet         4     80    144    141    143
+tet         4     79    143     77    140
+tet         4    359    299    301    298
+tet         4    245    249    306    305
+tet         4    133    137     74    134
+tet         4    134    137     74    135
+tet         4    133     70    131    134
+tet         4    135     74    134     72
+tet         4    117     53    114     51
+tet         4     91    154    157    155
+tet         4     57    116    120    117
+tet         4    100     40    101    103
+tet         4     23     85     89     86
+tet         4     34    100     38     97
+tet         4    112     52    116    115
+tet         4    163    100     97    101
+tet         4    164    103    167    163
+tet         4    111    174    108    170
+tet         4     38    104    102    101
+tet         4    121    187    124    125
+tet         4    180    182    120    183
+tet         4    163    103    100    101
+tet         4    119     58    122    121
+tet         4    183    182    120    186
+tet         4     57     53    117     54
+tet         4     83     80     21     17
+tet         4     83     84     21     80
+tet         4    199    262    203    259
+tet         4    254    256    196    257
+tet         4    263    323    260    262
+tet         4    286    284    287    227
+tet         4     80     83     84    146
+tet         4    211    271    274    215
+tet         4     80     82     83    146
+tet         4    281    225    282    284
+tet         4    224    221    284    225
+tet         4     20     83     21     17
+tet         4     23     82     83     19
+tet         4    308     65    249     66
+tet         4    308     65     66    369
+tet         4      2     65      6     66
+tet         4      2     65     66    249
+tet         4    215    218    278    274
+tet         4    211    212    209    272
+tet         4    215    274    272    271
+tet         4     46     49     50    112
+tet         4     49     52     50    112
+tet         4    325    329    269    266
+tet         4    357    297    300    296
+tet         4    346    283    343    344
+tet         4    346    283    344    286
+tet         4     14     78     15     12
+tet         4    252    254    191    194
+tet         4    296    239    300    299
+tet         4    247     65     64      2
+tet         4    128    127     65     67
+tet         4    262    322    320    259
+tet         4    271    331    332    328
+tet         4    281    222    279    282
+tet         4    351    293    350    290
+tet         4    329    332    273    330
+tet         4    266    324    264    323
+tet         4    350    293    349    290
+tet         4    350    293    352    349
+tet         4    332    276    273    333
+tet         4    244    242    301    305
+tet         4    319    259    320    322
+tet         4      4     68      5      2
+tet         4      5     68      6      2
+tet         4    247    308    245    305
+tet         4    313    256    317    316
+tet         4    300    303    359    299
+tet         4    300    242    303    299
+tet         4    300    239    242    299
+tet         4    100     40     38    101
+tet         4     38     40     41    104
+tet         4    299    239    242    298
+tet         4    102     38     36     39
+tet         4    298    299    296    359
+tet         4     52     55    116    118
+tet         4    299    300    296    359
+tet         4    281    222    282    225
+tet         4      4     67     70     68
+tet         4    284    285    282    345
+tet         4     68     72      6     69
+tet         4     11     75     12      9
+tet         4     26     27     24     90
+tet         4     26     29     90     92
+tet         4     23     85     86     82
+tet         4    284    282    342    345
+tet         4    284    225    282    285
+tet         4     96     92    156     93
+tet         4    107    170    108    104
+tet         4    265    325    269    266
+tet         4    268    271    269    328
+tet         4    247    308      2    248
+tet         4     55     58     56    119
+tet         4    121     62    122    125
+tet         4    248    308      2    249
+tet         4    247    308    248    245
+tet         4    248    308    249    245
+tet         4    228    230    231    291
+tet         4     82     79     80     16
+tet         4    143    141     77    140
+tet         4     80    141     77    143
+tet         4     80     84    144    143
+tet         4    211    274    214    215
+tet         4    274    218    277    214
+tet         4    114     50    111     48
+tet         4    201    257    198    197
+tet         4    197    199    257    196
+tet         4     53    119     57    116
+tet         4    113    175    179    176
+tet         4    179    182    120    180
+tet         4    252    255    194    192
+tet         4    145    148     86    146
+tet         4    143     82    146    145
+tet         4    315    257    254    255
+tet         4    146     86     84    147
+tet         4    217    218    277    220
+tet         4    131    132     69    129
+tet         4      2     66      3    249
+tet         4    113     52    116    112
+tet         4    113    112    116    175
+tet         4    175    112    109    113
+tet         4     94     35     95     97
+tet         4    286    223    283    227
+tet         4    279    222    218    219
+tet         4    281    342    339    341
+tet         4    238    242    298    301
+tet         4    359    362    361    301
+tet         4    245    308    249    305
+tet         4    161     95     98     97
+tet         4    245    246    243    306
+tet         4    317    259    320    316
+tet         4    316    259    320    319
+tet         4    196    259    257    256
+tet         4    254    196    193    194
+tet         4    203    262    260    259
+tet         4    317    256    257    259
+tet         4    201    260    257    197
+tet         4     32     29     31     95
+tet         4    215    276    213    273
+tet         4    218    222    278    281
+tet         4     95     31     35     94
+tet         4     96     92     95    158
+tet         4    277    218    278    220
+tet         4    230    291    290    293
+tet         4    292    293    352    295
+tet         4    351    293    354    350
+tet         4     95     92     94    158
+tet         4     80    141     78     77
+tet         4     14     77     18     78
+tet         4     14     77     78     75
+tet         4    250    254    193    191
+tet         4     75     77     78    141
+tet         4    251    254    250    191
+tet         4     47     50     48    110
+tet         4     41     44     45    107
+tet         4    266    324    323    327
+tet         4     41     43     44    107
+tet         4     26     29     27     90
+tet         4    271    329    269    328
+tet         4    362    301    305    361
+tet         4    359    303    363    362
+tet         4    344    284    342    345
+tet         4    347    286    344    284
+tet         4    344    284    345    288
+tet         4    346    286    344    347
+tet         4    296    239    299    298
+tet         4    362    305    363    366
+tet         4    365    305    366    368
+tet         4    242    303    240    300
+tet         4    344    347    288    345
+tet         4    368    307    128    367
+tet         4    271    328    332    329
+tet         4     14     77     75     11
+tet         4    208    271    269    268
+tet         4     29     28     31     92
+tet         4    208    265    269    266
+tet         4    312    254    311    252
+tet         4     84     87    150     86
+tet         4    311    310    254    313
+tet         4    233    237    294    293
+tet         4     43     46    107    109
+tet         4    354    294    296    293
+tet         4    119    121    185    182
+tet         4     53     55    119    116
+tet         4    294    237    296    293
+tet         4     76    139     77     73
+tet         4     11     77     75     74
+tet         4     72     74    138    135
+tet         4    291    294    354    293
+tet         4    220    281    278    277
+tet         4    277    281    278    337
+tet         4     79    143    140    142
+tet         4     26     28     92     89
+tet         4     33     95     35     96
+tet         4    215    218    274    214
+tet         4    274    278    337    277
+tet         4    264    320    261    260
+tet         4     44    110     48    107
+tet         4    194    257    197    198
+tet         4     44     46    110    107
+tet         4    109     50    110    113
+tet         4    128     65    129     69
+tet         4    128     65     69    131
+tet         4    367    307    128    127
+tet         4    128    127     67    130
+tet         4    368    307     65    128
+tet         4    129     65     66     69
+tet         4    260    259    257    199
+tet         4    260    262    320    259
+tet         4    203    264    261    260
+tet         4    146    149     86    150
+tet         4     91    154    152    151
+tet         4     14     75     78     12
+tet         4    201    261    257    260
+tet         4    245    243    303    306
+tet         4     96    155    158    159
+tet         4    227    230    286    226
+tet         4    159    158     96    162
+tet         4    209    269    210    206
+tet         4    175    178    116    179
+tet         4    209    213    210    269
+tet         4    118    121    182    184
+tet         4    211    209    269    272
+tet         4    209    213    269    272
+tet         4    116     57    120    119
+tet         4    317    257    261    320
+tet         4    197    260    257    199
+tet         4    254    256    193    196
+tet         4    252    194    191    192
+tet         4    218    216    215    276
+tet         4    212    213    209    272
+tet         4    208    269    206    266
+tet         4    113    175    116    179
+tet         4    113    175    176    172
+tet         4    162    158     96     99
+tet         4    288    227    230    290
+tet         4    233    237    293    296
+tet         4    347    284    344    288
+tet         4     33     92     95     96
+tet         4    160     94    158     97
+tet         4    244    307    247    305
+tet         4    245    249    246    306
+tet         4      4     70      7     68
+tet         4     95     29     31     92
+tet         4    157     94    158    160
+tet         4    206    269    210    266
+tet         4    264    323    320    260
+tet         4    138     77    141    140
+tet         4    203    266    264    263
+tet         4     80     78     18     77
+tet         4     58    124     62    121
+tet         4    278    218    215    276
+tet         4    281    282    279    342
+tet         4    287    284    288    227
+tet         4    283    281    341    344
+tet         4    335    278    338    337
+tet         4     46     50     47    110
+tet         4    232    235    293    295
+tet         4    235    239    296    295
+tet         4    107     48     45    111
+tet         4    303    305    363    362
+tet         4    107     45    108    111
+tet         4    109    170    107    106
+tet         4    317    256    259    316
+tet         4     43    109    107    106
+tet         4    317    259    257    320
+tet         4     26     92     90     89
+tet         4    269    330    267    327
+tet         4    269    273    210    270
+tet         4    269    273    270    330
+tet         4     26     90     24     87
+tet         4     73     77     14     11
+tet         4     11     73     77     74
+tet         4    283    284    281    344
+tet         4    136     73     70     74
+tet         4    308     65      2    249
+tet         4     76     77     14     73
+tet         4     16     77     14     76
+tet         4     56     60     57    119
+tet         4     53     56     57    119
+tet         4     70     73      7     11
+tet         4     21     84     18     80
+tet         4     53     55     56    119
+tet         4     33     32     35     95
+tet         4     32     31     35     95
+tet         4     25     28     26     88
+tet         4     19     22     23     85
+tet         4     11     70     74     71
+tet         4    110     50    111    113
+tet         4     51     53    114     50
+tet         4      5     68      9      6
+tet         4    127     64     65     67
+tet         4    233    294    291    293
+tet         4    305    307    308    368
+tet         4    308    307     65    368
+tet         4     66     65      6     69
+tet         4    300    357    356    359
+tet         4    233    234    231    294
+tet         4     21     80     18     17
+tet         4     23     86     83     82
+tet         4    294    237    297    296
+tet         4     38     41     39    102
+tet         4    343    283    340    344
+tet         4    280    281    340    283
+tet         4     33     29     32     95
+tet         4    116    119    120    182
+tet         4    116     52    113     50
+tet         4    116    118    119    182
+tet         4    220    281    280    283
+tet         4    114     50    113    111
+tet         4    337    281    341    340
+tet         4    340    281    341    283
+tet         4    266    323    322    325
+tet         4    264    324    320    323
+tet         4    128     65    131     67
+tet         4    308     66    309    369
+tet         4    134     74     71     72
+tet         4     89    153     87    150
+tet         4     87     89    150     86
+tet         4     84    150    147     86
+tet         4    308    369    306    366
+tet         4    143    146     84    147
+tet         4    146    150     86    147
+tet         4    106    167    104    103
+tet         4    101    167    164    103
+tet         4    305    366    306    303
+tet         4    289    349    290    286
+tet         4    286    230    290    289
+tet         4    106    167    170    104
+tet         4    226    223    286    227
+tet         4    185    125    123    186
+tet         4    116    182    120    179
+tet         4     58     61     62    124
+tet         4    289    230    290    232
+tet         4    122     62    123    125
+tet         4    122     62     60    123
+tet         4    203    266    263    262
+tet         4    260    323    320    262
+tet         4    232    293    290    289
+tet         4    203    264    201    261
+tet         4    211    215    272    271
+tet         4    289    293    290    349
+tet         4    279    219    218    216
+tet         4    276    332    273    272
+tet         4    215    216    213    276
+tet         4    276    335    332    272
+tet         4    215    278    275    274
+tet         4    272    335    332    274
+tet         4    215    276    273    272
+tet         4    272    213    269    273
+tet         4    232    233    293    235
+tet         4    255    257    254    198
+tet         4    312    315    254    252
+tet         4    255    254    194    198
+tet         4    258    257    198    261
+tet         4    238    239    298    242
+tet         4    269    332    273    329
+tet         4    271    332    269    329
+tet         4    271    332    272    269
+tet         4    266    329    269    327
+tet         4    221    222    218    281
+tet         4    265    266    262    325
+tet         4    223    281    284    221
+tet         4    286    283    344    284
+tet         4    221    281    284    225
+tet         4    221    222    281    225
+tet         4    203    204    264    206
+tet         4    184    187    125    188
+tet         4    108    174    171    170
+tet         4    359    301    361    298
+tet         4    172    169    170    106
+tet         4    242    245    303    305
+tet         4    305    308    249    306
+tet         4    106    169    170    167
+tet         4    153     92    155     89
+tet         4     73     10      7     11
+tet         4    203    261    201    260
+tet         4    136    140     77     74
+tet         4    118    121    119    182
+tet         4    323    322    320    262
+tet         4    274    334    332    271
+tet         4    340    283    341    344
+tet         4    215    278    276    275
+tet         4    185    125    186    188
+tet         4    119     57    120    123
+tet         4    235    296    293    295
+tet         4     55     58    119    121
+tet         4    354    296    357    356
+tet         4     56     59     60    122
+tet         4    357    296    300    356
+tet         4    305    308    306    366
+tet         4    247     64      1      2
+tet         4    105    108    168    104
+tet         4      2     66      6      3
+tet         4    165    101     99    102
+tet         4     99     38     36    102
+tet         4    102    105    168    104
+tet         4     22     25     26     88
+tet         4     99     35     36     38
+tet         4     99     35     33     36
+tet         4    101     38     99    102
+tet         4     58     62     59    122
+tet         4     82     85     86    148
+tet         4     54     53    117     51
+tet         4     89     88    151     85
+tet         4     71     72     11      9
+tet         4    131     70     67     68
+tet         4      8     71     11      9
+tet         4     23     85     26     89
+tet         4     89     91    151     88
+tet         4     26     90     87     89
+tet         4     23     89     87     86
+tet         4     23     87     84     86
+tet         4    111    170    108    107
+tet         4     40     43    104    106
+tet         4     40     43     41    104
+tet         4     38     40    104    101
+tet         4     68     70      7     71
+tet         4    104     45    108    107
+tet         4    119     60     57    123
+tet         4      4     68      7      5
+tet         4      4     65     68      2
+tet         4     74    137    138    135
+tet         4     74    137    140    138
+tet         4     68     65      6      2
+tet         4     68     72      9      6
+tet         4    117    116    120    180
+tet         4    155     91    152     89
+tet         4     89     91     28     92
+tet         4    150    153    152     89
+tet         4     51     50    114     48
+tet         4     37     40     38    100
+tet         4     31     34     35     97
+tet         4     16     80     17     14
+tet         4     17     80     18     14
+tet         4    113    173    111    110
+tet         4    170    167    108    104
+tet         4    109    112     50    113
+tet         4    131    129     69    128
+tet         4    131    135     69    132
+tet         4    295    356    293    352
+tet         4    125    189    186    188
+tet         4    335    278    276    339
+tet         4    332    333    273    330
+tet         4    205    206    203    266
+tet         4    205    208    206    266
+tet         4    315    317    257    318
+tet         4    318    317    257    261
+tet         4    252    254    194    255
+tet         4    318    257    258    261
+tet         4    345    347    288    348
+tet         4    220    221    218    281
+tet         4    220    281    283    223
+tet         4    302    301    305    362
+tet         4    235    239    298    238
+tet         4    239    242    240    300
+tet         4    239    240    237    300
+tet         4    304    305    364    307
+tet         4    317    256    254    257
+tet         4    252    315    254    255
+tet         4    315    317    254    257
+tet         4    199    259    257    196
+tet         4    348    347    288    351
+tet         4    272    332    273    269
+tet         4    287    347    290    288
+tet         4    272    274    332    271
+tet         4    286    283    284    227
+tet         4    281    284    342    344
+tet         4    284    225    288    227
+tet         4    334    278    335    337
+tet         4    203    204    201    264
+tet         4    351    291    354    293
+tet         4    351    288    291    290
+tet         4    206    267    264    266
+tet         4    225    227    228    288
+tet         4    247    307     65    308
+tet         4    247    307    308    305
+tet         4    242    243    303    245
+tet         4    369     65     66    129
+tet         4    305    368    308    366
+tet         4    328    329    269    325
+tet         4    206    210    267    266
+tet         4    266    327    267    264
+tet         4    205    265    208    266
+tet         4    253    256    193    254
+tet         4    230    231    291    233
+tet         4    230    233    291    293
+tet         4    297    237    300    296
+tet         4    298    359    356    358
+tet         4    295    239    296    298
+tet         4    294    297    357    296
+tet         4    232    233    230    293
+tet         4     16     79     77     76
+tet         4     70     73     11     74
+tet         4     23     26     24     87
+tet         4     23     24     21     87
+tet         4    279    216    218    276
+tet         4     43     46     44    107
+tet         4    194    257    196    197
+tet         4    257    259    260    320
+tet         4    116    118    182    179
+tet         4    184    187    121    125
+tet         4    184    125    185    188
+tet         4     59     62     60    122
+tet         4    118    184    182    181
+tet         4     14     78     18     15
+tet         4     75     74     77    138
+tet         4     72     74     75    138
+tet         4     11     74     75     72
+tet         4     11     72     75      9
+tet         4     92     91     94    155
+tet         4    116     55    119    118
+tet         4     28     91     94     92
+tet         4     92     28     31     94
+tet         4     33     29     95     92
+tet         4    182    123    120    186
+tet         4    101    165    167    104
+tet         4    113    172    173    109
+tet         4    178    118    116    179
+tet         4    112     52     50    113
+tet         4    113    173    110    109
+tet         4    126    189    123    125
+tet         4     62    126    123    125
+tet         4     38     41    102    104
+tet         4     39     41     42    105
+tet         4     89    148    152     86
+tet         4     19     83     20     17
+tet         4    133     74     70    134
+tet         4    135    134    131     72
+tet         4     79    145    143    142
+tet         4    114    116    117    180
+tet         4    114     53    116     50
+tet         4    111    113    114    177
+tet         4     35     97     38    101
+tet         4     34     37     38    100
+tet         4    158    155     94    157
+tet         4     57     53    116    117
+tet         4    117     53    116    114
+tet         4     79     82     80    143
+tet         4    114    116    180    177
+tet         4    113    116    177    179
+tet         4     82     86    146    145
+tet         4    254    257    194    198
+tet         4    254    257    196    194
+tet         4     16     79     80     77
+tet         4     19     80     83     17
+tet         4    191    254    193    194
+tet         4    272    215    213    273
+tet         4    208    211    269    271
+tet         4    208    269    209    206
+tet         4     19     82     80     16
+tet         4     62     63     60    126
+tet         4    233    231    291    294
+tet         4    228    230    291    288
+tet         4    228    227    230    288
+tet         4    242    244    245    305
+tet         4    115     52    116    118
+tet         4     93     92    156     90
+tet         4    242    243    240    303
+tet         4     89     90     87    153
+tet         4    244    247    245    305
+tet         4     26     28     29     92
+tet         4     14     75     12     11
+tet         4    303    302    305    362
+tet         4    205    266    203    262
+tet         4    257    260    261    320
+tet         4    262    266    322    325
+tet         4    278    279    276    339
+tet         4    278    222    279    281
+tet         4    278    279    339    281
+tet         4    364    305    365    368
+tet         4    281    279    339    342
+tet         4    359    299    362    301
+tet         4    296    239    237    300
+tet         4    295    296    293    356
+tet         4    354    294    357    296
+tet         4    354    296    356    293
+tet         4    295    296    356    298
+tet         4    359    303    362    299
+tet         4    313    256    253    254
+tet         4    312    311    254    315
+tet         4    313    256    254    317
+tet         4    244    305    304    307
+tet         4    128    307     65    127
+tet         4    286    346    290    347
+tet         4    286    347    290    287
+tet         4    288    347    290    351
+tet         4    329    330    269    327
+tet         4    329    273    269    330
+tet         4    349    347    290    346
+variable
+imt1                             1
+  1.00000E+000  2.00000E+000  2.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000
+  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  1.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000  2.00000E+000
+itp1                             1
+  1.20000E+001  1.20000E+001  1.20000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.20000E+001  2.00000E+000  1.20000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001
+  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000
+  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001
+  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001
+  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000
+  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001
+  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.20000E+001  1.20000E+001  1.20000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000
+  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001
+  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001
+  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000
+  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001
+  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001  0.00000E+000  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001  1.00000E+001
+icr1                             1
+  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001
+  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000
+  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001
+  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001
+  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000
+  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001
+  1.20000E+001  7.00000E+000  2.00000E+001  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000
+  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000
+  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000
+  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000
+  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000
+  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000
+  5.00000E+000  0.00000E+000  6.00000E+000  1.30000E+001  2.00000E+000  2.10000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001
+  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000
+  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001
+  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001
+  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000
+  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001
+  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.50000E+001  9.00000E+000  2.30000E+001  1.60000E+001
+  8.00000E+000  2.40000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000
+  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001
+  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001
+  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000
+  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001
+  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.10000E+001  3.00000E+000  1.90000E+001  1.70000E+001
+  1.00000E+000  2.50000E+001  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000
+  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000
+  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000
+  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000
+  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000
+  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  5.00000E+000  0.00000E+000  6.00000E+000  1.80000E+001
+  1.00000E+001  2.60000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000
+  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001
+  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001
+  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000
+  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001
+  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001  1.40000E+001  4.00000E+000  2.20000E+001
+isn1                             1
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000  0.00000E+000
+endvars
+flags    
+imt1            3         1
+mr_xpos 
+mr_xneg 
+ERROR
+    1    2    2    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2
+endflag
+material            3         0
+mr_xpos 
+mr_xneg 
+ERROR
+    2    2    1    1    2    1    1    2    2    2    2    1    2    2    1    2    1    1    1    2
+    1    2    2    2    1    1    2    2    2    2    1    1    1    2    1    1    2    1    1    1
+    1    1    1    1    1    1    1    2    2    1    1    1    1    2    1    1    2    2    1    2
+    1    2    2    2    2    1    2    1    2    2    1    1    1    1    2    1    1    1    2    2
+    2    2    1    2    2    2    2    2    2    1    1    2    1    1    1    2    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    2    2    2    1
+    1    2    1    1    1    1    1    1    2    2    2    1    1    2    1    2    1    1    2    1
+    2    2    1    2    2    1    2    2    2    1    1    2    2    2    2    2    2    2    2    1
+    2    2    2    1    2    1    1    2    1    1    1    2    1    2    2    2    1    2    1    1
+    2    2    1    2    2    1    1    1    2    1    1    2    2    2    1    1    2    1    1    2
+    1    2    2    1    2    2    1    1    1    1    1    2    2    1    1    1    1    2    1    1
+    2    2    2    2    1    1    1    2    2    1    2    2    2    1    2    2    1    1    1    2
+    2    2    1    1    2    2    2    2    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    2    1    1    1    1    1    1    2    2    1    2    1    2    2    1    1    1
+    2    1    1    1    1    1    2    1    2    1    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    1    2    1    2    2
+    2    2    1    2    1    1    1    1    1    1    2    1    1    2    2    2    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    2    2
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    2    2    2    2    1    2    1    2    2    1    1    2    2    1    2    2    2
+    2    1    1    2    2    2    2    1    2    2    2    1    2    2    2    2    2    2    2    2
+    2    2    2    1    1    2    2    2    2    2    1    1    2    1    2    1    2    2    1    2
+    1    1    1    1    1    2    2    1    1    2    2    2    1    1    2    2    2    2    1    1
+    1    1    2    2    1    2    2    1    1    1    2    1    1    2    1    2    1    2    1    1
+    1    1    2    2    2    2    2    2    1    2    2    2    2    2    2    2    2    1    2    2
+    1    1    2    2    2    2    1    1    1    1    2    1    2    1    1    2    1    1    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    1    2    1    2    2    1    2    2
+    1    2    1    1    2    1    1    1    2    2    2    1    1    1    2    2    2    1    2    1
+    1    1    1    2    1    2    1    2    2    2    1    1    1    2    2    1    2    1    2    1
+    2    1    2    2    1    2    2    2    2    2    2    2    1    1    1    2    2    2    1    1
+    2    2    1    1    1    2    2    1    2    1    1    2    2    2    2    2    1    2    2    1
+    2    1    1    2    1    2    1    2    2    2    1    1    1    2    1    2    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    2    2    2    1    2    2    1    1
+    2    1    2    2    1    1    1    1    2    1    2    2    2    2    1    2    1    1    1    1
+    2    1    1    1    1    2    2    2    1    2    1    1    1    2    1    1    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    1    1    2    1    2    2    1    1    1    2    1    1    2
+    2    2    2    1    1    2    1    2    1    2    2    2    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    2    1    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2    2    2    1    1    1    1    2    1    2    2
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
+    1    1    2    2    1    1    2    2    2    2    1    1    2    2    2    2    1    1    2    1
+    2    1    1    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2    2
+    2    2    2    2    2    2    2    2    2    2
+cycleno          0
+probtime   0.00000E+000
+endgmv

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.pset
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.pset	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/bar_tet4_200m.pset	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,44 @@
+pset ascii         4
+fault    1          9
+         1          2          3         64         65         66        127        128        129
+boundary_xm    2          9
+       190        191        192        250        251        252        310        311        312
+boundary_xp    3          9
+        61         62         63        124        125        126        187        188        189
+pts_notfault    4        360
+         4          5          6          7          8          9         10         11         12         13
+        14         15         16         17         18         19         20         21         22         23
+        24         25         26         27         28         29         30         31         32         33
+        34         35         36         37         38         39         40         41         42         43
+        44         45         46         47         48         49         50         51         52         53
+        54         55         56         57         58         59         60         61         62         63
+        67         68         69         70         71         72         73         74         75         76
+        77         78         79         80         81         82         83         84         85         86
+        87         88         89         90         91         92         93         94         95         96
+        97         98         99        100        101        102        103        104        105        106
+       107        108        109        110        111        112        113        114        115        116
+       117        118        119        120        121        122        123        124        125        126
+       130        131        132        133        134        135        136        137        138        139
+       140        141        142        143        144        145        146        147        148        149
+       150        151        152        153        154        155        156        157        158        159
+       160        161        162        163        164        165        166        167        168        169
+       170        171        172        173        174        175        176        177        178        179
+       180        181        182        183        184        185        186        187        188        189
+       190        191        192        193        194        195        196        197        198        199
+       200        201        202        203        204        205        206        207        208        209
+       210        211        212        213        214        215        216        217        218        219
+       220        221        222        223        224        225        226        227        228        229
+       230        231        232        233        234        235        236        237        238        239
+       240        241        242        243        244        245        246        247        248        249
+       250        251        252        253        254        255        256        257        258        259
+       260        261        262        263        264        265        266        267        268        269
+       270        271        272        273        274        275        276        277        278        279
+       280        281        282        283        284        285        286        287        288        289
+       290        291        292        293        294        295        296        297        298        299
+       300        301        302        303        304        305        306        307        308        309
+       310        311        312        313        314        315        316        317        318        319
+       320        321        322        323        324        325        326        327        328        329
+       330        331        332        333        334        335        336        337        338        339
+       340        341        342        343        344        345        346        347        348        349
+       350        351        352        353        354        355        356        357        358        359
+       360        361        362        363        364        365        366        367        368        369

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/matprops.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/matprops.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/matprops.spatialdb	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+  // number of physical properties
+  num-values = 3
+
+  // Names and units of physical properties
+  value-names =  density vs vp
+  value-units =  kg/m^3  m/s  m/s
+
+  // Number of locations where physical properties are given
+  // 1 location -> data is uniform -> data dimension is 0
+  num-locs = 1
+  data-dim = 0
+
+  // Problem is in 3-D
+  space-dim = 3
+
+  // Coordinate system (3-D with coordinates in meters)
+  cs-data = cartesian {
+    to-meters = 1.0
+    space-dim = 3
+  }
+}
+// Columns are x, y, z, density, vs, vp
+0.0  0.0  0.0   2500.0  1000.0  1732.0

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/mesh_tet4_200m.lagrit
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/mesh_tet4_200m.lagrit	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/mesh_tet4_200m.lagrit	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,186 @@
+* -*- LaGriT -*-
+*
+* Mesh of a bar with a through-going, vertical, strike-slip fault.
+*
+* Bar is 400m x 400m x 8000m
+* -4000 m <= x <= 4000 m
+* -200 m <= y <= 200 m
+* -200 m <= z <= 200 km
+* 
+* The domain is constructed with the region command after constructing the
+* boundaries of the domain using the surface command with the plane
+* attribute (construct the 3-D domain from boundary surfaces). The
+* fault plane is also a plane.
+*
+* ----------------------------------------------------------------------
+* PARAMETERS CONTROLLING MESH
+* ----------------------------------------------------------------------
+*
+* Boundaries of domain
+define / domain_xm / -4000.0
+define / domain_xp /  4000.0
+define / domain_ym / -200.0
+define / domain_yp / 200.0
+define / domain_zm / -200.0
+define / domain_zp /  200.0
+*
+* Mesh size and number of points along x & y
+define / dx / 200.0
+define / nx / 41
+define / ny / 3
+*
+* ----------------------------------------------------------------------
+* CREATE GEOMETRY
+* ----------------------------------------------------------------------
+*
+* Create a mesh object and name it box
+cmo / create / box /
+*
+* Create boundary surfaces
+surface / surf_xm / reflect / plane / &
+   domain_xm   0.0e+0   0.0e+0 / &
+   domain_xm   0.0e+0  10.0e+3 / &
+   domain_xm  10.0e+3   0.0e+0
+surface / surf_xp / reflect / plane / &
+    domain_xp   0.0e+0   0.0e+0 / &
+    domain_xp  10.0e+3   0.0e+0 / &
+    domain_xp   0.0e+0  10.0e+3
+surface / surf_ym / reflect / plane / &
+     0.0e+0  domain_ym   0.0e+0 / &
+     0.0e+0  domain_ym  10.0e+3 / &
+   -10.0e+3  domain_ym   0.0e+0
+surface / surf_yp / reflect / plane / &
+     0.0e+0   domain_yp   0.0e+0 / &
+   -10.0e+3   domain_yp   0.0e+0 / &
+     0.0e+0   domain_yp  10.0e+3
+surface / surf_zm / reflect / plane / &
+     0.0e+0    0.0e+0  domain_zm / &
+   -10.0e+3    0.0e+0  domain_zm / &
+     0.0e+0   10.0e+3  domain_zm
+surface / surf_zp / reflect / plane / &
+     0.0e+0    0.0e+0    domain_zp / &
+     0.0e+0   10.0e+3    domain_zp / &
+   -10.0e+3    0.0e+0    domain_zp
+*
+* Create fault plane
+surface / surf_fault / intrface / plane / &
+  0.0  0.0  0.0 / &
+  0.0  1.0  0.0 / &
+  0.0  0.0  1.0 /
+
+* ----------------------------------------------------------------------
+* Create regions
+* ----------------------------------------------------------------------
+
+** Material on +x side of fault
+region / r_xpos / &
+  le surf_xp and ge surf_fault  &
+  and le surf_yp and le surf_ym  &
+  and le surf_zp and le surf_zm /
+
+** Material on -x side of the fault
+region / r_xneg / &
+  le surf_xm and lt surf_fault  &
+  and le surf_yp and le surf_ym  &
+  and le surf_zp and le surf_zm /
+
+* ----------------------------------------------------------------------
+* Create material regions
+* ----------------------------------------------------------------------
+* Material regions do not include internal interfaces (regions do)
+*
+mregion / mr_xpos / &
+  le surf_xp and gt surf_fault  &
+  and le surf_yp and le surf_ym  &
+  and le surf_zp and le surf_zm /
+mregion / mr_xneg / &
+  le surf_xm and lt surf_fault  &
+  and le surf_yp and le surf_ym  &
+  and le surf_zp and le surf_zm /
+*
+* ----------------------------------------------------------------------
+* CREATE MESH
+* ----------------------------------------------------------------------
+*
+createpts / xyz / nx,ny,1 / &
+  domain_xm  domain_ym  domain_zp / &
+  domain_xp  domain_yp  domain_zp / &
+  1,1,0 /
+pset / rayend / seq / 0,0,0
+*
+
+loop / foreach / iregion / &
+  r_xpos r_xneg  &
+  / loop_end / &
+  regnpts / iregion / dx / pset,get,rayend / xyz / &
+    0.0e+3  0.0e+3  domain_zm / &
+    0.0e+3  1.0e+3  domain_zm / &
+    1.0e+3  0.0e+3  domain_zm / &
+    0,0
+cmo/setatt//itp/pset,get,rayend/dud
+pset / rayend / delete
+
+* ----------------------------------------------------------------------
+* Remove duplicate points
+* ----------------------------------------------------------------------
+filter / 1,0,0 / 10.0
+rmpoint / compress
+
+* ----------------------------------------------------------------------
+* Set node type and material
+* ----------------------------------------------------------------------
+setpts
+
+* ----------------------------------------------------------------------
+* Connect the points into a Delaunay tetrahedral mesh
+* ----------------------------------------------------------------------
+connect
+
+* ----------------------------------------------------------------------
+* Set elements to tetrahedra
+* ----------------------------------------------------------------------
+settets / geometry
+resetpts/cell_color
+
+* Remove parent/child hierarchy
+resetpts
+rmpoint / compress
+*
+* ----------------------------------------------------------------------
+* MESH QUALITY
+* ----------------------------------------------------------------------
+*
+quality
+*
+* ----------------------------------------------------------------------
+* IDENTIFY GROUPS
+* ----------------------------------------------------------------------
+*
+* Identify nodes on fault surface
+pset / fault / surface / surf_fault /
+*
+* Identify nodes on boundaries
+pset / boundary_xm / surface / surf_xm /
+pset / boundary_xp / surface / surf_xp /
+*
+* Identify all nodes except those on fault
+pset / tmp / seq / 1,0,0
+pset / pts_notfault / not / tmp / fault
+pset / tmp / delete
+*
+* ----------------------------------------------------------------------
+* WRITE MESH
+* ----------------------------------------------------------------------
+*
+* Export mesh and psets to files.
+cmo / setatt / box / ipolydat / no
+dump / gmv / bar_tet4_200m.gmv / / ascii
+
+pset / -all- / write / bar_tet4_200m.pset / ascii
+
+*
+* Terminate processing
+finish
+
+* ======================================================================
+* End of file

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave.cfg	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave.cfg	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,246 @@
+# -*- Python -*-
+
+# The settings in this file (shearwave.cfg) will be read if it is
+# specified on the command line:
+# 'pylith shearwave.cfg'
+#
+# If this file is given, the settings in this file will override all
+# others except for parameters specified directly on the command line.
+#
+# These settings define an explicit time-stepping problem using a fault with
+# kinematically-specified slip.
+
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# journal
+# ----------------------------------------------------------------------
+[pylithapp.journal.info]
+timedependent = 1
+#explicit = 1
+#petsc = 1
+#solverlinear = 1
+#meshiolagrit = 1
+#homogeneous = 1
+#explicitelasticity = 1
+#quadrature2d = 1
+#fiatlagrange = 1
+#faultcohesivekin = 1
+
+# ----------------------------------------------------------------------
+# mesh_generator
+# ----------------------------------------------------------------------
+[pylithapp.mesh_generator]
+#debug = 1
+importer = pylith.meshio.MeshIOLagrit
+
+[pylithapp.mesh_generator.importer]
+filename_gmv = bar_tet4_200m.gmv
+filename_pset = bar_tet4_200m.pset
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+total_time = 12.0*s
+default_dt = 0.0333333333333333333*s
+dimension = 3
+
+# Change to an explicit time stepping formulation
+formulation = pylith.problems.Explicit
+
+# Reuse solution from previous time step as initial guess for next time step
+formulation.solver.initial_guess_zero = False
+
+# Set bc to an array with 3 boundary conditions: 'x_neg', 'x_pos', and 'y_pos'.
+bc = [x_pos,x_neg,y_pos]
+
+# Change the DirichletPoints BC to the AbsorbingDampers BC
+bc.x_pos = pylith.bc.AbsorbingDampers
+bc.x_neg = pylith.bc.AbsorbingDampers
+
+# Set interfaces to an array with 1 fault: 'fault'.
+interfaces = [fault]
+
+# Set materials to an array with 2 materials 'pos' and 'neg'.
+materials = [pos,neg]
+
+# ----------------------------------------------------------------------
+# materials
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.materials.pos]
+
+# Label for material
+label = elastic material +x
+
+# Material id (itetclr) from LaGriT
+id = 1
+
+# Spatial database with physical properties for elastic material
+db.iohandler.filename = matprops.spatialdb
+
+# Set the basis functions and quadrature:
+# 3-D simplex cell with 2nd order quadrature
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.materials.neg]
+
+# Label for material
+label = elastic material -x
+
+# Material id (itetclr) from LaGriT
+id = 2
+
+# Spatial database with physical properties for elastic material
+db.iohandler.filename = matprops.spatialdb
+
+# Set the basis functions and quadrature:
+# 3-D simplex cell with 2nd order quadrature
+quadrature = pylith.feassemble.quadrature.Quadrature3D
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = tetrahedron
+quadrature.cell.quad_order = 2
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.bc.x_pos]
+# Absorbing boundary condition on +x face of bar
+
+# Pset name from LaGriT
+label = boundary_xp
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC +x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 2-D simplex cell in 3-D space with 2nd order quadrature
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.x_neg]
+# Absorbing boundary condition on +x face of bar
+
+# Pset name from LaGriT
+label = boundary_xm
+
+# Specify label and filename for spatial database with physical properties
+db.label = Absorbing BC -x
+db.iohandler.filename = matprops.spatialdb
+
+# Specify the basis functions and quadrature:
+# 2-D simplex cell in 3-D space with 2nd order quadrature
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+[pylithapp.timedependent.bc.y_pos]
+# Dirichlet boundary condition on all vertices except fault vertices
+
+fixed_dof = [0, 2]
+
+# Pset name from LaGriT
+label = pts_notfault
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+# NOTE: It is possible to assign an ID number to a fault (e.g.,
+# 'id = 10').  Care must be taken when doing this, however, because the
+# assigned ID will become the material ID for the cohesive element.
+# This ID must not conflict with any of the material ID numbers for
+# volume elements.  The default ID for a fault is 100.  If you have a
+# fault in your mesh you must:
+# 1.  If you create your own fault ID, make sure it does not conflict
+#     with any of you material ID's.
+# 2.  If you use the default fault ID, make sure that none of your
+#     material ID's are equal to 100.
+
+[pylithapp.timedependent.interfaces.fault]
+# Material id
+id = 100
+
+# Pset name from LaGriT
+label = fault
+
+# Specify the basis functions and quadrature:
+# 2-D simplex cell in 3-D space with 2nd order quadrature
+quadrature = pylith.feassemble.quadrature.Quadrature2Din3D
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = triangle
+quadrature.cell.quad_order = 2
+
+# Set the physical properties for the fault to properly condition the
+# system of equations.
+mat_db.iohandler.filename = matprops.spatialdb
+
+[pylithapp.timedependent.interfaces.fault.eq_src.slip_function]
+
+# Database specifying the final slip.
+slip.iohandler.filename = shearwave_slip.spatialdb
+
+# Database specifying peak slip rate.
+slip_rate.iohandler.filename = shearwave_sliprate.spatialdb
+
+# Database specifying time at which slip begins at each point.
+slip_time.iohandler.filename = shearwave_sliptime.spatialdb
+
+# ----------------------------------------------------------------------
+# PETSc
+# ----------------------------------------------------------------------
+[pylithapp.petsc]
+ksp_type = gmres
+pc_type = asm
+ksp_rtol = 1.0e-8
+ksp_max_it = 50
+ksp_gmres_restart = 10
+
+#ksp_monitor = true
+#ksp_view = true
+log_summary = true
+
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# We set three output manager parameters for each type of output.
+#   skip - Number of time steps to skip between output
+#   writer.filename - Root name for VTK file
+#   writer.time_format - Format of time stamp (decimal point is removed)
+
+# Give basename for VTK domain output of solution over domain.
+[pylithapp.problem.formulation.output.output]
+skip = 3
+writer.filename = output/shearwave.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK fault output.
+[pylithapp.timedependent.interfaces.fault.output]
+skip = 3
+writer.filename = output/shearwave-fault.vtk
+writer.time_format = %05.2f
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.pos.output]
+skip = 33
+writer.filename = output/shearwave-statevars-pos.vtk
+writer.time_format = %05.2f
+
+# Averate state variables over cell
+cell_filter = pylith.meshio.CellFilterAvg
+
+# Give basename for VTK output of state variables.
+[pylithapp.timedependent.materials.neg.output]
+skip = 33
+writer.filename = output/shearwave-statevars-neg.vtk
+writer.time_format = %05.2f
+
+# Averate state variables over cell
+cell_filter = pylith.meshio.CellFilterAvg

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_slip.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_slip.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_slip.spatialdb	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+  // number of values for final slip
+  num-values = 3
+
+  // Names and units of final slip values
+  value-names =  left-lateral-slip  reverse-slip  fault-opening
+  value-units =  m  m  m
+
+  // Number of locations where physical properties are given
+  // 1 location -> data is uniform -> data dimension is 0
+  num-locs = 1
+  data-dim = 0
+
+  // Problem is in 3-D
+  space-dim = 3
+
+  // Coordinate system (3-D with coordinates in meters)
+  cs-data = cartesian {
+    to-meters = 1.0
+    space-dim = 3
+  }
+}
+// Columns are x, y, z, left-lateral-slip, reverse-slip, fault-opening
+0.0  0.0  0.0   1.0  0.0  0.0

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliprate.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliprate.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliprate.spatialdb	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,27 @@
+#SPATIAL.ascii 1
+SimpleDB {
+  // number of values for peak slip rate
+  num-values = 1
+
+  // Name and units for peak slip rate
+  value-names =  slip-rate
+  value-units =  m/s
+
+  // Number of locations where physical properties are given
+  // 1 location -> data is uniform -> data dimension is 0
+  num-locs = 1
+
+  // Problem is in 3-D
+  data-dim = 0
+
+  // Problem is in 3-D
+  space-dim = 3
+
+  // Coordinate system (3-D with coordinates in meters)
+  cs-data = cartesian {
+    to-meters = 1.0
+    space-dim = 3
+  }
+}
+// Columns are x, y, z, slip-rate
+0.0  0.0  0.0    1.0

Added: short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliptime.spatialdb
===================================================================
--- short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliptime.spatialdb	                        (rev 0)
+++ short/3D/PyLith/trunk/examples/3d/bar_tet4/shearwave_sliptime.spatialdb	2008-03-31 07:04:29 UTC (rev 11654)
@@ -0,0 +1,25 @@
+#SPATIAL.ascii 1
+SimpleDB {
+  // number of values for slip initiation time
+  num-values = 1
+
+  // Names and units of slip initiation time
+  value-names =  slip-time
+  value-units =  s
+
+  // Number of locations where physical properties are given
+  // 1 location -> data is uniform -> data dimension is 0
+  num-locs = 1
+  data-dim = 0
+
+  // Problem is in 3-D
+  space-dim = 3
+
+  // Coordinate system (3-D with coordinates in meters)
+  cs-data = cartesian {
+    to-meters = 1.0
+    space-dim = 3
+  }
+}
+// Columns are x, y, z, slip-time
+0.0 0.0 0.0   0.1



More information about the cig-commits mailing list