cmo / create / cmo1 / / / triplane
*
* Define some parameters
*
define / NX / 51
define / NY / 51
define / NZ / 3
define / NZ / 1
define / XMAX / 200000.
define / YMAX / 200000.
define / ZMAX / 40000.
define / R1 / 1000.
define / R2 / 100000.
define / R3 / 200000.

define / THETA / 10.
define / N_TH_1 / 5
define / N_TH_2 / 9
define / N_TH_3 / 17

define / THETA / 90.
define / N_TH_1 / 17
define / N_TH_2 / 33
define / N_TH_3 / 65

*
* Define Geometry
*
surface / cyl_r1 / intrface / cylinder / 0. 0. 0. / 0. 0. ZMAX / R1
surface / cyl_r2 / intrface / cylinder / 0. 0. 0. / 0. 0. ZMAX / R2
surface / cyl_r3 / intrface / cylinder / 0. 0. 0. / 0. 0. ZMAX / R3

surface / b_box / reflect / box / 0. 0. 0. / XMAX YMAX ZMAX

region / r1 / le b_box and le cyl_r1
region / r2 / le b_box and ge cyl_r1 and le cyl_r2
region / r3 / le b_box and ge cyl_r2 and le cyl_r3

region / r4 / le b_box and ge cyl_r3

mregion / mr1 / le b_box and lt cyl_r1
mregion / mr2 / le b_box and gt cyl_r1 and lt cyl_r2
mregion / mr3 / le b_box and gt cyl_r2 and lt cyl_r3

mregion / mr4 / le b_box and gt cyl_r3
*
* Create Points
*
createpts / rtz / 5 N_TH_1 NZ / 0. 0. 0. / R1 THETA ZMAX / 1 1 1
createpts / rtz / 5 N_TH_2 NZ / R1 0. 0. / 5000. THETA ZMAX / 1 1 1
createpts / rtz / 10 N_TH_3 NZ / 5000. 0. 0. / R2 THETA ZMAX / &
   1 1 1 / 1 0 0 / 0.95 1. 1.
createpts / rtz / 10 N_TH_3 NZ / R2 0. 0. / R3 THETA ZMAX / &
   1 1 1 / 1 0 0 / 1.1 1. 1.

filter / 1 0 0
rmpoint / compress

setpts

dump / gmv / pts_cyl2.gmv / cmo1

connect

dump / gmv / connect_cyl2.gmv / cmo1

settets / geometry
resetpts/cell_color

dump / gmv / settets_cyl2.gmv / cmo1

rmmat / 1 / element
* resetpts / parents
rmpoint / compress
resetpts / itp


quality

dump / gmv / bm7.gmv / cmo1
dump / zone / bm7 / cmo1 / keepatt

resetpts / parents
rmpoint / compress
*
* Locate slip nodes
*
pset / p_r_minus / geom / rtz / 1 0 0 / 0. 0. -1. / 99999.0 THETA 40001.
pset / p_r_plus / geom / rtz / 1 0 0 / 0. 0. -1. / 100001.0 THETA 40001.
pset / p_fault / not p_r_plus p_r_minus

cmo/addatt/cmo1/f_offset/VDOUBLE/scalar/nnodes
cmo / setatt / cmo1 / f_offset / 1 0 0 / 0.0

pset / p_z1 / geom / xyz / 1 0 0 / -1. -1. 27990.9 / 1.e6 1.e6 40001.1
pset / p_z2 / geom / xyz / 1 0 0 / -1. -1. 26990.9 / 1.e6 1.e6 27000.1
pset / p_z3 / geom / xyz / 1 0 0 / -1. -1. 25990.9 / 1.e6 1.e6 26000.1
pset / p_z4 / geom / xyz / 1 0 0 / -1. -1. 24990.9 / 1.e6 1.e6 25000.1

pset / p_f_z1 / inter / p_fault p_z1
pset / p_f_z2 / inter / p_fault p_z2
pset / p_f_z3 / inter / p_fault p_z3
pset / p_f_z4 / inter / p_fault p_z4

cmo / setatt / cmo1 / f_offset / pset get p_f_z1 / 1000.
cmo / setatt / cmo1 / f_offset / pset get p_f_z2 / 750.
cmo / setatt / cmo1 / f_offset / pset get p_f_z3 / 500.
cmo / setatt / cmo1 / f_offset / pset get p_f_z4 / 250.


dump / gmv / id_fault.gmv / cmo1
dump / avs / id_fault.inp / cmo1

finish

massage / 500. / 0.1 / 0.001

dump / gmv / bm7_massage1.gmv / cmo1
dump / avs / bm7_massage1.inp / cmo1

pset / p_per / attribute / itp / 1 0 0 / 0 / eq
perturb/ pset get p_per /1.,1.,0

smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
recon 0

smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
recon 0

smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
recon 0

smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
recon 1

smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
smooth/position/esug/1 0 0 
recon 1

dump / gmv / smooth_recon.gmv / cmo1

dump / gmv / bm7_massage2.gmv / cmo1

dump / gmv / bm7_perturb.gmv / cmo1

finish