Particle size distribution in custom model

[Nata2340]

Hallo, I intend to incorporate the particle size distribution in a BaseModel (DNF-Model Half cell). I have followed the example for the MP-DFN model (as shown below), but I am not able to generate a plot. Am I missing any important steps or is this not the correct approach? I would appreciate any guidance or answer.

R_av_p_dim = params["Positive particle radius [m]"]

sd_p_dim = 0.6 * R_av_p_dim

R_min_p = 0
R_max_p = 3 * R_av_p_dim

def f_a_dist_p_dim(R):
return pybamm.lognormal(R, R_av_p_dim, sd_p_dim)

distribution_params = {
"Positive minimum particle radius [m]": R_min_p,
"Positive maximum particle radius [m]": R_max_p,
"Positive area-weighted "
+ "particle-size distribution [m-1]": f_a_dist_p_dim,
}
params.update(distribution_params, check_already_exists=False)
sim_custom = pybamm.Simulation(model, parameter_values=params, solver=solver)
sim_custom.solve(t_eval=[0, 3500])
output_variables = [
"X-averaged negative area-weighted particle-size distribution [m-1]",
"X-averaged positive area-weighted particle-size distribution [m-1]",
"Voltage [V]"
]
quickplot = pybamm.QuickPlot(
[sim, sim_custom], output_variables=output_variables, labels=["default lognormals", "custom"]
)
quickplot.plot(0)