"Throughput capacity [A.h]" is different for sol.cycles[0] and sol.cycles[0].steps[0]

[RuiheLi]

Dear PyBaMM,

I am using PyBaMM 22.9 and Python 3.9.6.

I find when using that set_initial_conditions_from, "Throughput capacity [A.h]" is different for sol.cycles[0] and sol.cycles[0].steps[0] for a later solution. Figure below:

Script below:

import pybamm 
import matplotlib as mpl
fs=17;
font = {'family' : 'DejaVu Sans','size'   : fs}
mpl.rc('font', **font)
V_max = 4.0;  V_min = 2.7; 
model_0 = pybamm.lithium_ion.DFN(options={
    "calculate discharge energy":"true", },)
exp_AGE_text = [(f"Discharge at 1 C until {V_min} V", 
    f"Charge at 1 C until {V_max} V", ),  ]
Experiment    = pybamm.Experiment( exp_AGE_text * 1)
Para_0 = pybamm.ParameterValues("OKane2022" )

Sim_0= pybamm.Simulation(
    model_0,
    experiment = Experiment, 
    parameter_values=Para_0, 
    solver = pybamm.CasadiSolver(),)
Sol_0 = Sim_0.solve(calc_esoh=False,) 

model_1 = model_0.set_initial_conditions_from(Sol_0, inplace=False)
Sim_1= pybamm.Simulation(
    model_1,
    experiment = Experiment, 
    parameter_values=Para_0, 
    solver = pybamm.CasadiSolver(),)
Sol_1 = Sim_1.solve(calc_esoh=False,) 

model_2 = model_0.set_initial_conditions_from(Sol_1, inplace=False)
Sim_2= pybamm.Simulation(
    model_2,
    experiment = Experiment, 
    parameter_values=Para_0, 
    solver = pybamm.CasadiSolver(),)
Sol_2 = Sim_2.solve(calc_esoh=False,) 

label = ["Sol_0","Sol_1","Sol_2",] 
output_variables3 = [
    "Throughput energy [W.h]",
    "Throughput capacity [A.h]", 
    "Terminal voltage [V]", 
]
quick_plot = pybamm.QuickPlot(
    [Sol_0,Sol_1,Sol_2,], output_variables3,label,variable_limits='fixed',
    time_unit='hours',n_rows=1) #figsize = (18,12),
quick_plot.dynamic_plot()
label = ["Sol_2 whole","Sol_2 cycles","Sol_2 steps",] 
quick_plot = pybamm.QuickPlot(
    [Sol_2,Sol_2.cycles[0],Sol_2.cycles[0].steps[0]], output_variables3,label,variable_limits='fixed',
    time_unit='hours',n_rows=1) #figsize = (18,12),
quick_plot.dynamic_plot()

Do you have any ideas on this? Big thanks!

[valentinsulzer]

I got something different running your code

[RuiheLi]

Sorry for the unclearness! The figure you showed is the first plot in the script. If you look at the second one, which corresponds to the last five lines, you may see the figure I post.

[valentinsulzer]

Ok, this is a bug, and gonna be a pain to fix. Converting this to issue.