[cig-commits] r11362 - in short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog: figs meshes parameters

brad at geodynamics.org brad at geodynamics.org
Thu Mar 6 17:20:22 PST 2008


Author: brad
Date: 2008-03-06 17:20:21 -0800 (Thu, 06 Mar 2008)
New Revision: 11362

Added:
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/figs/plot_mesh.py
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0250m.lagrit
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0500m.lagrit
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_1000m.lagrit
Removed:
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/debug.cfg
Modified:
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/createbc.jou
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/reverseslip_hex8_0250m.exo.gz
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/reverseslip_hex8_0500m.exo.gz
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/reverseslip_hex8_1000m.exo.gz
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg
   short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg
Log:
Added ground surface to output. Setup preliminary LaGriT mesh (want to try extruded mesh to see if mesh is better).

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/figs/plot_mesh.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/figs/plot_mesh.py	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/figs/plot_mesh.py	2008-03-07 01:20:21 UTC (rev 11362)
@@ -0,0 +1,80 @@
+#!/usr/bin/env python
+#
+# ======================================================================
+#
+#                           Brad T. Aagaard
+#                        U.S. Geological Survey
+#
+# {LicenseText}
+#
+# ======================================================================
+#
+
+from plot_geometry import PlotScene
+
+shape = "tet4"
+res = 1000
+
+class PlotMesh(PlotScene):
+
+  # PUBLIC METHODS /////////////////////////////////////////////////////
+
+  def run(self):    
+    from enthought.mayavi.sources.vtk_data_source import VTKDataSource
+    from enthought.mayavi.modules.surface import Surface
+
+    self._setupScene(showFault=False, showMaterials=False)
+    mesh = self._readMesh()
+
+    script = self.script
+    script.add_source(VTKDataSource(data=mesh))
+    script.engine.current_object.name = "Mesh"
+
+    surf = Surface()
+    script.add_module(surf)
+    
+    surf = Surface()
+    script.add_module(surf)
+    
+    import vtk_geometry
+    vtk_geometry.setCamera(script.engine.current_scene.scene.camera)
+    return
+
+
+  # PRIVATE METHODS ////////////////////////////////////////////////////
+
+  def _readMesh(self):
+    from enthought.tvtk.api import tvtk
+    import tables
+    import numpy
+
+    filename = "../meshes/reverseslip_%s_%04dm.h5" % (shape, res)
+    h5 = tables.openFile(filename, 'r')
+
+    cells = h5.root.topology.cells[:]
+    (ncells, ncorners) = cells.shape
+    vertices = h5.root.geometry.vertices[:] / 1e+3
+    (nvertices, spaceDim) = vertices.shape
+    h5.close()
+    
+    if shape == "tet4":
+        assert(spaceDim == 3)
+        assert(ncorners == 4)
+        cellType = tvtk.Tetra().cell_type
+    else:
+        raise ValueError("Unknown shape '%s'." % shape)
+
+    data = tvtk.UnstructuredGrid()
+    data.points = vertices
+    data.set_cells(cellType, cells)
+    return data
+
+
+# ----------------------------------------------------------------------
+if __name__ == "__main__":
+  app = PlotMesh()
+  app.main()
+
+
+# End of file
+


Property changes on: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/figs/plot_mesh.py
___________________________________________________________________
Name: svn:executable
   + *

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/createbc.jou	2008-03-07 01:05:39 UTC (rev 11361)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/createbc.jou	2008-03-07 01:20:21 UTC (rev 11362)
@@ -59,3 +59,12 @@
 group "symmetry_disp" remove node in fault
 nodeset 15 group symmetry_disp
 nodeset 15 name "symmetry_disp"
