From 64ba1317fc1c498c734a0ac737bd65e0a00cd439 Mon Sep 17 00:00:00 2001 From: anu1217 Date: Fri, 5 Apr 2024 18:23:06 -0500 Subject: [PATCH 1/6] Adding updated versions of Spherical Shell script and associated files --- Neutron_flux_vs_energy.png | Bin 0 -> 28755 bytes SS_h5_Reader.py | 28 +++++++++++++++ SphericalShell.py | 70 +++++++++++++++++++++++++++++++++++++ statepoint.10.h5 | Bin 0 -> 36808 bytes 4 files changed, 98 insertions(+) create mode 100644 Neutron_flux_vs_energy.png create mode 100644 SS_h5_Reader.py create mode 100644 SphericalShell.py create mode 100644 statepoint.10.h5 diff --git a/Neutron_flux_vs_energy.png b/Neutron_flux_vs_energy.png new file mode 100644 index 0000000000000000000000000000000000000000..1e364f811fb77737dbc7d68aaa74fb4c3df574ad GIT binary patch literal 28755 zcmeFZbyQXR_BXs~0Rd?QK?NkFRYF<>1t}FpL^`BPx-pPVNGc*N1|Z!C(y1s&Zn}}~ z?q@DM=iGCi-+i8Y-+$gQ-Z5UsaOi-&_gee=o!>db#Sn?7v*Xq2qEAw= z=mqqhS({U`ar&}RA;XQ=duWYwsBj6W$o;GmE3;o(clWbyZeCY1jqd-&CZS;cNVlb- z#&4rP!rpr4*HoW-Bo}8uk`EgG9A`P(D=aMwf5Pf;L*WOFdQm6}a&q#M$0??=N;y{$q5PTJ!tB2vkBovkDKGfuI~!LzH#6-1T)r9-_s(W z??V%@`6!3R#vVtB*bbv^tEmZkUW+ug?9Rse{P{C~%#X7ewNfUf&t19l>4Y?mFzB?fwG;Q zKtN#Nh%H>95{p~s3!QqJU~iEea;qI(T?YNdcEz?M*tc)rPVyAHRRwo@o{OvX(_OBY zuV25j>MM+vWNEt4m45#zImT6f~A2S~?~R5FIXeY^4g!-tciSeJf%yN-l0&KPCo30tBlSP3UPq9{VW z2wNfzZWW7ZKHh7wU87XC?X&jjbcX3EuNq{uMDV_ASIR{F_?37lZLzPE=w^<2%8AjW z;Y`{gEsV2mc%Nh87NVGMK=Iu~wn5G!T|`|bt09NE5Lu22I|~OT6ke$D$y{!W%m;^y zZ@d2EYErmK1XJhC^YhMN4-|#j5sBbx;lQDzDlNtjMedX&Q9}#!zM`1UPCKGY!!F%6 z-E?1N#M;m9R&o~~2K9(Y`P~6cw*XUv`fzoSfcg@z6iOIe*`#MHi z1dsa}6YXsMjEr`er@elu$e(X~iy6}(?^6-4t4r5%8=u@!YHy?V9_yMaCVckqlF&wGRGy~_h8hZ_8m`teSo2_q1>!9k+aw`9W4_MrB^IA zjdV!gy5agUyfDvIQ)0Uycr~&G&vZfz z#BSuu&yfaBJ*+hsJ(F~%r1QBpOx9e%DV zdlrE>R98pG85S0ty1F{AEWZ1>xN6BNb>H7y%(1}m8W|a#k&-%{teOzK6o=HN>FY`r z#Knn;C_Fs8PngN(xIlkOc7wgGWt2~NI7PbFX9{9k-rBY4_NNP3nVF=Hj*jsOTwGjO z#X3SfO}QV1`q!~L$?mAmyYWwGHZjZH8LBFzT!^~lq`YY7udeiDJ7}ud^aXXpjj44J z_1Z#M)PcZC;AE-mIx|2!ug_47sFXc-i1h3cjp*lp@d>c)YGV&B6(I{ zVoB@kpFMy6e6~Ryk!nx<3x*uMDs1ith3T!54}V1a_T-uDzT`W~iFX?}VGS>oh!XH; zbgT!>OjwGAQ`8meK6nslH&#DV?Mp1Lpg`W|5`>;+_jz&I{@3wSqy)VM7L#`eGquI$ zhAM_uSC>X=a!gv7IyyRr1Gq9P+de1?=vAUBtR&Yv)S^TkKTS?THgs6%KMC30w68GZ z<8HX|NsF%ZDl|EZo>0uy%w#8ts7D#bVa%Tv_4){%mFP=9ZWhaOu~m6obI|(dp?Iws z!`OZUgM(UGyHe@;6f|B1R}b~~b%_xyR8Teryn zs0zpFFC9xwO^uecDC*16ufZSx^{d+G_Tf6gr~ha7j>9o<9M#y=^t6*g6F;SkFWJX% zk5JFiQIM2L{1|GgJ(=tm1XF*%*87hilN4X`YUCOb5)l(GcRjE~YT1f)~)B}+g(*HvQoVGy>KwYJVwb%o8R;nBoJX=LkR>9|f%$BMgBym|8` zO3GUj&d9U1MK5d(9i3+7GuXEJrAX@Wo$tqOTiZFsRZqudK9MB!*hqY`7dLlu`b)MY zwt}MKMF(uUcC||6kR|Rpw!K$Sj$b*496fe{eO!J`5j-U%Tf}0MNc5JUdJe-DJ|z z2z}*}A$B*^Hx3$%AX0pxu>G@u*pldN*%4`=jw;&v#eY}P;#>lVF{tOKkM6?)B{y91 zarhf^6?!%riPm##Zl|~mkXrQnqyDU>l3F-)UP`DK9}867XbIbXUELpV_*>4;#5xWk zOPttHK%Z0MbrjuABbTRKp+L%xL~dV}Nye9 zU^|f(>5rmbYp7LB$h!BpdfJamadaCo(zQ@}k=r2dEBI&WvBoC+eH(QpIVC=ta2wsc z!Q`}s7XQ2rw$P<9wi6Ht|LW3_f))9hSOa&#@^~gUMNo0xt$&snXDq+d!jGT%lp^P_ z+7H}^%I zW>ovIGi$ADqB~8)jLaQwm2-2--@sMhE0mg(Uy4?7Iv4-IlcHuv_1LTw<9NiCq$NFo zb5tf&a#u2byyfg!-IQ&K%NA5BH&VF5|2c|PA^YGZN0(H+VV%vc;Wm2$vy?Zh5?fx2 zXP39=l#0xcvrM>3vIhcy{Q~UP6V}XJzd4(g+}-!rNtw(vJF7quPND4?^@A-v*2A;R z%xp3F=>q%Hfj)La4_b^{9aV(CDxQCI;}YR8dxFQ4q5TDCThTa|m9Ol}l{jRwyOC%# z{ax?X=U+Jmqn89@>+yK4ibE-%xF^x_hO@d-sBE@Zg|S>W=-2%wJ~Z1GaW`tC-K~y^ z@%GYJ7R>e4NPmO)x0KSQ1Hq|k$p$sgNlZIGQSWTb2SOKffAxU`8qNGwq9Dz9^|{z` zo~D5XN!;$VULY$b;)=Ig&(OY15wx^2w)tHg<7p@@ix zvfb9r`&}pFZk1nTW`4a>>h#pLjg~BiomF9l^~tLODTB^Woi1a3RHq~#1YR{U50&yJ z_S(IOSwy3sM4XssFb((-BW$=f)jHLk!@{aiy)x0{IYUKD%OWAc_#;~Is)`Ef4Mjx> z44}3Vg#%O~^gkkAdoHZjap;s=0ZEA;n}Zo9y&p5Rq|3CtN&6CBp`wmgYBD|&W3syZ zBX3*d^0r&P{DmrWyM5}}uezQKoTKk$^*YBNPepcNqinY`&_D2JSnTr?LXf zACrBEiHL6AyqW0cxzPVPO2m%2KIAkRii4B0>f0->9NQ5Yr36VvMlq)~w^8T4EjuF< zleEHn_wMG5NWW1_ zVaO6pUGRwfkTe2m{w7X@U=cYCU1Q_4TluaP&g!w&dYw+pi5pe}f^aABa-aIF4Cyz-d`Mb>u*a4%SNaet4O=u7rk9<9sP%oy3yl%?{c z!ppv_eYn>`z31``L{g|iUaBOj!Xiz%mVJ)vg?eNVk5x{hL|;`*QRx>Q7!1;z%C8nl{h4+%vDMRmRBD{ zbV_2pd&SAV^|o2HI6dp_XF z(7{F#tSQ#4&|OVUdJ0NPyO|EotQ9z>Z#{Q6l+8<1QjW>V$)%*H%ciMkgvP|AZO=*H zkIZ@agFc8(0Q;f+7@wL=`hV!S(&~GywjPpt&+gJ^&JsV&pl`S_oRd4^uNCC>j>@m4 zYqUb@MwGg_%;v^djRRY6)@S-E-DpF(46WLV8gnbBmdmfVw^^=wU3{Tg;epLzO0*R$ zwop0oS~oLwfh9zoIjPBqn}!(4iH@^fuJJ^oaVqH#ibx~almi+iPK4QYNM)%uJrC9p zOT660Do>%V%!JOdvEHm5g=E{9P=VKUdZA@V>ZFd64R&f~W_{ID0xXp0&RSLfXZx_2 zn81#XI}X!rc#ul2KnpzRHp>Gi*>LcS>(){YSq{WtLc!yzVJd!dtX^SzRP3)8G~Q~} zC}H#zW3o>T-$%;I7YxG8!ory3FWK2A0%>@#QQR8YsH`t5L@M#({@=eVeNT`&efjcb znH^_UtkQ!4>g!+#Xlq>Uy#ZR*ucDVGTmBbdk6KOuIBn-o}D@Ng zmFl3ft+&tQ9F|a>^fe|D;o*T)pSU#tY_Z}{^2UuDwo^aPTJ{x^XqP(GeQ{m=nejk4 z=;Xy~rn=C4lenzxzE_A4nNXCK;UKzcm^U_fm^GIyTPJ(+A@Q?>(i}>nsCeO87;)i>y-KEZanO|0JaBJot zTP)w@g=2T-Q(DR>6@OR`uBe6e4L;0sGL^-8{ABFCWs60P=Z4W?G7X%=OO(;5jj2`~ zzAsDI^}&pRXN7EQE0`gRpwloqFREF`E*v=+BEeB zI?l?=hn{K59dbe#!sO%}gB;YXw=mT^okQ%Qm%46>waxFSQ}Q`Y`>NB;N`{UZo};i6 z;Ck(a$!xjAqrEI_{X%6+fZut*k#?bEmKS)_q>>WNRR$S@H{4C}ha#q_|A}5?=5-We zus%Qdt6(t8!rb0tqtVe;WX@K^qi+QGzWDftZ!3QGwzr&ylYWf!oCy~+v8ogyjQ8!tZ;6&Os{mUrjc;am&06bFuPSxUX%0*8SK=jC8m zz}*M}w|T1I$(*+C0^2%Y3sj9N+~Yl6Pv6Hl1UoM`7l))Whw9%R;kYiPBOIslCFAG{ z9#)j3ajDB%i>+tz62u=#oFWxtRd7EdacPMS>kP0CY`&er?EJ)3RKC3^DB@8_Q{m>B zX5wioveVh5(dwq`poM%f3RxjG1&U9fK6!lO(<(;!)tt@cVoplr40sYPD)#(G*~>j` zs3{94U6Ix-e*c=6LijTp1-A4UQCtWP6bkO4D*JQHO*Ip;tsA$QB`9`qZI_0tj-q+9 z;*inbKaD79nQU!3%F*HY!EJBxh9O3_!8b;_TD6<)L|#u&R>gtl z9eY{XD`^#9-+8TI0Fprm4ZwBs%w$x|dNsNY93teD=NwOI#H*HPw#jJm}pA47_Mj!K;%^;aHx9oNu|B{?O z^4UvvgiT@T6s?-KX7~Q3rHJd9WUf|c!c1i-E6dzBM|&-@@K8XwRnBUPA3lC<%!I>& z9a4A<0_RA^Q8URKd6h+J=dsJC+Oxs#cc#D>kJL}%FWxhVUp-PPcC0|*{Xv)jl>Mw7 zSUm)f zf+L^-$cPWC52#=<1RU-A{qk&BX8}p1+8#03arYqbD?fgKF!wD+f$ZWrO~Ax3Y-1c>pOZ+Z8?4wYDt*LWtp@qtbbm}({V`f zT>jEq!r_F2HTH^$yoWzz932Z9GtrIGu_|~qzpPi>W*W51sRH5c0FEV^yShUKO?Xge zESfcq9Ic+#QRXIj_!&{jv|!z_-yCo@1nM_zDgZsIR5PCSK$CbbCNAVv#8vMm>|6Pe zM+$kk8)tg!)zqEF#hW|3D_?C>5t1vQimdvc1Hmkp#Kjyy37D}NYrJJ>^wnNQjOa6) zwZ!5x3Ka1DQN|aayMNoM$R~vExAEya&benQSVNDRcyy$=0-mg_JbL-+)yQ;m0!31ax8ORZ!GCg z6=%7?MTW-QiZsa8kRdp>6kTIt5w8o^jw`&x_Z$|qrR?SV8xK&@E%6eq(7d#2PhtP1rG%mVd(|p* z3vK)f+POldF616Dp47qA8AF`iPU;ytPm4&4>JK7btF%D>Soi5G&b|s^!|&@8#$ii( z?p#_Z9wW8vFDB1;P=t@VbLY+(VPV?S9!prfrE?d;2Bx1gkjyXrfdqcd=%yLl5*EG^p*NUKxu&#d#=PAQR)CS=rQ=?wKwiFO0q zi3mfGA5OcDli#^}SAVG7Q^a}U6sWs?aM+U|OGmJ)9_^RN<}&YkLW*pmv_k+ifCiyS9K{7j3~tR@b$R>fODnH& zHz1`|@b=WRhRU#A!usCSRP&Nk&33ewP_<`ytYNv^FhsvQTOSTr65TtF%Fpx+P6rNM zCvDBXy;j{v(@&~kyE_!;>%X0AbcIv)!XrRF?kBEo)e}4fSGl=~dTrz-jF?F6o-A>k zeGan57w967K5W$CyvxL_?3%Y(L~^=iy{*bN1>V`SSE|2##kiMUex{1Wc}K<>%NbM< zz{Ax591%Q%lHObx5b)SqLM5;!xlA~m2ar=Y>rs)OY4x+v*H^7v_R6D*%kQuIQF#D| zXY-9~`#_J#|0nJV+x?MK0tUF(g6Oa+9ghda4|uHp=q|Qnc}~Lc+_))HZr2Gt8sMIs zwd+x_R;ZpHE6y<-LxPr%+hwelFv-LocJsC^8kg<2#mSRyQ-Qukk+1C8ImJj9uk!iA zzMq9w4;*GX2;Dd4$a`M#|4nTR5QDb0nIA{XW?e(j^^CCl3BsSY#~43_d%00GvWDi9 zN24NNzD%^AY)1JAT6F0FuY?5Jt(!M}*yCMfjP6w*XB4FZ-#}|QNmB8MTb>`Ippf5W z5OR}ysvn=u?xS(dtW{*yF+w`V`?|f)fN`p_$5x<(#P&6UUE()y*JoWy!g4*hHs7p0 z+pr-xv9IsHkoxjdb)TYT#vRTV+a1tOEq_11wqxc-go>2vAG*0HaE?fL^U_jz&vtR4 zqJ`?V@Wd-|?&wSKbqdYe@xi!&Gb~EgTeb~stjMVmpF09@ii3j#h4NVMytlR8;?3mX4vJ&E z(B2y{XI?gySTuw*j;->QtbZwvSg2CxMbd-7l$5lztJ>OeAimwecuyJKf6m<$-D`^> zz9V&0qF}b_r3FtQ)h$ngK9yp_wr;_)206uwOio<{I1c3X)MU8t8hNbUB%O0|uXK`` zexO+PSBktiqd|C2%&C-7!M~Ed;zxao@C~~^+PlFth8*)1>6W%YjetPzzA{NJ#kc43 z5ILV!ZaCl;v{BF&jDT8M3^)@(bErzBUa_*Xk3A+}I<|l99N}rtHF5xx#Hj2ZOQ5JU z3N7CxA4h=(b_@6rYA$tCYj+@=K-s!dsa*+Wt{T?J^BXOn&hfO)PoJfdT;6h2TI^ET zrOdSLEg1ixg0RBLrOamC;#Fr)TCzJ|NnA>~Jak$cr$vIv(2DR{98Q$_jN7&KM}=Hf z)T#f2&XSkM-uiyAMelo*z*JvR7TSF!O&{9@(g|`v{vF}d(XsF^J@kyAAT_9X>*iUK z-l`}u=Y=4sZuwSy>T2G(+=g&E=mZR&4PkznW&B$i@lsb-5{Jk=!&TA#&#)gDw+j5j z0>X)dchvuCUaALFYFdh01|6ir4vk^YAflqqv3F%^Gf>*47W9NI_q(Ag;(BggY^vM%!4Q`qsCJB1O`k*U1I!xR&f#_c~=JUM;(X(vl z&Yk7t$gIw%eblH{ckvCC@5Idqs4zl;4VlV^*DA`lj_CE zls7Hx&(zh9MFpvJQ)RbL*GBR@Sc0oja;c|FySiS6#PwLfuwJht<(9m%vQ_Vk=U_U~ z0Fjo1m-l(i2oPliASym@vkf3bm7L>>YHP2W)?sE$m3lGL_4U)RO9oY|KpKY+Pvg0> zi0b&(0yVn}2wdXpFqEXQxVRTF)+L0_gRbnNjx$sa?j!&Y8oz$M;^b5~DN&=6OsQSs zPz6kiIJDlWP`olU3#xMr8!oXILs&BiSxQ@5^D42e6J-A%s9us?$MZzi+y@lbo3(LK z4$4-+#uZNXUFXG_5PeI~jX z@k|8@Sq^?V`ILa{TGKC@{?%q)kM)Eb4vOkan*;>RpFb}uJ%$p;A%=ze@q-Zl4O{p7 zf%3_netK0>L3+L%Oa-lPuM)8Nj4X_U-C%1{>d-$#>VWPkK1uahh*ajZeP?HHSL%y(|i@_URh*=et>vm1*;Ax`tXzVo|F zglA6P`ohi2et}}My$G|x=2jNcml1Kd>ZML_sE}o+cP>Hk!S)XuPhsEA1fC^~NzD(V zfpBW^{(&lXAAOoIT>d-%r6kS#$a{g^2TVQqN2VS)8|yQIbRw5%D-KpwvYdCiq;zBy z#o+g(4FQr30XjmDx5gaK)2;pznv=Rezv+h1k)k<#U8$B7X5D_UU}N))-H!R#ap^>w zZ0uvDY6pGl?L~?fd;5>QQVK4IUN$Mi6wIrVE?e@Cfrfcv$78SIw0G@nRz)xggU~Rn zBO1ZT@wTN?qY%NS56btIAue)3m0MID-E%Hhy39-i2cMJ0yT?I$`x~DzdR-p>V zsQA#FAbxP#?V#gZ((se{w%qWYO@Vhoz;3Z#!mR|EkfDDI#l&bT6CP?65Pm`14 zT~BtTH~%KU2g+`s;kg~db5V;_3r^AWucJjvTc|AeFG^6$T`{&|?8?#&@x5M_FteYe zt*A7m`S{BhaaLi~>JTRar6mcQpAYknmPWu^4`C2#7us z1}>#zz8T3tH^CwHXi#lq{_FW5o^1-n)AT-lHwq8ef^8^Ki<|@JDD!g62F&x#TaQ1? zJE@@XS|_h7JK?=-r9n=mrI6vo_!YYK9ol$jKNA2-R2COiwIf>SEPh)NpgUaZyr)|i zp(^PQ0*J$W?_Hzz61FA}+U#S0GwTHeCUK;s=eWLF|ICc%>?H2TuzFEqag{`lWAk*z zpI#pW&-5)y#QgB4eugJLbi(!9JrDkqkYZ|!CSFYm+P^Y)_bknSe*e)oUS{VoY`J$`%*ufnG#d zSZ|wgiE8GtB#wZgg2&*3V3JApDY7020Hftu7M80>-;9y&w-^EG0GMt;7b zSG^xX%BivQ6Q#kEDePW+$j%G44Lo~^sG)+20$K3gbq2o5p$KxJrnfQ zC@%?yo-aO>I&!_Cm9EH*#ev9WeWo)?!d)od`@mJp{)eFCNvAL34=3oQ1MUJlE5t5t zV!8avTJ}9H3Fh>3!~N5hf<8XOxcq0U)@7u<`CKDe2}%8~-nmx-KR2+`R&yCka;hH! z8qYLr-21|r_%1;^62105fo~dwzT^!7nyzZ2;3a@mL*T@ev6L=n{cSBvvjY4rE<`RM zLNmt%o^;c~P^Hb9$Zq^_R30>+@o3pmKf4bS9wlk+ZV0>Sx3mV{3xsjZIrfaE}IkRL>GQD!akT0P+MX z^Pczz4c4q5Oy;9@`tFGqe|T5QUQm-C-HYK_-rp-rG7Kl}MbV%eF&N2oCSruFs%^Yh z89aclK$4<0_SrJ~^OdJ7Hazc(345F%x1&6EMg+^UEBz!=)7I7;vQ{82-vQvXbCx~+ zn-DY@1>lCDuGQL>jDx0uf(w7RQT_LGn8T0Om%A1iI$TvOU9atlt)1B{Pl`&Q7)@0| z4Z@BbrP$_So(y*SrMtE!&MO$63_MBKaW{F2=={N!Ka4eC2x=X)JVAQx6u7J0f)(@Y z)8mgRqhPgEl;Xut0v@6hci~52))o=0+Q8_Hh#T~_qSVJ@+^2sq{XyklkyvJ7W;TR# zDY@BiFXFMq4`>fR;%jF<^jBq@gMqX z#jhVRW6>6wR{*8RIlU9h1&=Yat%;$nvGsw@>_!g@rfMYnU8+eI!ygu$XF4tNrOGU= z;ArOk>B;csr}CL~tDaoE`}gn5-?%|^ef3tIcIA880PH(`+nU16NqV-lwse%0V7{9h zEMou}$XFyAgSP->LyX~RU|_|)Qo#+)qMLs@d}^CD{BaTSV?_f^?2JVmkDbY<k)Evd`g|==8ltUoXnq{M_~KD$&t0-GCQ? z)6qB0j2X}!9uS-&cP3(n*QY{P+OM=I9iY+jVs8ocl$A*imi(mjz#k|A#7nG7GG-kT z?GSS*OMrVo(r{O})o)@HpuYx>KPq~|i41&|yZ~Hr$i#gRH8|dV;ffzAK+CFq8c7>9(pk=Y`CO+hDs#>`^$=8MPZde zX_OPXPUU=LELw2%2u!M+j4SuCXvo?vrUipqSoNb5YXAY%FUx8_sboY@GwoAe^{(Nl$v`emepo>wa`s#seG_m}@~| z4<+`k_e#DvFoa;=k zxG=8JU1W0_AU`#c_b1mN&Y_{r{* z;JI-UU|k@XI^4gxtxY`!UqR5R%eUdg3)pn*wEzaUIm!CRf$x6ghM%*9!2W9H7 zyawVqI_C2~(yn#P_BiC(x z#ut&;m`f?DjV|be_}=*ado7SRtwlDv{kAnEAS%e^+q@V>x2NiTXXsm=R~Q=EXxjMd z!oM%B$@bZhh1kf zyS-B6y7~}7zhM7)#;TxE?(tbQ4wm+Pa`LXQ9nUd76W_3I~C6!99 zs)s%X6C*)CdKFS_0KUE|>c_*2B@WZrpu&Ap0C^g*l5Y-qOId%Oph@`C+;r>q?I*FZ zrx2vIy}gZSHn!UnuW?WpE?(3x>hqK=1&0?hL4!iRP(L`_w^kDlx#tGI7DH4k$1!g5 za@RGpOiD}ut_j9m7-L~hnu4x*Z+GaB>2Ou*@F`sJG}x4tt+W0MZ~75aWSt)ZCZ`}q zF`D*d)dm0@wAt2wvPKb%_ao0X7isVHMh&5 zF|Pv&Xb*`?gVhslUU%_VO}BrHA3bvOzPZ^D&?}f>h(^4zRrRtNAmtRd@bi{WfS)NK z32mMZ=oEpeM)$8lb@J?LgI(tuFFdDEef_lOGl)rSl%)_ux;hu|2){USd_c(sx?8U< zfEugs4iz3^LE`gTCjK5R$fdUr+or!U*CRO*NXz%U)4X%%E{6&y(M>BqOjYxG8PILF zlhWdU5Nc2$zGW(sf`S4%Ve3=8TE#vuU%oUQQn_=?6INcxQWk)cI4oIHtTl7wcn5&^Tn#24 z*i~YslU3rf9cR^#)A3iT+4Sa{b?2K=mhY`JWe+I+V-89Zi?!pjkJ*QV!w*77C_r+a zg}IGb5{qgNcj+5fRxLOV5cW4~B_M!E%YD9J3xr^v^emW^U}a+)1%IB-1YAo7PLc2W zm@4?=0|EkW-oNh;s>@k0&Ph5cpi-zyiCv-l^~sAibNk!gQ`8hWR*<7D=t&w7Rc z&_l%Up$DauBM2u@n?vzAD&>`ZI_~Dl`ia|!kw^qf++!DLlcF^t{N}<@FCl;mV}Vl z*w~-t@}d1|l@k*afXC0!d@NAU(s?@T_p8e2!dd%fw**l>lRkl&d&^Ru6#h;kc~Aeif#ir?`GAI6 zk9aBG^$EXo-yJs{B8C2=(cmF%cDyZ7cT|VdPRYi0i=KOZK|cp$=OWO(q)O6n$i0my{1R^lQIr3HFALO#ESQ_s10oTv;kSN4@*A&5}*`X zWI_l)fnRAFdXy?#!!4%k6GEECZmJAT?({CNsbHS2W<4B#>pRhQ+51|wq`$0KEVg=b zQ4WRWpCfa~A~W-im%=z9sGYU|G~zrXC6 z{`%!H58=LbzVK(bhavm)rH2|X6pIa+M=!pGmKZ<~b&KJL)OERrzmg@oH>wj6LiV;- z%e>U7@E`&;RGRJxH-~%P!R*@_1hVJ>Z>25YMg>VIehyA}=84#xIUy^PS>ilY*txw+GIa4c`U z$l6ryckL=>KCNc$0(Yk&^vz|WR1F64f0#M4oV&0E~`J;Jw7)2 zz}iTh02%9}co0jp)a`^hoc;3Bsv%;%dBwh>fiGz)B|x zFhovY6-4gdo!q`Zf7j+Qa)-#U3oI_sMRz={h4||=ZgGY97TZ#yhdGDyOjX*q zGk{WEet6)nZyS)Ah@1C}h8CQZ9!`F^5V5$K-(%Ft?mf$N zI$)P3G6~jcZN7~@^4Uc#X~q$!q=Ia1i}X>mbaZA*4n+?#;CN7<(RvPTIR|~$(;mzE z)_2#-KFZou>^pgYU=4C2MRh95xxB1qb+zK*i`oKhO&nBiN#PY@MP@`k9}c$6q>B`vdr=z8!eknS=YFF=ESB2%toV`r#kVPGR=JK+d& zSMv20$=-hQgp6rw@-kb!wX+qGfT9m~A48cG3r}J&H0_$mSmm;ax@_m5gfHN`2! zi22xoMDNj~3kW@oTm^t@cb*9u;vEFZ{w$2hnY6@HBZe$WdCS!B$0v<54>7JELl|4U zCJqnXi>scZ8|jJ>K;cqDBAf^g6PJb0lQ0O37+>fS$0MV2KUQ+PKS-GVGdF1ArRg6x z<|T)Fj`GE<*`$&qBYuh{c50U8jvnxJ+~d0!+puTk^C>fqcD2dyT!&L}k96EWfO`d0 zLOLN!x^LgU>E-qCTlHR=>nkFIsl6KLCS;vHaJgxvGo%qaUY+kd0$qGd%C6^-`cx{0 zoy+|8eiGr>!K3Mxc+Q-hp%4DVDR=1r?RRO^^=go}k_t1VTpWcU!q{si2|iqxHO1~u zyT=E{7ZfoKw=^u)CnzCJJVl9M&t(*KIMG8L=9i>(;{c^M*-8{E9)WHs%(@Qrl!(K} zdbQ22@Ql_HWq>OKCk}Xrp{kBR0WO<*=zmB1C1d<}%MZ_V#hR6Nskmt%8hl}7;OnQv z&u@s4NfEkd>^{`g)EqoKb>NyGaS$!b9Do)wsk9yM{x4oVa+&$QGO>GUh|;=v_43q3 zld@#X8+-&H@Vz65v?tIu>hteHW@-ZJa@t@K`|kX?i0A0X)NP(wV{aqkLg_*YF-Qt9 z6u7(-lFF(8E`ZpS%%p%)pY;#TSW3Io>(;LY)s}`%<5V=06tA6hKIx32BtMA*_+AAC zway%Y0Emg|c-RR$nE83Tmcc<(|0RbJjGEBnxj#aiG-O4CHx>(ZA-;4CH#_Iv7c+>; z`NTko%HEz7LG>IVb=mxyfsZqPvP1oN+%800Tgzin;;y-8!z*^yjwXf7@%o%Iek>>) zxH-sQ3AnpY_j+eu_=E03VxZ-3P+Ip*gS-b61<@l1rbNTJBW|#qh!+Cc*w~mkhPBNS zhHJJy3BiVtJd?!wvRu!wTuZ4s3=l_4c4goIR;ob(K+7{gdInI1g=&%9 zO>YrY2p_f;k#&OfD=T{zP^6IfhOa!A+^w4yc~|=*k;fO@X4yL>oumhLhd*-nRAelg z!CH?7nZ8z0RzXcd)gRjXC2ZQoy|kcsG_s)WJm3rCs~?>*7;?WqO^iwEde{J{01M^! ziSK1Q8i5#EPDr8{sCqwO)hv>ob3G zH{XCc7X~j!9@qS}qzK{GK5mbq1e^KFxPq2bk8$)&^85FC;8!SK`NkE~W_4%zjA_oy z>G9(!r@y8Dp}F3WDKf%Sa}(9G*V#$0z=_|tE6U>VL0awN$IsA>@-8RxeSqmai&f^f z{^aZXdHqp&XRbDh{KK|(dyuL6_(A#0{C>KKD zRNtl__pokvIi|;Q#e0!t=Ebh!Nyw;)gOm9KFtD+G&^mec(@c;{aYmPh92Xn&E@>fS zFty{~c{^EgELI|4i`Utk0x~ePGn#a2%k(WW1ElO~Tx|UUhL3;Q((+I1>}yw-W9Sbp zdxAzx!O;N5$7%;>q(iwrl{So)A?NTfmN!Gt;Yg5B=LWrvUWJ{HkMpD(vAzRU)=3e0%cm?6wO2L-93I6v`z z3r@sLu1^h@8d^Z;(ksugDm>{lFhE)W+$zFh&<0c7GQND#GP*~_{zT$g3eZ(qaTV6z zXW=9A!OyJ&O@DWpdy(q4)WI^@+QPsCqk9O{I>W_93?rmdooTcOdy`VvcO#FV`?;1? zv5E75$6sj!-W@GJMPCUf?IiERgCqLeV+8@iNZvzaMT89x(@7_Fs~aJwYxGjjZRP&)!{6gP9N!LH zxO4dCSUC8{E0J+C80650KB9Q~V+>M*8goxMxc^rZP!f*)FpqKI2(AN^GyooC*!q;H z!!dxj^)PBCv%B*CeoeA>CgYc0sYS^yU)T&8&>=)FUZg};4XTJQa8`!3X!7Eb=M1)s zWt#?wF=r+%>xub)MHEynex&jS%aUJE&~cB0?aAWxt_M#`REbpf-|6>2KHP?q`%&BQde08@g^ zLtaIt9%hKLfr}A`UhmWT?99x~d-o<`EU4{OZDVfwHevGLN3Y(Pf<1h1b(HXacHg{? z@%mxUX5IxOh2>PrlrPJ>gf_?sd2;a<%cr56z5ePjLw;UUfV>BuCGxj4A#2z_^WuGN zV92SZ~=Gxd*=3H!C%CvB%E0g{4m5j1}O z56L*;r0f-qK!W+SJ_KB`$Xy^?w?gN+p4z zFA%U18z5W6PS#&eAB*!)5H`h1mpZ*YoL7HMOjMaG$Y%V&dTEW1;udcsS{!obxw&&0 zrMEp^6?*3F=gU>AY#}@AoxQ2LgHcD18)*#+q)FN#`4D$xZumDzx7cJUNC}2Zgsgg9 z_PLubm|Ixr=V<5^XxhQZU1*>67A^?D|87u?{E_Zsn2&k~_CHcT(4J4BOJSNS*%76u zaO}%qGpvqvP4c7l&tS+UTKb0Z7txtDAI^P9bg=tw{BE*-!f(UBl`1}Di+*^t$P#7c z>SW9AOHaT?E?XmhNWn(W5rI4iV^Vcs+Y={D1i=LU!*ABc`qonrvfHzyB<4Hn>Pc`! z8-M)h^l-`NZqmqmNP_t5VDRG65sR_QmDgiQCuoC1diQg97I61ti)3qRo_?8s`3+6*t^`I3fCYpL${f7Jz>0-e~WilGCW;wuP4@#~jti=*7F zr7NzSPK>2KP?)Q`Z}nge9Wop}48lSyfAF2@|8>=u2or~KeHE~&@&>(04+Gx_JJ62ssZpbGB*C>c>-6Jm0y_L z;_-eXRSp03VCb>-b*@XDB(hIzlDafjI2FQGb>Vpv`@1`Oh{ZvmLbVcUdK%Hpxr9d0 z7XUu9{dsZ|o;CrqShH$68X9Q1KpHF*;>rIXJjN;VqKxzL$v*qa7$Kn-pNzlVZKAT} z@T3H?`-xiR6Us`cRNb>C&@1Wa1mh9YAf7_RN+{>}fp~#)cx7d!_ziaX%~h7%gY2HP zENB>0B|FoO{7yU9W`)w%cP|!XhIdPT}61qO#rpLFZ#3x{%w% zdlkG-I-cuy!+J^hIih&^2aFal)68yGUKre63M9J9R7%z)DZX+aS~}% z!quOK;eeOVTCMaQKM*dUVnxS19kOa&@;;5sUguF3=N#5;-EwlFkgjwLPJ|{Q?{Y4| zEIjq$q(VdjYNzkRY!Qn>_0IZ10(M=rOEk8%OX5n+03%ki&um54e8@VfO{VT?TgKxg zE~l5o*_%k^u;K^JR`y~~_4u*FL$z#Uv42t7M##DmY24yv^sp@I^Hsv~Ug@FFdbYSa zL{gi0ugGPf3G8D(5tl(O&JQL3;#I!PrwQIYk|%dn)v-`*eVT~-59RW07Z6iLzdO%` z$zlrc&9x=))*>EnRQCSYv!;C)I8%m{&%}07G{6j4EVi|qXE3#4&~NcC*#~bP`v8#6 z3+sWBlN&#^kzuFPrWpx#7&R)qqDadc=x-78JdJ9jd48N;lLTh}%ez*i3uDa{FAae4 zm>N8bHxyagZ12DhA-wi>_UGbufzE_jh5KakKaWjxdHwZ*(dDF)U{#l%4KM9kSWr5 zK_AfA^>uMHEbj>&NV9w74%srIe7t>V^s%7#G{-4J4fGpw=5X?#HyS< z*z*)RN_&*e>KNB=&&e3-1$_jHz1S$`Wx3#~UDWVRc2q_z{S4*)!j|3l+Ur@)5+<@L zNX{Mx+BILNlZiv>%h(F&@qQrB4uLyA{kO&z&{XJQ1<0w zn$h1W32|};s4HQ6CfgR)wPpd zR~bp5(RgnHW5`$o8(oELo`dEJYc~#VXV6p}!b4nS!@#Mm-IUUckM34;CETu3*Y>P8j-h;N^!;kydWF0lZlxKeTQ(&=`bJw^|{v>!WB2hN@@mH7T-KbYE`y=UkiIyhNg1$gM4yr(v zxj^d1tR((NuhEU|MSO#wX_;$l$XrT{;Pea}@;m+`ozSqo={#qse`q<0mvS+ID#%o{dS=rB24Rcj9p>F5f^ zbb(i>5}tV7e3NKRVsGI3@Q!$|Ov`l+*<3~@oOi*Xa3Rk&`%MkG90$M2O`QK#+L;GZ zy{&KjyG_{c4r{&ZUGMumpLI-X!iC~#r{aYt1r}ui-m=LP z8p}f&00o@1PP0gl#B!?erklP!G_ZVqF&;LJ0=r^FvZe-V1z}gb&|q!@%FZWFVbCpz zlR9{1$9h4|*h=l_nR%%9=03ZE6sVCB%bkYJzo)~an@J^C`knn{1^O`(*M)F%EF6zm0=crE3a7sVke{uA)(@R#xs=J@?$RIP`dT<{ut@dx~dV z!}-?j7p_xm7Kdq_T6z6|Pr`F2?Xm0k_wt4KtNIKQ%gyn51Q7S9hql;}D~Ohs)$Bdt z9zTa2?wX2ux8@m-aR7=m$>3TnvTw17I@D%hm|KFPzQtvob)2hr*PkgFX}Wj}A05ue zN9FWahscRVwIwc$b%hieZoj)~_I`cZ@fWT=Z)JD?M9`V+#u>TG6^CKxxt+48Dy_jU zLmKG|1LI3lk@5_#_&LwbZ51Q zExWwq+ED=#d0hRCs2mj5K@O9Km7y$qiu8G&7&Mq-wCtL_G;q}uSF-5!kSZy8_3b6i z@CdmbLp{w}ZZKPfM@5~2y-ZeKenNbF{Hp%$!KTwqaCD0wEVR{+us5_`P+tHVQQX_u z4ns)q30x$X_F$*yiXZ$tZymKhzpH3c?a+K=ANv37>)lOf zh&e92Kjx^LCU(SNlx6?Rqbg-@LYR}aP@=aijUy(vvf!YiY&Fw7B@&rXEc>zxR#@M} z(PPoS+Xz#ai}~(U1fvpv=mix69s+91=E=5t#+{3S34j(OUz#(xwqnH}EEqF3ds|eH z7W0Lkz6+^V&$qBshFf-4DqZ{O{VjXfr@g(1{tM{s(QYh#l{*7h^DbIpJh^w|jL(jT;7Jv>OhcOS5MQjSO@NZ;1Zm?spBkvQ-Zq?%bO69ql4jZF?5<$1w*} z6Eu*MP}?m&c;!S(EZNC#{4n2SAg&nAi3?FbsO&GPYxQt=cLJEtm8Nn@j4xQZI_c?& zITB@^JQ89LyiQt~Wws;%{VDewk`Cnrdf(NjzgCa2_=g6y+g-Vm zIt}x>t2Ws8w;Rz)o~P(+@=5UUSK$@3Y%h$+9T*Xv4BxxBv!FKZzWD<9Qx>Dn)2dHh z-6)Sg0-#v;d(ibIqQ-u8Ief#br2fw6V3U(sTmR(vvDX{a&yraAK+cQFX-eZR7C`3p_QMav!qqJ-zYIYVq4+(>HQV#vPtSc83ECOd4%ql0g zp^W8QDND^S&9EAtxKrSA&N&OkN)6AC(U|q9bIvN|enFn_(HE-cmo1f#*_Ps~vn5;L zb4!E>HF@uvc|~zyIaDs*<`B@nS&_1YFW=j{&rnp2fb+UfndFkx zMw8bPJJ5&o>W}I4g!Kn(FnN49>Qyh<(q*CBwZHOyZNlSN2P2^;l(v)D5NsrByppz> zvkjwl<7RH8eYbdtJ6l*>0*T7Xh*mi9+g3{Hvf(6#%HeZE)~Z7b2Ut0E2hKcq)qKs_ z)LU~r2K2qPhYe;VvEI5pupc8tzSaL3F56s`C18Hus;g;Bj^n_#+Aynl&Ttpksbybg zjkH+y+NJ*~MTD5Rv2}`MQr5v18_<&q534J35<2=z zc-|A_V6X?iz;jQ(?;hX1@oMNya5>UfesFaRks)zF{-v4Xp z<7!Gj$S0n07*FTD*BWzevFGmTe8w{i*N>9{0XJf{0!t_p1|nEYQJ5TEKZu!ws!Mn> zp$K}tmy@k<2os!$8v zwOH=lvXQwWi2`m(S&saCSuZf5R(8G|*F*mwu@jsHKvM#p7rik|o-2zD%-NZ0r(uE( zvvp{BaB6WRZ^pjnk0+)3ITr=0b&~d+;4Kf|FK_WZArdI!)-Pnpt7pa&mij|KPf&Fs zIkt%sEtD1^_hjFWn`A}`?w-O`bYQs$Z9>|%`PZ*n)bQTDrR7no_wZY3Qw)IMyMcBD z1SKWl5>$|0Apsr2Dmf3@S}qTo*yQkVGL3xYnw|m(8rI>_d`eYp6cye5#jK?I} zrtX*eYu^Y%^COBZ_e8vE5~q}JT=g;E{5a zu-qy)t=oCqzuR^RZRhtCS%wGkOXI3i1ey~QR)yq0bFga_sUVH>FR%7b*DS?TJa4-! zq52vJxb?`i#kBW|R1`fN|NO6HTF3%=0f{;Fx4^#NFy%_4YCcdCBMn-G0tQi7pmaSB z;e$%`1vW<1{&$b2EZp`m^s@s2UFz|Sqgw+7*Zk!K5kNqH!RfvX_Hcuf_sHcK_0tO- z?tb0rcI&fDa{rek5rWY?uu+bI&lmgnu_Up)peKC1HtTQzJM)A7YTpyyInLZvZo?-w zsVHtbQs%~G4ktAFPZb>~XXc|*ti%}gq*m(d(@mL{L1=*34)=YaeBNbxs8+3o(2k-Y z8OWWaTH}LU`_^0^$lqY=?H#q;x|L3SA$$*S{OEP8{V#B6zfwt|zvcy)8UW!NLsr?6x&ed=r*EmAola2Zx2I0|$VDECXmtRhrn8M+~mFFe1Fl zZE%jk%~&@>;rR@5Ipqh0V^#zIIt7HyF>sV`qbd5S+{9={d-<(#Y;swK>z;1_-r@Vz zo$&8-!^XH{$t}(>_jd}bzk-VYI9+3QNaf#FO(KrVVte|`i2C{U>#)ekK(vSsor_-K zHn9^@68k!o%oG3&hVawk!9G7Q>xjws=RLc&7F<%cSG3LJhn+2~<8LA4ZAea|q$uMGDobn?HcGR4ZgkA6|NkaYqT`dvc4~r;TNLF8h6J)Nt!BJ+^ZR=zQ`O*eIe)fLyuK&L?sTb3&Z+@D_mgLx5Jh&OisVI zb-2GH0HI=8#yEOU;`|o$nK1j8yZ>UoSt`Vy{Vve)Pvm1!pwn+u5HS6?!+{Nb$Zv-$ z)Cj+n$MFeaVN^)IXdg!t%G`Z3sIN}tL`N{>HjiDq(eE+)(H~C=m=j5CUZ&&Vcyd#R zmDYj;=<<6uXE{_-_STaA}qgYKp!ybDbLz5fOKE{$}B(IPK6NEu}@ndFgF{ zRuhu?bir@`{@2q!iizjq+&MpN)ce{i1iZ(3qhU?xJ5yC<0o&yiyk^h)PJWCC6bTT6 zXQTU#AzhX|lZznrurG*%l`*~LlD)BsiCE{hC66Hy6@sK8Du)2qnqOX6JMvKQSYdA1 zZnPm0Xwj9v3Csh?Q^uyI{tq3h*NnUo>_Yii6Rhq)m^YeEdym@PcM7OqPoeu20xj8^ z5z{t#0vo%c=Bl{2ROX7hhY;NJiWdjhY|g@_e+>X4rzF*U0q52oac`80rXqF+1qI=M zqI{sjR|w`B;2;e$TA8k2A?-b? zWbN6&dxYK3`ub^8!2FSAdWs4h?&>NY;4 zQ{qeHs26Y)yx(&E>eWP2Xk@rAms}6P$2qLzWlj>8yu6dy7naMexVL;HsL*N!~0SEFImi{gf`XymR!9Pqp|3dDRjX-{wx8dkv9aJoT~X6Sb1n+wl@FsHOLFf& zoLJ16@Jk9)%MZNgXY^jK^HlNH=qZ{_rcRtorhm+wx#EW`8V`J1wJkGSVL=E)j-@wh zeRA~%1B1C@5rHCM#4{R=&V>Uh{aymmAq8_AL-5hjekFs*h8QnE&f(MxuM$ zN@KB~L%9^uXmzy~3=;{Rh?TJpV&Ta+EGh102A-}4-H@Hb=)Ys5H}dAv&c?hTqNPHN zTUH^h>2<%|V{Ejy1X(lSnD*qv258I&-*qwMd)$MHd<@F6X|t9*v|+cK5pd??7l`Xm zs}hQLAJ4a?ZyUkk6G(@e zJFr2w0JUaa@hw;rpPx4Dl(RMxQdJOGrO_{^dksU_tK>q*CAvC zFWLI&LYr(!bOl#p2{*38D3ny*Cr;Q|AfXm@EQQY1{L#UAa}6RzHrwdx9wR4Cy}Jx& zw_|3f3W-7Eaw?XUM*~1#ln*R3@0*)r;bc;nKmXaD=%fA^Jq^&4_wH+g#nLa|_6gPq zV|c(olxV<-4*3A#kan%cVEyqyc9(b8t-3fpff|>>nOa8jb%$DK)_8RgWz?G^oMRu+ z?QTFA$wyHuMj^P>73OL`fKm{BruTJ)0?|Ej@}{$>nf(I`%6hPCACGwv9(-N_Ta+IS=AKCbO0%MC+0G)7T`$hgv)*O# zvAfzSXA7kOl**^_>C5(p9YBM1q>sS8=fLmZFBKTbg_Cq<1>Z)h{CKG4vh^;+FKg{g zRt}9R?}{~(m~3%`Pn0H-{o{3RF507?xUnMaqgDbvmc=}Z-)mzfMVoB}(xs#X@h%T< zI0P5pl;q8x|H!>F7;xYT{WBlbG%fZX+F*(3Yo6{b>n}7j81uZbw?w-zwM0vTyjP3y z$+532j+cRyLe-*hzsh0u0)0`042rcP1LFLq_zTuOLzXb_J|3p$6F)Buf&>{(hIIi@ z;`$N3`*b=qH%2PbPHXh=qgB!UO1-O|aKfSVhBqb%Hs1O)s~VD#uNR;@$vU2*?_XSG9`ER?fE^aojYz(`k) zcH7P~f$Ict0Gjv~BrMAh^;~D!^m}~b`w@yJq?ir32IA;~cS3u6#qBNAv>6(bLILfq zQb#_mA`KWLtsyxf$>ly+^-5C((ipERZWu}sZKP*okh~&RvS>RC-C8?HEnP?*^;SaE zYOs4?Tg;9bnKkRRFAf{fs~`g%tGFlI4Ku=i>aZXdB7RSwzS3BoaGf+JggcklhJ*S4d8tJ^|Tg_>L^-xH0!8jA7C>p+QuL<)-mLK@#XsWSkMBcRgA) z<)fd&i1d_XK18P+J8-zKZPWezm7*XX`<5akdJf{$Il&HMpG9|xxHbuYngpgOm;oaV znB1b|yd83_7PL=< kKoMyCys&?NA|6Nkc7Si5+L3@b9Pmf!>uuC!uHNqZFG9&yw*UYD literal 0 HcmV?d00001 diff --git a/SS_h5_Reader.py b/SS_h5_Reader.py new file mode 100644 index 0000000..4c36ab3 --- /dev/null +++ b/SS_h5_Reader.py @@ -0,0 +1,28 @@ +import matplotlib.pyplot as plt +import openmc + + +# Get results from statepoint +with openmc.StatePoint('statepoint.10.h5') as sp: + t = sp.get_tally(name="Flux spectrum") + k = sp.get_tally(name="Neutron tally") + + # Get the energies from the energy filter + energy_filter = t.filters[0] + energies = energy_filter.bins[:, 0] + + # Get the flux values + mean = t.get_values(value='mean').ravel() + + #Flux/elastic/absorption tallies: + tal = k.get_values(value='mean').ravel() + print(tal) + +# Plot flux spectrum +fix, ax = plt.subplots() +ax.loglog(energies, mean, drawstyle='steps-post') +ax.set_xlabel('Energy [eV]') +ax.set_ylabel('Flux [neutron-cm/source]') +ax.grid(True, which='both') +plt.savefig('Neutron_flux_vs_energy.png') +plt.show() diff --git a/SphericalShell.py b/SphericalShell.py new file mode 100644 index 0000000..21dcb00 --- /dev/null +++ b/SphericalShell.py @@ -0,0 +1,70 @@ +# -*- coding: utf-8 -*- +""" +Created on Thu Feb 8 08:16:12 2024 + +@author: Anupama Rajendra +""" +import openmc +import numpy as np +#import tkinter as tk + +# Create materials & export to XML: +#Simulating tungsten shell: +W = openmc.Material(name='W_Shell') +W.set_density('g/cm3', 19.28) +W.add_element('W', 1.0) +materials = openmc.Materials([W]) +materials.export_to_xml() + +# Create geometry +#Spherical shell: +R_1= openmc.Sphere(r=1000) #sphere of radius 1000cm +inside_sphere_1 = -R_1 +outside_sphere_1 = +R_1 +R_2 = openmc.Sphere(r=1005, boundary_type='vacuum') +inside_sphere_2 = -R_2 +outside_sphere_2 = +R_2 +R_3 = outside_sphere_1 & inside_sphere_2 #filled with tungsten + +# Mapping materials to geometry: +Void = openmc.Cell(fill=None, region = inside_sphere_1) +Shell = openmc.Cell(fill=W, region=R_3) +geometry = openmc.Geometry([Void, Shell]) +geometry.export_to_xml() + + +# # Source distribution: +PointSource = openmc.stats.Point(xyz=(0.0, 0.0, 0.0)) +Prob = openmc.stats.Discrete(14E+06, 1.0) + +# Assign simulation settings +settings = openmc.Settings() +settings.batches = 10 +settings.inactive = 1 +settings.particles = 100000 +settings.source = openmc.Source(space=PointSource, energy=Prob, strength = 10.0, particle = 'neutron') +settings.run_mode = 'fixed source' +settings.export_to_xml() + +# Define tallies +neutron_tally = openmc.Tally(name="Neutron tally") +neutron_tally.scores = ['flux', 'elastic', 'absorption'] +# Implementing filter for neutron tally through W shell +cell_filter = openmc.CellFilter([Shell]) +neutron_tally.filters = [cell_filter] + +# Creating a tally to get the flux energy spectrum. +# An energy filter is created to assign to the flux tally. +e_min, e_max = 5e2, 14.001e6 +groups = 500 +energies = np.logspace(np.log10(e_min), np.log10(e_max), groups + 1) +energy_filter = openmc.EnergyFilter(energies) + +spectrum_tally = openmc.Tally(name="Flux spectrum") +# Implementing energy and cell filters for flux spectrum tally +spectrum_tally.filters = [energy_filter, cell_filter] +spectrum_tally.scores = ['flux'] + +# Collecting and exporting tallies to .xml +tallies = openmc.Tallies([neutron_tally, spectrum_tally]) +tallies.export_to_xml() \ No newline at end of file diff --git a/statepoint.10.h5 b/statepoint.10.h5 new file mode 100644 index 0000000000000000000000000000000000000000..ba20132544ee0085ea0f5dcc7298d7e0fc5b92bc GIT binary patch literal 36808 zcmeHw2|$fa+y7}_q|%;aG69G(K`VJ2`(m4Jjd(eB;LT zdU`kgML`h`Cxm&wAdOJIMlW)hdGvbIcKkM*cuViuj8lJsSFNdpp}IQ|PnxlOq06>gQyuqN=K-tfcDTsG)4DrJ-i8t*+&uY^$lF zuH>Yw?5OVGpslQ?p*2@;*01?G5Mi;lb$77#adDdkZ#LiAKT5tTN-ApdO2q%_coo$p zDoRU~)rr(uIz5Sp%s7i9DgCQQtQY%N`fU4>DDag)M~uzhLQPHm^-e}|&qPDXfXq$G z1&8DKhkO;o*qv}VisVr^{uurSBjXiRL`fz|0}eMr){7sNj-5jI7N8{j)^{dHZMOV2 z{r52-Cn(J}-k;b}^nJvbR0%V3UFPI=0a4Ewah}c|J|6DYK3=x&-ku&_WUVIl7erwr zi6lwt$m018k-2_jR$m!$-X6YQ_Kwz`UXI?5vwg#_?w@Vr**wYJ+_(No3y~g?Kz&*eTFnsxt}J1OPg za1GemJhLNR7EZBN3YUne0z~B58Tcf@p^e6aS`3S@wHlZ5OP(GDzoRnJql}{vi~B56 zL?;Cg=~QK2Kaqo3Kva2dtUlkrnG?- zVpVjEQ)_OZggw}-w6DT!2d}Eem#yRjhAPK zbqk1pKOR35jR#-dEZ3a-e*7FeNBjeZ{uMk}h~@t=9{glHX1V6vUk{Ieb)F?o#)a&W zq4^${Kju8^*YYuk$3NN?ek>>QWFSb~6iNBB;6W*`x1*!OY?w41uy*rsAelnEe2Eqp zN`y#vxAwI4A`b9f$tDQDQz` zj<&8Y{g^u9GTB!l+eC0UrQhE6!x%{}aX*yHAKv&k z;r|2zzn4qOXrhg5_oskBc^{Zi68?ui#m-7GoB9coGNe2wfgCRpCZ&HX7kT}F!uR`i z3W_?lku5tEJ{*wp!7nHKBo7DVIE63AkVL2_+bzj#e^~E_c1#gckrT)z|JC~`R+Pj@ z!H;~~-`@Wb1b*KSP~J~nWc{PSKzc6`Ur1y}j;m1O6g~Q{zMm-Tjlb9TIr~2ISNp!Z zHRX*+deUsaLFyCWYXarD)Fd+T9!febG0cZ^aUc$EiSD~~${(x?&gM6p7uDB)bL#|@ z^5XDBBGKW=dh;v&8WAJ$Q~EJd?CiH*10T zv7JJl+y(r~P7(WuhLI#!G~d=w4J3bFT*^34OXBBzl>V`6@z44GD%4-g)xT;7)BZ2Y zSL7cO!te1Q^+;6Y_cWzHqWJ%~KZT@G_^D&{MESZ&-vb*AS)cjpB=f) zCN03DZ}Vuic{sHh@tUoKNv_h_{!zZZ*F(R5B>#U8f8FJclU)r*pMhB=syt%A98f_jEsYTu0$6p(MiSf~WZj*g|)D`oPQy zpj~;v@4&`LU=X=(_9ZhC1az9V8PG&QWbJ#r)7B`6idwI1TMz}ipKaQs%Nq^)!qxM6 zcSpmb;U=xVwrDv2L3Ho!WhWsxlO`fD=_I6eg{Lv(L_JP*&O793MKc>#_sOsjvQ7!R#|=U$f7#Y0}B?Z|}JMUb+) za<_sd0WP_8->kMx0NI%;s>u8VcuHR|vS>jfyxQSFe`H4@^qlB?|EMw%0!5RuKZ_^9 z!hoz5BB4o8a`_M&VIT=!3Ok)~Sep!k8#iCLkedurewzfw7pA}t-cUX7kQ8`TTi8JcIVJTY954b+Vm-1=GB;IX>rg`2?@h~(UeZymS-51w!{M;^Nh z&^vIpl|KjM2hu7uoT% z3;tyD`8LGd;Y&-7DuVkB`Ip{H-GKs?rCM&~cR*p{=s1_{T{sr!Ry9La3<*0nP7_j# zfxA$EZ&acLbfxcH58PJ*d0}pW%}phsfA+?+Nx4$+zEo}VF02$9A|q>RdrE=RIfpk! zw+xE*BzSF1DFY{gQpCwn4tm;q_&ax$1Lq_2#-he@_++83b4h|iuhyIkAZDB9!VK1h*X>>^g=pG6k(1>Y(kA$v4tFZ>i4#(O0 ztZ#&a6>0l}t~Wxra#nUEPZJpXah|%ow+RXvFqW?FCa^oqBdf8t89w<-L*cDvU=;t@ zx?8*jbcbZ~{ZF@mxNg(yM+~hnX>jvK;J#K!!IBfskF`P$W0kL>+hcf^vOrh=)nkwr zJ-KCgTO0VdjhU&xYJ-#dy@Jwi?O>FV(D-<)9U6DLju`KM0tenoYbP*0g{ijJcOJ(+ zg~{9#&*@}4;D)=q$&!)|I5En3~CjNVVs16~8t8 zY5Q|v*5<;m!n)vrhZbgkrV9!W`R>wt+XYqK@roY}UchF*L#&5WUcklu)BV#_-B5gE z=$@5rH;Aotepy)74eqOoZ!MH~36>JcZw`mQ1fe&pUD(H7g2}CTgCffw7)LHWXR3N& z^5KbwD;mAP5Mk2IlGY1R;bD?@`TM}dYU0x4ao~FCXz%3VKA0{ zf$`NrSW_saI*tv&>#WX&NsoshOe^b&v*#N~jMaW(#xe|>oOs%TvxebWZ{Xm2!?#e~ zj6$o&-ojx&xP8cI1lDTUj`yUDKs3C}WBoh=>$44f=52lls+C$DJ-P27KI=$z3i~KL z*k){K>NyHI*G95&jiVs91#`ZpJ_dE?la##9jRBjP^td6-I6U@v&~EHD4&G)Fx(D0G zVff3vy0MiLz#wC6k)1yQeY>7~zAZZmpZ8YD_mPTV%@!QSXK0l+<&FQej2M=W}AJAdVg;(}e>Cj`GJi)TD*Xgl&E2{b*%Q9f0FKA;= zr7~bH*YMk)E@s49dz%^xk{L0h0p6oJ5=_{^s)QJyETZpkrRS!~jLFTXSJJr4jIq%s z>pWXC4=enll;YDd4|~^TSHIPX^9}~73so9<|AKSg5kjGqy6)RAGloZ;@ zighclPNofF!!DJpZW!fc#{y!ntPyl!$Fh_qY?w;eF|*48ePw(LuszFy9M*d;z&^2; zoGPtbfb}Mo+b>k$z>JSEk2xLZz~cwLMm=T|>t74cSrsy`zyjEBUd#<9(cO+j36RxKg zdv7X?S&Wu&r%eiDp@j1anQ zHbGC{M*|am)*h!ZNc0Qc#?CuvVwJap%tE^~F|NA?d^l?@%z?|ZlDA6>3+_0z{HeV* zrpsJ<)$Fx4w#0NKAl`ckCaHSC^Y+XVEP-)T#NFddF{&a9nUNNaaU7h8Q$C%$}87u!Xb(P&`00u#_FnEX_-0*gJ!b8b*U57ST2iojpg!w$N= zev-C8A7i}3Qx_4ckBQujegEvUK2~37Dwe)`CAKGl-DPrYCFU#BApF?N05cz7-@AF- z0An8Cn-c1`3fpu0NyN^XRhW9RoUr7P)mSpaknm0JHCTRmUz+@-HCO{zC;o!wTC9&# zT9>PKEfz59E|BbKh$TvGnrA^}gk30KQ?n|;2rHvEay_`*7@LkLHRl^J#`-wYhQi*ydfW3$L53!+hH=_iewt4!d)(BZQH8Jti6afG);mJ+_`k zAVa8TJ!a%wGmxRM0UHrp@2hue1GX{dp62xQ1}yXq^@di*jTjie7VB=1T&MP1Ikx{fQmc+sq%el3nNjp5E+4%pi!QLdJ5 z2TWeO@AIy8j#$w-zd^e)M@%^T-KP|FCoI@JIeJx=6NY_uGupa%8+Ii{fcwpbZ5Uy1 z85_NTGsY6BkbNW087t!1kjy0Hf_+|7KEjsdf_;%teby?y9h+TScgx$3eOpiGA(;s9 zBy&G;&0gzk0u}K?InwvW-?prNO`o-kkqkCK=KJb>`>yxjxBmanzY`o;L3&;S$-|NC z??3+TK3$i@h$r(0ZY2er$??CvzcvCC`+E+#nt^10WDP0sqkUWn>q+zkGXIa(gYtt( z3~h(VjNkZw^?mo8{e3?1o*&!&){yO^vs`oTd%Iuk9Ptkr`u)DcNWGsI63FBdgP-#$ zk;yqLf6j*|lN0oR&i6;Be(Se^R8q>m%DH}A&%e@&{?q|8?+@DSEe~rh# zx?W;W#^tZFUP7u4MSi~RDv20cN8)K4Zni!iUSz`WzpMNG z{$-e1>lXCFlwcTmdnBmevM z`>i~amwmq-CqhGP#o4yyB>u4w|9;$#LF5J$`^}@jZ^y&Gq5p43el=e+TSqD(o?Os2 zYf`|O9RJ(;=k(2uf7!p*Zf- zz6Kk9E_XmCdv5wUAEn>4G5t9oo=i5`_;bELI+Y@?zyH2%PG0{i^R-;E-2Xjao0ZAA zbNTL&arxEw`;R%^`hNUEs?L9Z{8>uIgY1zN`1^6op9PQKj^}3W>$8sH|LgO$>i;4> zfnT5VzkT^$K7NY_MNTOGKaOvIK3{7fT9W}NXO5TZBV#x^H<7SZf;)t-`(%rL&nSI zHz<6Mm-BDv|LMs0cxZinBo*-h-dQ&%7hj{_-UB~M|L^YyU%^Cr0%b%^ihYld8Hw;8 zx{}s)=x~&$bglmRgUF4Fi(US6J_9nje$~(U{K#Za%b)X6q$FYO&-wnhsK1h{*#l~4 z9Yu~PTE8oX)mHil0Rn$74l=y5U8#H$PjU7i>bb$M@VL9H{Tm z4`{PMBVM8ByHQ61#y!i1cLgp$P$HO>VB<-{MhFf2p8kc+BSG@L$SdL&z|3_Bjir3?fpPyfSyK4Flf|<6ZqcOFH5<% z8r@CzroVRdFu}x#>W)4(XQeXLYf`y z+8?OY2x`q*_GyNj5kozF-ZK9r!t)1$m4>ui(EgQL3TG`62xmVuPYw*)A+LVt<=)=K z1RGn)&Y`+MBzR2o-n88af#dikiR1@?sIIWti(t@CnAhZ0f2z?L$!H%ENIqLa*rc9O z&?s(%x>}}t^?1Sw$7?R}Z@z4g8idqajR$WLBux|&g#COGmA2Tvea%gTE1$MMVwZ43 z{?(2JE&RoVl9to)w__bqCC@p*Q2D!r#zLC)F%mxLfb31Ftb=U?M!Qp;-1uF{-znH> z$=hcH%j8U^IHk>K^D5c~slqb^tK#PiEVX@+WvKkMLfTG31E0$1YFa;ZdbtA)_rn(i zexvR^_nrizm}E}8WaL}Ild=Q-de`@$9)t3ey4%_ai_$y;SFG58Qk!3Gzo(E$2zkhx z*q&m6;(8pXn=MiZE2wHu7U{VlR+gNQ7 zb6#O_hd~x0F=c^d>7XleTvehIwzY(?<=jhMtB5{VfC3{VBmnPe}XYcy}*B zd}uHI!Kgim!8UMj%UCBt$mW1;tn+r1o6Y<4s9FL+ZRPpy%mH>E+1t<+k#g-)Y= zgDr2z*_cq@S%KYmID*lArl90Si(U~f;spkt*X~7@jKi+ByzK0KPtVHNtqDQeL`iIf98Z!0!2BCpzc+#KT3$SSBw)Xa(om36+_-$Qd4JvHp- zwACCT7^}T(F9`}pYouf^C>-h~3^ed;FHAatrYcnytlm6DIA9;YJ7YW+v1wTJmx*(s zJT1GryrOexb<(=SBH28MJH+~guXh4EP_!zdHC7nWX5w3=0}|1Tz^jKHK8d3Ftf9iI z*y7QYVSa|+K0y@2OC_N5Gzf_u!_4OM_7cjiKv;Mp0cr1Wyad5`bWdl<(M&%Y`Qe=- zH8#>9o^v}I?(I8=*tRxKa;xzm>5|$_d?|^j(9tYP|MVgxd3G~HRZJSvk8CKEmQ_GI zWRARTwuwaH6|3=5d0z-Ah*c;4!9^t9e?wX9l?V#p*`705kbv$CuX|~I3XdKLN(DE# zW+KbmEGw@ItD%>7o9OC!Q;^8ccLr@~63B2uL5th(G~zp{B;hB?jHm-Xt2J{aqMBu? zSEp$fA*b$BqqsLIh;yk9(w~w&w1`et@OScfWl zdfL4Cib)!hZQJzp-g0@AxpUt8wcJWb!rq5B$~6@Q;2Go!(`As%N_s!T zf;4m#F@!$-pn$gdOtIvy%|Ovnql~-#lu%9F^m$MJD=0C2kEdU`HoDawem!?v8frWP z)rJ-d2$yagZOw5N<%yiTCZnZ|;t2c>0_|7OiEPsZCg&w6wpQs#wo)!SpQL=aTopr^ z$(0#1HQ7k(in`;fE1HPit?_k7WDc_Lus$C9VHtWT7JXS>=PK$xTfaqYjW){JVE$HT z<~lm?X2;;MkQIo!Q29w_z;y&qq>)%SB7BXt|}W#gI}S-y)iz8;Il1 z-o>hh2IzjtiQ95pZXw69(@gw&Ymr))`_la2o5*RqeOhJwD#W4CR~W)kfR4=jTzW;> z1o4`$d+ZsLj|R7_lo(lTjO6&@y+t%`qIViidO3XtXsOkkA;;&1D8Fq&GWzIxM3*`A za&YDr+QV9CyUyPb-BsJ{Fx7Ymah;v7xvy|DlIhS&5dL@^B<$XkM}!`Q=b zKeIr~b)(ayT1(LhHo@bErESpr%js{!_)5{dr&=DFty>WTtJ8jOsxlOolZGpHv_bQB z86UdnRf_5^iR!t0u|oM&4W_iy_s~Whec+W{PG~sX3uwmfp@VW8UzTonLOgX+6$ud) zXj1$3)}1cfP_fWF4&8?ps3!GTu*FhmGWKSG=Ct8Ue6+KsqI7)rn z!s+yu$<3%ehF{UzF9f}?wOJA{(~SIsxLa?kg(7-~62{}|t!Q`6C$YPeZv`z z?2bLkn&Iz6Ya3RV+;cyTj*MQa3^49Q+)LNksXjc7+Q*LCB$aleMfP5AKWvCWx23`j z$1I;CsmK;l=Gs`a+`K|_e@+*AT-6%OI&}_t?NiU&8~XyW$i@!6fs3fLTH%xPo)^f< zV*bVD@$sm;)^%D(yBobIvbrtvFaepEh|V)~>PCXjMLTFeBp};KdZ7~29z?aV-`lM* z1-)D&>B!b=MGYsg_ve;z6|jGOigy za*=(5Yj8mGAX0MJ97^YP9rZ&@qxrER#F?oby?JLIVvajy_7e97>2wb6yvlt8-5R}P z6~*utU15K9B3daQt#^KtA~Z9M9D-IKo~MwH66dMqT;(4@1=p9;mT(p#BTO=JZ^sDA z@LD@_`)Co;*=+Mdr1BjypJLH0s4Yf)Pd>;oZhwb*RxDocz4b0?+4t$iZh=wc8)TFR zUPQfLvc5J%W(;k5nc^{ISB4zVIPA+B9YYT~c?SyLl%x5heA%7*$B|Q}|I=>v3e>wg z=;o!n6R1|K?@LrpC9?mxDo>>%-)F$Kjj79q+3QJ0T%KoMfSx;;X)|Z|_-j?t0 z_$f3X1Fcir0-~o-;SHl(gM|%MkuS)+A}9K2St}A|+C2J7ZU&Wh90741C6NRBj+-oqWkzfDx1$Zc;@TXR`8+&RUeVP_|cULGK-Rs zk^M9Dp#u3&(o;igjN?p^UKd(Ac&g0z6E#$XO1yo*(~Ux2UYA<(&`G%(X1UlJSl5*bw34ay18Le2gW#*?}|NOd=NCp6H)mm@N?9LEX+)s17tZtu~E-qj0z{A@sL@2Nhh z`T=b!tI;d%W`k$CHaE3)eney2_o~x`vqL-HreXcIX*5(^`DKV|0rVL}Q`wAvLf*Xl z6O^?UfcoY2Vvn;wqY;L)x%6)qK)1^?#k;&X*lFDy9c#n^=dh}(t@m-z&-uW#;3)@) z-x?75@{kI~9vN)s+Q$hCs<#J^<7hytNYUoK2^WZ;5vY~VpoPsxGuTZ}Z~=W^e^itl z9R#g7{Yo0Rq1osPR(_rycI>Hjk0IfX@IimEn8=|X7Es+Ne6ZjGAMEEnc7uud`#KgC)e;|;@EsUAt>2iwP@VctA==u`GRw7r)fijI%db}z+4X@6H5uvl%@Q`;p1i$s}}&&rBI>TYK7Cx`HmGN0y* zT8t>@yQw^hT)hac%>%caiK4Ll0H>LP%_8tIx6l|E6$R&XpWR<3iSc`J!U{%?&;WqfT24yWbsC^xzdEo|Ew;e!sdH zT&&KO%vKSD+fAz#GF`-AmqM6YV1gK^H^jd-S|bilcd6xXQi+36OS$rk1rp%b=;gH4 zPaN7r-!$XdC85J{p!HOPIHa^}mp8AI1Ra&4_~03F*z*SQ9np}2)U6hqDvc%JMewy( zjzdyVpW&c#u*z!;bdNV(AEHIXV&CB)j<{K-+o_8jU z2ef2C%@o_85vBr7FSZPFvdhA?K#PwP#wu`6-?>-UK^6uhZW?U3rves-BVAID%R-to z?dH#4RN(pWsS7McvVdrO4We~a;Yo~yDyy6v$Xvr6{nV}syxY!aTl&j^wpBqrGC=Z+wJq6sx>cFm7z*nI! z4>H?6c8q$cL)rzJ(j&q0P{*jS*}YPo`2Ds2*$@5l;Msj`l-o)J4lI?wy>3Jv#xB#9 zJ#^6kY({?Pbch1*LEzNj8cnDYVHK3xssO#|6*EWXYXa_QlUCfg0*FuCDr)$k2}VM< zBqXI3L0q2Qj>20WqX6~SZUk|x zU}7(t7^kiT%T?DOxn-~f1Q$E|@03-7K5NMk56LCKxZq)yyuK1_lv{eUasLuvxltNn zqpAdAgU4#$V@u%Ge!87zQc7@YXxvp;WC=vn7|JKEQ-aL0l%Z3}OJH>6$aDg$5}f?F z=Ht;CZFqD~eIbjj5{R~Xy3w^Q0bFI4j~zh??uDdJe0Ep@Q@)=BMbwm_<09W}FTEv@ zpJh5Q#YU7HI<=B_kF>$Zq<_71wjx|nRepDtPa8Ik=rPZBnJ z6Rk_+Au~qa$!052epGhSpIj~vt2OJMmpiM2${lsd#6Y6FT~arnFj0fdnzWAfTjhX1 z>A91xh#JJ&$&1^G%Yo=Ne*F)1st|s45{irV^;+eXM?{Dh&%h*A80WSA<;5K2f3$Gv#XV0zf?UyVhgk;5FX?%&JhM5=)9ARmY9$pOCfmf5Uor|Gle7&i<;$oO`zZ7(6 zyC@8+7vsY!MPVxBiXJ2r`}2`bL6x(jaBAxKxxL1VfJWr(?N_sO0{?h8LK+89utI3#?pCF+wibAKQrv$&_XD7da#a*1rL#7 zT&(ZP`QZgkjODuB!Z3P}+b3`tFZi92yuErt2#WKJx5n@0fxP=vW{$&x@aO>d&V`d) z5LK2msL3b@r>&bh^s+cX!p3WL^f>{Ddl_WCmxB}J&+M@lp%H+E{R0|_%Q+ym=h|?0 z89%5<74$g7vBTp%ebpAv`Cz;IlXz7+Hn>rJEztNGFP!benC3^az%5>NHLGG?cq3o$ zXs5>l88N9QB9=TrZ#IA576vA;Pq}ZwqQe8hGjHT<${Atj7l#d(UvPsKTjs&F;|!4a z@l3CGB^TJ8-S8^#JRP`+gy#(Sa)HCzRC9l7I^YX^{!E053ohuU96ZlV3mOfaRX#qP zz`L`9AUH}5Jr#RcRxIHJr}~Qi{4i=TsVU#ohH-$@s+}VCDyk;YvIx_nJ+iDo%3b-eY79l^hU;-gvA}(^w?`uU$C0?A3g_@0X1KbsBVyIi zJ47uJ=~*?*1dP@#`bCdN&~o1tAGS|S;4h?Abh2jzt*+Sd^oAcJm>B6OwNAf5Q{)F`eW=#ReC5d zJTY98H-Km)PY4uh(F3)(=doRl18CP92HUmqbdWIu8qC}Kk*svfIIN$WWToyP8>!)UVN$4>6A;HuI(hX!$+V zugSzC`DhefNqm$uzo`^42<)?ej2lA|tUFW5sA$+YynY zQty!F%2P&Ok$2Dv$6(yn%n?)@v_g9YYY__985Imxd5fyriz7Yb@=?8G$fbqa!zgsz zDEhJBEo9Zpnm)a47{%LPe(`e6Eo5AH<7txEAY!@QbH+M67mX~E-=r)%gxce_9=5{e zp+u8o&z_%ujjGI40_>-+A$&g`aXuW6Y4adGK$f0ad;nJ}|w3p+g5wrMpG+t>u z!t`YT<>cQZo*i@*p<)fmj?n>>DAZsomwgp^oN4`ZXJ0>xyeYVs@og63SZ0JVYYw0? z3m3Upr?L@u6!sYf^&$6gxyPR(Gtk*xGz?Qxedy`68+@^M(~+Nzg!S6%Jt#9k)^%X` z5<2|aFNATr2U)kWs(nVOXm^27^5v2qB=GX_Q^OaR(DTxF(QQR95h2(>`rYwlw05%r z&jbH%6f^U_vs*3^J(SXM^EZ8obOws)WjK-%?{H^<)u(PW$}zL#kzo=#khj35o;Ytd zW*Z(TuS-G0r!_83IrX5@zTw-e`clxyyOhe+mtUgG!_>{$3CXB&op%L4O%EFKG|}iW zAlkhh;uhO^_aN2{`#mCrFQL)v7mGd_w8h__a;3ryo~4@pIqPfs1J?NQw@f(W}+O`=$f6YD9!=s2ZEjDzuh6|h@-n&RjKNotiDN=!JNp1SKM}$0BAAV0 z=jT_M?HE86du_Yhd9NVxDsKYsnOA5H{=WOJ^c?h>o9k(yZa;cKvqfVAWT6vj=^{aC z1H^Vcf6p=X3fj(HoFQ|$AKeN$w^5e+GMW}WUcXjm0Ns7)K9+tU8|`S5s3>^Vi)4zn z#7WqsqlG=Tw8BgK(0+~~wH}5H6eD#zLFPpdDtp`)WBBe8nhN5>`CRTn?$W+TWTEP54HnHo%ouYH9Kyx2WG@i{2?2Cu-h zbT?`dP>KsrPCz1OpX}SX?-{Zvnb&t_;~BKcugJ>q@(Uzy-AOPwdJ&!P$c)PxX9tTN zIUzcaJ|SiQ+SLa6CgpwXWsqg>raev?bCqv-(-$5mi%RdnQpRLHR z`@8?{eojAfzemwbQa!%!_x@NCPwMX&an5AeV~_nDe={<@MLm&*w6X2$mINA zz!FHwFXHFqml!K3#t-H6UhnEjILt7%_E=T;gAMv%dxA!I8LCowQjz0VPT;e^uy@*dHlDvIv teH=YKh{x=Bf6a+9Ze%Ct`=8VAgg<5b{pavV{~HgY%+2+8Ie$t#{vUAN5S#!2 literal 0 HcmV?d00001 From 566d393118ea1c0178bb6b57ca1ae8bf209e54b4 Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Fri, 5 Apr 2024 18:25:56 -0500 Subject: [PATCH 2/6] Rename Neutron_flux_vs_energy.png to SphericalShell/Neutron_flux_vs_energy.png --- .../Neutron_flux_vs_energy.png | Bin 1 file changed, 0 insertions(+), 0 deletions(-) rename Neutron_flux_vs_energy.png => SphericalShell/Neutron_flux_vs_energy.png (100%) diff --git a/Neutron_flux_vs_energy.png b/SphericalShell/Neutron_flux_vs_energy.png similarity index 100% rename from Neutron_flux_vs_energy.png rename to SphericalShell/Neutron_flux_vs_energy.png From f973c848b1f835c6e1e2fb16b81f00227f927c2a Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Fri, 5 Apr 2024 18:26:16 -0500 Subject: [PATCH 3/6] Rename SS_h5_Reader.py to SphericalShell/SS_h5_Reader.py --- SS_h5_Reader.py => SphericalShell/SS_h5_Reader.py | 0 1 file changed, 0 insertions(+), 0 deletions(-) rename SS_h5_Reader.py => SphericalShell/SS_h5_Reader.py (100%) diff --git a/SS_h5_Reader.py b/SphericalShell/SS_h5_Reader.py similarity index 100% rename from SS_h5_Reader.py rename to SphericalShell/SS_h5_Reader.py From 83a0a104d22dbe7770f5385cd44a435d91e731f5 Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Fri, 5 Apr 2024 18:26:46 -0500 Subject: [PATCH 4/6] Rename statepoint.10.h5 to SphericalShell/statepoint.10.h5 --- statepoint.10.h5 => SphericalShell/statepoint.10.h5 | Bin 1 file changed, 0 insertions(+), 0 deletions(-) rename statepoint.10.h5 => SphericalShell/statepoint.10.h5 (100%) diff --git a/statepoint.10.h5 b/SphericalShell/statepoint.10.h5 similarity index 100% rename from statepoint.10.h5 rename to SphericalShell/statepoint.10.h5 From e015c7dce92a87e94ab1fec965ffb9c12f2dd1eb Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Fri, 5 Apr 2024 18:27:16 -0500 Subject: [PATCH 5/6] Rename SphericalShell.py to SphericalShell/SphericalShell.py --- SphericalShell.py => SphericalShell/SphericalShell.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) rename SphericalShell.py => SphericalShell/SphericalShell.py (98%) diff --git a/SphericalShell.py b/SphericalShell/SphericalShell.py similarity index 98% rename from SphericalShell.py rename to SphericalShell/SphericalShell.py index 21dcb00..4dba3e9 100644 --- a/SphericalShell.py +++ b/SphericalShell/SphericalShell.py @@ -67,4 +67,4 @@ # Collecting and exporting tallies to .xml tallies = openmc.Tallies([neutron_tally, spectrum_tally]) -tallies.export_to_xml() \ No newline at end of file +tallies.export_to_xml() From d0efcb25b7117b16e5b231d398e92bd9043ccea3 Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Tue, 14 May 2024 16:07:03 -0500 Subject: [PATCH 6/6] Update SphericalShell/SphericalShell.py --- SphericalShell/SphericalShell.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/SphericalShell/SphericalShell.py b/SphericalShell/SphericalShell.py index 4dba3e9..d0860b5 100644 --- a/SphericalShell/SphericalShell.py +++ b/SphericalShell/SphericalShell.py @@ -42,7 +42,7 @@ settings.batches = 10 settings.inactive = 1 settings.particles = 100000 -settings.source = openmc.Source(space=PointSource, energy=Prob, strength = 10.0, particle = 'neutron') +settings.source = openmc.Source(space=PointSource, energy=Prob, strength = 1.0, particle = 'neutron') settings.run_mode = 'fixed source' settings.export_to_xml()