****************************** 
* 
* Carl W. Gable
* gable@lanl.gov
* LaGriT Control File
* http://lagrit.lanl.gov
* http://meshing.lanl.gov
* 
****************************** 
*
* LaGriT input control file to create mesh for SCEC Strike Slip benchmark.
* This mesh has increased refinement around the fault and fault tip.
*
* Carl Gable
* gable@lanl.gov
*
****************************************************
*
cmo / create / cmo1 / / / hex
cmo/setatt//ipolydat/no
*
* Working units are meters.
*
* Create mesh with 1000 meter node spacing in x,y,z directions.
*
createpts / brick / xyz / 25 25 25 / 0. 0. -24000. / 24000. 24000. 0. / 1 1 1
*
* Refine elements around the fault from 1000m to 500m
*
define / ZBOT / -17000
*
* Identify nodes in a box defined by:
*
define / XMIN /  9000.9
define / XMAX / 14000.1
define / YMIN /     0.0
define / YMAX / 17000.0
define / ZMIN /-17000.0
define / ZMAX /     1.0
pset / prefine / geom/xyz/ 1 0 0 / XMIN YMIN ZMIN / XMAX YMAX ZMAX
*
* Create a set of elements from the set of nodes
*
eltset / erefine / inclusive / pset get prefine
*
* Refine the element set.
*
refine/ eltset / eltset get erefine
eltset / erefine / delete
*
dump / gmv / hex_octree_l1.gmv / cmo1 / ascii
*
* Refine to 250m
*
define / XMIN / 11000.4
define / XMAX / 12000.6
define / YMIN /     0.0
define / YMAX / 17000.0
define / ZMIN /-17000.0
define / ZMAX /     1.0
pset / prefine1 / geom/xyz/ 1 0 0 / XMIN YMIN ZMIN / XMAX YMAX ZMAX
*
*
* Identify nodes in a cylinder,
* radius = 0 - 2000.
* theta  = 0 - 360
* z      = ZMIN - ZMAX
* Cylinder is centered on x,y,z = 12000.,16000.,0.
*
pset / prefine2 / geom/rtz/ 1 0 0 / &
    0. 0. ZBOT / 2000. 360. 1. / 12000. 16000. 0.
pset / prefine / union / prefine1 prefine2
eltset / erefine / inclusive / pset get prefine
refine/ eltset / eltset get erefine
eltset / erefine / delete
*
dump / gmv / hex_octree_l2.gmv / cmo1 / ascii
*
* Refine to 125m
*
define / XMIN / 11999.0
define / XMAX / 12000.1
define / YMIN /     0.0
define / YMAX / 17000.0
define / ZMIN /-17000.0
define / ZMAX /     1.0
pset / prefine1 / geom/xyz/ 1 0 0 / XMIN YMIN ZMIN / XMAX YMAX ZMAX
*
pset / prefine2 / geom/rtz/ 1 0 0 / &
    0. 0. ZBOT / 1000. 360. 1. / 12000. 16000. 0.
pset / prefine / union / prefine1 prefine2
eltset / erefine / inclusive / pset get prefine
refine/ eltset / eltset get erefine
eltset / erefine / delete
*
dump / gmv / hex_octree_l3.gmv / cmo1 / ascii
*
* Create a tet mesh from the points of the hex mesh.
*
cmo / create / cmo_tet / / / tet
cmo/setatt//ipolydat/no
copypts / cmo_tet / cmo1
cmo / delete / cmo1
cmo / select / cmo_tet
filter / 1 0 0
rmpoint / compress
*
cmo / setatt / cmo_tet / imt / 1 0 0 / 1
cmo / setatt / cmo_tet / itp / 1 0 0 / 0
*
* Connect point distribution to create a tet mesh.
*
connect / noadd
resetpts / itp
*
* Identify the fault nodes.
* Output files.
*
define / OUTPUT_zone  / tet_variable_refine
define / OUTPUT_table / tet_variable_refine.table
define / OUTPUT_gmv   / tet_variable_refine.gmv
define / OUTPUT_avs   / tet_variable_refine.inp
define / CMO_NAME     / cmo_tet
*
infile identify_fault_nodes.mlgi
*
infile output_files.mlgi
*
finish