+
+# ----------------------------------------------------------------------
+# Create nodeset for ground surface
+# ----------------------------------------------------------------------
+group "ground_surf" add node in surface 20
+group "ground_surf" add node in surface 27
+nodeset 16 group ground_surf
+nodeset 16 name "ground_surf"
+

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit	2008-03-07 01:20:21 UTC (rev 11362)
@@ -0,0 +1,205 @@
+* -*- LaGriT -*-
+*
+* ----------------------------------------------------------------------
+* PARAMETERS CONTROLLING MESH
+* ----------------------------------------------------------------------
+*
+* Boundaries of domain
+define / domain_xmin /   0.0e+3
+define / domain_xmax /  24.0e+3
+define / domain_ymin /   0.0e+3
+define / domain_ymax /  24.0e+3
+define / domain_zmin / -24.0e+3
+define / domain_zmax /   0.0e+3
+*
+* ----------------------------------------------------------------------
+* CREATE GEOMETRY
+* ----------------------------------------------------------------------
+*
+* Create a mesh object and name it box
+cmo / create / box /
+*
+* Create boundary surfaces
+surface / surf_xneg / reflect / plane / &
+   domain_xmin   0.0e+0   0.0e+0 / &
+   domain_xmin   0.0e+0  10.0e+3 / &
+   domain_xmin  10.0e+3   0.0e+0
+surface / surf_xpos / reflect / plane / &
+    domain_xmax   0.0e+0   0.0e+0 / &
+    domain_xmax  10.0e+3   0.0e+0 / &
+    domain_xmax   0.0e+0  10.0e+3
+surface / surf_yneg / reflect / plane / &
+     0.0e+0  domain_ymin   0.0e+0 / &
+     0.0e+0  domain_ymin  10.0e+3 / &
+   -10.0e+3  domain_ymin   0.0e+0
+surface / surf_ypos / reflect / plane / &
+     0.0e+0   domain_ymax   0.0e+0 / &
+   -10.0e+3   domain_ymax   0.0e+0 / &
+     0.0e+0   domain_ymax  10.0e+3
+surface / surf_zneg / reflect / plane / &
+     0.0e+0    0.0e+0  domain_zmin / &
+   -10.0e+3    0.0e+0  domain_zmin / &
+     0.0e+0   10.0e+3  domain_zmin
+surface / surf_zpos / reflect / plane / &
+     0.0e+0    0.0e+0    domain_zmax / &
+     0.0e+0   10.0e+3    domain_zmax / &
+   -10.0e+3    0.0e+0    domain_zmax
+*
+* Create fault plane
+surface / surf_fault / intrface / plane / &
+  4.0e+3        domain_ymax   0.0e+0 / &
+  4.0e+3        0.0e+0        0.0e+0 / &
+  domain_xmax   0.0e+0      -20.0e+3 /
+* Create material interface
+surface / surf_mantle / intrface / plane / &
+  0.0     0.0     -12.0e+3 / &
+  1.0e+3  0.0     -12.0e+3 / &
+  0.0     1.0e+3  -12.0e+3 /
+
+* ----------------------------------------------------------------------
+* Create regions
+* ----------------------------------------------------------------------
+
+** Elastic on hanging wall side of the fault
+region / r_elastic_hang / &
+  le surf_xpos and ge surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and le surf_zpos and ge surf_mantle /
+
+** Elastic on footwell side of the fault
+region / r_elastic_foot / &
+  le surf_xneg and lt surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and le surf_zpos and ge surf_mantle /
+
+** Viscoelastic on hanging wall side of the fault
+region / r_visco_hang / &
+  le surf_xpos and ge surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and lt surf_mantle and le surf_zneg /
+
+** Viscoelastic on footwall side of the fault
+region / r_visco_foot / &
+  le surf_xneg and lt surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and lt surf_mantle and le surf_zneg /
+
+* ----------------------------------------------------------------------
+* Create material regions
+* ----------------------------------------------------------------------
+* Material regions do not include internal interfaces (regions do)
+*
+
+* Material 1
+mregion / mr_elastic_hang / &
+  le surf_xpos and gt surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and le surf_zpos and gt surf_mantle /
+
+* Material 2
+mregion / mr_elastic_foot / &
+  le surf_xneg and lt surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and le surf_zpos and gt surf_mantle /
+
+* Material 3
+mregion / mr_visco_hang / &
+  le surf_xpos and gt surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and lt surf_mantle and le surf_zneg /
+
+* Material 4
+mregion / mr_visco_foot / &
+  le surf_xneg and lt surf_fault  &
+  and le surf_ypos and le surf_yneg  &
+  and lt surf_mantle and le surf_zneg /
+*
+* ----------------------------------------------------------------------
+* CREATE MESH
+* ----------------------------------------------------------------------
+*
+createpts / xyz / nx,ny,1 / &
+  domain_xmin  domain_ymin  domain_zmax / &
+  domain_xmax  domain_ymax  domain_zmax / &
+  1,1,0 /
+pset / rayend / seq / 0,0,0
+*
+
+loop / foreach / iregion / &
+  r_elastic_hang r_elastic_foot r_visco_hang r_visco_foot &
+  / loop_end / &
+  regnpts / iregion / dx / pset,get,rayend / xyz / &
+    0.0e+3  0.0e+3  domain_zmin / &
+    0.0e+3  1.0e+3  domain_zmin / &
+    1.0e+3  0.0e+3  domain_zmin / &
+    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 / tmp / surface / surf_fault /
+pset / fault / geom / xyz / pset,get,tmp / &
+   domain_xmin -0.001e+3 -16.001e+3 / &
+   domain_xmax 16.001e+3   0.001e+3 /
+pset / tmp / delete
+*
+* Identify nodes on boundaries
+pset / boundary_xneg / surface / surf_xneg /
+pset / boundary_xpos / surface / surf_xpos /
+
+pset / tmp / surface / surf_ypos /
+pset / boundary_ypos / not / tmp, boundary_xneg, boundary_xpos
+pset / tmp / delete
+
+pset / tmp / surface / surf_zneg /
+pset / boundary_zneg / not / tmp, &
+  boundary_xneg, boundary_xpos, boundary_ypos
+pset / tmp / delete
+
+pset / tmp / surface / surf_yneg /
+pset / boundary_yneg / not / tmp, boundary_xneg, boundary_xpos, &
+  boundary_zneg, fault
+pset / tmp / delete
+
+pset / boundary_zpos / surface / surf_zpos /
+
+* ----------------------------------------------------------------------
+finish
+
+* End of file

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0250m.lagrit
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0250m.lagrit	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0250m.lagrit	2008-03-07 01:20:21 UTC (rev 11362)
@@ -0,0 +1,28 @@
+* -*- LaGriT -*-
+*
+* ----------------------------------------------------------------------
+* PARAMETERS CONTROLLING MESH
+* ----------------------------------------------------------------------
+*
+* Mesh size and number of points along x & y
+define / dx / 250.0
+define / nx / 97
+define / ny / 97
+
+* ----------------------------------------------------------------------
+* GENERATE MESH AND MARK BOUNDARIES
+* ----------------------------------------------------------------------
+input / mesh_tet4.lagrit
+
+* ----------------------------------------------------------------------
+* WRITE MESH
+* ----------------------------------------------------------------------
+cmo / setatt / box / ipolydat / no
+dump / gmv / reverseslip_tet4_0250m.gmv / / binary
+pset / -all- / write / reverseslip_tet4_0250m.pset / ascii
+
+* ----------------------------------------------------------------------
+finish
+
+* ======================================================================
+* End of file

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0500m.lagrit
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0500m.lagrit	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_0500m.lagrit	2008-03-07 01:20:21 UTC (rev 11362)
@@ -0,0 +1,28 @@
+* -*- LaGriT -*-
+*
+* ----------------------------------------------------------------------
+* PARAMETERS CONTROLLING MESH
+* ----------------------------------------------------------------------
+*
+* Mesh size and number of points along x & y
+define / dx / 500.0
+define / nx / 49
+define / ny / 49
+
+* ----------------------------------------------------------------------
+* GENERATE MESH AND MARK BOUNDARIES
+* ----------------------------------------------------------------------
+input / mesh_tet4.lagrit
+
+* ----------------------------------------------------------------------
+* WRITE MESH
+* ----------------------------------------------------------------------
+cmo / setatt / box / ipolydat / no
+dump / gmv / reverseslip_tet4_0500m.gmv / / binary
+pset / -all- / write / reverseslip_tet4_0500m.pset / ascii
+
+* ----------------------------------------------------------------------
+finish
+
+* ======================================================================
+* End of file

