reset
reset error

# Create Domain
brick x 10 y 10 z 50

# Electrodes
# xplane surfaces
create surface rectangle width 1 xplane
move surface 7 location x 5 y 0 z -15 include_merged

create surface rectangle width 1 xplane
move surface 8 location x 5 y 0 z -5 include_merged

create surface rectangle width 1 xplane
move surface 9 location x 5 y 0 z 5 include_merged

create surface rectangle width 1 xplane
move surface 10 location x 5 y 0 z 15 include_merged

# yplane surfaces
create surface rectangle width 1 yplane
move surface 11 location x 0 y 5 z -15 include_merged

create surface rectangle width 1 yplane
move surface 12 location x 0 y 5 z -5 include_merged

create surface rectangle width 1 yplane
move surface 13 location x 0 y 5 z 5 include_merged

create surface rectangle width 1 yplane
move surface 14 location x 0 y 5 z 15 include_merged

# xplane surfaces
create surface rectangle width 1 xplane
move surface 15 location x -5 y 0 z -15 include_merged

create surface rectangle width 1 xplane
move surface 16 location x -5 y 0 z -5 include_merged

create surface rectangle width 1 xplane
move surface 17 location x -5 y 0 z 5 include_merged

create surface rectangle width 1 xplane
move surface 18 location x -5 y 0 z 15 include_merged

# yplane surfaces
create surface rectangle width 1 yplane
move surface 19 location x 0 y -5 z -15 include_merged

create surface rectangle width 1 yplane
move surface 20 location x 0 y -5 z -5 include_merged

create surface rectangle width 1 yplane
move surface 21 location x 0 y -5 z 5 include_merged

create surface rectangle width 1 yplane
move surface 22 location x 0 y -5 z 15 include_merged

#group 'geom' equals surface 1 to 6
#group 'elec' equals surface 7 to 22

surface 7 to 22 size 0.5
surface 7 to 22 scheme TriMesh
mesh surface 7 to 22


volume 1 size 0.5
volume 1 scheme tetmesh
mesh volume 1