lc=10.0;
//lc=1.40; gives around 25600 elements //
//// Points ////

// y=0.0 //

Point( 1)={  0.0,  0.0,0.0,lc};
Point( 2)={  5.0,  0.0,0.0,lc};
Point( 3)={ 35.0,  0.0,0.0,lc};
Point( 4)={ 45.0,  0.0,0.0,lc};
Point( 5)={ 80.0,  0.0,0.0,lc};
Point( 6)={120.0,  0.0,0.0,lc};
Point( 7)={145.0,  0.0,0.0,lc};
Point( 8)={160.0,  0.0,0.0,lc};

// y=5.0 //

Point( 9)={  0.0,  5.0,0.0,lc};
Point(10)={  5.0,  5.0,0.0,lc};
Point(11)={ 35.0,  5.0,0.0,lc};
Point(12)={ 45.0,  5.0,0.0,lc};

// y=15.0 //

Point(13)={ 75.0, 15.0,0.0,lc};
Point(14)={ 80.0, 15.0,0.0,lc};
Point(15)={ 85.0, 15.0,0.0,lc};

// y=20.0 //

Point(16)={115.0, 20.0,0.0,lc};
Point(17)={120.0, 20.0,0.0,lc};

// y=25.0 //

Point(18)={ 25.0, 25.0,0.0,lc};
Point(19)={ 35.0, 25.0,0.0,lc};
Point(20)={ 75.0, 25.0,0.0,lc};
Point(21)={ 85.0, 25.0,0.0,lc};

// y=30.0 //

Point( 22)={ 70.0, 30.0,0.0,lc};
Point( 23)={ 75.0, 30.0,0.0,lc};
Point( 24)={140.0, 30.0,0.0,lc};
Point( 25)={145.0, 30.0,0.0,lc};
// y=35.0 //

Point( 26)={  0.0, 35.0,0.0,lc};
Point( 27)={  5.0, 35.0,0.0,lc};
Point( 28)={ 25.0, 35.0,0.0,lc};
Point( 29)={ 35.0, 35.0,0.0,lc};
Point( 30)={110.0, 35.0,0.0,lc};
Point( 31)={115.0, 35.0,0.0,lc};

// y=40.0 //

Point( 32)={ 65.0, 40.0,0.0,lc};
Point( 33)={ 70.0, 40.0,0.0,lc};

// y=45.0 //

Point( 34)={  0.0, 45.0,0.0,lc};
Point( 35)={  5.0, 45.0,0.0,lc};
Point( 36)={155.0, 45.0,0.0,lc};
Point( 37)={160.0, 45.0,0.0,lc};

// y=50.0 //

Point( 38)={ 60.0, 50.0,0.0,lc};
Point( 39)={ 65.0, 50.0,0.0,lc};
Point( 40)={105.0, 50.0,0.0,lc};
Point( 41)={110.0, 50.0,0.0,lc};
Point( 42)={135.0, 50.0,0.0,lc};
Point( 43)={140.0, 50.0,0.0,lc};

// y=55.0 //

Point( 44)={ 55.0, 55.0,0.0,lc};
Point( 45)={ 60.0, 55.0,0.0,lc};
Point( 46)={ 65.0, 55.0,0.0,lc};
Point( 47)={150.0, 55.0,0.0,lc};
Point( 48)={155.0, 55.0,0.0,lc};

// y=60.0 //

Point( 49)={ 50.0, 60.0,0.0,lc};
Point( 50)={ 55.0, 60.0,0.0,lc};
Point( 51)={100.0, 60.0,0.0,lc};
Point( 52)={105.0, 60.0,0.0,lc};
Point( 53)={130.0, 60.0,0.0,lc};
Point( 54)={135.0, 60.0,0.0,lc};

// y=65.0 //

Point( 55)={ 40.0, 65.0,0.0,lc};
Point( 56)={ 50.0, 65.0,0.0,lc};
Point( 57)={ 55.0, 65.0,0.0,lc};
Point( 58)={ 65.0, 65.0,0.0,lc};
Point( 59)={145.0, 65.0,0.0,lc};
Point( 60)={150.0, 65.0,0.0,lc};

// y=70.0 //

Point(61)={ 30.0, 70.0,0.0,lc};
Point(62)={ 40.0, 70.0,0.0,lc};
Point(63)={ 95.0, 70.0,0.0,lc};
Point(64)={100.0, 70.0,0.0,lc};
Point(65)={125.0, 70.0,0.0,lc};
Point(66)={130.0, 70.0,0.0,lc};

