# refine base box shape mesh in area of repo and fault
# Background spacing is 50m 
# refine halves are 25.0 12.5 6.25 
#
# Materials are Background 5, fault 4, repo 3
# Zones are Top left 1, Top right 2
#

read avs hex_base.inp mo
cmo/addatt/mo/itetlev/vint/scalar/nelements
cmo/setatt/mo/ipolydat no
cmo status mo
cmo select mo

# PLANES for fault 
surface / sleft / intrface / plane / &
680.0 0. 0. / 800.0 1000.0 550.0 / 800.0 0.0 550.0

#surface / sright / intrface / plane / &
#700.0 0. 0. / 820.0 1000.0 550.0 / 820.0 0.0 550.0 
# Expand fault width to capture elements of 50m spacing
surface / sright / intrface / plane / &
720.0 0. 0. / 840.0 1000.0 550.0 / 840.0 0.0 550.0 

# BOX for repository middle at 300. (200x200x20)
# choose height to include 20m at resolution 6.25m
# this sets refinement slightly larger than final set
define XMIN 100.
define XMAX 300.
define YMIN 400.
define YMAX 600.
define ZMIN 287.
define ZMAX 313.

# Refine Level 1 (25 m)
# Select element sets
# Refine octree based on element sets
define LEVEL 1
infile get_eltsets.mlgi
  cmo/setatt/mo/itetclr/ eltset,get,eflt 4
  cmo/setatt/mo/itetclr/ eltset,get,ebox 3
  dump gmv tmp_lvl0_clrs.gmv mo

refine/eltset/ eltset,get,erefine
  dump gmv tmp_lvl1.gmv mo


# Refine Level 2 (12.5 m)
# Select element sets
# Refine octree based on element sets
define LEVEL 2
infile get_eltsets.mlgi
  cmo/setatt/mo/itetclr/ eltset,get,eflt 4
  cmo/setatt/mo/itetclr/ eltset,get,ebox 3
  dump gmv tmp_lvl1_clrs.gmv mo

refine/eltset/ eltset,get,erefine
  dump gmv tmp_lvl2.gmv mo

# Refine Level 3 (6.25 m)
# Select element sets

# Refine octree based on element sets
# reset surface back to thinner width
cmo select mo
surface /sright / remove
surface / sright / intrface / plane / &
710.0 0. 0. / 830.0 1000.0 550.0 / 830.0 0.0 550.0 

define LEVEL 3
infile get_eltsets.mlgi
  cmo/setatt/mo/itetclr/ eltset,get,eflt 4
  cmo/setatt/mo/itetclr/ eltset,get,ebox 3
  dump gmv tmp_lvl2_clrs.gmv mo

refine/eltset/ eltset,get,erefine
  dump gmv hex_lvl3_octree.gmv mo

# Clean up attributes for final hex mesh 
grid2grid / tree_to_fe / mo_clean / mo
  cmo / select / mo_clean
  rmpoint / compress
  resetpts / itp

  infile get_eltsets.mlgi
  cmo/setatt/mo_clean/itetclr/ eltset,get,eflt 4
  cmo/setatt/mo_clean/itetclr/ eltset,get,ebox 3

## WRITE FINAL LEVEL 3 (6.25 m) BASE MESH 
  cmo/setatt/mo_clean/ipolydat no
  dump / avs / hex_lvl3.inp / mo_clean / 1 1 0 0
  dump / hex_lvl3.gmv / mo_clean
  cmo printatt mo_clean -all- minmax

# check point distribution of top to fault refine area
cmo select mo_clean
pset/psub/attribute zic/1,0,0 lt 400. 
rmpoint pset,get,psub
rmpoint compress

read avs tet_top.inp motop
addmesh append mov mo_clean motop
dump gmv merge_octree_top.gmv mov

finish