diff --git a/crun.i b/crun.i index b14b92f..a6a4f06 100644 --- a/crun.i +++ b/crun.i @@ -1,2 +1,2 @@ CONTINUE -nps 1e7 +nps 1e10 diff --git a/e2-estia.i b/e2-estia.i index 4ede172..3f4631a 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -2,9 +2,8 @@ Input File ESTIA c ------------------------------------------------------- c --------------- CELL CARDS -------------------------- c ------------------------------------------------------- -read file=e2-ess_cells.i echo +read file=geometry/ess_cells.i echo c tallies -700 0 -700 701 9 2.45e-05 -701 702 0 -702 703 0 -703 @@ -14,453 +13,19 @@ c tallies 707 0 -707 708 0 -708 709 0 -709 -c elliptic guide -c NBOA beam -1001 9 2.45e-05 1100 -1101 -1301 -1302 -1303 -1304 -1002 9 2.45e-05 1101 -1102 (-1321:-1322) -1311 -1501 -1502 -1003 9 2.45e-05 1102 -1103 (-1331:-1332) -1311 -1501 -1502 -1004 9 2.45e-05 1103 -1104 (-1341:-1342) -1311 -1501 -1502 -1005 9 2.45e-05 1104 -1105 (-1351:-1352) -1311 -1501 -1502 -1006 9 2.45e-05 1105 -1106 (-1361:-1362) -1311 -1502 -1501 -1007 9 2.45e-05 1106 -1107 (-1371:-1372) -1311 -1502 -1501 -1008 9 2.45e-05 1107 -1108 (-1381:-1382) -1311 -1502 -1501 -1011 9 2.45e-05 1100 -1101 #(1201 1202 1203 1204) 1305 1306 1307 1308 -1012 9 2.45e-05 ((1101 -1115 1213 -1507 -1508) - :(1115 -1105 1214 -1509 -1510)) #(1211 -1503 -1504) -1013 9 2.45e-05 1105 -1107 #(1212 -1505 -1506) 1214 -1509 -1510 -1014 9 2.45e-05 ((1117 -1108 1213 -1507 -1508) - :(1107 -1117 1214 -1509 -1510) - ) #(1211 -1503 -1504) -c NBOA outside -1021 14 -8.92 (1100 -1101 1201 1202 1203 1204 - (1301 : 1302 : 1303 : 1304)) -1022 14 -8.92 (1101 -1103 1211 -1503 -1504) - (1311:1501:1502) -1023 14 -8.92 (1103 -1104 1211 -1503 -1504) - (1311:1501:1502) -1024 14 -8.92 (1104 -1105 1211 -1503 -1504) - (1311:1501:1502) -1025 14 -8.92 (1105 -1107 1212 -1505 -1506) - (1311:1501:1502) -1026 14 -8.92 (1107 -1108 1211 -1504 -1503) - (1311:1501:1502) -c NBOA wedges -1032 14 -8.92 ((1101 -1102 1321 1322) - :(1102 -1103 1331 1332) - ) -1502 -1311 -1033 14 -8.92 1103 -1104 1341 1342 -1502 -1311 -1034 14 -8.92 1104 -1105 1351 1352 -1502 -1311 -1035 14 -8.92 ((1105 -1106 1361 1362) - :(1106 -1107 1371 1372) - ) -1502 -1311 -1036 14 -8.92 1107 -1108 1381 1382 -1502 -1311 -c Al window and He in front of NBOA -1041 5 0.0582256 61 -1110 1305 1306 1307 1308 -1042 9 2.45e-05 1110 -1100 1305 1306 1307 1308 -c Al window and He behind NBOA -1051 9 2.45e-05 1108 -1118 1213 -1508 -1507 -1052 5 0.0582256 1118 -64 1213 -1508 -1507 -c extended cut-out in reflector -1135 0 -61 3 58 1205 1206 1207 1208 #700 -1136 5 0.0582256 -61 3 #(41 -42) 58 1201 1202 1203 1204 #700 #1135 -c split of wedge of monolith components for increased importance -1236 5 0.0582256 ( 142 -332 -556 325 ) ( -1201 : -1202 : -1203 : -1204 : -58) - ( -39 : 40 : -41 : 42 ) -61 3 21 -46 45 -44 43 - dxc1:n=0.001 -1240 121 0.0862739 ( 142 -332 -556 325 ) - ( -43 : -45 : 46 : 44 ) -58 ( -53 : 54 : 55 ) - dxc1:n=0.001 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1 -1241 121 0.0862739 ( 142 -332 -556 325 ) ( -1201 : -1202 : -1203 : -1204 ) - -61 58 ( -43 : -45 : 46 : 44 ) - dxc1:n=0.1 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1 -c partitioned air in bunker -1255 4 -1.2E-3 3 -68 -2 -101 66 -67 64 112 - ( 329 : 4615 : -342 : 353 : 4100 ) - ( -66 : 4615 : -342 : 353 : -3111 ) - #(-1605 -1623) #(-1622 1623 -1624) - #(-1639 1624 -1637) #(-1658 1637 -1660) - #3137 -1256 4 -1.2E-3 ( -329 -4615 342 -353 ) 64 -3110 - #( 2100 -2101 2201 2202 2203 2204 ) -1257 4 -1.2E-3 (( -329 -4615 342 -353 ) 3110 -4805 - #((3110 -3111 -3541) - :(3111 -4100 -3531))) - #(-3440 -3404 3446) - #(-3441 -3326 3446) - #(-3442 -3326 3446) -1258 4 -1.2E-3 (( 66 3111 329 #(4811 -4812 -4820) ) - :( -329 -4615 4805 - #( 5101 4604 5501 5502 -5505 5504 ) ) - #((4811 -4812 -4820) - :(-4814 -4815 -4602 4813 4805 -4101) - :(4812 -4602 ((4601 -4603 -4818 -4819):(4603 -4604 -4818))) - :(4602 4611 4612 4613 4614 -4615)) - ) 342 -353 -68 #4201 #4202 #4203 #4204 #4131 -c part around inserts E01,E02,E03 -1489 54 0.0833854 109 556 -334 325 -326 -64 330 3 - ( -153 : 154 : 159 ) -329 -1490 54 0.0833854 -556 -1107 - ( -336 : -330 : -341 : -459 : 64 : 460 : 342 ) - ( -348 : -330 : -353 : -459 : 64 : 460 : 354 ) - -334 -329 -326 325 -64 330 3 ( -153 : 154 : 159 ) - ( -330 : -569 : -459 : 64 : 460 : 570 : 564 ) -1491 54 0.0833854 -556 1107 - ( -336 : -330 : -341 : -459 : 64 : 460 : 342 ) - ( -348 : -330 : -353 : -459 : 64 : 460 : 354 ) - -334 -329 -326 325 -64 330 3 ( -153 : 154 : 159 ) - ( -330 : -569 : -459 : 64 : 460 : 570 : 564 ) -1494 54 0.0833854 -556 -1103 - ( -336 : -61 : -337 : -575 : 457 : 330 : 338 : ( 575 565 -457 ) ) - ( -348 : -61 : -349 : -575 : 457 : 330 : 350 : ( 575 554 -457 ) ) - ( -153 : 154 : 159 ) 3 61 -330 -326 325 335 -332 - ( -61 : -565 : -1001 : 457 : 330 : 566 : 349 : 564 : ( -457 349 575 ) ) -1495 54 0.0833854 -556 1103 - ( -336 : -61 : -337 : -575 : 457 : 330 : 338 : ( 575 565 -457 ) ) - ( -348 : -61 : -349 : -575 : 457 : 330 : 350 : ( 575 554 -457 ) ) - ( -153 : 154 : 159 ) 3 61 -330 -326 325 335 -332 - ( -61 : -565 : -1001 : 457 : 330 : 566 : 564 : ( -457 349 575 ) ) -1586 3 0.0876276 ( -554 : -575 : 457 ) -330 #589 - -458 -345 351 -352 1103 348 -1588 3 0.0876276 -64 -332 459 355 -356 330 - 348 1107 #589 -1591 3 0.0876276 ( -349 : -1001 : 457 ) -330 - -458 1000 567 -568 -564 -2 1103 - (-1101 : -1213 : 1508 : 1507) -1593 3 0.0876276 -64 -332 459 571 -572 330 -564 -2 1107 - ( -1213 : 1508 : 1507) - #(1107 -1117 1214 -1509 -1510) -1596 3 0.0876276 ( -565 : -575 : 457 ) -330 - -458 -345 339 -340 1103 336 -1598 3 0.0876276 -64 -332 459 343 -344 330 336 1107 -c SKADI in-bunker components -c collimators -1901 0 64 3 ((-1604 -1623):(-1622 1623 -1624) - :(-1639 1624 -1637):(-1658 1637 -1660)) - #1940 #1945 #1957 #1958 #1963 #1967 #1982 #1983 -1940 5 0.0582256 1604 -1605 -1623 64 3 -1945 73 0.0844385 ( -1594 : 1595 : -1615 : 1616 ) -1607 1606 -1618 1617 -1599 - 1598 1614 -1957 0 ( -1690 : 1691 ) ( -1647 : 1646 : 1645 ) ( -1643 : 1644 : 1645 ) -1637 - 1636 -1638 ( -1648 : -1653 : -1598 : -1650 : 1599 : 1654 : 1649 ) -1958 5 0.0582256 1638 -1639 -1641 1640 -1963 73 0.0844385 ( 1595 : -1594 : -1651 : 1652 ) -1649 1648 -1654 1653 1598 - -1599 1650 -1967 5 0.0582256 1657 -1658 -1660 1659 -1982 73 0.0844385 -1622 ( -1631 : -1598 : -1632 : 1633 : 1599 ) - -1624 1623 -1983 73 0.0844385 -1638 ( -1650 : -1598 : -1653 : 1654 : 1599 ) - -1691 1690 -c -c BBG beam top/beam bottom/outside -2001 9 2.45e-05 2100 -2101 -2301 -2302 -2303 -2304 3 #701 -2011 9 2.45e-05 2100 -2101 -2301 -2302 -2305 -2306 3 -2021 14 -8.92 2100 -2101 2201 2202 2203 2204 3 - (2301:2302:2303:2306:(2304 2305)) -c NFGA beam top -3001 0 3100 -3101 -3301 -3311 -3501 -3502 -3002 0 3101 -3102 -3302 -3312 -3501 -3502 -3003 0 3102 -3103 -3302 -3313 -3501 -3502 -3004 0 3103 -3104 -3302 -3314 -3501 -3502 -3005 0 3104 -3105 -3302 -3315 -3501 -3502 -3006 0 3105 -3106 -3302 -3316 -3501 -3502 -3007 0 3106 -3107 -3302 -3317 -3501 -3502 -c NFGA beam bottom -3011 0 3100 -3101 -3301 -3321 -3501 -3502 -3012 0 3101 -3102 -3302 -3322 -3501 -3502 -3013 0 3102 -3103 -3302 -3323 -3501 -3502 -3014 0 3103 -3104 -3302 -3324 -3501 -3502 -3015 0 3104 -3105 -3302 -3325 -3501 -3502 -3016 0 3105 -3106 -3302 -3326 -3501 -3502 -3017 0 3106 -3107 -3302 -3327 -3501 -3502 -c NFGA outside -3021 5 0.0582256 3100 -3101 -3503 -3504 3201 - #(-3301 #(3311 3321) -3501 -3502) -3022 5 0.0582256 3101 -3102 -3503 -3504 3202 - #(-3302 #(3312 3322) -3501 -3502) -3023 5 0.0582256 3102 -3103 -3503 -3504 3202 - #(-3302 #(3313 3323) -3501 -3502) -3024 5 0.0582256 3103 -3104 -3503 -3504 3202 - #(-3302 #(3314 3324) -3501 -3502) -3025 5 0.0582256 3104 -3105 -3503 -3504 3202 - #(-3302 #(3315 3325) -3501 -3502) -3026 5 0.0582256 3105 -3106 -3503 -3504 3202 - #(-3302 #(3316 3326) -3501 -3502) -3027 5 0.0582256 3106 -3107 -3503 -3504 3202 - #(-3302 #(3317 3327) -3501 -3502) -c NFSA Vaccum tank and connection to collimator -3122 0 ((((3112 -3131 -3522 3) - :(3131 -4805 -3532 3)) #( - (3100 -3101 3201 -3504 -3503) - :(3101 -3107 3202 -3504 -3503))) - :(-3434 -3404 -3541) - :(-3435 -3326 -3541) - :(-3436 -3326 -3541)) - #702 #3151 -3123 5 0.0582256 ((3110 -3111 3522 -3521) - :(3110 -3112 -3521) - :(3113 -3111 3532 -3521) - :(3111 -4805 3532 -3531)) 3 - #(-3431 -3326) - #(-3432 -3326) - #(-3433 -3326) -3124 0 3446 3541 3437 -3440 -3326 #3131 #3134 -3125 0 3446 3541 3438 -3441 -3326 #3132 #3135 -3126 0 3446 3541 3439 -3442 -3326 #3133 #3136 -3131 5 0.0582256 3434 -3431 3522 3443 -3326 -3132 5 0.0582256 3435 -3432 3522 3443 -3326 -3133 5 0.0582256 3436 -3433 3522 3443 -3326 -3134 71 0.0843223 ((-3431 3446 -3443) - :(-3440 3444 -3443) - :(-3440 3446 -3445)) 3434 -3135 71 0.0843223 ((-3432 3446 -3443) - :(-3441 3444 -3443) - :(-3441 3446 -3445)) 3435 -3136 71 0.0843223 ((-3433 3446 -3443) - :(-3442 3444 -3443) - :(-3442 3446 -3445)) 3436 -3137 49 0.074981 3421 3422 -3423 3424 3425 3426 -c B-Al plate -3141 10 -2.7 3131 -3113 3532 -3522 3 -c MirroBor outside vessel -3142 17 -1.36 3110 -3111 3521 -3541 3 - #(-3431 -3326) - #(-3432 -3326) - #(-3433 -3326) -c NFSA Guide support -3151 71 0.0843223 (3120 -3121 (3401 3402 3403 3404 3) - #(3411 3412 -3413 3414)) - :(3446 -3404 -3437) - :(3446 -1000 -3438) - :(3446 -1000 -3439) -c -----Collimators---- -c CPC1 - Collimator before VS -4001 0 -4301 -4302 -4303 -4304 4805 -4101 -2 -4010 10 -2.7 -4814 -4815 -4602 4813 4805 -4807 #4001 -4011 14 -8.92 -4814 -4815 -4602 4813 4807 4819 #4001 -4012 14 -8.92 4201 4202 -4705 4704 -4708 -4101 #4001 -c Vacuum tank around VS -4101 0 (4705 4616 4617 4610 4619 -4618) #4121 -4102 0 (4704 -4705 ((-4706 -4707 -4708 -4603):(4603 -4702))) - #(-4101 4704 4201 4202) - #(5604 -5601 -5611) - #5701 #5702 #5711 #5712 #5751 #5761 #4121 - #5021 #703 #704 -4103 0 (4703 -4704 -4701) - #(-4101 4704 4201 4202) - #5711 #5712 #5721 #5722 #5731 -4104 0 (4603 -4604 4702 -5301 -5302 -5303 -5304) -4105 0 (5604 -5601 -5612) - #(5604 -5601 (5612:(5613 5616):(5614 5615) - :(5614 5616):(-5614 -5618):(-5616 5617)) -5611) - #5701 #5702 -4111 5 0.0582256 ((4602 4611 4612 4613 4614 -4615): - (4705 -4602 ((4601 -4603 -4818 -4819) - :(4603 -4604 -4818)))) - #(4616 4617 4610 4619 -4618) -4112 5 0.0582256 (4704 -4705 ((4601 -4603 -4818 -4819):(4603 -4604 -4818))) - #((-4706 -4707 -4708 -4603):(4603 -4702)) - #(-5301 -5302 -5303 -5304) - #(-4301 -4302 -4303 -4304) -4113 5 0.0582256 (4812 -4704 ((4601 -4603 -4818 -4819):(4603 -4604 -4818))) - #(4703 -4704 -4701) -4114 5 0.0582256 4811 -4812 -4820 #(4703 -4704 -4701) -c chopper disk -4121 5 0.0582256 -4621 4622 -4623 -c chopper spindle -4131 3 0.0876276 -4612 4631 -4632 -4604 -c Concrete base (5% B4C) -4201 16 -2.4 (4801 4802 4804 -4811 4805 4806) -4202 16 -2.4 (4801 4802 4811 -4812 4805 4806) - #(-4820) - #(-4817 -4814 -4815) -4203 16 -2.4 (4801 4802 4812 -4813 4805 4806) - #(-4818 -4819) - #(-4816 -4814 -4815) -4204 16 -2.4 (4801 4802 4813 4803 4805 4806) - #(-4818 -4819) - #(-4814 -4815) -c CPC2 with Cu nose -5001 0 ((5101 4604 -5102 -5301 -5302 -5303 -5304 3) - :(5102 -7100 -5532 3)) #705 -5002 5 0.0582256 (5102 -7100 5532 -5531 3) -5021 14 -8.92 5401 5402 -5403 5404 5405 5406 - #(-5301 -5302 -5303 -5304) -5022 14 -8.92 5101 4604 -5102 5501 5502 -4602 5504 3 - #5001 -5023 71 0.0843223 5101 4604 -5102 5501 5502 4602 -5505 3 - #5001 -c VS Lead shielding (FIX MATERIALS) -c boral part top and bottom -5601 10 -2.7 5602 -5601 (5612:(5613 5616):(5614 5615) - :(5614 5616):(-5614 -5618):(-5616 5617)) - -5611 -5602 10 -2.7 5604 -5603 (5612:(5613 5616):(5614 5615) - :(5614 5616):(-5614 -5618):(-5616 5617)) - -5611 -c lead part -5603 13 -11.34 5603 -5602 (5612 - :(5613 5616):(5614 5615) - :(5614 5616) - :(-5614 -5618):(-5616 5617)) - -5611 -c VS holder stands -5701 5 0.0582256 -5701 5704 -5703 -5702 5 0.0582256 -5702 5704 -5703 -c simplified model of stages (hollow cylinders) -c x-stages -5711 5 0.0582256 5732 -5731 -5711 - #(5742 -5741 -5721) -5712 5 0.0582256 5732 -5731 -5712 - #(5742 -5741 -5722) -c z-stages -5721 5 0.0582256 5733 -5732 -5713 - #(5744 -5743 -5723) -5722 5 0.0582256 5733 -5732 -5714 - #(5744 -5743 -5724) -c omega-stage -5731 5 0.0582256 5734 -5733 -5715 - #(5746 -5745 -5725) -c B4C blades -5751 20 -2.3 ((-5751 -5754 5755) - :(-5752 -5755 5703)) -5753 5756 -5757 -5761 20 -2.3 ((-5762 -5764 5765) - :(-5763 -5765 5703)) -5761 -5766 5767 -c -----Wall Cut-Out--- -6001 4 -1.2E-3 ((68 -6100 -6201 -6202 -6203 -6204) - :(6100 -6101 -6301 -6302 -6303 -6304)) 3 - #(-5102 5501 5502 -5505 5504) - #(5102 -7100 -5531) -6002 4 -1.2E-3 ((6101 -6102 -6401 -6402 -6403 -6404) - :(6102 -6103 -6501 -6502 -6503 -6504)) 3 - #(5102 -7100 -5531) - #((7100 -7102 7201 7202 7203 7204) - :(7102 -7103 7205 7206 7207 7208)) - #(7103 -6103 -8921) -6003 4 -1.2E-3 (6103 -72 -6601 -6602 -6603 -6604 -6605) 3 - #(7103 8405 -8921) - #(7102 -7103 7205 7206 7207 7208) - #(-8405 8401 8402 8403 8404) -c Vacuum before and inside Selene Guide 1 -6101 0 ((7103 8405 -8922 3) - :(-8405 8406 -8411 8412 8413 8402 3)) - #706 #8311 #8411 #8412 #8413 -6102 5 0.0582256 (7103 8405 8922 -8921 3) -c Vaccum after Selene 1, until and inside Selen 2 -6111 0 ((-8406 -8100 -8932 3) - :(8100 -8101 -8501 -8502 -8503 -8504 3) - :(8101 -8121 -8932 3) - :(8121 -8122 8601 -8611 -8612 -8613 3)) - #707 #8511 -6112 5 0.0582256 (-8406 -8121 8932 -8931 3 #(8100 -8101)) -c -----Bunker Wall Collimator (BWC)----- -7001 0 (7100 -7103 -7301 -7302 -7303 -7304 3) - #(7101 -7401 7402) -7021 14 -8.92 7100 -7101 7201 7202 7203 7204 3 - #(7100 -7101 -7301 -7302 -7303 -7304) -7022 14 -8.92 7101 -7111 7201 7202 7203 7204 3 - #(7101 -7111 -7301 -7302 -7303 7401) - #(7101 -7111 -7301 -7302 -7402 -7304) -7023 14 -8.92 7111 -7112 7201 7202 7203 7204 3 7110 -7024 14 -8.92 7111 -7112 -7110 3 - #(7111 -7112 -7301 -7302 -7303 7401) - #(7111 -7112 -7301 -7302 -7402 -7304) -7025 14 -8.92 7112 -7102 7201 7202 7203 7204 3 - #(7112 -7103 -7301 -7302 -7303 7401) - #(7112 -7103 -7301 -7302 -7402 -7304) -7026 14 -8.92 7102 -7103 7205 7206 7207 7208 3 - #(7102 -7103 -7301 -7302 -7303 7401) - #(7102 -7103 -7301 -7302 -7402 -7304) -c ---- Outside of Bunker -c -- Shielding Walls -- -c concrete floor -8001 49 0.074981 ((68 -72 8209 -70 -104 -2) - :(72 -8101 8209 -59 8201 8202) - :(8101 -8709 8209 -59 8301 8302)) 3 -c ceiling outside -8002 49 0.074981 ( 72 -8101 8207 8203 8201 8202) -c left outer wall (E01) -8003 49 0.074981 ( 72 -8100 8201 -8205 59 -8207 3) -c Selene 1/2 separator -8004 200 -3.9 (8100 -8101 8201 8202 59 -8207 3) - #(8100 -8101 -8221 -8222 -8223 -8224) -c Optical cave -8005 200 -3.9 (8101 -8709 59 8301 8302 8303 3) - #(8101 -8709 -8311 -8312 -8313 -8314) -c air -8008 4 -1.2E-3 (72 -8100 -8211 -8212 -8213 -8214 3) - #(-8406 -8931) - #(-8405 8406 8401 8402 8403 8404) -8009 4 -1.2E-3 (8101 -8709 -8311 -8312 -8313 -8314 3) - #(8121 -8122 8601 8602 8603 8604) - #(-8121 -8931) - #707 -c B4C concrete inner layers -8011 16 -2.4 (8211 8215 -8213 -8214 72 -8100 3) -8012 16 -2.4 (8205 8202 8213 -8218 72 -8100 3) -8013 16 -2.4 (8212 -8216 -8213 -8214 8111 -8100 3) -8014 16 -2.4 (8202 -8206 8217 -8213 -8214 8111 -8100 3) -c Steel middle layers -8020 71 0.0843223 (8212 8202 -8213 -8214 72 -8111 3) -8021 71 0.0843223 (-8215 8205 -8213 -8214 72 -8100 3) -8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3) -8023 71 0.0843223 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3) -c -c Selene 1 carrier -8301 18 -2.4 -8405 8406 ((8401 8402 8403 8404) - #(-8411 8402 8412 8413)) -c Selene 1 glass (change material!!) -8311 12 -2.51 -8405 8406 -8903 -8904 -8518 - #(-8901 -8902) #(-8515 -8517) 3 -c Copper in wall -8401 14 -8.92 (8100 -8101 -8221 -8222 -8223 -8224 3) - #(8100 -8101 -8501 -8502 -8503 -8504) -c Copper in Selene 1 -8411 14 -8.92 (8111 -8112 8518 8402 8412 8413 3) - :(8111 -8112 8511 8402 -8901 3) - :(8111 -8112 8519 8402 -8515 -8517 3) -8412 14 -8.92 (8113 -8114 8518 8402 8412 8413 3) - :(8113 -8114 8512 8402 -8901 3) - :(8113 -8114 8519 8402 -8515 -8517 3) -8413 14 -8.92 (8115 -8116 8518 8402 8412 8413 3) - :(8115 -8116 8513 8402 -8901 3) - :(8115 -8116 8519 8402 -8515 -8517 3) -c Selene 2 carrier -8501 18 -2.4 (8121 -8122 8601 8602 8603 8604 3) - #(8121 -8122 8601 -8611 -8612 -8613 3) -c Selene 2 glass (change material!!) -8511 12 -2.51 8121 -8122 -8913 -8914 8915 - #(-8911 -8912) #(-8515 -8517) 3 -c Air outside the guide shielding -8601 4 -1.2E-3 72 -8709 -104 -8201 8203 8204 - #708 -8602 4 -1.2E-3 ((112 -8202) - :(8202 -8302 8101)) 72 -8709 8203 8204 - #709 -c Experimental cave -c First largest area: up-stream/down-stream/roof/floor/walls -c then cut corner and add extended roof -8701 49 0.074981 ((8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705)) - :(8709 8710 8702 8703 8711 8712)) - #((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810 - (-8804 : -8805)) - :(-8809 -8810 -8802 -8803 -8811 -8812) - :(-8311 -8312 8709 8809 -8313 -8803)) - #(( 72 -8101 8201 8202 8203 8204 3) - :(8101 -8709 8301 8302 8303 8304 3)) -8702 4 -1.2E-3 ((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810 - (-8804 : -8805)) - :(-8809 -8810 -8802 -8803 -8811 -8812) - :(-8311 -8312 8709 8809 -8313 -8803)) 3 - #6111 #8501 #8511 +710 0 -710 +read file=geometry/estia_inbunker.i echo +read file=geometry/estia_guides.i echo +read file=geometry/estia_cave.i echo c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ c ------------------------------------------------------- c --------------- SURFACE CARDS ------------------------- c ------------------------------------------------------- c modified ESS geometry -read file=e2-ess_surfaces.i echo +read file=geometry/ess_surfaces.i echo c tallies -700 s 117.0 -155.00 13.5 3.0 +c 700 s 117.0 -155.00 13.5 3.0 701 1 s 542.0 -6.3 6.3 3.0 702 1 s 1010 -0.8 0 5.0 703 1 s 1082 -25.0 0 5.0 @@ -470,209 +35,11 @@ c tallies 707 1 s 2190 3.5 0 12.0 708 s 730 -1380 0 30.0 709 s 1020 -1170 0 30.0 +710 s 1328 -1756 -13 12.0 c Estia geometry -read file=e2-estia_surfaces.i echo +read file=geometry/estia_surfaces.i echo -c ------------------------------------------------------- -c --------------- MATERIAL CARDS ------------------------ -c ------------------------------------------------------- -c -c Material : Stainless304 rho=0.0876276 -m3 6000.70c 0.00031864 14028.70c 0.00170336 15031.70c 6.95038e-05 - 16032.70c 4.4752e-05 24000.50c 0.0174813 25055.70c 0.00174159 - 26054.70c 0.00338 26056.70c 0.053455 26057.70c 0.00128218 - 28000.50c 0.00815128 pnlib=70u -mx3:h 6012.70h j j j j j j j j j -c -C NIST STANDARD DRY AIR rho=1.2E-3 -m4 06000.70C 0.000150190 - 07014.70C 0.781575139 - 07015.70C 0.002855327 - 08016.70C 0.210748199 - 18036.70C 0.000015718 - 18038.70C 0.000002952 - 18040.70C 0.004652475 -c Material : Aluminium rho=0.0582256 -m5 13027.70c 0.054381 12024.70c 0.00207 12025.70c 0.0002621 - 12026.70c 0.0002885 14028.70c 0.000217 22048.70c 8.7e-05 - 24052.24c 7.8e-05 24053.24c 1e-05 25055.70c 0.000408 26056.70c - 0.000214 26057.70c 5e-06 29063.70c 4e-05 29065.70c 2e-05 - 30000.70c 0.000145 pnlib=70u -mt5 al27.10t -c -c polyethylene PSI -0.94 g/cm3 -c m6 1001.70c -0.1437 6000.70c -0.8563 -c mt6 poly.01t -c mx6:h j 6012.70h -c borated polyethylene PSI (5% B4C) -0.95 g/cm3 -c m7 1001.70c -0.1183 5010.70c -0.0268 5011.70c -0.1077 6000.70c -0.7472 -c mx7:h j j j 6012.70h -c mt7 poly.01t -c Material : helium rho=2.45e-05 -c PSI Borax for HE-shutter -c EP + B4C 50%-50% (1.8 g/cc) -c m8 5010.70c -0.0782 5011.70c -0.3133 6000.70c -0.4481 7015.70c -0.132 & -c 1001.70c -0.0284 -c mx8:h j j 6012.70h j j -c mt8 lwtr.01t -c Helium -m9 2004.70c 2.45e-05 pnlib=70u -c -c BorAL 1100 alloy (2.7 g/cc) -m10 5010.70c -0.04 5011.70c -0.01 13027.70c -0.9288 & - 14000.70c -0.01 25055.70c -5.0E-3 26000.50c -5.0E-3 & - 29000.50c -1.0E-3 30000.70c -2.0E-4 -c -c Material : H2O rho=0.100283 -m11 1001.70c 0.0668456 1002.70c 1.003e-05 8016.70c 0.0334278 - pnlib=70u -mt11 lwtr.01t -c -c bor-float glass 13% B2O3 2.51 g/cm3 -m12 8016.70c -0.4078 11023.70c -0.096441 14000.70c -0.336563 & - 20000.60c -0.107205 5010.70c -0.0104 5011.70c -0.0416 -c lead (11.34 g/cm3) -m13 82000.50c 1.0 -c Copper density=8.92 g/cm3 -m14 29000.50c 1.0 -c PSI Beton Ba/Fe/Bor 4.9 (4.9% B4C) (4.99 g/cm3) -c m15 1001.70c -3.62e-3 -c 6000.70c -5.23e-3 -c 8016.70c -0.08511 -c 12000.60c -7e-4 -c 13027.70c -4.71e-3 -c 14000.70c -1.983e-2 -c 16000.70c -1.368e-2 -c 20000.60c -2.425e-2 -c 24000.50c -0.13896 -c 25055.70c -1.544e-2 -c 26000.50c -0.52608 -c 27059.50c -0.007 -c 28000.50c -7.72e-2 -c 56138.50c -3.985e-2 -c 5010.70c -7.66e-3 -c 5011.70c -3.068e-2 -c mx15:h j 6012.70h j j j j j j j j j j j j j j -c mt15 poly.01t -c Boron Concrete (5% B4C)(2.4 g/cc) -m16 1001.70c -6.09e-3 - 6000.70c -0.1726 - 8016.70c -0.404 - 12000.60c -0.0322 - 13027.70c -0.0115 - 14000.70c -3.3906e-2 - 20000.60c -0.3222 - 26054.70c -4.22e-4 - 26056.70c -6.92e-3 - 26057.70c -1.62e-4 - 5010.70c -2e-3 - 5011.70c -8e-3 -mx16:h j 6012.70h j j j j j j j j j j -mt16 lwtr.01t -c Mirrobor (1.36 g/cm3) -m17 1001.70c -0.01 - 6000.70c -0.3 - 8016.70c -0.05 - 5010.70c -0.127 - 5011.70c -0.513 -mx17:h j 6012.70h j j j -mt17 poly.01t -c Mineral cast 2% Boron (2.4 g/cm3) -m18 1001.70c -0.05 - 6000.70c -0.26 - 8016.70c -0.2325 - 12000.60c -0.007 - 13027.70c -0.023 - 14000.70c -0.38 - 26000.50c -0.02 - 28000.50c -0.004 - 5010.70c -0.0047 - 5011.70c -0.0188 -mx18:h j 6012.70h j j j j j j j j -mt18 poly.01t -c -c B4C sintered (2.3 g/cc) -m20 5010.70c -0.1841 5011.70c -0.74345 8016.70c -0.07245 -c -c Material : Concrete rho=0.074981 -m49 1001.70c 0.00776555 1002.70c 1.16501e-06 8016.70c 0.0438499 - 11023.70c 0.00104778 12000.60c 0.000148662 13027.70c 0.00238815 - 14028.70c 0.0158026 16032.70c 5.63433e-05 19000.60c 0.000693104 - 20000.60c 0.00291501 26054.70c 1.84504e-05 26056.70c - 0.000286826 26057.70c 6.5671e-06 26058.70c 8.75613e-07 - pnlib=70u -mt49 lwtr.01t -c -c Material : CastIron rho=0.0833854 -m54 6000.70c 0.00636 14028.70c 0.0051404 25055.70c 0.000386 - 26054.70c 0.00419 26056.70c 0.065789 26057.70c 0.00152 - pnlib=70u -mx54:h 6012.70h j j j j j -c -c Material : Steel71 rho=0.0843223 -m71 6000.70c 0.000101187 14028.70c 0.000295128 24000.50c - 0.000168645 25055.70c 0.000505934 26054.70c 0.0048414 26056.70c - 0.0759996 26057.70c 0.00175516 26058.70c 0.00023358 28058.70c - 0.000172218 28060.70c 6.6335e-05 28061.70c 2.88126e-06 - 28062.70c 9.19272e-06 28064.70c 2.33992e-06 29063.70c - 0.000116651 29065.70c 5.19931e-05 pnlib=70u -mx71:h 6012.70h j j j j j j j j j j j j j j -c -c Material : Copper rho=0.0844385 -c (rho=8.91001 g/cc) -m73 29063.70c 0.0583892 29065.70c 0.0260493 pnlib=70u -c -c Material : Be5H2O rho=0.122325 -m118 1001.70c 0.00334228 4009.70c 0.117311 8016.70c 0.00167114 - pnlib=70u -mx118:h j model j -mt118 be.60t lwtr.01t -c -c Material : Be10H2O rho=0.121164 -m119 1001.70c 0.00668456 4009.70c 0.111137 8016.70c 0.00334278 - pnlib=70u -c -mx119:h j model j -mt119 be.60t lwtr.01t -c -c Material : Tungsten_15.1g rho=0.0494621 -m120 74182.70c 0.0131232 74183.70c 0.00708653 74184.70c 0.0151734 - 74186.70c 0.014079 pnlib=70u -c -c Material : Iron_10H2O rho=0.0862739 -m121 1001.70c 0.00668924 8016.70c 0.00334462 26054.70c 0.00445623 - 26056.70c 0.0699533 26057.70c 0.00161553 26058.70c 0.000214997 - pnlib=70u -mt121 lwtr.01t -c -c Material : CarstonConc rho=0.0824808 -m122 1001.70c 0.0229186 1002.70c 1.23628e-06 5010.70c 0.0001187 - 5011.70c 0.00047781 6000.70c 0.00623398 8016.70c 0.0406466 - 14028.70c 0.0049816 14029.70c 0.0002529 14030.70c 0.00016675 - 20040.70c 0.0062902 20042.70c 4.1981e-05 20043.70c 8.759e-06 - 20044.70c 0.0001353 20046.70c 2.5956e-07 20048.70c 1.21339e-05 - 26054.70c 1.1341e-05 26056.70c 0.00017803 26057.70c 4.1115e-06 - 26058.70c 5.47166e-07 pnlib=70u -mx122:h j j j j 6012.70h j j j j j j j j j - j j j j j -mt122 poly.01t -c -c Material : ParaOrtho%0.5 rho=0.041975 -m123 1001.70c 0.000209875 1005.60c 0.0417651 pnlib=70u -mx123:h j 1001.70h -mt123 hpara.10t hortho.10t -c -c Material : ESS Magnadense Heavy Concrete (density -3.9g/cm³) -m200 01001.70c 1.3431e-01 01002.70c 1.5448e-05 08016.70c 5.2962e-01 - 08017.70c 2.0174e-04 11023.70c 5.1116e-03 13027.70c 6.0596e-03 - 14028.70c 3.9343e-02 14029.70c 1.9987e-03 14030.70c 1.3190e-03 - 15031.70c 3.6290e-03 16032.70c 2.2703e-04 16033.70c 1.7929e-06 - 16034.70c 1.0158e-05 16036.70c 2.4026e-08 19039.70c 9.1400e-04 - 19040.70c 1.1512e-07 19041.70c 6.5961e-05 20040.70c 1.2235e-02 - 20042.70c 8.1657e-05 20043.70c 1.7038e-05 20044.70c 2.6327e-04 - 20046.70c 5.0454e-07 20048.70c 2.3602e-05 26054.70c 1.5462e-02 - 26056.70c 2.4273e-01 26057.70c 5.6057e-03 26058.70c 7.4601e-04 -mt200 lwtr.01t +read file=geometry/materials.i echo c c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ c ------------------------------------------------------- @@ -722,68 +89,22 @@ c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ c ----------------------------------------------------------- c ------------------- TALLY CARDS --------------------------- c -fc24 neutron flux at VS -f24:n 704 -e24 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 - 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 - 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 - 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 - 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 - 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 - 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 - 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 - 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 - 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 - 1e-0 1.58 2.51 3.98 6.31 - 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 - 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 - 1e+3 1.58e+3 2.51e+3 -fq24 e f +c read file=tallies/bunker_mesh.i echo +read file=tallies/inbunker_f4.i echo +c read file=tallies/chopper_pit_doses.i echo +c read file=tallies/bunker_wall_f4.i echo +c read file=tallies/e1_e2_f4.i echo +read file=tallies/e1_e2_f5.i echo +c read file=tallies/shutter_f4.i echo +c read file=tallies/shutter_f5.i echo +read file=tallies/optical_cave_f5.i echo c -fc106 deposited energy at chopper disk and spindle -+f106 4121 4131 -e106 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 - 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 - 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 - 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 - 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 - 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 - 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 - 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 - 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 - 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 - 1e-0 1.58 2.51 3.98 6.31 - 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 - 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 - 1e+3 1.58e+3 2.51e+3 -fq106 e f -c -fc6 deposited energy at VS stages x/z/rot -+f6 (5711 5712) (5721 5722) (5731) -e6 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 - 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 - 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 - 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 - 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 - 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 - 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 - 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 - 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 - 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 - 1e-0 1.58 2.51 3.98 6.31 - 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 - 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 - 1e+3 1.58e+3 2.51e+3 -fq6 e f c ========================================== c c ========================================== c dxt:n 730 -1380 0 31 120 1020 -1170 0 31 120 c 779.4 -1093.7 -2.6 13 80 0.1 1e-15 dxc0:n 0 508r 0.25 64r 1 115r -c has to be for specific detector as pd0 does not work properly -c 1-255 256 257 268 259-265 266-268 269-508 509-end -c pd5 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 70r c c 1 2 3 4 5 6-14 15-22 23 24 25 26-34 35 36 37 38 39 c 40-42 43-266 267 268-488 489-494 // @@ -804,7 +125,7 @@ c ------------------------------------------------------- c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 -nps 1e8 +nps 1e9 c histp mode n p n h / d t s a c volume calculated for tally cells (spheres) @@ -815,7 +136,6 @@ imp:n,p 0 0 1 687r cut:h,/,d,t,s,a 1e+08 1e-3 c cut:n j 0.99E-6 cut:n j 1e-11 -c esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000 cut:p 1e+08 0.5 phys:n 3000 0 phys:p 100 @@ -824,11 +144,11 @@ phys:h 3000 lca 2 1 1 0023 1 1 0 1 1 0 lea 1 4 1 0 1 0 0 1 c * Source for simulation from proton beam (needs different WW and ESplit) -c read file=source_proton.i echo +c read file=sources/proton.i echo c * Source Term from Valentina for 2m adapted to size -c read file=e2-source_2m.i echo -c * McStas source terms for gamma doeses (needs no WW) -read file=mcstas_source_vs.i echo +c read file=sources/e2_2m.i echo +c * McStas source terms for gamma doeses +read file=sources/mcstas_source_vs.i echo prdmp 5e7 -60 1 2 5e6 print 30 32 $ 120 128 150 160 161 162 c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ diff --git a/e2-ess_cells.i b/geometry/ess_cells.i similarity index 99% rename from e2-ess_cells.i rename to geometry/ess_cells.i index fc11dee..2bbe7b5 100644 --- a/e2-ess_cells.i +++ b/geometry/ess_cells.i @@ -45,7 +45,7 @@ c Main ESS plant geometry without the instrument 32 5 0.0582256 ( -34 : 35 : 36 ) -37 34 -38 33 119 0.121164 ( -34 : 37 : 38 ) -37 34 -19 34 5 0.0582256 ( -34 : 37 : 19 ) -37 320 -21 -35 0 -61 3 21 -42 41 -40 39 #700 #1135 #1136 +35 0 -61 3 21 -42 41 -40 39 #1135 #1136 36 5 0.0582256 -61 3 21 ( -39 : 40 : -41 : 42 ) -46 45 -44 43 ( -142 : 332 : 556 : -325 ) 37 0 -61 -3 21 -320 33 -48 47 diff --git a/e2-ess_surfaces.i b/geometry/ess_surfaces.i similarity index 100% rename from e2-ess_surfaces.i rename to geometry/ess_surfaces.i diff --git a/geometry/estia_cave.i b/geometry/estia_cave.i new file mode 100644 index 0000000..c30179c --- /dev/null +++ b/geometry/estia_cave.i @@ -0,0 +1,16 @@ +c Experimental cave +c First largest area: up-stream/down-stream/roof/floor/walls +c then cut corner and add extended roof +8701 49 0.074981 ((8709 8710 8701 8703 8706 8707 8708 8709 8710 (8704 : 8705)) + :(8709 8710 8702 8703 8711 8712)) + #((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810 + (-8804 : -8805)) + :(-8809 -8810 -8802 -8803 -8811 -8812) + :(-8311 -8312 8709 8809 -8313 -8803)) + #(( 72 -8101 8201 8202 8203 8204 3) + :(8101 -8709 8301 8302 8303 8304 3)) +8702 4 -1.2E-3 ((-8809 -8810 -8801 -8803 -8806 -8807 -8808 -8809 -8810 + (-8804 : -8805)) + :(-8809 -8810 -8802 -8803 -8811 -8812) + :(-8311 -8312 8709 8809 -8313 -8803)) 3 + #6111 #8501 #8511 diff --git a/geometry/estia_guides.i b/geometry/estia_guides.i new file mode 100644 index 0000000..d0b87a6 --- /dev/null +++ b/geometry/estia_guides.i @@ -0,0 +1,112 @@ +c geometry from bunker wall until experimental cave entrance +c -----Wall Cut-Out--- +6001 4 -1.2E-3 ((68 -6100 -6201 -6202 -6203 -6204) + :(6100 -6101 -6301 -6302 -6303 -6304)) 3 + #(-5102 5501 5502 -5505 5504) + #(5102 -7100 -5531) +6002 4 -1.2E-3 ((6101 -6102 -6401 -6402 -6403 -6404) + :(6102 -6103 -6501 -6502 -6503 -6504)) 3 + #(5102 -7100 -5531) + #((7100 -7102 7201 7202 7203 7204) + :(7102 -7103 7205 7206 7207 7208)) + #(7103 -6103 -8921) +6003 4 -1.2E-3 (6103 -72 -6601 -6602 -6603 -6604 -6605) 3 + #(7103 8405 -8921) + #(7102 -7103 7205 7206 7207 7208) + #(-8405 8401 8402 8403 8404) +c Vacuum before and inside Selene Guide 1 +6101 0 ((7103 8405 -8922 3) + :(-8405 8406 -8411 8412 8413 8402 3)) + #706 #8311 #8411 #8412 #8413 +6102 5 0.0582256 (7103 8405 8922 -8921 3) +c Vaccum after Selene 1, until and inside Selen 2 +6111 0 ((-8406 -8100 -8932 3) + :(8100 -8101 -8501 -8502 -8503 -8504 3) + :(8101 -8121 -8932 3) + :(8121 -8122 8601 -8611 -8612 -8613 3)) + #707 #8511 +6112 5 0.0582256 (-8406 -8121 8932 -8931 3 #(8100 -8101)) +c -----Bunker Wall Collimator (BWC)----- +7001 0 (7100 -7103 -7301 -7302 -7303 -7304 3) + #(7101 -7401 7402) +7021 14 -8.92 7100 -7101 7201 7202 7203 7204 3 + #(7100 -7101 -7301 -7302 -7303 -7304) +7022 14 -8.92 7101 -7111 7201 7202 7203 7204 3 + #(7101 -7111 -7301 -7302 -7303 7401) + #(7101 -7111 -7301 -7302 -7402 -7304) +7023 14 -8.92 7111 -7112 7201 7202 7203 7204 3 7110 +7024 14 -8.92 7111 -7112 -7110 3 + #(7111 -7112 -7301 -7302 -7303 7401) + #(7111 -7112 -7301 -7302 -7402 -7304) +7025 14 -8.92 7112 -7102 7201 7202 7203 7204 3 + #(7112 -7103 -7301 -7302 -7303 7401) + #(7112 -7103 -7301 -7302 -7402 -7304) +7026 14 -8.92 7102 -7103 7205 7206 7207 7208 3 + #(7102 -7103 -7301 -7302 -7303 7401) + #(7102 -7103 -7301 -7302 -7402 -7304) +c ---- Outside of Bunker +c -- Shielding Walls -- +c concrete floor +8001 49 0.074981 ((68 -72 8209 -70 -104 -2) + :(72 -8101 8209 -59 8201 8202) + :(8101 -8709 8209 -59 8301 8302)) 3 +c ceiling outside +8002 49 0.074981 ( 72 -8101 8207 8203 8201 8202) +c left outer wall (E01) +8003 49 0.074981 ( 72 -8100 8201 -8205 59 -8207 3) +c Selene 1/2 separator +8004 200 -3.9 (8100 -8101 8201 8202 59 -8207 3) + #(8100 -8101 -8221 -8222 -8223 -8224) +c Optical cave +8005 200 -3.9 (8101 -8709 59 8301 8302 8303 3) + #(8101 -8709 -8311 -8312 -8313 -8314) +c air +8008 4 -1.2E-3 (72 -8100 -8211 -8212 -8213 -8214 3) + #(-8406 -8931) + #(-8405 8406 8401 8402 8403 8404) +8009 4 -1.2E-3 (8101 -8709 -8311 -8312 -8313 -8314 3) + #(8121 -8122 8601 8602 8603 8604) + #(-8121 -8931) + #707 #710 +c B4C concrete inner layers +8011 16 -2.4 (8211 8215 -8213 -8214 72 -8100 3) +8012 16 -2.4 (8205 8202 8213 -8218 72 -8100 3) +8013 16 -2.4 (8212 -8216 -8213 -8214 8111 -8100 3) +8014 16 -2.4 (8202 -8206 8217 -8213 -8214 8111 -8100 3) +c Steel middle layers +8020 71 0.0843223 (8212 8202 -8213 -8214 72 -8111 3) +8021 71 0.0843223 (-8215 8205 -8213 -8214 72 -8100 3) +8022 71 0.0843223 (8205 8202 8218 -8207 72 -8100 3) +8023 71 0.0843223 (8216 8202 (-8217:8206) -8213 -8214 8111 -8100 3) +c +c Selene 1 carrier +8301 18 -2.4 -8405 8406 ((8401 8402 8403 8404) + #(-8411 8402 8412 8413)) +c Selene 1 glass (change material!!) +8311 12 -2.51 -8405 8406 -8903 -8904 -8518 + #(-8901 -8902) #(-8515 -8517) 3 +c Copper in wall +8401 14 -8.92 (8100 -8101 -8221 -8222 -8223 -8224 3) + #(8100 -8101 -8501 -8502 -8503 -8504) +c Copper in Selene 1 +8411 14 -8.92 (8111 -8112 8518 8402 8412 8413 3) + :(8111 -8112 8511 8402 -8901 3) + :(8111 -8112 8519 8402 -8515 -8517 3) +8412 14 -8.92 (8113 -8114 8518 8402 8412 8413 3) + :(8113 -8114 8512 8402 -8901 3) + :(8113 -8114 8519 8402 -8515 -8517 3) +8413 14 -8.92 (8115 -8116 8518 8402 8412 8413 3) + :(8115 -8116 8513 8402 -8901 3) + :(8115 -8116 8519 8402 -8515 -8517 3) +c Selene 2 carrier +8501 18 -2.4 (8121 -8122 8601 8602 8603 8604 3) + #(8121 -8122 8601 -8611 -8612 -8613 3) +c Selene 2 glass (change material!!) +8511 12 -2.51 8121 -8122 -8913 -8914 8915 + #(-8911 -8912) #(-8515 -8517) 3 +c Air outside the guide shielding +8601 4 -1.2E-3 72 -8709 -104 -8201 8203 8204 + #708 +8602 4 -1.2E-3 ((112 -8202) + :(8202 -8302 8101)) 72 -8709 8203 8204 + #709 diff --git a/geometry/estia_inbunker.i b/geometry/estia_inbunker.i new file mode 100644 index 0000000..f7ac4b3 --- /dev/null +++ b/geometry/estia_inbunker.i @@ -0,0 +1,312 @@ +c Estia geometry inside monolith and bunker +c elliptic guide +c NBOA beam +1001 9 2.45e-05 1100 -1101 -1301 -1302 -1303 -1304 +1002 9 2.45e-05 1101 -1102 (-1321:-1322) -1311 -1501 -1502 +1003 9 2.45e-05 1102 -1103 (-1331:-1332) -1311 -1501 -1502 +1004 9 2.45e-05 1103 -1104 (-1341:-1342) -1311 -1501 -1502 +1005 9 2.45e-05 1104 -1105 (-1351:-1352) -1311 -1501 -1502 +1006 9 2.45e-05 1105 -1106 (-1361:-1362) -1311 -1502 -1501 +1007 9 2.45e-05 1106 -1107 (-1371:-1372) -1311 -1502 -1501 +1008 9 2.45e-05 1107 -1108 (-1381:-1382) -1311 -1502 -1501 +1011 9 2.45e-05 1100 -1101 #(1201 1202 1203 1204) 1305 1306 1307 1308 +1012 9 2.45e-05 ((1101 -1115 1213 -1507 -1508) + :(1115 -1105 1214 -1509 -1510)) #(1211 -1503 -1504) +1013 9 2.45e-05 1105 -1107 #(1212 -1505 -1506) 1214 -1509 -1510 +1014 9 2.45e-05 ((1117 -1108 1213 -1507 -1508) + :(1107 -1117 1214 -1509 -1510) + ) #(1211 -1503 -1504) +c NBOA outside +1021 14 -8.92 (1100 -1101 1201 1202 1203 1204 + (1301 : 1302 : 1303 : 1304)) +1022 14 -8.92 (1101 -1103 1211 -1503 -1504) + (1311:1501:1502) +1023 14 -8.92 (1103 -1104 1211 -1503 -1504) + (1311:1501:1502) +1024 14 -8.92 (1104 -1105 1211 -1503 -1504) + (1311:1501:1502) +1025 14 -8.92 (1105 -1107 1212 -1505 -1506) + (1311:1501:1502) +1026 14 -8.92 (1107 -1108 1211 -1504 -1503) + (1311:1501:1502) +c NBOA wedges +1032 14 -8.92 ((1101 -1102 1321 1322) + :(1102 -1103 1331 1332) + ) -1502 -1311 +1033 14 -8.92 1103 -1104 1341 1342 -1502 -1311 +1034 14 -8.92 1104 -1105 1351 1352 -1502 -1311 +1035 14 -8.92 ((1105 -1106 1361 1362) + :(1106 -1107 1371 1372) + ) -1502 -1311 +1036 14 -8.92 1107 -1108 1381 1382 -1502 -1311 +c Al window and He in front of NBOA +1041 5 0.0582256 61 -1110 1305 1306 1307 1308 +1042 9 2.45e-05 1110 -1100 1305 1306 1307 1308 +c Al window and He behind NBOA +1051 9 2.45e-05 1108 -1118 1213 -1508 -1507 +1052 5 0.0582256 1118 -64 1213 -1508 -1507 +c extended cut-out in reflector +1135 0 -61 3 58 1205 1206 1207 1208 +1136 5 0.0582256 -61 3 #(41 -42) 58 1201 1202 1203 1204 #1135 +c split of wedge of monolith components for increased importance +1236 5 0.0582256 ( 142 -332 -556 325 ) ( -1201 : -1202 : -1203 : -1204 : -58) + ( -39 : 40 : -41 : 42 ) -61 3 21 -46 45 -44 43 + dxc1:n=0.001 +1240 121 0.0862739 ( 142 -332 -556 325 ) + ( -43 : -45 : 46 : 44 ) -58 ( -53 : 54 : 55 ) + dxc1:n=0.001 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1 +1241 121 0.0862739 ( 142 -332 -556 325 ) ( -1201 : -1202 : -1203 : -1204 ) + -61 58 ( -43 : -45 : 46 : 44 ) + dxc1:n=0.1 dxc2:n=0.1 dxc3:n=0.1 dxc4:n=0.1 +c partitioned air in bunker +1255 4 -1.2E-3 3 -68 -2 -101 66 -67 64 112 + ( 329 : 4615 : -342 : 353 : 4100 ) + ( -66 : 4615 : -342 : 353 : -3111 ) + #(-1605 -1623) #(-1622 1623 -1624) + #(-1639 1624 -1637) #(-1658 1637 -1660) + #3137 +1256 4 -1.2E-3 ( -329 -4615 342 -353 ) 64 -3110 + #( 2100 -2101 2201 2202 2203 2204 ) +1257 4 -1.2E-3 (( -329 -4615 342 -353 ) 3110 -4805 + #((3110 -3111 -3541) + :(3111 -4100 -3531))) + #(-3440 -3404 3446) + #(-3441 -3326 3446) + #(-3442 -3326 3446) +1258 4 -1.2E-3 (( 66 3111 329 #(4811 -4812 -4820) ) + :( -329 -4615 4805 + #( 5101 4604 5501 5502 -5505 5504 ) ) + #((4811 -4812 -4820) + :(-4814 -4815 -4602 4813 4805 -4101) + :(4812 -4602 ((4601 -4603 -4818 -4819):(4603 -4604 -4818))) + :(4602 4611 4612 4613 4614 -4615)) + ) 342 -353 -68 #4201 #4202 #4203 #4204 #4131 +c part around inserts E01,E02,E03 +1489 54 0.0833854 109 556 -334 325 -326 -64 330 3 + ( -153 : 154 : 159 ) -329 +1490 54 0.0833854 -556 -1107 + ( -336 : -330 : -341 : -459 : 64 : 460 : 342 ) + ( -348 : -330 : -353 : -459 : 64 : 460 : 354 ) + -334 -329 -326 325 -64 330 3 ( -153 : 154 : 159 ) + ( -330 : -569 : -459 : 64 : 460 : 570 : 564 ) +1491 54 0.0833854 -556 1107 + ( -336 : -330 : -341 : -459 : 64 : 460 : 342 ) + ( -348 : -330 : -353 : -459 : 64 : 460 : 354 ) + -334 -329 -326 325 -64 330 3 ( -153 : 154 : 159 ) + ( -330 : -569 : -459 : 64 : 460 : 570 : 564 ) +1494 54 0.0833854 -556 -1103 + ( -336 : -61 : -337 : -575 : 457 : 330 : 338 : ( 575 565 -457 ) ) + ( -348 : -61 : -349 : -575 : 457 : 330 : 350 : ( 575 554 -457 ) ) + ( -153 : 154 : 159 ) 3 61 -330 -326 325 335 -332 + ( -61 : -565 : -1001 : 457 : 330 : 566 : 349 : 564 : ( -457 349 575 ) ) +1495 54 0.0833854 -556 1103 + ( -336 : -61 : -337 : -575 : 457 : 330 : 338 : ( 575 565 -457 ) ) + ( -348 : -61 : -349 : -575 : 457 : 330 : 350 : ( 575 554 -457 ) ) + ( -153 : 154 : 159 ) 3 61 -330 -326 325 335 -332 + ( -61 : -565 : -1001 : 457 : 330 : 566 : 564 : ( -457 349 575 ) ) +1586 3 0.0876276 ( -554 : -575 : 457 ) -330 #589 + -458 -345 351 -352 1103 348 +1588 3 0.0876276 -64 -332 459 355 -356 330 + 348 1107 #589 +1591 3 0.0876276 ( -349 : -1001 : 457 ) -330 + -458 1000 567 -568 -564 -2 1103 + (-1101 : -1213 : 1508 : 1507) +1593 3 0.0876276 -64 -332 459 571 -572 330 -564 -2 1107 + ( -1213 : 1508 : 1507) + #(1107 -1117 1214 -1509 -1510) +1596 3 0.0876276 ( -565 : -575 : 457 ) -330 + -458 -345 339 -340 1103 336 +1598 3 0.0876276 -64 -332 459 343 -344 330 336 1107 +c SKADI in-bunker components +c collimators +1901 0 64 3 ((-1604 -1623):(-1622 1623 -1624) + :(-1639 1624 -1637):(-1658 1637 -1660)) + #1940 #1945 #1957 #1958 #1963 #1967 #1982 #1983 +1940 5 0.0582256 1604 -1605 -1623 64 3 +1945 73 0.0844385 ( -1594 : 1595 : -1615 : 1616 ) -1607 1606 -1618 1617 -1599 + 1598 1614 +1957 0 ( -1690 : 1691 ) ( -1647 : 1646 : 1645 ) ( -1643 : 1644 : 1645 ) -1637 + 1636 -1638 ( -1648 : -1653 : -1598 : -1650 : 1599 : 1654 : 1649 ) +1958 5 0.0582256 1638 -1639 -1641 1640 +1963 73 0.0844385 ( 1595 : -1594 : -1651 : 1652 ) -1649 1648 -1654 1653 1598 + -1599 1650 +1967 5 0.0582256 1657 -1658 -1660 1659 +1982 73 0.0844385 -1622 ( -1631 : -1598 : -1632 : 1633 : 1599 ) + -1624 1623 +1983 73 0.0844385 -1638 ( -1650 : -1598 : -1653 : 1654 : 1599 ) + -1691 1690 +c +c BBG beam top/beam bottom/outside +2001 9 2.45e-05 2100 -2101 -2301 -2302 -2303 -2304 3 #701 +2011 9 2.45e-05 2100 -2101 -2301 -2302 -2305 -2306 3 +2021 14 -8.92 2100 -2101 2201 2202 2203 2204 3 + (2301:2302:2303:2306:(2304 2305)) +c NFGA beam top +3001 0 3100 -3101 -3301 -3311 -3501 -3502 +3002 0 3101 -3102 -3302 -3312 -3501 -3502 +3003 0 3102 -3103 -3302 -3313 -3501 -3502 +3004 0 3103 -3104 -3302 -3314 -3501 -3502 +3005 0 3104 -3105 -3302 -3315 -3501 -3502 +3006 0 3105 -3106 -3302 -3316 -3501 -3502 +3007 0 3106 -3107 -3302 -3317 -3501 -3502 +c NFGA beam bottom +3011 0 3100 -3101 -3301 -3321 -3501 -3502 +3012 0 3101 -3102 -3302 -3322 -3501 -3502 +3013 0 3102 -3103 -3302 -3323 -3501 -3502 +3014 0 3103 -3104 -3302 -3324 -3501 -3502 +3015 0 3104 -3105 -3302 -3325 -3501 -3502 +3016 0 3105 -3106 -3302 -3326 -3501 -3502 +3017 0 3106 -3107 -3302 -3327 -3501 -3502 +c NFGA outside +3021 5 0.0582256 3100 -3101 -3503 -3504 3201 + #(-3301 #(3311 3321) -3501 -3502) +3022 5 0.0582256 3101 -3102 -3503 -3504 3202 + #(-3302 #(3312 3322) -3501 -3502) +3023 5 0.0582256 3102 -3103 -3503 -3504 3202 + #(-3302 #(3313 3323) -3501 -3502) +3024 5 0.0582256 3103 -3104 -3503 -3504 3202 + #(-3302 #(3314 3324) -3501 -3502) +3025 5 0.0582256 3104 -3105 -3503 -3504 3202 + #(-3302 #(3315 3325) -3501 -3502) +3026 5 0.0582256 3105 -3106 -3503 -3504 3202 + #(-3302 #(3316 3326) -3501 -3502) +3027 5 0.0582256 3106 -3107 -3503 -3504 3202 + #(-3302 #(3317 3327) -3501 -3502) +c NFSA Vaccum tank and connection to collimator +3122 0 ((((3112 -3131 -3522 3) + :(3131 -4805 -3532 3)) #( + (3100 -3101 3201 -3504 -3503) + :(3101 -3107 3202 -3504 -3503))) + :(-3434 -3404 -3541) + :(-3435 -3326 -3541) + :(-3436 -3326 -3541)) + #702 #3151 +3123 5 0.0582256 ((3110 -3111 3522 -3521) + :(3110 -3112 -3521) + :(3113 -3111 3532 -3521) + :(3111 -4805 3532 -3531)) 3 + #(-3431 -3326) + #(-3432 -3326) + #(-3433 -3326) +3124 0 3446 3541 3437 -3440 -3326 #3131 #3134 +3125 0 3446 3541 3438 -3441 -3326 #3132 #3135 +3126 0 3446 3541 3439 -3442 -3326 #3133 #3136 +3131 5 0.0582256 3434 -3431 3522 3443 -3326 +3132 5 0.0582256 3435 -3432 3522 3443 -3326 +3133 5 0.0582256 3436 -3433 3522 3443 -3326 +3134 71 0.0843223 ((-3431 3446 -3443) + :(-3440 3444 -3443) + :(-3440 3446 -3445)) 3434 +3135 71 0.0843223 ((-3432 3446 -3443) + :(-3441 3444 -3443) + :(-3441 3446 -3445)) 3435 +3136 71 0.0843223 ((-3433 3446 -3443) + :(-3442 3444 -3443) + :(-3442 3446 -3445)) 3436 +3137 49 0.074981 3421 3422 -3423 3424 3425 3426 +c B-Al plate +3141 10 -2.7 3131 -3113 3532 -3522 3 +c MirroBor outside vessel +3142 17 -1.36 3110 -3111 3521 -3541 3 + #(-3431 -3326) + #(-3432 -3326) + #(-3433 -3326) +c NFSA Guide support +3151 71 0.0843223 (3120 -3121 (3401 3402 3403 3404 3) + #(3411 3412 -3413 3414)) + :(3446 -3404 -3437) + :(3446 -1000 -3438) + :(3446 -1000 -3439) +c -----Collimators---- +c CPC1 - Collimator before VS +4001 0 -4301 -4302 -4303 -4304 4805 -4101 -2 +4010 10 -2.7 -4814 -4815 -4602 4813 4805 -4807 #4001 +4011 14 -8.92 -4814 -4815 -4602 4813 4807 4819 #4001 +4012 14 -8.92 4201 4202 -4705 4704 -4708 -4101 #4001 +c Vacuum tank around VS +4101 0 (4705 4616 4617 4610 4619 -4618) #4121 +4102 0 (4704 -4705 ((-4706 -4707 -4708 -4603):(4603 -4702))) + #(-4101 4704 4201 4202) + #(5604 -5601 -5611) + #5701 #5702 #5711 #5712 #5751 #5761 #4121 + #5021 #703 #704 +4103 0 (4703 -4704 -4701) + #(-4101 4704 4201 4202) + #5711 #5712 #5721 #5722 #5731 +4104 0 (4603 -4604 4702 -5301 -5302 -5303 -5304) +4105 0 (5604 -5601 -5612) + #(5604 -5601 (5612:(5613 5616):(5614 5615) + :(5614 5616):(-5614 -5618):(-5616 5617)) -5611) + #5701 #5702 +4111 5 0.0582256 ((4602 4611 4612 4613 4614 -4615): + (4705 -4602 ((4601 -4603 -4818 -4819) + :(4603 -4604 -4818)))) + #(4616 4617 4610 4619 -4618) +4112 5 0.0582256 (4704 -4705 ((4601 -4603 -4818 -4819):(4603 -4604 -4818))) + #((-4706 -4707 -4708 -4603):(4603 -4702)) + #(-5301 -5302 -5303 -5304) + #(-4301 -4302 -4303 -4304) +4113 5 0.0582256 (4812 -4704 ((4601 -4603 -4818 -4819):(4603 -4604 -4818))) + #(4703 -4704 -4701) +4114 5 0.0582256 4811 -4812 -4820 #(4703 -4704 -4701) +c chopper disk +4121 5 0.0582256 -4621 4622 -4623 +c chopper spindle +4131 3 0.0876276 -4612 4631 -4632 -4604 +c Concrete base (5% B4C) +4201 16 -2.4 (4801 4802 4804 -4811 4805 4806) +4202 16 -2.4 (4801 4802 4811 -4812 4805 4806) + #(-4820) + #(-4817 -4814 -4815) +4203 16 -2.4 (4801 4802 4812 -4813 4805 4806) + #(-4818 -4819) + #(-4816 -4814 -4815) +4204 16 -2.4 (4801 4802 4813 4803 4805 4806) + #(-4818 -4819) + #(-4814 -4815) +c CPC2 with Cu nose +5001 0 ((5101 4604 -5102 -5301 -5302 -5303 -5304 3) + :(5102 -7100 -5532 3)) #705 +5002 5 0.0582256 (5102 -7100 5532 -5531 3) +5021 14 -8.92 5401 5402 -5403 5404 5405 5406 + #(-5301 -5302 -5303 -5304) +5022 14 -8.92 5101 4604 -5102 5501 5502 -4602 5504 3 + #5001 +5023 71 0.0843223 5101 4604 -5102 5501 5502 4602 -5505 3 + #5001 +c VS Lead shielding (FIX MATERIALS) +c boral part top and bottom +5601 10 -2.7 5602 -5601 (5612:(5613 5616):(5614 5615) + :(5614 5616):(-5614 -5618):(-5616 5617)) + -5611 +5602 10 -2.7 5604 -5603 (5612:(5613 5616):(5614 5615) + :(5614 5616):(-5614 -5618):(-5616 5617)) + -5611 +c lead part +5603 13 -11.34 5603 -5602 (5612 + :(5613 5616):(5614 5615) + :(5614 5616) + :(-5614 -5618):(-5616 5617)) + -5611 +c VS holder stands +5701 5 0.0582256 -5701 5704 -5703 +5702 5 0.0582256 -5702 5704 -5703 +c simplified model of stages (hollow cylinders) +c x-stages +5711 5 0.0582256 5732 -5731 -5711 + #(5742 -5741 -5721) +5712 5 0.0582256 5732 -5731 -5712 + #(5742 -5741 -5722) +c z-stages +5721 5 0.0582256 5733 -5732 -5713 + #(5744 -5743 -5723) +5722 5 0.0582256 5733 -5732 -5714 + #(5744 -5743 -5724) +c omega-stage +5731 5 0.0582256 5734 -5733 -5715 + #(5746 -5745 -5725) +c B4C blades +5751 20 -2.3 ((-5751 -5754 5755) + :(-5752 -5755 5703)) -5753 5756 -5757 +5761 20 -2.3 ((-5762 -5764 5765) + :(-5763 -5765 5703)) -5761 -5766 5767 diff --git a/e2-estia_surfaces.i b/geometry/estia_surfaces.i similarity index 100% rename from e2-estia_surfaces.i rename to geometry/estia_surfaces.i diff --git a/geometry/materials.i b/geometry/materials.i new file mode 100644 index 0000000..1927c59 --- /dev/null +++ b/geometry/materials.i @@ -0,0 +1,200 @@ +c ------------------------------------------------------- +c --------------- MATERIAL CARDS ------------------------ +c ------------------------------------------------------- +c +c Material : Stainless304 rho=0.0876276 +m3 6000.70c 0.00031864 14028.70c 0.00170336 15031.70c 6.95038e-05 + 16032.70c 4.4752e-05 24000.50c 0.0174813 25055.70c 0.00174159 + 26054.70c 0.00338 26056.70c 0.053455 26057.70c 0.00128218 + 28000.50c 0.00815128 pnlib=70u +mx3:h 6012.70h j j j j j j j j j +c +C NIST STANDARD DRY AIR rho=1.2E-3 +m4 06000.70C 0.000150190 + 07014.70C 0.781575139 + 07015.70C 0.002855327 + 08016.70C 0.210748199 + 18036.70C 0.000015718 + 18038.70C 0.000002952 + 18040.70C 0.004652475 +c Material : Aluminium rho=0.0582256 +m5 13027.70c 0.054381 12024.70c 0.00207 12025.70c 0.0002621 + 12026.70c 0.0002885 14028.70c 0.000217 22048.70c 8.7e-05 + 24052.24c 7.8e-05 24053.24c 1e-05 25055.70c 0.000408 26056.70c + 0.000214 26057.70c 5e-06 29063.70c 4e-05 29065.70c 2e-05 + 30000.70c 0.000145 pnlib=70u +mt5 al27.10t +c +c polyethylene PSI -0.94 g/cm3 +c m6 1001.70c -0.1437 6000.70c -0.8563 +c mt6 poly.01t +c mx6:h j 6012.70h +c borated polyethylene PSI (5% B4C) -0.95 g/cm3 +c m7 1001.70c -0.1183 5010.70c -0.0268 5011.70c -0.1077 6000.70c -0.7472 +c mx7:h j j j 6012.70h +c mt7 poly.01t +c Material : helium rho=2.45e-05 +c PSI Borax for HE-shutter +c EP + B4C 50%-50% (1.8 g/cc) +c m8 5010.70c -0.0782 5011.70c -0.3133 6000.70c -0.4481 7015.70c -0.132 & +c 1001.70c -0.0284 +c mx8:h j j 6012.70h j j +c mt8 lwtr.01t +c Helium +m9 2004.70c 2.45e-05 pnlib=70u +c +c BorAL 1100 alloy (2.7 g/cc) +m10 5010.70c -0.04 5011.70c -0.01 13027.70c -0.9288 & + 14000.70c -0.01 25055.70c -5.0E-3 26000.50c -5.0E-3 & + 29000.50c -1.0E-3 30000.70c -2.0E-4 +c +c Material : H2O rho=0.100283 +m11 1001.70c 0.0668456 1002.70c 1.003e-05 8016.70c 0.0334278 + pnlib=70u +mt11 lwtr.01t +c +c bor-float glass 13% B2O3 2.51 g/cm3 +m12 8016.70c -0.4078 11023.70c -0.096441 14000.70c -0.336563 & + 20000.60c -0.107205 5010.70c -0.0104 5011.70c -0.0416 +c lead (11.34 g/cm3) +m13 82000.50c 1.0 +c Copper density=8.92 g/cm3 +m14 29000.50c 1.0 +c PSI Beton Ba/Fe/Bor 4.9 (4.9% B4C) (4.99 g/cm3) +c m15 1001.70c -3.62e-3 +c 6000.70c -5.23e-3 +c 8016.70c -0.08511 +c 12000.60c -7e-4 +c 13027.70c -4.71e-3 +c 14000.70c -1.983e-2 +c 16000.70c -1.368e-2 +c 20000.60c -2.425e-2 +c 24000.50c -0.13896 +c 25055.70c -1.544e-2 +c 26000.50c -0.52608 +c 27059.50c -0.007 +c 28000.50c -7.72e-2 +c 56138.50c -3.985e-2 +c 5010.70c -7.66e-3 +c 5011.70c -3.068e-2 +c mx15:h j 6012.70h j j j j j j j j j j j j j j +c mt15 poly.01t +c Boron Concrete (5% B4C)(2.4 g/cc) +m16 1001.70c -6.09e-3 + 6000.70c -0.1726 + 8016.70c -0.404 + 12000.60c -0.0322 + 13027.70c -0.0115 + 14000.70c -3.3906e-2 + 20000.60c -0.3222 + 26054.70c -4.22e-4 + 26056.70c -6.92e-3 + 26057.70c -1.62e-4 + 5010.70c -2e-3 + 5011.70c -8e-3 +mx16:h j 6012.70h j j j j j j j j j j +mt16 lwtr.01t +c Mirrobor (1.36 g/cm3) +m17 1001.70c -0.01 + 6000.70c -0.3 + 8016.70c -0.05 + 5010.70c -0.127 + 5011.70c -0.513 +mx17:h j 6012.70h j j j +mt17 poly.01t +c Mineral cast 2% Boron (2.4 g/cm3) +m18 1001.70c -0.05 + 6000.70c -0.26 + 8016.70c -0.2325 + 12000.60c -0.007 + 13027.70c -0.023 + 14000.70c -0.38 + 26000.50c -0.02 + 28000.50c -0.004 + 5010.70c -0.0047 + 5011.70c -0.0188 +mx18:h j 6012.70h j j j j j j j j +mt18 poly.01t +c +c B4C sintered (2.3 g/cc) +m20 5010.70c -0.1841 5011.70c -0.74345 8016.70c -0.07245 +c +c Material : Concrete rho=0.074981 +m49 1001.70c 0.00776555 1002.70c 1.16501e-06 8016.70c 0.0438499 + 11023.70c 0.00104778 12000.60c 0.000148662 13027.70c 0.00238815 + 14028.70c 0.0158026 16032.70c 5.63433e-05 19000.60c 0.000693104 + 20000.60c 0.00291501 26054.70c 1.84504e-05 26056.70c + 0.000286826 26057.70c 6.5671e-06 26058.70c 8.75613e-07 + pnlib=70u +mt49 lwtr.01t +c +c Material : CastIron rho=0.0833854 +m54 6000.70c 0.00636 14028.70c 0.0051404 25055.70c 0.000386 + 26054.70c 0.00419 26056.70c 0.065789 26057.70c 0.00152 + pnlib=70u +mx54:h 6012.70h j j j j j +c +c Material : Steel71 rho=0.0843223 +m71 6000.70c 0.000101187 14028.70c 0.000295128 24000.50c + 0.000168645 25055.70c 0.000505934 26054.70c 0.0048414 26056.70c + 0.0759996 26057.70c 0.00175516 26058.70c 0.00023358 28058.70c + 0.000172218 28060.70c 6.6335e-05 28061.70c 2.88126e-06 + 28062.70c 9.19272e-06 28064.70c 2.33992e-06 29063.70c + 0.000116651 29065.70c 5.19931e-05 pnlib=70u +mx71:h 6012.70h j j j j j j j j j j j j j j +c +c Material : Copper rho=0.0844385 +c (rho=8.91001 g/cc) +m73 29063.70c 0.0583892 29065.70c 0.0260493 pnlib=70u +c +c Material : Be5H2O rho=0.122325 +m118 1001.70c 0.00334228 4009.70c 0.117311 8016.70c 0.00167114 + pnlib=70u +mx118:h j model j +mt118 be.60t lwtr.01t +c +c Material : Be10H2O rho=0.121164 +m119 1001.70c 0.00668456 4009.70c 0.111137 8016.70c 0.00334278 + pnlib=70u +c +mx119:h j model j +mt119 be.60t lwtr.01t +c +c Material : Tungsten_15.1g rho=0.0494621 +m120 74182.70c 0.0131232 74183.70c 0.00708653 74184.70c 0.0151734 + 74186.70c 0.014079 pnlib=70u +c +c Material : Iron_10H2O rho=0.0862739 +m121 1001.70c 0.00668924 8016.70c 0.00334462 26054.70c 0.00445623 + 26056.70c 0.0699533 26057.70c 0.00161553 26058.70c 0.000214997 + pnlib=70u +mt121 lwtr.01t +c +c Material : CarstonConc rho=0.0824808 +m122 1001.70c 0.0229186 1002.70c 1.23628e-06 5010.70c 0.0001187 + 5011.70c 0.00047781 6000.70c 0.00623398 8016.70c 0.0406466 + 14028.70c 0.0049816 14029.70c 0.0002529 14030.70c 0.00016675 + 20040.70c 0.0062902 20042.70c 4.1981e-05 20043.70c 8.759e-06 + 20044.70c 0.0001353 20046.70c 2.5956e-07 20048.70c 1.21339e-05 + 26054.70c 1.1341e-05 26056.70c 0.00017803 26057.70c 4.1115e-06 + 26058.70c 5.47166e-07 pnlib=70u +mx122:h j j j j 6012.70h j j j j j j j j j + j j j j j +mt122 poly.01t +c +c Material : ParaOrtho%0.5 rho=0.041975 +m123 1001.70c 0.000209875 1005.60c 0.0417651 pnlib=70u +mx123:h j 1001.70h +mt123 hpara.10t hortho.10t +c +c Material : ESS Magnadense Heavy Concrete (density -3.9g/cm³) +m200 01001.70c 1.3431e-01 01002.70c 1.5448e-05 08016.70c 5.2962e-01 + 08017.70c 2.0174e-04 11023.70c 5.1116e-03 13027.70c 6.0596e-03 + 14028.70c 3.9343e-02 14029.70c 1.9987e-03 14030.70c 1.3190e-03 + 15031.70c 3.6290e-03 16032.70c 2.2703e-04 16033.70c 1.7929e-06 + 16034.70c 1.0158e-05 16036.70c 2.4026e-08 19039.70c 9.1400e-04 + 19040.70c 1.1512e-07 19041.70c 6.5961e-05 20040.70c 1.2235e-02 + 20042.70c 8.1657e-05 20043.70c 1.7038e-05 20044.70c 2.6327e-04 + 20046.70c 5.0454e-07 20048.70c 2.3602e-05 26054.70c 1.5462e-02 + 26056.70c 2.4273e-01 26057.70c 5.6057e-03 26058.70c 7.4601e-04 +mt200 lwtr.01t diff --git a/e2-source_2m.i b/sources/e2_2m.i similarity index 99% rename from e2-source_2m.i rename to sources/e2_2m.i index 3fafe20..e9cdc76 100644 --- a/e2-source_2m.i +++ b/sources/e2_2m.i @@ -1,7 +1,10 @@ c source definition based on tally_181.txt c but adapted to model geometry and actual size of Estia beam extraction +c make sure to supply the correct WWINP file +WWP:n,p j j 4 j -1 j 6.55305E14 $ WW was generated for other source scaling +c the source definition with neutrons/s sdef dir=d1 erg=fdir d2 x=199.98 y=d23 z=d24 vec=1 0 0 tr=101 - par=n ara=165.0 wgt=4.125e-3 + par=n ara=165.0 wgt=3.6041775E13 si1 -1 0 0.5 0.866 0.9659 0.9848 0.99619 0.9975640 0.998629 0.999390 0.999847 1 sp1 0 5.01272E+12 1.04029E+12 2.01432E+12 1.39792E+12 diff --git a/mcstas_source_vs.i b/sources/mcstas_source_vs.i similarity index 100% rename from mcstas_source_vs.i rename to sources/mcstas_source_vs.i diff --git a/source_proton.i b/sources/proton.i similarity index 59% rename from source_proton.i rename to sources/proton.i index e0aec2f..1137f5f 100644 --- a/source_proton.i +++ b/sources/proton.i @@ -1,5 +1,11 @@ c Original definition of source from proton beam -sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2 +c weight windows from external file generated by WWG +c make sure to supply the correct WWINP file +WWP:n,p j j 10 j -1 j 8.7374e15 $ WW was generated for source scaling 1 +c energy splitting for higher neutron energies and roulette below 1MeV +esplt:n 0.05 1e-6 0.2 1e-3 1 10 3 100 10 3000 +c the source definition with protons/s +sdef dir=1 erg=2000 par=9 vec=0.0 1 0.0 x=d1 y=-30 z=d2 wgt=8.7374e15 si1 A -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 sp1 D 0 0.12109 0.23438 0.33984 0.4375 0.52734 0.60938 0.68359 0.75 diff --git a/tallies/bunker_mesh.i b/tallies/bunker_mesh.i new file mode 100644 index 0000000..55f0707 --- /dev/null +++ b/tallies/bunker_mesh.i @@ -0,0 +1,8 @@ +c large mesh tally covering the beamline from monolith +tmesh + rmesh1:n dose 31 1 2 1 trans 1 + cora1 500 10i 1000 20i 1500 16i 2300 + corb1 -305 5i -55 4i 45 5i 295 + corc1 -200 1i -100 5i 100 2i 300 +endmd +fm1 8.7374e15 diff --git a/tallies/bunker_wall_f4.i b/tallies/bunker_wall_f4.i new file mode 100644 index 0000000..f331ee5 --- /dev/null +++ b/tallies/bunker_wall_f4.i @@ -0,0 +1,18 @@ +c f4 tally spheres behind collimators in the bunker wall +fc34 neutron flux behind CPC2 and BWC +f34:n 705 706 +e34 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq34 e f diff --git a/tallies/chopper_pit_doses.i b/tallies/chopper_pit_doses.i new file mode 100644 index 0000000..767f0d5 --- /dev/null +++ b/tallies/chopper_pit_doses.i @@ -0,0 +1,36 @@ +c energy deposition dose (MeV/g) for chopper pit components +fc106 deposited energy at chopper disk and spindle ++f106 4121 4131 +e106 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq106 e f +c +fc116 deposited energy at VS stages x/z/rot ++f116 (5711 5712) (5721 5722) (5731) +e116 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq116 e f diff --git a/tallies/e1_e2_f4.i b/tallies/e1_e2_f4.i new file mode 100644 index 0000000..fb829b8 --- /dev/null +++ b/tallies/e1_e2_f4.i @@ -0,0 +1,18 @@ +c f4 tally for large spheres just outside the bunker wall at E01 and E03 +fc44 neutron flux outside bunker wall E01 and E03 +f44:n 708 709 +e44 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq44 e f diff --git a/tallies/e1_e2_f5.i b/tallies/e1_e2_f5.i new file mode 100644 index 0000000..12dbf3e --- /dev/null +++ b/tallies/e1_e2_f5.i @@ -0,0 +1,20 @@ +c f5 detector tally just outside the bunker wall at E01 and E03 +fc45 neutron flux outside bunker wall E01 and E03 +f45:n 730 -1380 0 0.0 1020 -1170 0 0.0 +e45 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq45 e f +c 1-255 256 257 268 259-265 266-268 269-508 509-end +pd45 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 70r diff --git a/tallies/inbunker_f4.i b/tallies/inbunker_f4.i new file mode 100644 index 0000000..fa79786 --- /dev/null +++ b/tallies/inbunker_f4.i @@ -0,0 +1,35 @@ +c f4 tallies often used for reference at selected placed inside the bunker +fc14 neutron flux at BBG +f14:n 701 +e14 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq14 e f +fc24 neutron flux before CPC1 and behind CPC1 in beam and next to beam +f24:n 702 704 703 +e24 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq24 e f diff --git a/tallies/optical_cave_f5.i b/tallies/optical_cave_f5.i new file mode 100644 index 0000000..375e071 --- /dev/null +++ b/tallies/optical_cave_f5.i @@ -0,0 +1,20 @@ +c f5 detector tally besides beam just outside the selene 1/2 separator wall +fc65 neutron flux behind selene1/2 separator +f65:n 1328 -1756 -13 0.0 +e65 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq65 e f +c 1-255 256 257 268 259-265 266-268 269-508 509-end +pd65 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 70r diff --git a/tallies/shutter_f4.i b/tallies/shutter_f4.i new file mode 100644 index 0000000..586bab8 --- /dev/null +++ b/tallies/shutter_f4.i @@ -0,0 +1,18 @@ +c f4 tally sphere within the beam just outside the selene 1/2 separator wall +fc54 neutron flux behind instrument shutter +f54:n 707 +e54 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq54 e f diff --git a/tallies/shutter_f5.i b/tallies/shutter_f5.i new file mode 100644 index 0000000..0c854eb --- /dev/null +++ b/tallies/shutter_f5.i @@ -0,0 +1,20 @@ +c f5 detector tally within the beam just outside the selene 1/2 separator wall +fc55 neutron flux behind instrument shutter +f55:n 1282 -1794 -13 0.0 +e55 1e-10 1.58e-10 2.51e-10 3.98e-10 6.31e-10 + 1e-9 1.58e-9 2.51e-9 3.98e-9 6.31e-9 + 1e-8 1.58e-8 2.51e-8 3.98e-8 6.31e-8 + 1e-7 1.58e-7 2.51e-7 3.98e-7 6.31e-7 + 1e-6 1.58e-6 2.51e-6 3.98e-6 6.31e-6 + 1e-5 1.58e-5 2.51e-5 3.98e-5 6.31e-5 + 1e-4 1.58e-4 2.51e-4 3.98e-4 6.31e-4 + 1e-3 1.58e-3 2.51e-3 3.98e-3 6.31e-3 + 1e-2 1.58e-2 2.51e-2 3.98e-2 6.31e-2 + 1e-1 1.58e-1 2.51e-1 3.98e-1 6.31e-1 + 1e-0 1.58 2.51 3.98 6.31 + 1e+1 1.58e+1 2.51e+1 3.98e+1 6.31e+1 + 1e+2 1.58e+2 2.51e+2 3.98e+2 6.31e+2 + 1e+3 1.58e+3 2.51e+3 +fq55 e f +c 1-255 256 257 268 259-265 266-268 269-508 509-end +pd55 0.001 254r 0.1 0.001 1 0.001 6r 1 2r 0.001 239r 0.1 110r 1 70r