How do I read the data in XDMF and use it for the next calculation

Hello, everyone, I want to use the results of the last calculation as the initial value of the next calculation, and the following methods are used to save the file:

V = FunctionSpace(mesh, 'P', 1)
t_total = 600    
num_steps = 13000   
dt = t_total / num_steps   
xdmffile_c = XDMFFile('10_10_0_600s_slab/diffusion.xdmf')
xdmffile_c.parameters["flush_output"] = True
for n in range(num_steps):
    t += dt
    if abs(t - t_total) < 100:
        xdmffile_c.write(c4, t)

How can I extract the value at 600s in the xdmf for the next calculation and calculate its initial value?

See Read mesh from XDMF file (write_checkpoint) - #3 by dokken

I am currently trying to use time series, I don’t know if this method is feasible?

timeseries_c = TimeSeries('10_10_0_10s_slab/diffusion_series')
for n in range(num_steps):
    timeseries_c.store(c4.vector(), t)

I’ve never used the TimeSeries class, so I cannot say anything about the feasibility of the method.

As far as I can tell my quickly glancing at the code, is that it does not store the DOFMAP, thus this can only be used on the same number of processors, i.e. if the original code is executed on 1 process, you cannot execute the second code in parallel. Feel free to test and validate/invalidate my claim.