// y=75.0 //

Point( 67)={ 15.0, 75.0,0.0,lc};
Point( 68)={ 25.0, 75.0,0.0,lc};
Point( 69)={ 30.0, 75.0,0.0,lc};
Point( 70)={ 90.0, 75.0,0.0,lc};
Point( 71)={ 95.0, 75.0,0.0,lc};
Point( 72)={140.0, 75.0,0.0,lc};
Point( 73)={145.0, 75.0,0.0,lc};

// y=80.0 //

Point( 74)={  0.0, 80.0,0.0,lc};
Point( 75)={ 15.0, 80.0,0.0,lc};
Point( 76)={ 85.0, 80.0,0.0,lc};
Point( 77)={ 90.0, 80.0,0.0,lc};
Point( 78)={120.0, 80.0,0.0,lc};
Point( 79)={125.0, 80.0,0.0,lc};

// y=85.0 //

Point( 80)={ 15.0, 85.0,0.0,lc};
Point( 81)={ 25.0, 85.0,0.0,lc};
Point( 82)={ 80.0, 85.0,0.0,lc};
Point( 83)={ 85.0, 85.0,0.0,lc};
Point( 84)={115.0, 85.0,0.0,lc};
Point( 85)={120.0, 85.0,0.0,lc};
Point( 86)={135.0, 85.0,0.0,lc};
Point( 87)={140.0, 85.0,0.0,lc};

// y=90.0 //

Point( 88)={ 75.0, 90.0,0.0,lc};
Point( 89)={ 80.0, 90.0,0.0,lc};
Point( 90)={110.0, 90.0,0.0,lc};
Point( 91)={115.0, 90.0,0.0,lc};

// y=95.0 //

Point( 92)={ 70.0, 95.0,0.0,lc};
Point( 93)={ 75.0, 95.0,0.0,lc};
Point( 94)={105.0, 95.0,0.0,lc};
Point( 95)={110.0, 95.0,0.0,lc};
Point( 96)={130.0, 95.0,0.0,lc};
Point( 97)={135.0, 95.0,0.0,lc};

// y=100.0 //

Point( 98)={ 60.0,100.0,0.0,lc};
Point( 99)={ 70.0,100.0,0.0,lc};
Point(100)={100.0,100.0,0.0,lc};
Point(101)={105.0,100.0,0.0,lc};
Point(102)={125.0,100.0,0.0,lc};
Point(103)={130.0,100.0,0.0,lc};

// y=105.0 //

Point(104)={ 50.0,105.0,0.0,lc};
Point(105)={ 60.0,105.0,0.0,lc};
Point(106)={ 95.0,105.0,0.0,lc};
Point(107)={100.0,105.0,0.0,lc};
Point(108)={120.0,105.0,0.0,lc};
Point(109)={125.0,105.0,0.0,lc};

// y=110.0 //

Point(110)={ 35.0,110.0,0.0,lc};
Point(111)={ 50.0,110.0,0.0,lc};
Point(112)={ 90.0,110.0,0.0,lc};
Point(113)={ 95.0,110.0,0.0,lc};
Point(114)={115.0,110.0,0.0,lc};
Point(115)={120.0,110.0,0.0,lc};

// y=115.0 //

Point(116)={ 20.0,115.0,0.0,lc};
Point(117)={ 35.0,115.0,0.0,lc};
Point(118)={ 85.0,115.0,0.0,lc};
Point(119)={ 90.0,115.0,0.0,lc};
Point(120)={110.0,115.0,0.0,lc};
Point(121)={115.0,115.0,0.0,lc};

// y=120.0 //

Point(122)={  0.0,120.0,0.0,lc};
Point(123)={ 20.0,120.0,0.0,lc};
Point(124)={ 80.0,120.0,0.0,lc};
Point(125)={ 85.0,120.0,0.0,lc};
Point(126)={105.0,120.0,0.0,lc};
Point(127)={110.0,120.0,0.0,lc};

// y=125.0 //

Point(128)={ 70.0,125.0,0.0,lc};
Point(129)={ 80.0,125.0,0.0,lc};
Point(130)={100.0,125.0,0.0,lc};
Point(131)={105.0,125.0,0.0,lc};

// y=130.0 //

