[cig-commits] r16943 - in short/3D/PyLith/trunk/examples/3d/hex8: . spatialdb/powerlaw

[willic3 at geodynamics.org]

Author: willic3
Date: 2010-06-08 16:16:32 -0700 (Tue, 08 Jun 2010)
New Revision: 16943

Modified:
   short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step05.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg
Log:
More minor tweaks to examples, and explicitly included reference_strain_rate
in powerlaw_gendb.cfg.



Modified: short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg	2010-06-08 22:58:33 UTC (rev 16942)
+++ short/3D/PyLith/trunk/examples/3d/hex8/spatialdb/powerlaw/powerlaw_gendb.cfg	2010-06-08 23:16:32 UTC (rev 16943)
@@ -2,6 +2,7 @@
 [powerlaw_gendb]
 
 reference_value = strain_rate
+reference_strain_rate = 1.0e-6/s
 
 [powerlaw_gendb.db_exponent]
 label = Power-law exponent

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step05.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step05.cfg	2010-06-08 22:58:33 UTC (rev 16942)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step05.cfg	2010-06-08 23:16:32 UTC (rev 16943)
@@ -129,7 +129,7 @@
 db_initial.data = [0.0*cm,-0.5*cm]
 # Velocities.
 db_rate = spatialdata.spatialdb.UniformDB
-db_rate.label = Dirichlet rate BC on +x
+db_rate.label = Dirichlet rate BC on -x
 db_rate.values = [displacement-rate-x,displacement-rate-y,rate-start-time]
 db_rate.data = [0.0*cm/year,-1.0*cm/year,0.0*year]
 

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg	2010-06-08 22:58:33 UTC (rev 16942)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg	2010-06-08 23:16:32 UTC (rev 16943)
@@ -146,10 +146,6 @@
 quadrature.cell = pylith.feassemble.FIATLagrange
 quadrature.cell.dimension = 2
 
-# Due to the complex time history used for the fault, a number of additional
-# info fields are available.
-output.vertex_info_fields = [normal_dir,strike_dir,dip_dir,final_slip_creep,final_slip_one,final_slip_two,final_slip_three,slip_time_creep,slip_time_one,slip_time_two,slip_time_three]
-
 # Set earthquake sources to an array consisting of creep and 3 ruptures.
 eq_srcs = [creep,one,two,three]
 eq_srcs.creep.origin_time = 00.0*year
@@ -175,10 +171,13 @@
 # Define slip rate and origin time for fault creep.
 [pylithapp.timedependent.interfaces.fault.eq_srcs.creep]
 slip_function = pylith.faults.ConstRateSlipFn
-origin_time = 0.0*s
 slip_function.slip_rate.iohandler.filename = spatialdb/sliprate_creep.spatialdb
 slip_function.slip_time.iohandler.filename = spatialdb/sliptime.spatialdb
 
+# Due to the complex time history used for the fault, a number of additional
+# info fields are available.
+output.vertex_info_fields = [normal_dir,strike_dir,dip_dir,final_slip_creep,final_slip_one,final_slip_two,final_slip_three,slip_time_creep,slip_time_one,slip_time_two,slip_time_three]
+
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
@@ -188,6 +187,8 @@
 # ask for output every 50 years. The time stamps of the output files are
 # in years (rather than the default of seconds), and we give a format for
 # the time stamp.
+# We also request velocity output in addition to displacements.
+vertex_data_fields = [displacement,velocity]
 output_freq = time_step
 time_step = 50.0*year
 writer.filename = output/step07.vtk
@@ -198,6 +199,8 @@
 [pylithapp.problem.formulation.output.subdomain]
 # Name of nodeset for ground surface.
 label = face_zpos
+# We also request velocity output in addition to displacements.
+vertex_data_fields = [displacement,velocity]
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
 skip = 0