analysis.plot(
        model=submdl,
        time=time,
        volume=volume,
        speciesCounts=speciesCounts,