import openmc

mats = openmc.Materials()

mat = openmc.Material(1)
mat.name = "Core"
mat.set_density('sum')
mat.add_nuclide('Pu240', 1.06605E-04)
mat.add_nuclide('Pu241', 8.37189E-06)
mat.add_nuclide('U235', 6.02719E-06)
mat.add_nuclide('U238', 2.60196E-03)
mat.add_nuclide('Pu239', 1.65490E-03)
mat.add_nuclide('Pu242', 3.42379E-07)
mat.add_nuclide('Am241', 2.91589E-06)
mat.add_element('Cr', 1.84075E-03)
mat.add_element('Ni', 7.90457E-04)
mat.add_element('Fe', 7.27891E-03)
mat.add_element('Al', 1.10408E-04)
mat.add_element('C', 5.56078E-02)
mat.add_element('Mo', 2.09495E-04)
mat.add_nuclide('Mn55', 9.83153E-05)
mat.add_element('Cu', 2.79163E-06)
mat.add_element('Si', 9.24362E-05)
mats.append(mat)

mat = openmc.Material(2)
mat.name = "Axial blanket"
mat.set_density('sum')
mat.add_nuclide('U235', 8.28476E-05)
mat.add_nuclide('U238', 3.96350E-02)
mat.add_element('Cr', 1.44924E-03)
mat.add_element('Ni', 6.18487E-04)
mat.add_element('Fe', 5.82611E-03)
mat.add_element('C', 1.37151E-05)
mat.add_nuclide('Mn55', 7.06576E-05)
mat.add_nuclide('H1', 3.97871E-06)
mat.add_element('Si', 7.39858E-05)
mat.add_element('Cl', 6.87973E-06)
mat.add_nuclide('F19', 2.03717E-05)
mats.append(mat)

mat = openmc.Material(3)
mat.name = "Radial blanket"
mat.set_density('sum')
mat.add_nuclide('U235', 8.34523E-05)
mat.add_nuclide('U238', 3.99253E-02)
mat.add_element('Cr', 1.10716E-03)
mat.add_element('Ni', 4.52773E-04)
mat.add_element('Fe', 4.51901E-03)
mat.add_element('C', 8.40260E-06)
mat.add_nuclide('Mn55', 4.35132E-05)
mat.add_nuclide('H1', 2.43087E-06)
mat.add_element('Si', 6.05928E-05)
mat.add_element('Cl', 4.21476E-06)
mat.add_nuclide('F19', 1.24811E-05)
mats.append(mat)

mat = openmc.Material(4)
mat.name = "Drawer gap"
mat.set_density('sum')
mat.add_element('Cr', 3.45718E-03)
mat.add_element('Ni', 1.59962E-03)
mat.add_element('Fe', 1.92112E-02)
mat.add_element('C', 2.62585E-04)
mat.add_nuclide('Mn55', 2.34487E-04)
mat.add_element('Si', 1.50898E-04)
mats.append(mat)

mat = openmc.Material(5)
mat.name = "Matrix"
mat.set_density('sum')
mat.add_element('Cr', 1.10960E-03)
mat.add_element('Ni', 4.53554E-04)
mat.add_element('Fe', 4.52968E-03)
mat.add_nuclide('Mn55', 4.34957E-05)
mat.add_element('Si', 6.07700E-05)
mats.append(mat)

mats.export_to_xml()
