[cig-commits] r19832 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210

[brad at geodynamics.org]

Author: brad
Date: 2012-03-21 09:19:41 -0700 (Wed, 21 Mar 2012)
New Revision: 19832

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
Log:
Use desired total time.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou	2012-03-21 15:26:38 UTC (rev 19831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/hex8_200m.jou	2012-03-21 16:19:41 UTC (rev 19832)
@@ -1,4 +1,9 @@
+# -*- Python -*- (syntax highlighting)
 # ----------------------------------------------------------------------
+#
+# Brad T. Aagaard, U.S. Geological Survey
+#
+# ----------------------------------------------------------------------
 # Create hex8 mesh at 200m resolution.
 # ----------------------------------------------------------------------
 
@@ -8,16 +13,44 @@
 playback 'geometry.jou'
 
 # ----------------------------------------------------------------------
-# Set discretization size
+# Set discretization size and scheme
 # ----------------------------------------------------------------------
-volume all size 200
+#{dx=200*m}
+surface all scheme map
+surface ?? scheme pave
+playback 'gradient.jou'
 
 # ----------------------------------------------------------------------
 # Generate the mesh
 # ----------------------------------------------------------------------
+mesh surface all
+surface all smooth scheme condition number beta 1.3 cpu 10
+smooth surface all
+
 mesh volume all
 
 # ----------------------------------------------------------------------
+# Refine around fault
+# ----------------------------------------------------------------------
+#{refineWidth=6.001*km}
+#{refineLength=42.001*km}
+#{refineDepth=21.001*km}
+group "refinecells" add Hex with x_min >= {-0.5*refineWidth} and x_max <= {+0.5*refineWidth} and y_min >= {-0.5*refineLength} and y_max <= {+0.5*refineLength} and z_min >= {-refineDepth}
+refine hex in group refinecells depth 1
+delete group refinecells
+
+# ----------------------------------------------------------------------
+# Smooth to improve quality
+# ----------------------------------------------------------------------
+#{condnum=3.0}
+#{loop(5)}
+volume all smooth scheme condition number beta {condnum} cpu 10
+smooth volume all
+#{condnum=condnum-0.2}
+#{endloop}
+
+
+# ----------------------------------------------------------------------
 # Setup boundary conditions.
 # ----------------------------------------------------------------------
 playback 'createbc.jou'
@@ -28,3 +61,4 @@
 export mesh "hex8_200m.exo" dimension 3 overwrite
 
 
+# End of file

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg	2012-03-21 15:26:38 UTC (rev 19831)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210/pylithapp.cfg	2012-03-21 16:19:41 UTC (rev 19832)
@@ -49,8 +49,8 @@
 normalizer.mass_density = 2700.0*kg/m**3
  
 [pylithapp.timedependent.formulation.time_step]
-#total_time = 15.01*s
-total_time = 0.0*s
+total_time = 15.01*s
+#total_time = 0.0*s
 dt = 0.01*s
 
 
@@ -163,4 +163,4 @@
 log_summary = true
 
 
-# End of file
\ No newline at end of file
+# End of file