Oxygen Diffusion submodel

hello,
can someone help me with my code i tried several ways to plot the Oxygen concentration but have’nt succed
`import pybamm

model = pybamm.lithium_ion.DFN()

Set up the oxygen diffusion submodel

oxygen_diffusion = pybamm.oxygen_diffusion.Full(model.param)
particle_mechanics = pybamm.particle_mechanics.CrackPropagation(model.param, domain =“positive”, x_average = “full”,
options ={“operating mode”: “CCCV”} )
thermal = pybamm.thermal.lumped.Lumped(model.param)

Define experiments (e.g., cycling with degradation)

experiment = pybamm.Experiment([
“Discharge at 1C until 3.0 V”,
“Charge at 1C until 4.2 V”,
“Hold at 4.2 V for 1 hour”
] * 3) # 3 cycles

Solve

sim = pybamm.Simulation(model, parameter_values=param, experiment=experiment)
sol = sim.solve()

Plot results

pybamm.dynamic_plot(sol, [
“Terminal voltage [V]”,
“Positive total SEI thickness [m]”,
“Cell temperature [C]”,
“Positive oxygen concentration [mol.m-3]”
])
`

That submodel is only for lead-acid batteries

thank you for the clarification