*
* 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 blocks
*
createpts / brick / xyz / 25 25 25 / 0. 0. -24000. / 24000. 24000. 0. / 1 1 1
*
* Refine to 500m
*
define / ZBOT / -17000
pset / prefine1 / geom/xyz/ 1 0 0 /9000.9 0. ZBOT/14000.1 17000. 1.
pset / prefine2 / geom/rtz/ 1 0 0 / &
    0. 0. ZBOT / 3. 360. 1. / 12000. 16000. 0.
pset / prefine / union / prefine1 prefine2
eltset / erefine / inclusive / pset get prefine
refine/ eltset / eltset get erefine
dump / gmv / hex_octree_l1.gmv / cmo1 / ascii

eltset / erefine / delete
*
* Refine to 250m
*
define / ZBOT / -17000
pset / prefine1 / geom/xyz/ 1 0 0 /11000.4 0. ZBOT/12000.6 17000. 1.
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
dump / gmv / hex_octree_l2.gmv / cmo1 / ascii

eltset / erefine / delete
*
* Refine to 125m
*
define / ZBOT / -17000
pset / prefine1 / geom/xyz/ 1 0 0 /11999.0 0.0 ZBOT/12000.1 17000. 1.
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
dump / gmv / hex_octree_l3.gmv / cmo1 / ascii

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
*
finish