Plot_Dynamic_Data.py

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

"""
Copyright (c) 2021 by Erhan YILMAZ (https://orcid.org/0000-0003-4788-9022)
This file is licensed under CC BY-SA (Attribution-ShareAlike)
It is provided "as is", without express or implied warranty, for educational and informational purposes only.

This python script plots dynamic data collected for
ThunderSky-Winston-LIFEPO4-100Ah-WIDE cell according to test descriptions below in the files.

http://mocha-java.uccs.edu/BMS1/index.html
http://mocha-java.uccs.edu/ECE5710/index.html
http://mocha-java.uccs.edu/ECE5710/ECE5710-Notes02.pdf
"""

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt


temps=[-5, 5, 15, 25, 35, 45]
mags= [20, 20, 40, 50, 50, 50]

for erhan in range(len(temps)):

    for erhan2 in range(3):
        if(temps[erhan]>0):
            script  = "dyn_data/THUN100_DYN_%02d_P%02d_S%1d.csv" %(mags[erhan], temps[erhan], erhan2+1)
        else:
            script  = "dyn_data/THUN100_DYN_%02d_N%02d_S%1d.csv" %(mags[erhan], np.abs(temps[erhan]), erhan2+1)
            
        data = pd.read_csv(script)
        current = np.asarray(data['current'])
        voltage = np.asarray(data['voltage'])
        time    =  np.asarray(data['time']/3600.0)
        
        fig,a =  plt.subplots(2,1)
        x = np.arange(1,5)
        a[0].plot(time, current)
        a[0].set_title('Dynamic Script %1d T=%02d' % (erhan2+1, temps[erhan]))
        #a[0].set_xlabel('Time(h)')
        a[0].set_ylabel('Current(A)')
        a[1].plot(time, voltage)
        #a[1].set_title('Dynamic Script %1d T=%02d' % (erhan2+1, temps[erhan]))
        a[1].set_xlabel('Time(h)')
        a[1].set_ylabel('Voltage(V)')
        plt.show()