import openmc

mats = openmc.Materials()

mat = openmc.Material(1)
mat.name = "Inner Core 1"
mat.set_density('sum')
mat.add_nuclide('Pu240', 3.2178E-04)
mat.add_nuclide('Pu241', 5.3569E-05)
mat.add_nuclide('U235', 1.1872E-05)
mat.add_nuclide('U238', 5.5620E-03)
mat.add_nuclide('Pu239', 8.3844E-04)
mat.add_nuclide('Pu238', 1.1083E-06)
mat.add_nuclide('Pu242', 1.7447E-05)
mat.add_nuclide('Am241', 1.4681E-05)
mat.add_element('Cr', 2.6792E-03)
mat.add_element('Ni', 1.1908E-03)
mat.add_element('Fe', 1.2805E-02)
mat.add_element('Al', 4.1057E-06)
mat.add_element('Na', 9.1676E-03)
mat.add_element('O', 1.3684E-02)
mat.add_element('C', 3.6621E-05)
mat.add_element('Mo', 2.3662E-04)
mat.add_element('Mn', 2.2447E-04)
mat.add_element('Cu', 2.4580E-05)
mat.add_element('Si', 1.6182E-04)
mat.add_element('Ca', 2.1044E-06)
mat.add_element('Cl', 2.9484E-07)
mat.add_element('Co', 8.2889E-07)
mats.append(mat)

mat = openmc.Material(2)
mat.name = "Outer Core 1"
mat.set_density('sum')
mat.add_nuclide('Pu240', 1.1749E-04)
mat.add_nuclide('Pu241', 1.4641E-05)
mat.add_nuclide('U235', 1.2575E-05)
mat.add_nuclide('U238', 5.7820E-03)
mat.add_nuclide('Pu239', 8.8583E-04)
mat.add_nuclide('Pu238', 4.6444E-07)
mat.add_nuclide('Pu242', 1.7620E-06)
mat.add_nuclide('Am241', 2.9643E-06)
mat.add_element('Cr', 2.6909E-03)
mat.add_element('Ni', 1.1885E-03)
mat.add_element('Fe', 1.3230E-02)
mat.add_element('Al', 1.8593E-05)
mat.add_element('Na', 9.0718E-03)
mat.add_element('O', 1.4192E-02)
mat.add_element('C', 3.3605E-05)
mat.add_element('Mo', 2.3875E-04)
mat.add_element('Mn', 2.2395E-04)
mat.add_element('Cu', 2.4882E-05)
mat.add_element('Si', 1.5377E-04)
mat.add_element('Ca', 2.0828E-06)
mat.add_element('Cl', 2.9546E-07)
mat.add_element('Co', 5.4095E-08)
mats.append(mat)

mat = openmc.Material(3)
mat.name = "Pu-Al ring"
mat.set_density('sum')
mat.add_nuclide('Pu240', 5.0694E-05)
mat.add_nuclide('Pu241', 2.8433E-06)
mat.add_nuclide('U235', 1.2151E-05)
mat.add_nuclide('U238', 5.7042E-03)
mat.add_nuclide('Pu239', 1.0697E-03)
mat.add_nuclide('Pu238', 1.6248E-09)
mat.add_nuclide('Pu242', 1.1005E-07)
mat.add_nuclide('Am241', 1.9758E-06)
mat.add_element('Cr', 2.6230E-03)
mat.add_element('Ni', 1.1648E-03)
mat.add_element('Fe', 1.3160E-02)
mat.add_element('Al', 1.1477E-04)
mat.add_element('Na', 9.3420E-03)
mat.add_element('O', 1.4386E-02)
mat.add_element('C', 4.0021E-05)
mat.add_element('Mo', 1.3306E-05)
mat.add_element('Mn', 2.2020E-04)
mat.add_element('Cu', 2.5731E-05)
mat.add_element('Si', 1.5879E-04)
mat.add_element('Ca', 2.1452E-06)
mat.add_element('Cl', 3.0352E-07)
mat.add_element('Co', 8.2412E-07)
mats.append(mat)