Point(132)={ 60.0,130.0,0.0,lc};
Point(133)={ 70.0,130.0,0.0,lc};
Point(134)={ 95.0,130.0,0.0,lc};
Point(135)={100.0,130.0,0.0,lc};

// y=135.0 //

Point(136)={ 50.0,135.0,0.0,lc};
Point(137)={ 60.0,135.0,0.0,lc};
Point(138)={ 85.0,135.0,0.0,lc};
Point(139)={ 95.0,135.0,0.0,lc};

// y=140.0 //

Point(140)={ 30.0,140.0,0.0,lc};
Point(141)={ 50.0,140.0,0.0,lc};
Point(142)={ 75.0,140.0,0.0,lc};
Point(143)={ 85.0,140.0,0.0,lc};

// y=145.0 //

Point(144)={  0.0,145.0,0.0,lc};
Point(145)={ 30.0,145.0,0.0,lc};
Point(146)={ 65.0,145.0,0.0,lc};
Point(147)={ 75.0,145.0,0.0,lc};

// y=150.0 //

Point(148)={ 55.0,150.0,0.0,lc};
Point(149)={ 65.0,150.0,0.0,lc};

// y=155.0 //

Point(150)={ 45.0,155.0,0.0,lc};
Point(151)={ 55.0,155.0,0.0,lc};

// y=160.0 //

Point(152)={  0.0,160.0,0.0,lc};
Point(153)={ 45.0,160.0,0.0,lc};
Point(154)={160.0,160.0,0.0,lc};

//// Lines ////

/// x-direction ///

// y=0.0 //

Line(1)={1,2};
Line(2)={2,3};
Line(3)={3,4};
Line(4)={4,5};
Line(5)={5,6};
Line(6)={6,7};
Line(7)={7,8};

// y=5.0 //

Line(8)={9,10};
Line(9)={11,12};

// y=15.0 //

Line(10)={13,14};
Line(11)={14,15};

// y=20.0 //

Line(12)={16,17};

// y=25.0 //

Line(13)={18,19};
Line(14)={20,21};

// y=30.0 //

Line(15)={22,23};
Line(16)={24,25};

// y=35.0 //

Line(17)={26,27};
Line(18)={28,29};
Line(19)={30,31};

// y=40.0 //

Line(20)={32,33};

// y=45.0 //

Line(21)={34,35};
Line(22)={36,37};

// y=50.0 //

Line(23)={38,39};
Line(24)={40,41};
Line(25)={42,43};

// y=55.0 //

Line(26)={44,45};
Line(27)={45,46};
Line(28)={47,48};

// y=60.0 //

Line(29)={49,50};
Line(30)={51,52};
Line(31)={53,54};

// y=65.0 //

Line(32)={55,56};
Line(33)={57,58};
Line(34)={59,60};

// y=70.0 //

Line(35)={61,62};
Line(36)={63,64};
Line(37)={65,66};

// y=75.0 //

Line(38)={67,68};
Line(39)={68,69};
Line(40)={70,71};
Line(41)={72,73};

// y=80.0 //

Line(42)={74,75};
Line(43)={76,77};
Line(44)={78,79};

// y=85.0 //

Line(45)={80,81};
Line(46)={82,83};
Line(47)={84,85};
Line(48)={86,87};

// y=90.0 //

Line(49)={88,89};
Line(50)={90,91};

// y=95.0 //

Line(51)={92,93};
Line(52)={94,95};
Line(53)={96,97};

// y=100.0 //

Line(54)={98,99};
Line(55)={100,101};
Line(56)={102,103};

// y=105.0 //

Line(57)={104,105};
Line(58)={106,107};
Line(59)={108,109};

// y=110.0 //

Line(60)={110,111};
Line(61)={112,113};
Line(62)={114,115};

// y=115.0 //

Line(63)={116,117};
Line(64)={118,119};
Line(65)={120,121};

// y=120.0 //

Line(66)={122,123};
Line(67)={124,125};
Line(68)={126,127};

// y=125.0 //

Line(69)={128,129};
Line(70)={130,131};

// y=130.0 //

Line(71)={132,133};
Line(72)={134,135};

// y=135.0 //

Line(73)={136,137};
Line(74)={138,139};

// y=140.0 //

Line(75)={140,141};
Line(76)={142,143};

// y=145.0 //

Line(77)={144,145};
Line(78)={146,147};

// y=150.0 //

