[cig-commits] r16423 - short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4

brad at geodynamics.org brad at geodynamics.org
Mon Mar 15 16:37:21 PDT 2010


Author: brad
Date: 2010-03-15 16:37:21 -0700 (Mon, 15 Mar 2010)
New Revision: 16423

Modified:
   short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg
Log:
Fixed errors in settings for lumped, dynamic example.

Modified: short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg
===================================================================
--- short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg	2010-03-15 22:05:22 UTC (rev 16422)
+++ short/3D/PyLith/trunk/playpen/friction/bar_shearwave/quad4/shearwave.cfg	2010-03-15 23:37:21 UTC (rev 16423)
@@ -5,12 +5,17 @@
 # ----------------------------------------------------------------------
 # problem
 # ----------------------------------------------------------------------
-formulation = pylith.problems.Explicit
-#formulation = pylith.problems.ExplicitLumped
-formulation.solver = pylith.problems.SolverNonlinear
+[pylithapp.timedependent]
 
+#formulation = pylith.problems.Explicit
+#formulation.solver = pylith.problems.SolverNonlinear
+formulation = pylith.problems.ExplicitLumped
+
 # Nondimensionalize problem using wave propagation parameters.
 normalizer = spatialdata.units.NondimElasticDynamic
+normalizer.mass_density = 1000.0*kg/m**3
+normalizer.shear_wave_speed = 1.0*m/s
+normalizer.wave_period = 1.0*s
 
 bc = [x_pos,x_neg,all_nofault]
 
@@ -57,7 +62,7 @@
 quadrature.cell.dimension = 1
 quadrature.cell.quad_order = 2
 
-[pylithapp.timedependent.bc.y_pos]
+[pylithapp.timedependent.bc.all_nofault]
 # Dirichlet BC on all vertices except fault vertices
 
 # Pin x and z degrees of freedom
@@ -92,7 +97,7 @@
 db_initial_tractions = spatialdata.spatialdb.UniformDB
 db_initial_tractions.label = Initial fault tractions
 db_initial_tractions.values = [traction-shear,traction-normal]
-db_initial_tractions.data = [7000.0*MPa, -10000.0*MPa]
+db_initial_tractions.data = [-6001.0*MPa, -10000.0*MPa]
 
 friction = pylith.friction.StaticFriction
 



More information about the CIG-COMMITS mailing list