import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import sys
#import plot_func
#reload (plot_func)
import os


data = pd.read_csv('Figure1.csv')

X = data[data.keys()[0]]*1000
Y1 = data[data.keys()[1]]
Y2 = data[data.keys()[2]]*1e6


plt.figure(figsize=(10,7))

plt.subplot(2,1,1)
plt.plot(X, Y1, 'navy' ,linestyle='-',  label = '',linewidth=3)


plt.xlim(0, 1600, 1000)
plt.yticks(np.arange(-70, 40, 20))
plt.ylabel ('V$_m$ (mV)',fontsize=12)
# plt.xlabel ('Time (ms)',fontsize=12)
plt.title('A', loc='left' ,fontsize=20)
#plt.legend(loc = 'best',fontsize=10)

plt.grid()
plt.subplots_adjust(bottom=0.5, right=0.8, top=1)
plt.tight_layout(pad=0.4, w_pad=0.5, h_pad=1.0)

plt.subplot(2,1,2)
plt.plot(X, Y2, 'navy' ,linestyle='-',  label = '',linewidth=3)

plt.xlim(0, 1600, 1000)
plt.yticks(np.arange(50,300,50))
plt.ylabel ('[Ca$^{2+}$]$_i$ (nM)',fontsize=12)
plt.xlabel ('Time (ms)',fontsize=12)
plt.title('B', loc='left',fontsize=20)
plt.grid()
plt.subplots_adjust(bottom=0.5, right=0.8, top=1)
plt.tight_layout(pad=0.4, w_pad=0.5, h_pad=1.0)

plt.savefig('C:/Nima/ABI/Physiome Journal/sinus/Python_codes/figure1.png')
plt.show()