[cig-commits] r14092 - in cs/cigma/trunk/examples: . citcomcu gale pylith
luis at geodynamics.org
luis at geodynamics.org
Wed Feb 18 08:14:53 PST 2009
Author: luis
Date: 2009-02-18 08:14:52 -0800 (Wed, 18 Feb 2009)
New Revision: 14092
Added:
cs/cigma/trunk/examples/citcomcu/03-compare-solutions.sh
cs/cigma/trunk/examples/gale/03-compare-solutions.sh
cs/cigma/trunk/examples/pylith/
cs/cigma/trunk/examples/pylith/01-create-meshes.sh
cs/cigma/trunk/examples/pylith/02-run-pylith.sh
cs/cigma/trunk/examples/pylith/03-compare-solutions.sh
cs/cigma/trunk/examples/pylith/04-visualize-residuals.sh
cs/cigma/trunk/examples/pylith/README
cs/cigma/trunk/examples/pylith/bottom-layer.vtk
cs/cigma/trunk/examples/pylith/reverse-slip-fault.vtk
cs/cigma/trunk/examples/pylith/strike-slip-fault.vtk
cs/cigma/trunk/examples/pylith/top-layer.vtk
Log:
Updated examples directory
Added: cs/cigma/trunk/examples/citcomcu/03-compare-solutions.sh
===================================================================
--- cs/cigma/trunk/examples/citcomcu/03-compare-solutions.sh (rev 0)
+++ cs/cigma/trunk/examples/citcomcu/03-compare-solutions.sh 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,11 @@
+#!/bin/bash
+set -x
+
+var=velocity
+
+cigma compare citcomcu.case32.0.0.vtr:${var} citcomcu.case16.0.0.vtr:${var} -o citcomcu.h5:/case_32_16 -v
+cigma compare citcomcu.case32.0.0.vtr:${var} citcomcu.case8.0.0.vtr:${var} -o citcomcu.h5:/case_32_08 -v
+
+vtk-residuals -m citcomcu.case32.0.0.vtr:${var} -i citcomcu.h5:/case_32_16 -o res-citcomcu-${var}-32-16.vtk
+vtk-residuals -m citcomcu.case32.0.0.vtr:${var} -i citcomcu.h5:/case_32_08 -o res-citcomcu-${var}-32-08.vtk
+
Property changes on: cs/cigma/trunk/examples/citcomcu/03-compare-solutions.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/gale/03-compare-solutions.sh
===================================================================
--- cs/cigma/trunk/examples/gale/03-compare-solutions.sh (rev 0)
+++ cs/cigma/trunk/examples/gale/03-compare-solutions.sh 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,16 @@
+#!/bin/bash
+set -x
+
+var=PressureField
+bm=bm.gale.circular_inclusion
+
+cigma compare p256.vts:${var} p128.vts:${var} -o circular_inclusion.h5:/pressure_256_128 -v
+cigma compare p256.vts:${var} p64.vts:${var} -o circular_inclusion.h5:/pressure_256_064 -v
+
+cigma compare ${bm}.pressure p256.vts:${var} -o circular_inclusion.h5:/pressure_256 -v
+cigma compare ${bm}.pressure p128.vts:${var} -o circular_inclusion.h5:/pressure_128 -v
+cigma compare ${bm}.pressure p64.vts:${var} -o circular_inclusion.h5:/pressure_064 -v
+
+h5attr circular_inclusion.h5:/pressure_256_128
+h5attr circular_inclusion.h5:/pressure_256_064
+
Property changes on: cs/cigma/trunk/examples/gale/03-compare-solutions.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/pylith/01-create-meshes.sh
===================================================================
--- cs/cigma/trunk/examples/pylith/01-create-meshes.sh (rev 0)
+++ cs/cigma/trunk/examples/pylith/01-create-meshes.sh 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,18 @@
+#!/bin/bash
+
+#
+# Run LaGriT to generate the tet4 meshes
+# The lagrit command is a symlink to my ~/opt/lagrit/lagrit_lin_o executable.
+#
+for r in 100m 0500m 0250m; do
+ lagrit < mesh_tet4_${r}.lagrit
+done
+
+#
+# Use CUBIT mesh generator to create the hex8 meshes
+# The cubit command is a symlink to my ~/opt/cubit/cubitclaro-11.0.0-Linux/bin/cubit executable.
+# We augment the number of resolutions to use a refinement factor of 1.2
+for r in 1000m 0833m 0694m 0578m 0500m 0482m 0250m; do
+ cubit -batch -nographics mesh_hex8_${r}.jou
+done
+
Property changes on: cs/cigma/trunk/examples/pylith/01-create-meshes.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/pylith/02-run-pylith.sh
===================================================================
--- cs/cigma/trunk/examples/pylith/02-run-pylith.sh (rev 0)
+++ cs/cigma/trunk/examples/pylith/02-run-pylith.sh 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,20 @@
+#!/bin/bash
+
+
+time pylith strikeslip_cubit.cfg strikeslip_hex8_1000m.cfg timedep.cfg 2>&1 |tee hex8_1000m.log
+time pylith strikeslip_cubit.cfg strikeslip_hex8_0833m.cfg timedep.cfg 2>&1 |tee hex8_0833m.log
+time pylith strikeslip_cubit.cfg strikeslip_hex8_0694m.cfg timedep.cfg 2>&1 |tee hex8_0694m.log
+time pylith strikeslip_cubit.cfg strikeslip_hex8_0578m.cfg timedep.cfg 2>&1 |tee hex8_0578m.log
+time pylith strikeslip_cubit.cfg strikeslip_hex8_0500m.cfg timedep.cfg 2>&1 |tee hex8_0500m.log
+time pylith strikeslip_lagrit.cfg strikeslip_tet4_1000m.cfg timedep.cfg 2>&1 |tee tet4_1000m.log
+time pylith strikeslip_lagrit.cfg strikeslip_tet4_0500m.cfg timedep.cfg 2>&1 |tee tet4_0500m.log
+
+
+time pylith reverseslip_cubit.cfg reverseslip_hex8_1000m.cfg timedep.cfg 2>&1 |tee hex8_1000m.log
+time pylith reverseslip_cubit.cfg reverseslip_hex8_0833m.cfg timedep.cfg 2>&1 |tee hex8_0833m.log
+time pylith reverseslip_cubit.cfg reverseslip_hex8_0694m.cfg timedep.cfg 2>&1 |tee hex8_0694m.log
+time pylith reverseslip_cubit.cfg reverseslip_hex8_0500m.cfg timedep.cfg 2>&1 |tee hex8_0500m.log
+time pylith reverseslip_cubit.cfg reverseslip_hex8_0578m.cfg timedep.cfg 2>&1 |tee hex8_0578m.log
+time pylith reverseslip_lagrit.cfg reverseslip_tet4_1000m.cfg timedep.cfg 2>&1 |tee tet4_1000m.log
+time pylith reverseslip_lagrit.cfg reverseslip_tet4_0500m.cfg timedep.cfg 2>&1 |tee tet4_0500m.log
+
Property changes on: cs/cigma/trunk/examples/pylith/02-run-pylith.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/pylith/03-compare-solutions.sh
===================================================================
--- cs/cigma/trunk/examples/pylith/03-compare-solutions.sh (rev 0)
+++ cs/cigma/trunk/examples/pylith/03-compare-solutions.sh 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,16 @@
+#!/bin/bash
+#
+source compare_strikeslip.sh
+
+# alternatively
+#n=0000
+#cigma compare strikeslip_hex8_0500m_t${n}.vtk:displacements strikeslip_hex8_0578m_t${n}.vtk:displacements -m strikeslip_hex8_0833m_t${n}.vtk -o strikeslipnog.h5:/displacements-hex8_0500m-hex8_0578m-t${n} -v
+#cigma compare strikeslip_hex8_0500m_t${n}.vtk:displacements strikeslip_hex8_0694m_t${n}.vtk:displacements -m strikeslip_hex8_0833m_t${n}.vtk -o strikeslipnog.h5:/displacements-hex8_0500m-hex8_0694m-t${n} -v
+#cigma compare strikeslip_hex8_0500m_t${n}.vtk:displacements strikeslip_hex8_0833m_t${n}.vtk:displacements -m strikeslip_hex8_0833m_t${n}.vtk -o strikeslipnog.h5:/displacements-hex8_0500m-hex8_0833m-t${n} -v
+#cigma compare strikeslip_hex8_0500m_t${n}.vtk:displacements strikeslip_hex8_1000m_t${n}.vtk:displacements -m strikeslip_hex8_0833m_t${n}.vtk -o strikeslipnog.h5:/displacements-hex8_0500m-hex8_1000m-t${n} -v
+
+compare_strikeslip hex8_0500m hex8_0578m -t 100
+compare_strikeslip hex8_0500m hex8_0694m -t 100
+compare_strikeslip hex8_0500m hex8_0833m -t 100
+compare_strikeslip hex8_0500m hex8_1000m -t 100
+
Property changes on: cs/cigma/trunk/examples/pylith/03-compare-solutions.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/pylith/04-visualize-residuals.sh
===================================================================
--- cs/cigma/trunk/examples/pylith/04-visualize-residuals.sh (rev 0)
+++ cs/cigma/trunk/examples/pylith/04-visualize-residuals.sh 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,34 @@
+#!/bin/bash
+set -x
+
+#
+#h5attr strikeslipnog.h5:/displacements-hex8_0500m-hex8_1000m-t0000
+#
+#vtk-residuals -m strikeslip_hex8_0833m_t0000.vtk
+# -i strikeslipnog.h5:/displacements-hex8_0500m-hex8_1000m-t0000
+# -o strikeslipnog-displacements-hex8_0500m-hex8_1000m-t0000.vtk
+#
+
+function make_vtk
+{
+ local a=$1
+ local b=$2
+ local var=displacements
+ local steps=`seq -f '%04g' 0 10 100`
+ shift 2
+ for n in $steps; do
+ vtk-residuals \
+ -i "strikeslipnog.h5:/$var-${a}-${b}-t${n}" \
+ -o "ssnog-${var}-${a}-${b}-t${n}.vtk" \
+ $*
+ done
+
+}
+
+m="results/strikeslip_hex8_0500m_t0000.vtk"
+opts="-m $m --divide-by-cell-volumes --output-log-values"
+make_vtk hex8_0500m hex8_0578m $opts
+make_vtk hex8_0500m hex8_0694m $opts
+make_vtk hex8_0500m hex8_0833m $opts
+make_vtk hex8_0500m hex8_1000m $opts
+
Property changes on: cs/cigma/trunk/examples/pylith/04-visualize-residuals.sh
___________________________________________________________________
Name: svn:executable
+ *
Added: cs/cigma/trunk/examples/pylith/README
===================================================================
--- cs/cigma/trunk/examples/pylith/README (rev 0)
+++ cs/cigma/trunk/examples/pylith/README 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,7 @@
+
+The original files for this benchmark can be obtained from the following repository location
+
+ * http://geodynamics.org/svn/cig/short/3D/PyLith/benchmarks/trunk/quasistatic/strikeslipnog/
+
+
+
Added: cs/cigma/trunk/examples/pylith/bottom-layer.vtk
===================================================================
--- cs/cigma/trunk/examples/pylith/bottom-layer.vtk (rev 0)
+++ cs/cigma/trunk/examples/pylith/bottom-layer.vtk 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,20 @@
+# vtk DataFile Version 1.0
+This line is a comment
+ASCII
+DATASET POLYDATA
+POINTS 8 float
+ 0.0 0.0 -12000.0
+24000.0 0.0 -12000.0
+ 0.0 24000.0 -12000.0
+24000.0 24000.0 -12000.0
+ 0.0 0.0 -24000.0
+24000.0 0.0 -24000.0
+ 0.0 24000.0 -24000.0
+24000.0 24000.0 -24000.0
+POLYGONS 6 30
+4 0 2 6 4
+4 1 5 7 3
+4 0 4 5 1
+4 2 3 7 6
+4 0 1 3 2
+4 4 5 7 6
Added: cs/cigma/trunk/examples/pylith/reverse-slip-fault.vtk
===================================================================
--- cs/cigma/trunk/examples/pylith/reverse-slip-fault.vtk (rev 0)
+++ cs/cigma/trunk/examples/pylith/reverse-slip-fault.vtk 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,11 @@
+# vtk DataFile Version 1.0
+This line is a comment
+ASCII
+DATASET POLYDATA
+POINTS 4 float
+ 4000.0 0.0 0.0
+ 4000.0 16000.0 0.0
+20000.0 16000.0 -16000.0
+20000.0 0.0 -16000.0
+POLYGONS 1 5
+4 0 1 2 3
Added: cs/cigma/trunk/examples/pylith/strike-slip-fault.vtk
===================================================================
--- cs/cigma/trunk/examples/pylith/strike-slip-fault.vtk (rev 0)
+++ cs/cigma/trunk/examples/pylith/strike-slip-fault.vtk 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,11 @@
+# vtk DataFile Version 1.0
+This line is a comment
+ASCII
+DATASET POLYDATA
+POINTS 4 float
+12000.0 0.0 0.0
+12000.0 16000.0 0.0
+12000.0 16000.0 -16000.0
+12000.0 0.0 -16000.0
+POLYGONS 1 5
+4 0 1 2 3
Added: cs/cigma/trunk/examples/pylith/top-layer.vtk
===================================================================
--- cs/cigma/trunk/examples/pylith/top-layer.vtk (rev 0)
+++ cs/cigma/trunk/examples/pylith/top-layer.vtk 2009-02-18 16:14:52 UTC (rev 14092)
@@ -0,0 +1,20 @@
+# vtk DataFile Version 1.0
+This line is a comment
+ASCII
+DATASET POLYDATA
+POINTS 8 float
+ 0.0 0.0 0.0
+24000.0 0.0 0.0
+ 0.0 24000.0 0.0
+24000.0 24000.0 0.0
+ 0.0 0.0 -12000.0
+24000.0 0.0 -12000.0
+ 0.0 24000.0 -12000.0
+24000.0 24000.0 -12000.0
+POLYGONS 6 30
+4 0 2 6 4
+4 1 5 7 3
+4 0 4 5 1
+4 2 3 7 6
+4 0 1 3 2
+4 4 5 7 6
More information about the CIG-COMMITS
mailing list