Added: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_1000m.lagrit
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_1000m.lagrit	                        (rev 0)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4_1000m.lagrit	2008-03-07 01:20:21 UTC (rev 11362)
@@ -0,0 +1,28 @@
+* -*- LaGriT -*-
+*
+* ----------------------------------------------------------------------
+* PARAMETERS CONTROLLING MESH
+* ----------------------------------------------------------------------
+*
+* Mesh size and number of points along x & y
+define / dx / 1000.0
+define / nx / 25
+define / ny / 25
+
+* ----------------------------------------------------------------------
+* GENERATE MESH AND MARK BOUNDARIES
+* ----------------------------------------------------------------------
+input / mesh_tet4.lagrit
+
+* ----------------------------------------------------------------------
+* WRITE MESH
+* ----------------------------------------------------------------------
+cmo / setatt / box / ipolydat / no
+dump / gmv / reverseslip_tet4_1000m.gmv / / binary
+pset / -all- / write / reverseslip_tet4_1000m.pset / ascii
+
+* ----------------------------------------------------------------------
+finish
+
+* ======================================================================
+* End of file

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/reverseslip_hex8_0250m.exo.gz
===================================================================
(Binary files differ)

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/reverseslip_hex8_0500m.exo.gz
===================================================================
(Binary files differ)

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/reverseslip_hex8_1000m.exo.gz
===================================================================
(Binary files differ)

