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

# read and add attribute for use with LEVEL values
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 

# BOX for repository middle at 300. (200x200x20)
# choose height to include 20m at resolution 3.125m
# 300 +- 12.5m

define XMIN 100.
define XMAX 300.
define YMIN 400.
define YMAX 600.
define ZMIN 287.
define ZMAX 313.

# Refine Level 1 (12.5 m)
# Select element sets
# Refine octree based on element sets
define LEVEL 1

infile get_eltsets_max.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 (6.25 m)
# Select element sets
# Refine octree based on element sets

# move surfaces out just a tad to capture coarse elements
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 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 (3.125 m)
# Select element sets
# Refine octree based on element sets
define LEVEL 3

infile get_eltsets_min.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 / motmp / mo
  cmo / select / motmp
  rmpoint / compress
  resetpts / itp

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

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


# Refine once more in case higher resolution needed
cmo select mo

# Refine Level 4 (1.5625)
# Select element sets
# Refine octree based on element sets
define LEVEL 4

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

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

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

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

## FINAL LEVEL Level 4 (1.5625) BASE MESH 
  cmo/setatt/motmp/ipolydat no
  dump / avs / hex_lvl4.inp / motmp / 1 1 0 0
  dump / hex_lvl4.gmv / motmp
  cmo printatt motmp -all- minmax
  cmo/delete/motmp

finish