[cig-commits] r18367 - short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction

[brad at geodynamics.org]

Author: brad
Date: 2011-05-13 22:05:40 -0700 (Fri, 13 May 2011)
New Revision: 18367

Added:
   short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step03.cfg
Modified:
   short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg
Log:
Added step03.

Modified: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg	2011-05-14 03:08:08 UTC (rev 18366)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step02.cfg	2011-05-14 05:05:40 UTC (rev 18367)
@@ -26,7 +26,7 @@
 # problem
 # ----------------------------------------------------------------------
 [pylithapp.timedependent.formulation.time_step]
-total_time = 10.0*year
+total_time = 50.0*year
 dt = 1.0*year
 
 # ----------------------------------------------------------------------

Added: short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step03.cfg
===================================================================
--- short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step03.cfg	                        (rev 0)
+++ short/3D/PyLith/branches/v1.5-stable/examples/2d/subduction/step03.cfg	2011-05-14 05:05:40 UTC (rev 18367)
@@ -0,0 +1,186 @@
+# -*- 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 step03.cfg
+#
+# Output will be directed to directory output.
+
+# ----------------------------------------------------------------------
+# problem
+# ----------------------------------------------------------------------
+[pylithapp.timedependent.formulation.time_step]
+total_time = 100.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
+
+# Set earthquake sources to an array consisting of creep and 3 ruptures.
+eq_srcs = [creep,rupture]
+eq_srcs.creep.origin_time = 00.0*year
+eq_srcs.rupture.origin_time = 25.0*year
+
+# Coseismic rupture ---------
+[pylithapp.timedependent.interfaces.fault_slabtop.eq_srcs.rupture.slip_function]
+slip.iohandler.filename = fault_slip_coseismic.spatialdb
+slip.query_type = linear
+slip.label = Final slip
+
+# Slip time (relative to origin 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]
+
+# Creep ---------
+[pylithapp.timedependent.interfaces.fault_slabtop.eq_srcs.creep]
+slip_function = pylith.faults.ConstRateSlipFn
+
+# The slip time and final slip are defined in spatial databases.
+[pylithapp.timedependent.interfaces.fault_slabtop.eq_srcs.creep.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/step03.vtk
+
+# Ground surface
+[pylithapp.problem.formulation.output.subdomain]
+writer.filename = output/step03-groundsurf.vtk
+
+# Faults
+[pylithapp.problem.interfaces.fault_slabtop.output]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+writer.filename = output/step03-fault-slabtop.vtk
+
+[pylithapp.problem.interfaces.fault_slabbot.output]
+writer.time_format = %04.0f
+writer.time_constant = 1.0*year
+writer.filename = output/step03-fault-slabbot.vtk
+
+# Materials
+[pylithapp.timedependent.materials.continent_crust.output]
+writer.filename = output/step03-concrust.vtk
+
+[pylithapp.timedependent.materials.continent_mantle.output]
+writer.filename = output/step03-conmantle.vtk
+
+[pylithapp.timedependent.materials.ocean_crust.output]
+writer.filename = output/step03-oceancrust.vtk
+
+[pylithapp.timedependent.materials.ocean_mantle.output]
+writer.filename = output/step03-oceanmantle.vtk
+
+
+# End of file