Deleted: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/debug.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/debug.cfg	2008-03-07 01:05:39 UTC (rev 11361)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/debug.cfg	2008-03-07 01:20:21 UTC (rev 11362)
@@ -1,9 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-mesh_generator.distributor.debug = 1
-
-[pylithapp.journal.debug]
-pylithapp = 1
-problem = 1
-implicit = 1
-

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg	2008-03-07 01:05:39 UTC (rev 11361)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/pylithapp.cfg	2008-03-07 01:20:21 UTC (rev 11362)
@@ -17,6 +17,11 @@
 elasticityimplicit = 1
 faultcohesivekin = 1
 
+[pylithapp.journal.debug]
+pylithapp = 1
+problem = 1
+implicit = 1
+
 # ----------------------------------------------------------------------
 # mesh_generator
 # ----------------------------------------------------------------------
@@ -34,9 +39,18 @@
 default_dt = 1.0*s
 dimension = 3
 formulation = pylith.problems.Implicit
+
+# Array of DirichletPoints (default) boundary conditions
 bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
+
+# Fault
 interfaces = [fault]
 
+# Output solution over domain + ground surface (subdomain)
+[pylithapp.timedependent.implicit]
+output = [domain,subdomain]
+output.subdomain = pylith.meshio.OutputSolnSubset
+
 # ----------------------------------------------------------------------
 # boundary conditions
 # ----------------------------------------------------------------------
@@ -114,3 +128,6 @@
 ksp_monitor = true
 ksp_view = true
 log_summary = true
+
+
+# End of file

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg	2008-03-07 01:05:39 UTC (rev 11361)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_cubit.cfg	2008-03-07 01:20:21 UTC (rev 11362)
@@ -60,3 +60,11 @@
 label = 10
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.subdomain]
+label = 16
+
+
+# End of file

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg	2008-03-07 01:05:39 UTC (rev 11361)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_hex8_1000m.cfg	2008-03-07 01:20:21 UTC (rev 11362)
@@ -10,9 +10,12 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-[pylithapp.problem.formulation.output.output.writer]
+[pylithapp.problem.formulation.output.domain.writer]
 filename = ../results/reverseslip_hex8_1000m.vtk
 
+[pylithapp.problem.formulation.output.subdomain.writer]
+filename = ../results/reverseslip_hex8_1000m-groundsurf.vtk
+
 [pylithapp.timedependent.interfaces.fault.output.writer]
 filename = ../results/reverseslip_hex8_1000m-fault.vtk
 

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg	2008-03-07 01:05:39 UTC (rev 11361)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_lagrit.cfg	2008-03-07 01:20:21 UTC (rev 11362)
@@ -77,3 +77,12 @@
 label = fault
 quadrature.cell = pylith.feassemble.FIATSimplex
 quadrature.cell.shape = triangle
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+[pylithapp.problem.formulation.output.subdomain]
+label = boundary_zpos
+
+
+# End of file

Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg	2008-03-07 01:05:39 UTC (rev 11361)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/parameters/reverseslip_tet4_1000m.cfg	2008-03-07 01:20:21 UTC (rev 11362)
@@ -5,11 +5,33 @@
 # mesh_generator
 # ----------------------------------------------------------------------
 [pylithapp.mesh_generator.importer]
-filename_gmv = ../meshes/strikeslip_tet4_1000m.gmv
-filename_pset = ../meshes/strikeslip_tet4_1000m.pset
+filename_gmv = ../meshes/reverseslip_tet4_1000m.gmv
+filename_pset = ../meshes/reverseslip_tet4_1000m.pset
+#record_header_32bit = False
 
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-[pylithapp.problem.formulation.output.output.writer]
-filename = ../results/strikeslip_tet4_1000m.vtk
+[pylithapp.problem.formulation.output.domain.writer]
+filename = ../results/reverseslip_tet4_1000m.vtk
+
+[pylithapp.problem.formulation.output.subdomain.writer]
+filename = ../results/reverseslip_tet4_1000m-groundsurf.vtk
+
+[pylithapp.timedependent.interfaces.fault.output.writer]
+filename = ../results/reverseslip_tet4_1000m-fault.vtk
+
+[pylithapp.timedependent.materials.elastic_hang.output.writer]
+filename = ../results/reverseslip_tet4_1000m-elastic-hang.vtk
+
+[pylithapp.timedependent.materials.elastic_foot.output.writer]
+filename = ../results/reverseslip_tet4_1000m-elastic-foot.vtk
+
+[pylithapp.timedependent.materials.viscoelastic_hang.output.writer]
+filename = ../results/reverseslip_tet4_1000m-viscoelastic-hang.vtk
+
+[pylithapp.timedependent.materials.viscoelastic_foot.output.writer]
+filename = ../results/reverseslip_tet4_1000m-viscoelastic-foot.vtk
+
+
+# End of file



More information about the cig-commits mailing list