def METHOD_NAME(field, data):
    nucleus_densityIa = (
        data["gas", "metal_ia_density"] * SNIa_abundance[atom]
    )
    nucleus_densityII = (
        data["gas", "metal_ii_density"] * SNII_abundance[atom]
    )
    return nucleus_densityIa + nucleus_densityII