Click on image to enlarge.
Gallery written: Thu Jun 21 13:22:47 2007
| aspect_ratio
|
edge_ratio
|
fault_surface2
|
fault_surface3
|
fault_surface4
|
| fault_surface
|
read / gmv / tets_id_flt.gmv / cmo
cmo / setatt / cmo / imt / 1 0 0 / 1
cmo / setatt / cmo / itetclr / 1 0 0 / 1
cmo / copyatt / cmo / cmo / itetclr / id_block
read / gmv / fault01_ref1.gmv / cmo_tri1
*
cmo / DELATT / cmo / dfield
compute / distance_field / cmo / cmo_tri1 / dfield
pset / pmat3 / attribute / dfield / 1 0 0 / 2000. / gt
eltset / emat3 / inclusive / pset get pmat3
cmo / setatt / cmo / itetclr / eltset get emat3 / 3
resetpts / itp
settets / newtets
resetpts / cell_color
resetpts / itp
dump / gmv / recolor.gmv / cmo
define / MO_NAME / cmo
infile id_element_type.mlgi
dump / gmv / elem_type.gmv / cmo
crush_thin_tets / cmo / 0.05
rmpoint / compress
resetpts / itp
infile id_element_type.mlgi
dump / gmv / elem_type0.05.gmv / cmo
crush_thin_tets / cmo / 0.1
rmpoint / compress
resetpts / itp
infile id_element_type.mlgi
dump / gmv / elem_type0.10.gmv / cmo
crush_thin_tets / cmo / 0.15
rmpoint / compress
resetpts / itp
infile id_element_type.mlgi
dump / gmv / elem_type0.15.gmv / cmo
crush_thin_tets / cmo / 0.20
rmpoint / compress
resetpts / itp
infile id_element_type.mlgi
dump / gmv / elem_type0.20.gmv / cmo
crush_thin_tets / cmo / 0.25
rmpoint / compress
resetpts / itp
infile id_element_type.mlgi
dump / gmv / elem_type0.25.gmv / cmo
crush_thin_tets / cmo / 0.30
rmpoint / compress
resetpts / itp
infile id_element_type.mlgi
dump / gmv / elem_type0.30.gmv / cmo
dump / avs2 / elem_type0.30.inp / cmo
*
*******************************************
* Identify elements ABOVE and BELOW the fault plane.
*******************************************
*
read / gmv / fault01_ref1.gmv / cmo_tri1
*
cmo / select / cmo
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
dump / gmv / tmp_tri.gmv / cmo_tri1
cmo / DELATT / cmo / dfield
compute / distance_field / cmo / cmo_tri1 / dfield
cmo / DELATT / cmo / if_flt
cmo / addatt / cmo / if_flt / vint / scalar / nnodes
cmo / select / cmo
*
* The value here need to be larger than before because
* the fault mesh point distribution does not match the 3D tet
* mesh point distribution.
*
pset / pfault / attribute / dfield / 1 0 0 / 30.0 / lt
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
finish