mat = openmc.Material(4)
mat.name = "Inner Axial reflector 1"
mat.set_density('sum')
mat.add_nuclide('U235', 8.1353E-05)
mat.add_nuclide('U238', 3.7936E-02)
mat.add_element('Cr', 1.6464E-03)
mat.add_element('Ni', 6.8937E-04)
mat.add_element('Fe', 6.0603E-03)
mat.add_element('C', 3.6728E-05)
mat.add_element('Mo', 1.0620E-05)
mat.add_element('Mn', 1.4319E-04)
mat.add_element('Cu', 1.9742E-05)
mat.add_element('Si', 1.0623E-04)
mats.append(mat)

mat = openmc.Material(6)
mat.name = "Inner Axial reflector 2"
mat.set_density('sum')
mat.add_nuclide('U235', 8.5866E-05)
mat.add_nuclide('U238', 4.0055E-02)
mat.add_element('Cr', 1.4497E-03)
mat.add_element('Ni', 6.0752E-04)
mat.add_element('Fe', 5.2553E-03)
mat.add_element('C', 2.7168E-05)
mat.add_element('Mo', 1.0042E-05)
mat.add_element('Mn', 1.2717E-04)
mat.add_element('Cu', 1.9068E-05)
mat.add_element('Si', 9.5203E-05)
mats.append(mat)

mat = openmc.Material(5)
mat.name = "Outer Axial Reflector 1"
mat.set_density('sum')
mat.add_nuclide('U235', 8.1048E-05)
mat.add_nuclide('U238', 3.7795E-02)
mat.add_element('Cr', 1.6557E-03)
mat.add_element('Ni', 6.9310E-04)
mat.add_element('Fe', 6.0903E-03)
mat.add_element('C', 3.6634E-05)
mat.add_element('Mo', 1.0678E-05)
mat.add_element('Mn', 1.4383E-04)
mat.add_element('Cu', 1.9824E-05)
mat.add_element('Si', 1.0608E-04)
mats.append(mat)

mat = openmc.Material(7)
mat.name = "Outer Axial Reflector 2"
mat.set_density('sum')
mat.add_nuclide('U235', 8.4740E-05)
mat.add_nuclide('U238', 3.9527E-02)
mat.add_element('Cr', 1.4696E-03)
mat.add_element('Ni', 6.1564E-04)
mat.add_element('Fe', 5.3245E-03)
mat.add_element('C', 2.7315E-05)
mat.add_element('Mo', 1.0217E-05)
mat.add_element('Mn', 1.2863E-04)
mat.add_element('Cu', 1.9225E-05)
mat.add_element('Si', 9.5483E-05)
mats.append(mat)

mat = openmc.Material(8)
mat.name = "Radial Reflector 1"
mat.set_density('sum')
mat.add_nuclide('U235', 8.5846E-05)
mat.add_nuclide('U238', 4.0046E-02)
mat.add_element('Cr', 1.1782E-03)
mat.add_element('Ni', 4.7578E-04)
mat.add_element('Fe', 4.2400E-03)
mat.add_element('C', 1.8572E-05)
mat.add_element('Mo', 8.1814E-06)
mat.add_element('Mn', 1.0494E-04)
mat.add_element('Cu', 1.7075E-05)
mat.add_element('Si', 6.7685E-05)
mats.append(mat)

mat = openmc.Material(9)
mat.name = "Radial Reflector 2"
mat.set_density('sum')
mat.add_element('Cr', 1.1910E-03)
mat.add_element('Ni', 4.8172E-04)
mat.add_element('Fe', 4.3056E-03)
mat.add_element('C', 1.9752E-05)
mat.add_element('Mo', 8.2719E-06)
mat.add_element('Mn', 1.0599E-04)
mat.add_element('Cu', 1.7186E-05)
mat.add_element('Si', 6.8789E-05)
mats.append(mat)

mat = openmc.Material(10)
mat.name = "Empty Matrix"
mat.set_density('sum')
mat.add_element('Cr', 1.1910E-03)
mat.add_element('Ni', 4.8172E-04)
mat.add_element('Fe', 4.3056E-03)
mat.add_element('C', 1.9752E-05)
mat.add_element('Mo', 8.2719E-06)
mat.add_element('Mn', 1.0599E-04)
mat.add_element('Cu', 1.7186E-05)
mat.add_element('Si', 6.8789E-05)
mats.append(mat)

mats.export_to_xml()
