[cig-commits] [commit] master: More work setting up tpv28 simulation. (cd1dce4)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith_benchmarks

On branch  : master
Link       : https://github.com/geodynamics/pylith_benchmarks/compare/b18768a541005fcf437aa50bdb2cd57e8e8adc11...cd1dce4c9bdb81b261e07c7335dbe10ad958892c

>---------------------------------------------------------------

commit cd1dce4c9bdb81b261e07c7335dbe10ad958892c
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Wed Feb 19 17:57:29 2014 -0800

    More work setting up tpv28 simulation.


>---------------------------------------------------------------

cd1dce4c9bdb81b261e07c7335dbe10ad958892c
 dynamic/scecdynrup/{tpv205 => tpv28}/pylithapp.cfg | 47 ++++++++++++----------
 dynamic/scecdynrup/{tpv102 => tpv28}/tet4.cfg      |  0
 dynamic/scecdynrup/{tpv102 => tpv28}/tet4_200m.cfg |  0
 3 files changed, 25 insertions(+), 22 deletions(-)

diff --git a/dynamic/scecdynrup/tpv205/pylithapp.cfg b/dynamic/scecdynrup/tpv28/pylithapp.cfg
similarity index 82%
copy from dynamic/scecdynrup/tpv205/pylithapp.cfg
copy to dynamic/scecdynrup/tpv28/pylithapp.cfg
index 355faf1..bd20140 100644
--- a/dynamic/scecdynrup/tpv205/pylithapp.cfg
+++ b/dynamic/scecdynrup/tpv28/pylithapp.cfg
@@ -8,9 +8,12 @@ timedependent = 1
 explicitlumped = 1
 petsc = 1
 meshiocubit = 1
+distributor = 1
 elasticityexplicit = 1
 elasticityexplicittet4 = 1
 meshimporter = 1
+#quadrature3d = 1
+#fiatlagrange = 1
 faultcohesivedyn = 1
 
 [pylithapp.journal.debug]
@@ -41,13 +44,13 @@ dimension = 3
 
 elastic_prestep = False
 formulation = pylith.problems.Explicit
-formulation.norm_viscosity = 0.3
+formulation.norm_viscosity = 0.4
 normalizer = spatialdata.units.NondimElasticDynamic
 normalizer.shear_wave_speed = 3333*m/s
 
 [pylithapp.timedependent.formulation.time_step]
-total_time = 12.001*s
-dt = 0.01*s
+total_time = 13.001*s
+dt = 0.005*s
 
 # ----------------------------------------------------------------------
 # materials
@@ -62,6 +65,11 @@ id = 1
 db_properties.label = Elastic properties
 db_properties.iohandler.filename = matprops.spatialdb
 
+db_initial_stress = spatialdata.spatialdb.UniformDB
+db_initial_stress.label = Initial stress tensor
+db_initial_stress.values = [stress-xx,stress-yy,stress-zz,stress-xy,stress-yz,stress-yz]
+db_initial_stress.data = [-60.0*MPa,-60.0*MPa,0.0,-29.38*MPa,0.0*MPa,0.0*MPa]
+
 # ----------------------------------------------------------------------
 # boundary conditions
 # ----------------------------------------------------------------------
@@ -99,6 +107,7 @@ up_dir = [0, 1, 0]
 db.label = Absorbing BC -z
 db.iohandler.filename = matprops.spatialdb
 
+
 # ----------------------------------------------------------------------
 # faults
 # ----------------------------------------------------------------------
@@ -107,36 +116,31 @@ interfaces = [fault]
 interfaces.fault = pylith.faults.FaultCohesiveDyn
 
 [pylithapp.timedependent.interfaces.fault]
-id = 10
+id = 100
 label = fault
 
 friction = pylith.friction.SlipWeakening
-friction.label = Slip weakening
-
+friction.label = Slip Weakening friction
 friction.db_properties = spatialdata.spatialdb.SimpleDB
 friction.db_properties.iohandler.filename = friction.spatialdb
-friction.db_properties.label = Slip weakening
+friction.db_properties.label = Slip weakening on main fault
+friction.db_properties.query_type = nearest
 
 traction_perturbation = pylith.faults.TractPerturbation
 
 [pylithapp.timedependent.interfaces.fault.traction_perturbation]
 db_initial = spatialdata.spatialdb.SimpleDB
-db_initial.iohandler.filename = tractions.spatialdb
-db_initial.label = Initial fault tractions
+db_initial.label = Initial fault tractions on fault
+db_initial.query_type = nearest
+db_initial.iohandler.filename = initial_traction.spatialdb
 
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-[pylithapp.timedependent.formulation]
-output = [domain,subdomain]
+[pylithapp.problem.formulation]
+output = [subdomain]
 output.subdomain = pylith.meshio.OutputSolnSubset
 
-[pylithapp.problem.formulation.output.domain]
-vertex_data_fields=[displacement,velocity]
-output_freq = time_step
-time_step = 12.0*s
-writer = pylith.meshio.DataWriterHDF5Ext
-
 [pylithapp.problem.formulation.output.subdomain]
 label = face_zpos
 vertex_data_fields=[displacement,velocity]
@@ -145,23 +149,22 @@ time_step = 0.04999*s
 writer = pylith.meshio.DataWriterHDF5Ext
 
 [pylithapp.timedependent.interfaces.fault.output]
-vertex_info_fields = [strike_dir,normal_dir,traction_initial_value,static_coefficient,dynamic_coefficient,slip_weakening_parameter,cohesion]
+vertex_info_fields = [traction_initial_value,static_coefficient,dynamic_coefficient,slip_weakening_parameter,time_weakening_time,time_weakening_parameter]
 vertex_data_fields = [slip,slip_rate,traction]
 output_freq = time_step
 time_step = 0.04999*s
 writer = pylith.meshio.DataWriterHDF5Ext
 
 [pylithapp.timedependent.materials.elastic.output]
-cell_info_fields = []
 cell_data_fields = []
 cell_filter = pylith.meshio.CellFilterAvg
-output_freq = time_step
-time_step = 12.0*s
 writer = pylith.meshio.DataWriterHDF5Ext
 
-
 # ----------------------------------------------------------------------
 # PETSc
 # ----------------------------------------------------------------------
 [pylithapp.petsc]
 log_summary = true
+
+
+# End of file
diff --git a/dynamic/scecdynrup/tpv102/tet4.cfg b/dynamic/scecdynrup/tpv28/tet4.cfg
similarity index 100%
copy from dynamic/scecdynrup/tpv102/tet4.cfg
copy to dynamic/scecdynrup/tpv28/tet4.cfg
diff --git a/dynamic/scecdynrup/tpv102/tet4_200m.cfg b/dynamic/scecdynrup/tpv28/tet4_200m.cfg
similarity index 100%
copy from dynamic/scecdynrup/tpv102/tet4_200m.cfg
copy to dynamic/scecdynrup/tpv28/tet4_200m.cfg