import openmc

mats = openmc.Materials()

mat = openmc.Material(1)
mat.name = "Inner Core 1"
mat.set_density('sum')
mat.add_nuclide('Pu240', 1.17624E-04)
mat.add_nuclide('Pu241', 1.33836E-05)
mat.add_nuclide('U235', 1.26065E-05)
mat.add_nuclide('U238', 5.79290E-03)
mat.add_nuclide('Pu239', 8.86528E-04)
mat.add_nuclide('Pu238', 3.33696E-07)
mat.add_nuclide('Pu242', 1.40289E-06)
mat.add_nuclide('Am241', 2.95925E-06)
mat.add_element('Cr', 2.69344E-03)
mat.add_element('Ni', 1.19793E-03)
mat.add_element('Fe', 1.28611E-02)
mat.add_element('Al', 4.02385E-06)
mat.add_element('Na', 9.27752E-03)
mat.add_element('O', 1.37526E-02)
mat.add_element('C', 3.66446E-05)
mat.add_element('Mo', 2.36059E-04)
mat.add_element('Mn', 2.25651E-04)
mat.add_element('Cu', 2.46789E-05)
mat.add_element('Si', 1.62374E-04)
mat.add_element('Cl', 2.98410E-07)
mat.add_element('Co', 8.32631E-07)
mats.append(mat)

mat = openmc.Material(2)
mat.name = "Outer Core 1"
mat.set_density('sum')
mat.add_nuclide('Pu240', 1.17751E-04)
mat.add_nuclide('Pu241', 1.52377E-05)
mat.add_nuclide('U235', 1.26600E-05)
mat.add_nuclide('U238', 5.81701E-03)
mat.add_nuclide('Pu239', 8.87794E-04)
mat.add_nuclide('Pu238', 4.70056E-07)
mat.add_nuclide('Pu242', 1.76624E-06)
mat.add_nuclide('Am241', 2.40381E-06)
mat.add_element('Cr', 2.68925E-03)
mat.add_element('Ni', 1.18886E-03)
mat.add_element('Fe', 1.32494E-02)
mat.add_element('Al', 3.62390E-06)
mat.add_element('Na', 9.12725E-03)
mat.add_element('O', 1.42712E-02)
mat.add_element('C', 3.40434E-05)
mat.add_element('Mo', 2.39077E-04)
mat.add_element('Mn', 2.24037E-04)
mat.add_element('Cu', 2.48329E-05)
mat.add_element('Si', 1.54945E-04)
mat.add_element('Cl', 2.97156E-07)
mat.add_element('Co', 1.69138E-07)
mats.append(mat)

mat = openmc.Material(4)
mat.name = "Inner Axial reflector 1"
mat.set_density('sum')
mat.add_nuclide('U235', 8.29314E-05)
mat.add_nuclide('U238', 3.71138E-02)
mat.add_element('Cr', 1.64641E-03)
mat.add_element('Ni', 6.89374E-04)
mat.add_element('Fe', 6.06025E-03)
mat.add_element('C', 3.67281E-05)
mat.add_element('Mo', 1.06205E-05)
mat.add_element('Mn', 1.43193E-04)
mat.add_element('Cu', 1.97417E-05)
mat.add_element('Si', 1.06228E-04)
mats.append(mat)

mat = openmc.Material(6)
mat.name = "Inner Axial reflector 2"
mat.set_density('sum')
mat.add_nuclide('U235', 8.65264E-05)
mat.add_nuclide('U238', 3.87479E-02)
mat.add_element('Cr', 1.44970E-03)
mat.add_element('Ni', 6.07516E-04)
mat.add_element('Fe', 5.25530E-03)
mat.add_element('C', 2.71677E-05)
mat.add_element('Mo', 1.00416E-05)
mat.add_element('Mn', 1.27167E-04)
mat.add_element('Cu', 1.90679E-05)
mat.add_element('Si', 9.52031E-05)
mats.append(mat)

mat = openmc.Material(5)
mat.name = "Outer Axial Reflector 1"
mat.set_density('sum')
mat.add_nuclide('U235', 8.27657E-05)
mat.add_nuclide('U238', 3.70414E-02)
mat.add_element('Cr', 1.65506E-03)
mat.add_element('Ni', 6.92837E-04)
mat.add_element('Fe', 6.08815E-03)
mat.add_element('C', 3.66405E-05)
mat.add_element('Mo', 1.06735E-05)
mat.add_element('Mn', 1.43788E-04)
mat.add_element('Cu', 1.98180E-05)
mat.add_element('Si', 1.06087E-04)
mats.append(mat)

mat = openmc.Material(7)
mat.name = "Outer Axial Reflector 2"
mat.set_density('sum')
mat.add_nuclide('U235', 8.55655E-05)
mat.add_nuclide('U238', 3.83153E-02)
mat.add_element('Cr', 1.46821E-03)
mat.add_element('Ni', 6.15066E-04)
mat.add_element('Fe', 5.31966E-03)
mat.add_element('C', 2.73042E-05)
mat.add_element('Mo', 1.02046E-05)
mat.add_element('Mn', 1.28523E-04)
mat.add_element('Cu', 1.92143E-05)
mat.add_element('Si', 9.54632E-05)
mats.append(mat)

mat = openmc.Material(8)
mat.name = "Radial Reflector 1"
mat.set_density('sum')
mat.add_nuclide('U235', 8.66555E-05)
mat.add_nuclide('U238', 3.88076E-02)
mat.add_element('Cr', 1.17821E-03)
mat.add_element('Ni', 4.75780E-04)
mat.add_element('Fe', 4.23998E-03)
mat.add_element('C', 1.85721E-05)
mat.add_element('Mo', 8.18136E-06)
mat.add_element('Mn', 1.04936E-04)
mat.add_element('Cu', 1.70748E-05)
mat.add_element('Si', 6.76848E-05)
mats.append(mat)

mat = openmc.Material(9)
mat.name = "Radial Reflector 2"
mat.set_density('sum')
mat.add_element('Cr', 1.19183E-03)
mat.add_element('Ni', 4.82115E-04)
mat.add_element('Fe', 4.31023E-03)
mat.add_element('C', 1.98501E-05)
mat.add_element('Mo', 8.27748E-06)
mat.add_element('Mn', 1.06058E-04)
mat.add_element('Cu', 1.71921E-05)
mat.add_element('Si', 6.88696E-05)
mats.append(mat)

mat = openmc.Material(10)
mat.name = "Empty Matrix"
mat.set_density('sum')
mat.add_element('Cr', 1.19183E-03)
mat.add_element('Ni', 4.82115E-04)
mat.add_element('Fe', 4.31023E-03)
mat.add_element('C', 1.98501E-05)
mat.add_element('Mo', 8.27748E-06)
mat.add_element('Mn', 1.06058E-04)
mat.add_element('Cu', 1.71921E-05)
mat.add_element('Si', 6.88696E-05)
mats.append(mat)

mats.export_to_xml()
