****************************** 
* 
* Carl W. Gable
* gable@lanl.gov
* LaGriT Control File
* http://lagrit.lanl.gov
* http://meshing.lanl.gov
* 
****************************** 
read / gmv / elem_type0.30.gmv / cmo
read / gmv / faults_merge3.gmv / cmo_tri1
*
* Identify elements ABOVE and BELOW the fault plane.
* The algorithm tests 
*
surface / s_flt1 / reflect / sheet / cmo_tri1
cmo / select / cmo
region / r_above / ge s_flt1
region / r_below / lt s_flt1
eltset / e_below / region / r_below

cmo / setatt / cmo / id_block / 1 0 0 / 1
cmo / setatt / cmo / id_block / eltset get e_below / 2

cmo / select / cmo_tri1
refine2d
refine2d

cmo / DELATT / cmo / dfield
compute / distance_field / cmo / cmo_tri1 / dfield

cmo / addatt / cmo / if_flt / vint / scalar / nnodes
pset / pfault / attribute / dfield / 1 0 0 / lt / 3.

cmo / setatt / cmo / if_flt / 1 0 0 / 1
cmo / setatt / cmo / if_flt / pset get pfault / 2


dump / gmv / tets_id_flt_improved.gmv / cmo_tet2

finish