# ---------------------------------------------------------------------- # Set units to SI. # ---------------------------------------------------------------------- ${Units('si')} # # ---------------------------------------------------------------------- # Reset geometry. # ---------------------------------------------------------------------- reset # ---------------------------------------------------------------------- # SEGMENTS 1 & 2 # ---------------------------------------------------------------------- ${length=25.0*km} ${east=551580.17} ${north=3786913.11} ${strike=355} ${depth=200*km} ${x1=east-sind(strike)*(length/2)} ${y1=north-cosd(strike)*(length/2)} ${x2=east+sind(strike)*(length/2)} ${y2=north+cosd(strike)*(length/2)} create vertex x {x1} y {y1} z 0 ${idTbeg1=Id('vertex')} create vertex x {x2} y {y2} z 0 ${idTend1=Id('vertex')} create vertex x {x1} y {y1} z -{depth} ${idBbeg1=Id('vertex')} create vertex x {x2} y {y2} z -{depth} ${idBend1=Id('vertex')} create vertex x {x1} y 3580000 z 0 ${idTs=Id('vertex')} create vertex x {x1} y 3580000 z -{depth} ${idBs=Id('vertex')} create surface vertex {idTbeg1} {idTend1} {idBend1} {idBbeg1} surface 1 name "seg1" create surface vertex {idTbeg1} {idTs} {idBs} {idBbeg1} surface 2 name "seg2" # ---------------------------------------------------------------------- # SEGMENT 3 # ---------------------------------------------------------------------- ${length=25.0*km} ${east=545935.24} ${north=3810946.52} ${strike=339} ${x2=east+sind(strike)*(length/2)} ${y2=north+cosd(strike)*(length/2)} create vertex x {x2} y {y2} z 0 ${idTend2=Id('vertex')} create vertex x {x2} y {y2} z -{depth} ${idBend2=Id('vertex')} create surface vertex {idTend1} {idTend2} {idBend2} {idBend1} surface 3 name "seg3" # ---------------------------------------------------------------------- # SEGMENTS 3 & 4 # ---------------------------------------------------------------------- ${length=30.0*km} ${east=531159.77} ${north=3833506.84} ${strike=317} ${x2=east+sind(strike)*(length/2)} ${y2=north+cosd(strike)*(length/2)} create vertex x {x2} y {y2} z 0 ${idTend3=Id('vertex')} create vertex x {x2} y {y2} z -{depth} ${idBend3=Id('vertex')} create vertex x {x2} y 4080000 z 0 ${idTn=Id('vertex')} create vertex x {x2} y 4080000 z -{depth} ${idBn=Id('vertex')} create surface vertex {idTend2} {idBend2} {idBend3} {idTend3} surface 4 name "seg4" create surface vertex {idTend3} {idTn} {idBn} {idBend3} surface 5 name "seg5" merge curve all unite surface all # ---------------------------------------------------------------------- # BOX # ---------------------------------------------------------------------- brick x 100000 y 200000 z 150000 ${idV=Id("volume")} volume {idV} move x 537000 y 3810000 z -{100000} webcut volume {idV} with sheet body 1