analysis.plot(
        model=mdl,
        time=time,
        yDatas={'Protein': totalProt},
        fileName=os.path.join(output_directory, 'Total-protein.png')