SetFactory("OpenCASCADE");

Mesh.CharacteristicLengthMin = 0.003; 
Mesh.CharacteristicLengthMax = 0.015;

// ---------- Voxel ---------- 
dim = 1.2500; corner = -0.2500; v = newv; 
Box(v) = {corner,corner,corner,dim,dim,dim}; 
svox() = Boundary{Volume{v};}; 
//Physical Surface(0) = {svox()}; 

// -------- Vessel # 1 -------- 
id1 = 76; ox1 = 0.0000; oy1 = 0.0000; oz1 = 0.0000; 
vx1 = 0.0123; vy1 = 0.4787; vz1 = 0.4975; r1 = 0.0100; 
c1 = newv; Cylinder(c1) = {ox1,oy1,oz1,vx1,vy1,vz1,r1}; 
// ------ Vessel # 2 ------ 
id2 = 24; ox2 = 0.0123; oy2 = 0.4787; oz2 = 0.4975; 
vx2 = 0.0123; vy2 = -0.0163; vz2 = 0.0024; r2 = 0.0100; 
c2 = newv; Cylinder(c2) = {ox2,oy2,oz2,vx2,vy2,vz2,r2}; 
union2() = BooleanUnion{ Volume{c1}; Delete;}{ Volume{c2}; Delete;};
// ------ Vessel # 3 ------ 
id3 = 22; ox3 = 0.0246; oy3 = 0.4624; oz3 = 0.4999; 
vx3 = 0.0477; vy3 = -0.0721; vz3 = 0.0082; r3 = 0.0100; 
c3 = newv; Cylinder(c3) = {ox3,oy3,oz3,vx3,vy3,vz3,r3}; 
union3() = BooleanUnion{ Volume{union2()}; Delete;}{ Volume{c3}; Delete;};
// ------ Vessel # 4 ------ 
id4 = 360; ox4 = 0.0723; oy4 = 0.3903; oz4 = 0.5081; 
vx4 = 0.0312; vy4 = -0.0593; vz4 = 0.0072; r4 = 0.0100; 
c4 = newv; Cylinder(c4) = {ox4,oy4,oz4,vx4,vy4,vz4,r4}; 
union4() = BooleanUnion{ Volume{union3()}; Delete;}{ Volume{c4}; Delete;};
// ------ Vessel # 5 ------ 
id5 = 2; ox5 = 0.1035; oy5 = 0.3310; oz5 = 0.5153; 
vx5 = 0.0305; vy5 = -0.0561; vz5 = 0.0059; r5 = 0.0100; 
c5 = newv; Cylinder(c5) = {ox5,oy5,oz5,vx5,vy5,vz5,r5}; 
union5() = BooleanUnion{ Volume{union4()}; Delete;}{ Volume{c5}; Delete;};
// ------ Vessel # 6 ------ 
id6 = 8; ox6 = 0.1340; oy6 = 0.2749; oz6 = 0.5211; 
vx6 = 0.0343; vy6 = -0.0537; vz6 = -0.0628; r6 = 0.0100; 
c6 = newv; Cylinder(c6) = {ox6,oy6,oz6,vx6,vy6,vz6,r6}; 
union6() = BooleanUnion{ Volume{union5()}; Delete;}{ Volume{c6}; Delete;};
// ------ Vessel # 7 ------ 
id7 = 304; ox7 = 0.1683; oy7 = 0.2212; oz7 = 0.4584; 
vx7 = 0.0038; vy7 = -0.0123; vz7 = -0.0034; r7 = 0.0100; 
c7 = newv; Cylinder(c7) = {ox7,oy7,oz7,vx7,vy7,vz7,r7}; 
union7() = BooleanUnion{ Volume{union6()}; Delete;}{ Volume{c7}; Delete;};
// ------ Vessel # 8 ------ 
id8 = 318; ox8 = 0.1721; oy8 = 0.2089; oz8 = 0.4549; 
vx8 = 0.0019; vy8 = -0.0150; vz8 = -0.0105; r8 = 0.0100; 
c8 = newv; Cylinder(c8) = {ox8,oy8,oz8,vx8,vy8,vz8,r8}; 
union8() = BooleanUnion{ Volume{union7()}; Delete;}{ Volume{c8}; Delete;};
// ------ Vessel # 9 ------ 
id9 = 4; ox9 = 0.1740; oy9 = 0.1939; oz9 = 0.4444; 
vx9 = 0.0019; vy9 = -0.0150; vz9 = -0.0105; r9 = 0.0100; 
c9 = newv; Cylinder(c9) = {ox9,oy9,oz9,vx9,vy9,vz9,r9}; 
union9() = BooleanUnion{ Volume{union8()}; Delete;}{ Volume{c9}; Delete;};
// ------ Vessel # 10 ------ 
id10 = 1; ox10 = 0.1759; oy10 = 0.1789; oz10 = 0.4338; 
vx10 = 0.0419; vy10 = -0.1492; vz10 = -0.0873; r10 = 0.0100; 
c10 = newv; Cylinder(c10) = {ox10,oy10,oz10,vx10,vy10,vz10,r10}; 
union10() = BooleanUnion{ Volume{union9()}; Delete;}{ Volume{c10}; Delete;};
// ------ Vessel # 11 ------ 
id11 = 110; ox11 = 0.1759; oy11 = 0.1789; oz11 = 0.4338; 
vx11 = -0.0117; vy11 = 0.0257; vz11 = -0.0505; r11 = 0.0100; 
c11 = newv; Cylinder(c11) = {ox11,oy11,oz11,vx11,vy11,vz11,r11}; 
union11() = BooleanUnion{ Volume{union10()}; Delete;}{ Volume{c11}; Delete;};
// ------ Vessel # 12 ------ 
id12 = 88; ox12 = 0.1642; oy12 = 0.2047; oz12 = 0.3833; 
vx12 = -0.0037; vy12 = 0.0097; vz12 = -0.0211; r12 = 0.0100; 
c12 = newv; Cylinder(c12) = {ox12,oy12,oz12,vx12,vy12,vz12,r12}; 
union12() = BooleanUnion{ Volume{union11()}; Delete;}{ Volume{c12}; Delete;};
// ------ Vessel # 13 ------ 
id13 = 5; ox13 = 0.1605; oy13 = 0.2143; oz13 = 0.3622; 
vx13 = 0.0079; vy13 = 0.1520; vz13 = -0.0949; r13 = 0.0100; 
c13 = newv; Cylinder(c13) = {ox13,oy13,oz13,vx13,vy13,vz13,r13}; 
union13() = BooleanUnion{ Volume{union12()}; Delete;}{ Volume{c13}; Delete;};
// ------ Vessel # 14 ------ 
id14 = 230; ox14 = 0.1605; oy14 = 0.2143; oz14 = 0.3622; 
vx14 = -0.0520; vy14 = -0.0127; vz14 = -0.0787; r14 = 0.0100; 
c14 = newv; Cylinder(c14) = {ox14,oy14,oz14,vx14,vy14,vz14,r14}; 
union14() = BooleanUnion{ Volume{union13()}; Delete;}{ Volume{c14}; Delete;};
// ------ Vessel # 15 ------ 
id15 = 234; ox15 = 0.1085; oy15 = 0.2017; oz15 = 0.2836; 
vx15 = -0.0065; vy15 = -0.0172; vz15 = -0.0301; r15 = 0.0100; 
c15 = newv; Cylinder(c15) = {ox15,oy15,oz15,vx15,vy15,vz15,r15}; 
union15() = BooleanUnion{ Volume{union14()}; Delete;}{ Volume{c15}; Delete;};
// ------ Vessel # 16 ------ 
id16 = 122; ox16 = 0.1020; oy16 = 0.1845; oz16 = 0.2535; 
vx16 = -0.0147; vy16 = -0.0006; vz16 = -0.0383; r16 = 0.0100; 
c16 = newv; Cylinder(c16) = {ox16,oy16,oz16,vx16,vy16,vz16,r16}; 
union16() = BooleanUnion{ Volume{union15()}; Delete;}{ Volume{c16}; Delete;};
// ------ Vessel # 17 ------ 
id17 = 164; ox17 = 0.0873; oy17 = 0.1838; oz17 = 0.2152; 
vx17 = 0.0046; vy17 = -0.0078; vz17 = -0.0172; r17 = 0.0100; 
c17 = newv; Cylinder(c17) = {ox17,oy17,oz17,vx17,vy17,vz17,r17}; 
union17() = BooleanUnion{ Volume{union16()}; Delete;}{ Volume{c17}; Delete;};
// ------ Vessel # 18 ------ 
id18 = 89; ox18 = 0.0918; oy18 = 0.1761; oz18 = 0.1980; 
vx18 = -0.0027; vy18 = -0.0078; vz18 = -0.0099; r18 = 0.0100; 
c18 = newv; Cylinder(c18) = {ox18,oy18,oz18,vx18,vy18,vz18,r18}; 
union18() = BooleanUnion{ Volume{union17()}; Delete;}{ Volume{c18}; Delete;};
// ------ Vessel # 19 ------ 
id19 = 165; ox19 = 0.0918; oy19 = 0.1761; oz19 = 0.1980; 
vx19 = 0.0270; vy19 = -0.0474; vz19 = -0.0891; r19 = 0.0100; 
c19 = newv; Cylinder(c19) = {ox19,oy19,oz19,vx19,vy19,vz19,r19}; 
union19() = BooleanUnion{ Volume{union18()}; Delete;}{ Volume{c19}; Delete;};
// ------ Vessel # 20 ------ 
id20 = 123; ox20 = 0.1188; oy20 = 0.1287; oz20 = 0.1089; 
vx20 = -0.0891; vy20 = 0.0495; vz20 = 0.0495; r20 = 0.0100; 
c20 = newv; Cylinder(c20) = {ox20,oy20,oz20,vx20,vy20,vz20,r20}; 
union20() = BooleanUnion{ Volume{union19()}; Delete;}{ Volume{c20}; Delete;};
// ------ Vessel # 21 ------ 
id21 = 235; ox21 = 0.0297; oy21 = 0.1782; oz21 = 0.1584; 
vx21 = 0.0990; vy21 = -0.1188; vz21 = 0.0693; r21 = 0.0100; 
c21 = newv; Cylinder(c21) = {ox21,oy21,oz21,vx21,vy21,vz21,r21}; 
union21() = BooleanUnion{ Volume{union20()}; Delete;}{ Volume{c21}; Delete;};
// ------ Vessel # 22 ------ 
id22 = 231; ox22 = 0.1287; oy22 = 0.0594; oz22 = 0.2277; 
vx22 = -0.0891; vy22 = 0.1584; vz22 = 0.0495; r22 = 0.0100; 
c22 = newv; Cylinder(c22) = {ox22,oy22,oz22,vx22,vy22,vz22,r22}; 
union22() = BooleanUnion{ Volume{union21()}; Delete;}{ Volume{c22}; Delete;};
// ------ Vessel # 23 ------ 
id23 = 366; ox23 = 0.0396; oy23 = 0.2178; oz23 = 0.2772; 
vx23 = 0.1035; vy23 = 0.0216; vz23 = 0.0877; r23 = 0.0100; 
c23 = newv; Cylinder(c23) = {ox23,oy23,oz23,vx23,vy23,vz23,r23}; 
union23() = BooleanUnion{ Volume{union22()}; Delete;}{ Volume{c23}; Delete;};
// ------ Vessel # 24 ------ 
id24 = 111; ox24 = 0.1431; oy24 = 0.2394; oz24 = 0.3649; 
vx24 = -0.0144; vy24 = 0.0082; vz24 = -0.0184; r24 = 0.0100; 
c24 = newv; Cylinder(c24) = {ox24,oy24,oz24,vx24,vy24,vz24,r24}; 
union24() = BooleanUnion{ Volume{union23()}; Delete;}{ Volume{c24}; Delete;};
// ------ Vessel # 25 ------ 
id25 = 367; ox25 = 0.1431; oy25 = 0.2394; oz25 = 0.3649; 
vx25 = -0.0045; vy25 = 0.0181; vz25 = 0.0014; r25 = 0.0100; 
c25 = newv; Cylinder(c25) = {ox25,oy25,oz25,vx25,vy25,vz25,r25}; 
union25() = BooleanUnion{ Volume{union24()}; Delete;}{ Volume{c25}; Delete;};
// ------ Vessel # 26 ------ 
id26 = 319; ox26 = 0.1386; oy26 = 0.2574; oz26 = 0.3663; 
vx26 = -0.1287; vy26 = -0.1782; vz26 = 0.1188; r26 = 0.0100; 
c26 = newv; Cylinder(c26) = {ox26,oy26,oz26,vx26,vy26,vz26,r26}; 
union26() = BooleanUnion{ Volume{union25()}; Delete;}{ Volume{c26}; Delete;};
// ------ Vessel # 27 ------ 
id27 = 305; ox27 = 0.0099; oy27 = 0.0792; oz27 = 0.4851; 
vx27 = 0.2376; vy27 = 0.0990; vz27 = 0.0396; r27 = 0.0100; 
c27 = newv; Cylinder(c27) = {ox27,oy27,oz27,vx27,vy27,vz27,r27}; 
union27() = BooleanUnion{ Volume{union26()}; Delete;}{ Volume{c27}; Delete;};
// ------ Vessel # 28 ------ 
id28 = 26; ox28 = 0.2475; oy28 = 0.1782; oz28 = 0.5248; 
vx28 = -0.0341; vy28 = 0.0462; vz28 = -0.1043; r28 = 0.0100; 
c28 = newv; Cylinder(c28) = {ox28,oy28,oz28,vx28,vy28,vz28,r28}; 
union28() = BooleanUnion{ Volume{union27()}; Delete;}{ Volume{c28}; Delete;};
// ------ Vessel # 29 ------ 
id29 = 32; ox29 = 0.2134; oy29 = 0.2244; oz29 = 0.4204; 
vx29 = 0.0225; vy29 = 0.0016; vz29 = -0.0337; r29 = 0.0100; 
c29 = newv; Cylinder(c29) = {ox29,oy29,oz29,vx29,vy29,vz29,r29}; 
union29() = BooleanUnion{ Volume{union28()}; Delete;}{ Volume{c29}; Delete;};
// ------ Vessel # 30 ------ 
id30 = 10; ox30 = 0.2360; oy30 = 0.2260; oz30 = 0.3867; 
vx30 = 0.0225; vy30 = 0.0016; vz30 = -0.0337; r30 = 0.0100; 
c30 = newv; Cylinder(c30) = {ox30,oy30,oz30,vx30,vy30,vz30,r30}; 
union30() = BooleanUnion{ Volume{union29()}; Delete;}{ Volume{c30}; Delete;};
// ------ Vessel # 31 ------ 
id31 = 56; ox31 = 0.2585; oy31 = 0.2276; oz31 = 0.3529; 
vx31 = 0.0330; vy31 = 0.0069; vz31 = -0.0857; r31 = 0.0100; 
c31 = newv; Cylinder(c31) = {ox31,oy31,oz31,vx31,vy31,vz31,r31}; 
union31() = BooleanUnion{ Volume{union30()}; Delete;}{ Volume{c31}; Delete;};
// ------ Vessel # 32 ------ 
id32 = 9; ox32 = 0.2916; oy32 = 0.2346; oz32 = 0.2673; 
vx32 = 0.0253; vy32 = 0.0031; vz32 = -0.0197; r32 = 0.0100; 
c32 = newv; Cylinder(c32) = {ox32,oy32,oz32,vx32,vy32,vz32,r32}; 
union32() = BooleanUnion{ Volume{union31()}; Delete;}{ Volume{c32}; Delete;};
// ------ Vessel # 33 ------ 
id33 = 57; ox33 = 0.2916; oy33 = 0.2346; oz33 = 0.2673; 
vx33 = -0.0044; vy33 = 0.0031; vz33 = -0.0297; r33 = 0.0100; 
c33 = newv; Cylinder(c33) = {ox33,oy33,oz33,vx33,vy33,vz33,r33}; 
union33() = BooleanUnion{ Volume{union32()}; Delete;}{ Volume{c33}; Delete;};
// ------ Vessel # 34 ------ 
id34 = 298; ox34 = 0.2871; oy34 = 0.2376; oz34 = 0.2376; 
vx34 = 0.0431; vy34 = -0.0306; vz34 = 0.0328; r34 = 0.0100; 
c34 = newv; Cylinder(c34) = {ox34,oy34,oz34,vx34,vy34,vz34,r34}; 
union34() = BooleanUnion{ Volume{union33()}; Delete;}{ Volume{c34}; Delete;};
// ------ Vessel # 35 ------ 
id35 = 336; ox35 = 0.3303; oy35 = 0.2071; oz35 = 0.2704; 
vx35 = 0.0359; vy35 = -0.0103; vz35 = -0.0413; r35 = 0.0100; 
c35 = newv; Cylinder(c35) = {ox35,oy35,oz35,vx35,vy35,vz35,r35}; 
union35() = BooleanUnion{ Volume{union34()}; Delete;}{ Volume{c35}; Delete;};
// ------ Vessel # 36 ------ 
id36 = 14; ox36 = 0.3662; oy36 = 0.1968; oz36 = 0.2292; 
vx36 = 0.0359; vy36 = -0.0103; vz36 = -0.0412; r36 = 0.0100; 
c36 = newv; Cylinder(c36) = {ox36,oy36,oz36,vx36,vy36,vz36,r36}; 
union36() = BooleanUnion{ Volume{union35()}; Delete;}{ Volume{c36}; Delete;};
// ------ Vessel # 37 ------ 
id37 = 100; ox37 = 0.4020; oy37 = 0.1865; oz37 = 0.1879; 
vx37 = 0.0495; vy37 = 0.0040; vz37 = -0.0330; r37 = 0.0100; 
c37 = newv; Cylinder(c37) = {ox37,oy37,oz37,vx37,vy37,vz37,r37}; 
union37() = BooleanUnion{ Volume{union36()}; Delete;}{ Volume{c37}; Delete;};
// ------ Vessel # 38 ------ 
id38 = 236; ox38 = 0.4515; oy38 = 0.1905; oz38 = 0.1549; 
vx38 = 0.0257; vy38 = 0.0051; vz38 = -0.0164; r38 = 0.0100; 
c38 = newv; Cylinder(c38) = {ox38,oy38,oz38,vx38,vy38,vz38,r38}; 
union38() = BooleanUnion{ Volume{union37()}; Delete;}{ Volume{c38}; Delete;};
// ------ Vessel # 39 ------ 
id39 = 18; ox39 = 0.4772; oy39 = 0.1956; oz39 = 0.1385; 
vx39 = 0.0422; vy39 = 0.0051; vz39 = -0.0330; r39 = 0.0100; 
c39 = newv; Cylinder(c39) = {ox39,oy39,oz39,vx39,vy39,vz39,r39}; 
union39() = BooleanUnion{ Volume{union38()}; Delete;}{ Volume{c39}; Delete;};
// ------ Vessel # 40 ------ 
id40 = 28; ox40 = 0.5194; oy40 = 0.2006; oz40 = 0.1056; 
vx40 = 0.0819; vy40 = 0.0057; vz40 = 0.0129; r40 = 0.0100; 
c40 = newv; Cylinder(c40) = {ox40,oy40,oz40,vx40,vy40,vz40,r40}; 
union40() = BooleanUnion{ Volume{union39()}; Delete;}{ Volume{c40}; Delete;};
// ------ Vessel # 41 ------ 
id41 = 282; ox41 = 0.6013; oy41 = 0.2063; oz41 = 0.1184; 
vx41 = 0.0183; vy41 = -0.0117; vz41 = -0.0059; r41 = 0.0100; 
c41 = newv; Cylinder(c41) = {ox41,oy41,oz41,vx41,vy41,vz41,r41}; 
union41() = BooleanUnion{ Volume{union40()}; Delete;}{ Volume{c41}; Delete;};
// ------ Vessel # 42 ------ 
id42 = 70; ox42 = 0.6196; oy42 = 0.1947; oz42 = 0.1125; 
vx42 = 0.0183; vy42 = -0.0117; vz42 = -0.0096; r42 = 0.0100; 
c42 = newv; Cylinder(c42) = {ox42,oy42,oz42,vx42,vy42,vz42,r42}; 
union42() = BooleanUnion{ Volume{union41()}; Delete;}{ Volume{c42}; Delete;};
// ------ Vessel # 43 ------ 
id43 = 11; ox43 = 0.6379; oy43 = 0.1830; oz43 = 0.1029; 
vx43 = 0.0453; vy43 = -0.0146; vz43 = -0.1029; r43 = 0.0100; 
c43 = newv; Cylinder(c43) = {ox43,oy43,oz43,vx43,vy43,vz43,r43}; 
union43() = BooleanUnion{ Volume{union42()}; Delete;}{ Volume{c43}; Delete;};
// ------ Vessel # 44 ------ 
id44 = 74; ox44 = 0.6379; oy44 = 0.1830; oz44 = 0.1029; 
vx44 = 0.1342; vy44 = -0.0792; vz44 = 0.0552; r44 = 0.0100; 
c44 = newv; Cylinder(c44) = {ox44,oy44,oz44,vx44,vy44,vz44,r44}; 
union44() = BooleanUnion{ Volume{union43()}; Delete;}{ Volume{c44}; Delete;};
// ------ Vessel # 45 ------ 
id45 = 71; ox45 = 0.7721; oy45 = 0.1038; oz45 = 0.1581; 
vx45 = 0.0596; vy45 = -0.0543; vz45 = -0.0195; r45 = 0.0100; 
c45 = newv; Cylinder(c45) = {ox45,oy45,oz45,vx45,vy45,vz45,r45}; 
union45() = BooleanUnion{ Volume{union44()}; Delete;}{ Volume{c45}; Delete;};
// ------ Vessel # 46 ------ 
id46 = 330; ox46 = 0.7721; oy46 = 0.1038; oz46 = 0.1581; 
vx46 = 0.0132; vy46 = -0.0069; vz46 = 0.0292; r46 = 0.0100; 
c46 = newv; Cylinder(c46) = {ox46,oy46,oz46,vx46,vy46,vz46,r46}; 
union46() = BooleanUnion{ Volume{union45()}; Delete;}{ Volume{c46}; Delete;};
// ------ Vessel # 47 ------ 
id47 = 260; ox47 = 0.7853; oy47 = 0.0969; oz47 = 0.1873; 
vx47 = 0.0132; vy47 = -0.0127; vz47 = 0.0263; r47 = 0.0100; 
c47 = newv; Cylinder(c47) = {ox47,oy47,oz47,vx47,vy47,vz47,r47}; 
union47() = BooleanUnion{ Volume{union46()}; Delete;}{ Volume{c47}; Delete;};
// ------ Vessel # 48 ------ 
id48 = 75; ox48 = 0.7984; oy48 = 0.0842; oz48 = 0.2135; 
vx48 = 0.0333; vy48 = 0.0049; vz48 = 0.0241; r48 = 0.0100; 
c48 = newv; Cylinder(c48) = {ox48,oy48,oz48,vx48,vy48,vz48,r48}; 
union48() = BooleanUnion{ Volume{union47()}; Delete;}{ Volume{c48}; Delete;};
// ------ Vessel # 49 ------ 
id49 = 261; ox49 = 0.7984; oy49 = 0.0842; oz49 = 0.2135; 
vx49 = 0.0135; vy49 = -0.0644; vz49 = 0.1033; r49 = 0.0100; 
c49 = newv; Cylinder(c49) = {ox49,oy49,oz49,vx49,vy49,vz49,r49}; 
union49() = BooleanUnion{ Volume{union48()}; Delete;}{ Volume{c49}; Delete;};
// ------ Vessel # 50 ------ 
id50 = 331; ox50 = 0.8119; oy50 = 0.0198; oz50 = 0.3168; 
vx50 = -0.0198; vy50 = 0.0891; vz50 = -0.1089; r50 = 0.0100; 
c50 = newv; Cylinder(c50) = {ox50,oy50,oz50,vx50,vy50,vz50,r50}; 
union50() = BooleanUnion{ Volume{union49()}; Delete;}{ Volume{c50}; Delete;};

// ------- Definition of Entire Vessel Surface -------- 
sfinal() = Boundary{Volume{union50()};}; 
Physical Surface(1) = {sfinal()}; 

// ------- Definition of Final Volume -------- 
voxcyl() = BooleanDifference{Volume{v};}{Volume{union50()}; Delete;}; 
//Physical Volume(0) = {voxcyl()};