Line(79)={148,149};

// y=155.0 //

Line(80)={150,151};

// y=160.0 //

Line(81)={152,153};
Line(82)={153,154};

/// y-direction ///

// x=0.0 //

Line(83)={1,9};
Line(84)={9,26};
Line(85)={26,34};
Line(86)={34,74};
Line(87)={74,122};
Line(88)={122,144};
Line(89)={144,152};

// x=5.0 //

Line(90)={2,10};
Line(91)={27,35};

// x=15.0 //

Line(92)={67,75};
Line(93)={75,80};

// x=20.0 //

Line(94)={116,123};

// x=25.0 //

Line(95)={18,28};
Line(96)={68,81};

// x=30.0 //

Line(97)={61,69};
Line(98)={140,145};

// x=35.0 //

Line(99)={3,11};
Line(100)={19,29};
Line(101)={110,117};

// x=40.0 //

Line(102)={55,62};

// x=45.0 //

Line(103)={4,12};
Line(104)={150,153};

// x=50.0 //

Line(105)={49,56};
Line(106)={104,111};
Line(107)={136,141};

// x=55.0 //

Line(108)={44,50};
Line(109)={50,57};
Line(110)={148,151};

// x=60.0 //

Line(111)={38,45};
Line(112)={98,105};
Line(113)={132,137};

// x=65.0 //

Line(114)={32,39};
Line(115)={46,58};
Line(116)={146,149};

// x=70.0 //

Line(117)={22,33};
Line(118)={92,99};
Line(119)={128,133};      

// x=75.0 //

Line(120)={13,20};
Line(121)={20,23};
Line(122)={88,93};
Line(123)={142,147};

// x=80.0 //

Line(124)={5,14};
Line(125)={82,89};
Line(126)={124,129};

// x=85.0 //

Line(127)={15,21};
Line(128)={76,83};
Line(129)={118,125};
Line(130)={138,143};

// x=90.0 //

Line(131)={70,77};
Line(132)={112,119};

// x=95.0 //

Line(133)={63,71};
Line(134)={106,113};
Line(135)={134,139};

// x=100.0 //

Line(136)={51,64};
Line(137)={100,107};
Line(138)={130,135};

// x=105.0 //

Line(139)={40,52};
Line(140)={94,101};
Line(141)={126,131};

// x=110.0 //

Line(142)={30,41};
Line(143)={90,95};
Line(144)={120,127};

// x=115.0 //

Line(145)={16,31};
Line(146)={84,91};
Line(147)={114,121};

// x=120.0 //

Line(148)={6,17};
Line(149)={78,85};
Line(150)={108,115};

// x=125.0 //

Line(151)={65,79};
Line(152)={102,109};

// x=130.0 //

Line(153)={53,66};
Line(154)={96,103};

// x=135.0 //

Line(155)={42,54};
Line(156)={86,97};

// x=140.0 //

Line(157)={24,43};
Line(158)={72,87};

// x=145.0 //

Line(159)={7,25};
Line(160)={59,73};

// x=150.0 //

Line(161)={47,60};

// x=155.0 //

Line(162)={36,48};

// x=160.0 //

Line(163)={8,37};
Line(164)={37,154};

// additional line //

Line(165)={27,28};
Line(166)={19,11};

//// Line Loops ////

/// Control Rod Loop ///

Line Loop(167)={1,90,-8,-83};
Line Loop(168)={26,27,115,-33,-109,-108};

/// Na Filled CRP Loop ///

Line Loop(169)={13,100,-18,-95};
Line Loop(170)={3,103,-9,-99};
Line Loop(171)={17,91,-21,-85};
Line Loop(172)={10,11,127,-14,-120};
Line Loop(173)={38,96,-45,-93,-92};

/// Most inner Internal Blanket Loop ///

Line Loop(174)={2,99,-166,-13,95,-165,-17,-84,8,-90};

/// Outer Internal Blanket Loop ///

Line Loop(175)={9,-103,4,124,-10,120,121,-15,117,-20,114,-23,111,-26,108,-29,105,-32,102,-35,97,-39,-38,92,-42,-86,21,-91,165,18,-100,166};

/// Core Loop ///

Line Loop(176)={42,93,45,-96,39,-97,35,-102,32,-105,29,109,33,-115,-27,-111,23,-114,20,-117,15,-121,14,-127,-11,-124,5,148,-12,145,-19,142,-24,139,-30,136,-36,133,-40,131,-43,128,-46,125,-49,122,-51,118,-54,112,-57,106,-60,101,-63,94,-66,-87};

