# ______________________________________________________________________________________
reset

cd "myDir"

# create the top and base surfaces
create vertex 278715  3002020 0
create vertex 388715  3002020 0
create vertex 388715  3112020 0
create vertex 278715  3112020 0
create surface vertex 1 2 3 4

create vertex from vertex 1 2 3 4
move vert 5 6 7 8 z -26000
create surface vertex 5 6 7 8
create volume loft Surface 1 2 
delete surf 1 2

# ______________________________________________________________________________________
# # create a geological horizon to cut the volume

create surface from surf 3
surface 9 scale x 2 y 2 about vertex 9
move surf 9 x -15000 y -15000 include_merged

move surf 9 z -3000
webcut volume 3 with sheet body 4

move surf 9 z -3000
webcut volume 5 with sheet body 4

move surf 9 z -3000
webcut volume 6 with sheet body 4

del surf 9
# ______________________________________________________________________________________

imprint all
merge all
imprint all

# ______________________________________________________________________________________

vol 7 size 1800
mesh vol 7

draw vol 7
refine surf 30 numsplit 1

vol 6 size 600
mesh vol 6
draw vol 6

# # refine surf 20 numsplit 1
# draw vol all

vol 5 size 600
mesh vol 5
draw vol 5

refine surf 10 numsplit 1

vol 3 size 200
mesh vol 3 

imprint all
merge all
imprint all

#_______________________________________________________________________
# mesh quality assessment 

# shear metrics
# quality volume all shear global draw mesh

# scaled jacobian
# quality volume all scaled jacobian global draw mesh

# ______________________________________________________________________________________

block 1 Vol 3 5 6 7 
draw block 1

draw surf        11 23 31 36 37 34 22 14 13 21 33 38 12 24 32 39 8 
block 2 Surface  11 23 31 36 37 34 22 14 13 21 33 38 12 24 32 39 8
draw block 2  

# ______________________________________________________________________________________

export abaqus "myDir/abas.inp" nodefile "/Users/kt/Academic/6_UoE_RF/6_TomorrowVille_Sims/4_Cubit/nodes" elementfile "/Users/kt/Academic/6_UoE_RF/6_TomorrowVille_Sims/4_Cubit/block" dimension 3 overwrite everything

# ______________________________________________________________________________________

# set large exodus file off
# set exodus netcdf4 on
# export mesh "myDir/tb.e" block 1 2 overwrite