[cig-commits] r18365 - short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction
brad at geodynamics.org
brad at geodynamics.org
Fri May 13 18:48:31 PDT 2011
Author: brad
Date: 2011-05-13 18:48:30 -0700 (Fri, 13 May 2011)
New Revision: 18365
Added:
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_creep_slabtop.spatialdb
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip_coseismic.spatialdb
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step01.cfg
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg
Removed:
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip.spatialdb
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/postseismic.cfg
Modified:
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/mesh_tri3.exo
short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/pylithapp.cfg
Log:
Fixed bugs in setting fault slip. Added step01 and step02.
Added: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_creep_slabtop.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_creep_slabtop.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_creep_slabtop.spatialdb 2011-05-14 01:48:30 UTC (rev 18365)
@@ -0,0 +1,31 @@
+// -*- C++ -*- (syntax highlighting)
+//
+// This spatial database specifies the distribution of slip rate
+// associated with aseismic creep on the interface between the mantle
+// below the continental crust and the subducting oceanic crust.
+//
+// We specify a uniform creep rate below a depth of 75 km, tapering to
+// 0 at a depth of 60 km.
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = left-lateral-slip fault-opening
+ value-units = cm/year cm/year
+ num-locs = 4
+ data-dim = 1 // Data is specified along a line.
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // Specify coordinates in km for convenience.
+ space-dim = 2
+ } // cs-data
+} // SimpleDB
+// Columns are
+// (1) x coordinate (km)
+// (2) y coordinate (km)
+// (3) reverse-slip (cm)
+// (4) fault-opening (cm)
+0.0 99.0 0.0 0.0
+0.0 -60.0 0.0 0.0
+0.0 -75.0 -8.0 0.0
+0.0 -999.0 -8.0 0.0
Deleted: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip.spatialdb
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip.spatialdb 2011-05-14 01:25:45 UTC (rev 18364)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip.spatialdb 2011-05-14 01:48:30 UTC (rev 18365)
@@ -1,42 +0,0 @@
-// -*- C++ -*- (syntax highlighting)
-//
-// This spatial database specifies the distribution of slip on the
-// fault surface. This is a cross-section extracted from Gavin Hayes's
-// rupture model
-// (http://earthquake.usgs.gov/earthquakes/eqinthenews/2011/usc0001xgp/finite_fault.php).
-//
-#SPATIAL.ascii 1
-SimpleDB {
- num-values = 2
- value-names = left-lateral-slip fault-opening
- value-units = cm cm
- num-locs = 17
- data-dim = 1 // Data is specified along a line.
- space-dim = 2
- cs-data = cartesian {
- to-meters = 1.0e+3 // Specify coordinates in km for convenience.
- space-dim = 2
- } // cs-data
-} // SimpleDB
-// Columns are
-// (1) x coordinate (km)
-// (2) y coordinate (km)
-// (3) reverse-slip (cm)
-// (4) fault-opening (cm)
-0.0 99.0 0.0
-0.0 0.0 0.0
-0.0 -7.5475 -1143.00098
-0.0 -11.0926 -3202.27295
-0.0 -14.6378 -3050.15894
-0.0 -18.1829 -1995.05396
-0.0 -21.7280 -1167.55798
-0.0 -25.2732 -1326.13599
-0.0 -28.8183 -1260.37805
-0.0 -32.3634 -813.61987
-0.0 -35.9086 -255.54710
-0.0 -39.4537 -370.49261
-0.0 -42.9988 -394.57870
-0.0 -46.5439 -310.11920
-0.0 -50.0891 -349.24249
-0.0 -55.0 0.0
-0.0 -999.0 0.0
Copied: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip_coseismic.spatialdb (from rev 18363, short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip.spatialdb)
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip_coseismic.spatialdb (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/fault_slip_coseismic.spatialdb 2011-05-14 01:48:30 UTC (rev 18365)
@@ -0,0 +1,42 @@
+// -*- C++ -*- (syntax highlighting)
+//
+// This spatial database specifies the distribution of slip on the
+// fault surface. This is a cross-section extracted from Gavin Hayes's
+// rupture model
+// (http://earthquake.usgs.gov/earthquakes/eqinthenews/2011/usc0001xgp/finite_fault.php).
+//
+#SPATIAL.ascii 1
+SimpleDB {
+ num-values = 2
+ value-names = left-lateral-slip fault-opening
+ value-units = cm cm
+ num-locs = 17
+ data-dim = 1 // Data is specified along a line.
+ space-dim = 2
+ cs-data = cartesian {
+ to-meters = 1.0e+3 // Specify coordinates in km for convenience.
+ space-dim = 2
+ } // cs-data
+} // SimpleDB
+// Columns are
+// (1) x coordinate (km)
+// (2) y coordinate (km)
+// (3) reverse-slip (cm)
+// (4) fault-opening (cm)
+0.0 99.0 0.0 0.0
+0.0 0.0 0.0 0.0
+0.0 -7.5475 -1143.00098 0.0
+0.0 -11.0926 -3202.27295 0.0
+0.0 -14.6378 -3050.15894 0.0
+0.0 -18.1829 -1995.05396 0.0
+0.0 -21.7280 -1167.55798 0.0
+0.0 -25.2732 -1326.13599 0.0
+0.0 -28.8183 -1260.37805 0.0
+0.0 -32.3634 -813.61987 0.0
+0.0 -35.9086 -255.54710 0.0
+0.0 -39.4537 -370.49261 0.0
+0.0 -42.9988 -394.57870 0.0
+0.0 -46.5439 -310.11920 0.0
+0.0 -50.0891 -349.24249 0.0
+0.0 -55.0 0.0 0.0
+0.0 -999.0 0.0 0.0
Modified: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/mesh_tri3.exo
===================================================================
(Binary files differ)
Deleted: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/postseismic.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/postseismic.cfg 2011-05-14 01:25:45 UTC (rev 18364)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/postseismic.cfg 2011-05-14 01:48:30 UTC (rev 18365)
@@ -1,155 +0,0 @@
-# -*- Python -*-
-[pylithapp]
-
-# ----------------------------------------------------------------------
-# PROBLEM DESCRIPTION
-# ----------------------------------------------------------------------
-#
-#
-
-# ----------------------------------------------------------------------
-# RUNNING THE SIMULATON
-# ----------------------------------------------------------------------
-#
-# This is not a self-contained simulation configuration file. This
-# file specifies only the boundary conditions and earthquake
-# parameters for the simulation. The general quasi-static and mesh
-# parameters are specificed in the pylithapp.cfg file which PyLith
-# reads by default.
-#
-# To run the simulation:
-# pylith postseismic.cfg
-#
-# Output will be directed to directory output.
-
-# ----------------------------------------------------------------------
-# boundary conditions
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-# Set bc to an array of 3 boundary conditions:
-# 'bndry_east'
-# 'bndry_west'
-# 'bndry_bot'
-bc = [boundary_east,boundary_west,boundary_bottom]
-
-# For all boundaries, we fix the displacement normal to the boundary
-# (roller boundary condition) by retaining the default ZeroDispDB,
-# which specifies a zero value.
-#
-# The label corresponds to the name of the nodeset in CUBIT.
-
-# East boundary
-[pylithapp.timedependent.bc.boundary_east]
-bc_dof = [0]
-label = bndry_east
-db_initial.label = Dirichlet BC on east boundary
-
-# West boundary
-[pylithapp.timedependent.bc.boundary_west]
-bc_dof = [0]
-label = bndry_west
-db_initial.label = Dirichlet BC on west boundary
-
-# Bottom boundary
-[pylithapp.timedependent.bc.boundary_bottom]
-bc_dof = [1]
-label = bndry_bot
-db_initial.label = Dirichlet BC on bottom boundary
-
-# ----------------------------------------------------------------------
-# faults
-# ----------------------------------------------------------------------
-[pylithapp.timedependent]
-interfaces = [fault]
-
-# Set the type of fault interface condition.
-[pylithapp.timedependent.interfaces]
-fault = pylith.faults.FaultCohesiveKin
-
-# Set the parameters for the fault interface condition.
-[pylithapp.timedependent.interfaces.fault]
-# The label corresponds to the name of the nodeset in CUBIT.
-label = fault
-
-# We must define the quadrature information for fault cells.
-# The fault cells are 2D (surface).
-quadrature.cell = pylith.feassemble.FIATSimplex
-quadrature.cell.shape = line
-
-# The slip time and final slip are defined in spatial databases.
-[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
-slip.iohandler.filename = fault_slip.spatialdb
-slip.query_type = linear
-slip.label = Final slip
-
-slip_time = spatialdata.spatialdb.UniformDB
-slip_time.label = Slip time
-slip_time.values = [slip-time]
-slip_time.data = [0.0*year]
-
-# ----------------------------------------------------------------------
-# output
-# ----------------------------------------------------------------------
-[pylithapp.timedependent.implicit]
-# Set the output to an array of 2 output managers.
-# We will output the solution over the domain and the ground surface.
-output = [domain,subdomain]
-
-# Set subdomain component to OutputSolnSubset (subset of domain).
-output.subdomain = pylith.meshio.OutputSolnSubset
-
-# Domain
-[pylithapp.problem.formulation.output.domain]
-output_freq = time_step
-time_step = 4.99999*year
-writer.filename = output/postseismic.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Ground surface
-[pylithapp.problem.formulation.output.subdomain]
-label = groundsurf ; Name of CUBIT nodeset for ground surface.
-writer.filename = output/postseismic-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Fault
-[pylithapp.problem.interfaces.fault.output]
-writer.filename = output/postseismic-fault.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# Materials
-[pylithapp.timedependent.materials.continent_crust.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
-output_freq = time_step
-time_step = 9.99999*year
-writer.filename = output/postseismic-concrust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.continent_mantle.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
-output_freq = time_step
-time_step = 9.99999*year
-writer.filename = output/postseismic-conmantle.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.ocean_crust.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
-output_freq = time_step
-time_step = 9.99999*year
-writer.filename = output/postseismic-oceancrust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-[pylithapp.timedependent.materials.ocean_mantle.output]
-cell_filter = pylith.meshio.CellFilterAvgMesh
-output_freq = time_step
-time_step = 9.99999*year
-writer.filename = output/postseismic-oceanmantle.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
-
-# End of file
\ No newline at end of file
Modified: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/pylithapp.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/pylithapp.cfg 2011-05-14 01:25:45 UTC (rev 18364)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/pylithapp.cfg 2011-05-14 01:48:30 UTC (rev 18365)
@@ -121,6 +121,62 @@
# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Names of output files are set in stepXX.cfg. We consolidate all of the
+# output settings that are common to all of the simulations here.
+
+[pylithapp.timedependent.implicit]
+# Set the output to an array of 2 output managers.
+# We will output the solution over the domain and the ground surface.
+output = [domain,subdomain]
+
+# Set subdomain component to OutputSolnSubset (subset of domain).
+output.subdomain = pylith.meshio.OutputSolnSubset
+
+# Domain
+[pylithapp.problem.formulation.output.domain]
+output_freq = time_step
+time_step = 4.99999*year
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Ground surface
+[pylithapp.problem.formulation.output.subdomain]
+label = groundsurf ; Name of CUBIT nodeset for ground surface.
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# Materials
+[pylithapp.timedependent.materials.continent_crust.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+output_freq = time_step
+time_step = 9.99999*year
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.continent_mantle.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+output_freq = time_step
+time_step = 9.99999*year
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.ocean_crust.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+output_freq = time_step
+time_step = 9.99999*year
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+[pylithapp.timedependent.materials.ocean_mantle.output]
+cell_filter = pylith.meshio.CellFilterAvgMesh
+output_freq = time_step
+time_step = 9.99999*year
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+
+# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
# Set the solver options.
Copied: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step01.cfg (from rev 18362, short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/postseismic.cfg)
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step01.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step01.cfg 2011-05-14 01:48:30 UTC (rev 18365)
@@ -0,0 +1,141 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+#
+#
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+#
+# This is not a self-contained simulation configuration file. This
+# file specifies only the boundary conditions and earthquake
+# parameters for the simulation. The general quasi-static and mesh
+# parameters are specificed in the pylithapp.cfg file which PyLith
+# reads by default.
+#
+# To run the simulation:
+# pylith step01.cfg
+#
+# Output will be directed to directory output.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+total_time = 0.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 5 boundary conditions:
+# 'boundary_east_crust'
+# 'boundary_east_mantle'
+# 'boundary_west'
+# 'boundary_bottom_slab'
+# 'boundary_bottom_mantle'
+bc = [boundary_east_crust,boundary_east_mantle,boundary_west,boundary_bottom_slab,boundary_bottom_mantle]
+
+# For all boundaries, we fix the displacement normal to the boundary
+# (roller boundary condition) by retaining the default ZeroDispDB,
+# which specifies a zero value.
+#
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# East boundaries
+[pylithapp.timedependent.bc.boundary_east_crust]
+bc_dof = [0]
+label = bndry_east_crust
+db_initial.label = Dirichlet BC on east boundary (crust)
+
+[pylithapp.timedependent.bc.boundary_east_mantle]
+bc_dof = [0]
+label = bndry_east_mantle
+db_initial.label = Dirichlet BC on east boundary (mantle)
+
+# West boundary
+[pylithapp.timedependent.bc.boundary_west]
+bc_dof = [0]
+label = bndry_west
+db_initial.label = Dirichlet BC on west boundary
+
+# Bottom boundaries
+[pylithapp.timedependent.bc.boundary_bottom_slab]
+bc_dof = [1]
+label = bndry_bot_slab
+db_initial.label = Dirichlet BC on bottom boundary (slab)
+
+[pylithapp.timedependent.bc.boundary_bottom_mantle]
+bc_dof = [1]
+label = bndry_bot_mantle
+db_initial.label = Dirichlet BC on bottom boundary (mantle)
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+interfaces = [fault]
+
+# Set the type of fault interface condition.
+[pylithapp.timedependent.interfaces]
+fault = pylith.faults.FaultCohesiveKin
+
+# Set the parameters for the fault interface condition.
+[pylithapp.timedependent.interfaces.fault]
+# The label corresponds to the name of the nodeset in CUBIT.
+label = fault_coseismic
+
+# We must define the quadrature information for fault cells.
+# The fault cells are 1D (line).
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
+# The slip time and final slip are defined in spatial databases.
+[pylithapp.timedependent.interfaces.fault.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = fault_slip_coseismic.spatialdb
+slip.query_type = linear
+slip.label = Final slip
+
+# Slip time is uniform, so use UniformDB for convenience
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.label = Slip time
+slip_time.values = [slip-time]
+slip_time.data = [0.0*year]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Domain
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/step01.vtk
+
+# Ground surface
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/step01-groundsurf.vtk
+
+# Fault
+[pylithapp.problem.interfaces.fault.output]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+writer.filename = output/step01-fault.vtk
+
+# Materials
+[pylithapp.timedependent.materials.continent_crust.output]
+writer.filename = output/step01-concrust.vtk
+
+[pylithapp.timedependent.materials.continent_mantle.output]
+writer.filename = output/step01-conmantle.vtk
+
+[pylithapp.timedependent.materials.ocean_crust.output]
+writer.filename = output/step01-oceancrust.vtk
+
+[pylithapp.timedependent.materials.ocean_mantle.output]
+writer.filename = output/step01-oceanmantle.vtk
+
+
+# End of file
Added: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg 2011-05-14 01:48:30 UTC (rev 18365)
@@ -0,0 +1,169 @@
+# -*- Python -*-
+[pylithapp]
+
+# ----------------------------------------------------------------------
+# PROBLEM DESCRIPTION
+# ----------------------------------------------------------------------
+#
+#
+
+# ----------------------------------------------------------------------
+# RUNNING THE SIMULATON
+# ----------------------------------------------------------------------
+#
+# This is not a self-contained simulation configuration file. This
+# file specifies only the boundary conditions and earthquake
+# parameters for the simulation. The general quasi-static and mesh
+# parameters are specificed in the pylithapp.cfg file which PyLith
+# reads by default.
+#
+# To run the simulation:
+# pylith step02.cfg
+#
+# Output will be directed to directory output.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+total_time = 10.0*year
+dt = 1.0*year
+
+# ----------------------------------------------------------------------
+# boundary conditions
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+# Set bc to an array of 3 boundary conditions:
+# 'boundary_east_mantle'
+# 'boundary_west'
+# 'boundary_bottom_mantle'
+bc = [boundary_east_mantle,boundary_west,boundary_bottom_mantle]
+
+# For all boundaries, we fix the displacement normal to the boundary
+# (roller boundary condition) by retaining the default ZeroDispDB,
+# which specifies a zero value.
+#
+# The label corresponds to the name of the nodeset in CUBIT.
+
+# East boundary (mantle)
+[pylithapp.timedependent.bc.boundary_east_mantle]
+bc_dof = [0]
+label = bndry_east_mantle
+db_initial.label = Dirichlet BC on east boundary (mantle)
+
+# West boundary
+[pylithapp.timedependent.bc.boundary_west]
+bc_dof = [0]
+label = bndry_west
+db_initial.label = Dirichlet BC on west boundary
+
+# Bottom boundary (mantle)
+[pylithapp.timedependent.bc.boundary_bottom_mantle]
+bc_dof = [1]
+label = bndry_bot_mantle
+db_initial.label = Dirichlet BC on bottom boundary (mantle)
+
+# ----------------------------------------------------------------------
+# faults
+# ----------------------------------------------------------------------
+[pylithapp.timedependent]
+interfaces = [fault_slabtop,fault_slabbot]
+
+# Set the type of fault interface condition.
+[pylithapp.timedependent.interfaces]
+fault_slabtop = pylith.faults.FaultCohesiveKin
+fault_slabbot = pylith.faults.FaultCohesiveKin
+
+# Slab top --------------------
+[pylithapp.timedependent.interfaces.fault_slabtop]
+# The label corresponds to the name of the nodeset in CUBIT.
+label = fault_slabtop
+id = 100
+
+# We must define the quadrature information for fault cells.
+# The fault cells are 1D (line).
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
+# Switch to constant slip rate time function.
+[pylithapp.timedependent.interfaces.fault_slabtop.eq_srcs.rupture]
+slip_function = pylith.faults.ConstRateSlipFn
+
+# The slip time and final slip are defined in spatial databases.
+[pylithapp.timedependent.interfaces.fault_slabtop.eq_srcs.rupture.slip_function]
+slip_rate.iohandler.filename = fault_creep_slabtop.spatialdb
+slip_rate.query_type = linear
+slip_rate.label = Final slip
+
+# Slip time is uniform, so use UniformDB for convenience
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.label = Slip time
+slip_time.values = [slip-time]
+slip_time.data = [0.0*year]
+
+
+# Slab bottom --------------------
+[pylithapp.timedependent.interfaces.fault_slabbot]
+# The label corresponds to the name of the nodeset in CUBIT.
+label = fault_slabbot
+id = 101
+
+# We must define the quadrature information for fault cells.
+# The fault cells are 1D (line).
+quadrature.cell = pylith.feassemble.FIATSimplex
+quadrature.cell.shape = line
+
+# Switch to constant slip rate time function.
+[pylithapp.timedependent.interfaces.fault_slabbot.eq_srcs.rupture]
+slip_function = pylith.faults.ConstRateSlipFn
+
+# The creep rate and slip time are uniform, so use UniformDB for simplicity.
+[pylithapp.timedependent.interfaces.fault_slabbot.eq_srcs.rupture.slip_function]
+slip_rate = spatialdata.spatialdb.UniformDB
+slip_rate.label = Slip rate
+slip_rate.values = [left-lateral-slip, fault-opening]
+slip_rate.data = [8.0*cm/year, 0.0*cm/year]
+
+# Slip time is uniform, so use UniformDB for convenience
+slip_time = spatialdata.spatialdb.UniformDB
+slip_time.label = Slip time
+slip_time.values = [slip-time]
+slip_time.data = [0.0*year]
+
+# ----------------------------------------------------------------------
+# output
+# ----------------------------------------------------------------------
+# Domain
+[pylithapp.problem.formulation.output.domain]
+writer.filename = output/step02.vtk
+
+# Ground surface
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/step02-groundsurf.vtk
+
+# Faults
+[pylithapp.problem.interfaces.fault_slabtop.output]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+writer.filename = output/step02-fault-slabtop.vtk
+
+[pylithapp.problem.interfaces.fault_slabbot.output]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+writer.filename = output/step02-fault-slabbot.vtk
+
+# Materials
+[pylithapp.timedependent.materials.continent_crust.output]
+writer.filename = output/step02-concrust.vtk
+
+[pylithapp.timedependent.materials.continent_mantle.output]
+writer.filename = output/step02-conmantle.vtk
+
+[pylithapp.timedependent.materials.ocean_crust.output]
+writer.filename = output/step02-oceancrust.vtk
+
+[pylithapp.timedependent.materials.ocean_mantle.output]
+writer.filename = output/step02-oceanmantle.vtk
+
+
+# End of file
More information about the CIG-COMMITS
mailing list