/// Radial Blanket Loop ///

Line Loop(177)={66,-94,63,-101,60,-106,57,-112,54,-118,51,-122,49,-125,46,-128,43,-131,40,-133,36,-136,30,-139,24,-142,19,-145,12,-148,6,159,-16,157,-25,155,-31,153,-37,151,-44,149,-47,146,-50,143,-52,140,-55,137,-58,134,-61,132,-64,129,-67,126,-69,119,-71,113,-73,107,-75,98,-77,-88};

/// Radial Reflector Loop ///

Line Loop(178)={77,-98,75,-107,73,-113,71,-119,69,-126,67,-129,64,-132,61,-134,58,-137,55,-140,52,-143,50,-146,47,-149,44,-151,37,-153,31,-155,25,-157,16,-159,7,163,-22,162,-28,161,-34,160,-41,158,-48,156,-53,154,-56,152,-59,150,-62,147,-65,144,-68,141,-70,138,-72,135,-74,130,-76,123,-78,116,-79,110,-80,104,-81,-89}; 

/// Empty Matrix Loop ///

Line Loop(179)={80,-110,79,-116,78,-123,76,-130,74,-135,72,-138,70,-141,68,-144,65,-147,62,-150,59,-152,56,-154,53,-156,48,-158,41,-160,34,-161,28,-162,22,164,-82,-104};

Plane Surface(180)={167};
Plane Surface(181)={168};

Plane Surface(182)={169};
Plane Surface(183)={170};
Plane Surface(184)={171};
Plane Surface(185)={172};
Plane Surface(186)={173};

Plane Surface(187)={174};
Plane Surface(188)={175};
Plane Surface(189)={176};
Plane Surface(190)={177};
Plane Surface(191)={178};
Plane Surface(192)={179};

h1=10.0;

Extrude {0.0, 0.0, h1} { Surface{ 180}; } // 1  //
Extrude {0.0, 0.0, h1} { Surface{ 181}; } // 2  //
                                  
Extrude {0.0, 0.0, h1} { Surface{ 182}; } // 3  //
Extrude {0.0, 0.0, h1} { Surface{ 183}; } // 4  //
Extrude {0.0, 0.0, h1} { Surface{ 184}; } // 5  //
Extrude {0.0, 0.0, h1} { Surface{ 185}; } // 6  //
Extrude {0.0, 0.0, h1} { Surface{ 186}; } // 7  //
                                  
Extrude {0.0, 0.0, h1} { Surface{ 187}; } // 8  //
Extrude {0.0, 0.0, h1} { Surface{ 188}; } // 9  //
Extrude {0.0, 0.0, h1} { Surface{ 189}; } // 10 //
Extrude {0.0, 0.0, h1} { Surface{ 190}; } // 11 //
Extrude {0.0, 0.0, h1} { Surface{ 191}; } // 12 //
Extrude {0.0, 0.0, h1} { Surface{ 192}; } // 13 //

h2=45.0;

Extrude {0.0, 0.0, h2} { Surface{ 214}; } // 14 //
Extrude {0.0, 0.0, h2} { Surface{ 246}; } // 15 //
                                  
Extrude {0.0, 0.0, h2} { Surface{ 268}; } // 16 //
Extrude {0.0, 0.0, h2} { Surface{ 290}; } // 17 //
Extrude {0.0, 0.0, h2} { Surface{ 312}; } // 18 //
Extrude {0.0, 0.0, h2} { Surface{ 339}; } // 19 //
Extrude {0.0, 0.0, h2} { Surface{ 366}; } // 20 //
                                  
Extrude {0.0, 0.0, h2} { Surface{ 418}; } // 21 //
Extrude {0.0, 0.0, h2} { Surface{ 580}; } // 22 //
Extrude {0.0, 0.0, h2} { Surface{ 872}; } // 23 //
Extrude {0.0, 0.0, h2} { Surface{1214}; } // 24 //
Extrude {0.0, 0.0, h2} { Surface{1596}; } // 25 //
Extrude {0.0, 0.0, h2} { Surface{1788}; } // 26 //

h3=25.0;

Extrude {0.0, 0.0, h3} { Surface{1810}; } // 27 //
Extrude {0.0, 0.0, h3} { Surface{1842}; } // 28 //
                                  
