[cig-commits] r11380 - short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes
brad at geodynamics.org
brad at geodynamics.org
Fri Mar 7 18:15:59 PST 2008
Author: brad
Date: 2008-03-07 18:15:59 -0800 (Fri, 07 Mar 2008)
New Revision: 11380
Modified:
short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit
Log:
Fixed fault trimming. Merged materials into elastic and viscoelastic.
Modified: short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit
===================================================================
--- short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit 2008-03-08 02:12:46 UTC (rev 11379)
+++ short/3D/PyLith/benchmarks/trunk/quasistatic/reversenog/meshes/mesh_tet4.lagrit 2008-03-08 02:15:59 UTC (rev 11380)
@@ -50,6 +50,12 @@
4.0e+3 domain_ymax 0.0e+0 / &
4.0e+3 0.0e+0 0.0e+0 / &
domain_xmax 0.0e+0 -20.0e+3 /
+*
+* Create plane to trim fault
+surface / surf_ftrim / intrface / plane / &
+ domain_xmax 16.0e+3 0.0e+0 / &
+ 0.0e+0 16.0e+3 0.0e+0 / &
+ 0.0e+0 16.0e+3 domain_zmin /
* Create material interface
surface / surf_mantle / intrface / plane / &
0.0 0.0 -12.0e+3 / &
@@ -59,59 +65,84 @@
* ----------------------------------------------------------------------
* Create regions
* ----------------------------------------------------------------------
+* Regions do include internal interfaces (but no overlap).
+** Elastic region beyond end of fault
+region / r_elastic_end / &
+ le surf_xpos and le surf_xneg &
+ and lt surf_ftrim and le surf_ypos &
+ and le surf_zpos and ge surf_mantle /
+
+** Viscoelastic region beyond end of fault
+region / r_visco_end / &
+ le surf_xpos and le surf_xneg &
+ and lt surf_ftrim and le surf_ypos &
+ and lt surf_mantle and le surf_zneg /
+
** Elastic on hanging wall side of the fault
region / r_elastic_hang / &
le surf_xpos and ge surf_fault &
- and le surf_ypos and le surf_yneg &
+ and ge surf_ftrim and le surf_yneg &
and le surf_zpos and ge surf_mantle /
** Elastic on footwell side of the fault
region / r_elastic_foot / &
le surf_xneg and lt surf_fault &
- and le surf_ypos and le surf_yneg &
+ and ge surf_ftrim and le surf_yneg &
and le surf_zpos and ge surf_mantle /
** Viscoelastic on hanging wall side of the fault
region / r_visco_hang / &
le surf_xpos and ge surf_fault &
- and le surf_ypos and le surf_yneg &
+ and ge surf_ftrim and le surf_yneg &
and lt surf_mantle and le surf_zneg /
** Viscoelastic on footwall side of the fault
region / r_visco_foot / &
le surf_xneg and lt surf_fault &
- and le surf_ypos and le surf_yneg &
+ and ge surf_ftrim and le surf_yneg &
and lt surf_mantle and le surf_zneg /
* ----------------------------------------------------------------------
* Create material regions
* ----------------------------------------------------------------------
-* Material regions do not include internal interfaces (regions do)
+* Material regions do not include internal interfaces (regions do).
*
-* Material 1
+** Material 1: Elastic region beyond end of fault
+mregion / mr_elastic_end / &
+ le surf_xpos and le surf_xneg &
+ and lt surf_ftrim and le surf_ypos &
+ and le surf_zpos and gt surf_mantle /
+
+** Material 2: Viscoelastic region beyond end of fault
+mregion / mr_visco_end / &
+ le surf_xpos and le surf_xneg &
+ and lt surf_ftrim and le surf_ypos &
+ and lt surf_mantle and le surf_zneg /
+
+* Material 3
mregion / mr_elastic_hang / &
le surf_xpos and gt surf_fault &
- and le surf_ypos and le surf_yneg &
+ and gt surf_ftrim and le surf_yneg &
and le surf_zpos and gt surf_mantle /
-* Material 2
+* Material 4
mregion / mr_elastic_foot / &
le surf_xneg and lt surf_fault &
- and le surf_ypos and le surf_yneg &
+ and gt surf_ftrim and le surf_yneg &
and le surf_zpos and gt surf_mantle /
-* Material 3
+* Material 5
mregion / mr_visco_hang / &
le surf_xpos and gt surf_fault &
- and le surf_ypos and le surf_yneg &
+ and gt surf_ftrim and le surf_yneg &
and lt surf_mantle and le surf_zneg /
-* Material 4
+* Material 6
mregion / mr_visco_foot / &
le surf_xneg and lt surf_fault &
- and le surf_ypos and le surf_yneg &
+ and gt surf_ftrim and le surf_yneg &
and lt surf_mantle and le surf_zneg /
*
* ----------------------------------------------------------------------
@@ -126,7 +157,9 @@
*
loop / foreach / iregion / &
- r_elastic_hang r_elastic_foot r_visco_hang r_visco_foot &
+ r_elastic_end r_visco_end &
+ r_elastic_hang r_elastic_foot &
+ r_visco_hang r_visco_foot &
/ loop_end / &
regnpts / iregion / dx / pset,get,rayend / xyz / &
0.0e+3 0.0e+3 domain_zmin / &
@@ -150,7 +183,7 @@
* ----------------------------------------------------------------------
* Connect the points into a Delaunay tetrahedral mesh
* ----------------------------------------------------------------------
-connect
+connect / check_interface
* ----------------------------------------------------------------------
* Set elements to tetrahedra
@@ -175,8 +208,8 @@
* Identify nodes on fault surface
pset / tmp / surface / surf_fault /
pset / fault / geom / xyz / pset,get,tmp / &
- domain_xmin -0.001e+3 -16.001e+3 / &
- domain_xmax 16.001e+3 0.001e+3 /
+ domain_xmin -0.001e+3 -16.001e+3 / &
+ domain_xmax 16.001e+3 +0.001e+3 /
pset / tmp / delete
*
* Identify nodes on boundaries
@@ -200,6 +233,30 @@
pset / boundary_zpos / surface / surf_zpos /
* ----------------------------------------------------------------------
+* MERGE MATERIALS
+* ----------------------------------------------------------------------
+eltset / elastic1 / itetclr / eq / 1
+eltset / elastic2 / itetclr / eq / 3
+eltset / elastic3 / itetclr / eq / 4
+eltset / elastic / union / elastic1 elastic2 elastic3
+cmo / setatt / box / itetclr / eltset,get,elastic / 1
+eltset / elastic1 / delete
+eltset / elastic2 / delete
+eltset / elastic3 / delete
+eltset / elastic / delete
+
+eltset / visco1 / itetclr / eq / 2
+eltset / visco2 / itetclr / eq / 5
+eltset / visco3 / itetclr / eq / 6
+eltset / visco / union / visco1 visco2 visco3
+cmo / setatt / box / itetclr / eltset,get,visco / 2
+eltset / visco1 / delete
+eltset / visco2 / delete
+eltset / visco3 / delete
+eltset / visco / delete
+cmo / printatt / box / itetclr / minmax
+
+* ----------------------------------------------------------------------
finish
* End of file
More information about the cig-commits
mailing list