Extrude {0.0, 0.0, h3} { Surface{1864}; } // 29 //
Extrude {0.0, 0.0, h3} { Surface{1886}; } // 30 //
Extrude {0.0, 0.0, h3} { Surface{1908}; } // 31 //
Extrude {0.0, 0.0, h3} { Surface{1935}; } // 32 //
Extrude {0.0, 0.0, h3} { Surface{1962}; } // 33 //
                                  
Extrude {0.0, 0.0, h3} { Surface{2014}; } // 34 //
Extrude {0.0, 0.0, h3} { Surface{2176}; } // 35 //
Extrude {0.0, 0.0, h3} { Surface{2468}; } // 36 //
Extrude {0.0, 0.0, h3} { Surface{2810}; } // 37 //
Extrude {0.0, 0.0, h3} { Surface{3192}; } // 38 //
Extrude {0.0, 0.0, h3} { Surface{3384}; } // 39 //

h4=10.0;                                                                               
                                                                                       
Extrude {0.0, 0.0, h4} { Surface{3406}; } // 40 //                                     
Extrude {0.0, 0.0, h4} { Surface{3438}; } // 41 //
                                     
Extrude {0.0, 0.0, h4} { Surface{3460}; } // 42 //                                     
Extrude {0.0, 0.0, h4} { Surface{3482}; } // 43 //                                     
Extrude {0.0, 0.0, h4} { Surface{3504}; } // 44 //                                     
Extrude {0.0, 0.0, h4} { Surface{3531}; } // 45 //                                     
Extrude {0.0, 0.0, h4} { Surface{3558}; } // 46 //                                     
                                                                                       
Extrude {0.0, 0.0, h4} { Surface{3610}; } // 47 //                                     
Extrude {0.0, 0.0, h4} { Surface{3772}; } // 48 //                                     
Extrude {0.0, 0.0, h4} { Surface{4064}; } // 49 //                                     
Extrude {0.0, 0.0, h4} { Surface{4406}; } // 50 //                                     
Extrude {0.0, 0.0, h4} { Surface{4788}; } // 51 //  
Extrude {0.0, 0.0, h4} { Surface{4980}; } // 52 //

Color Red { Volume { 10 , 21 , 22 , 23 }; }
Color Black { Volume { 47, 48, 49, 50, 51 }; }
Color Yellow { Volume { 8 , 9 , 11 , 24 , 37 }; }
Color Blue { Volume { 12, 25, 38 }; }
Color Orange { Volume { 34, 35, 36}; }
Color Green { Volume { 1,2,3,4,5,6,7,14,15,16,17,18,19,20,27,28,29,30,31,32,33,40,41,42,43,44,45,46}; }
Color White { Volume { 13, 26, 39, 52 }; }

// 11 for reflective, 12 for vacuum , 13 for zero flux boundary condition //

// 1 for core, 2 for axial reflector, 3 for radial & inner blanket, 4 for radial reflector //
// 5 for Axial Blanket, 6 for Control Rod, 7 for Na-filled CRP, 8 for Empty Matrix //

Physical Surface(11)={180,181,182,183,184,185,186,187,188,189,190,191,192,201,381,277,463,747,1065,1439,1797,1977,1873,2059,2343,2661,3035,3393,3573,3469,3655,3939,4257,4631,4989,5169,5065,5251,5535,5853,6227,1595,1213,871,555,311,409,213,3191,2809,2467,2151,1907,2005,1809,4787,4405,4063,3747,3503,3601,3405,6383,6001,5659,5343,5099,5197,5001};
Physical Surface(12)={5002,5034,5056,5078,5100,5127,5154,5206,5368,5660,6002,6384,6576,1443,1779,3039,3375,4635,4971,6231,6567,1591,1783,3187,3379,4783,4975,6379,6571};

Physical Volume(1)={10 , 21 , 22 , 23};
Physical Volume(2)={47, 48, 49, 50, 51};
Physical Volume(3)={8 , 9 , 11 , 24 , 37};
Physical Volume(4)={12, 25, 38};
Physical Volume(5)={34, 35, 36};
Physical Volume(7)={1,2,3,4,5,6,7,14,15,16,17,18,19,20,27,28,29,30,31,32,33,40,41,42,43,44,45,46};
Physical Volume(8)={13, 26, 39, 52};

                                     
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      
                                      