notebook.community

Edit and run 





In [5]:



    


%matplotlib inline


import numpy as np

import matplotlib.pyplot as plt




width = 2
t=np.linspace(0,100,100/0.01+2)

i5hh=np.genfromtxt('/home/jovyan/3/HH/i5.csv', delimiter=',')
#plt.plot(i5hh)
#plt.show()
i15hh=np.genfromtxt('/home/jovyan/3/HH/i15.csv', delimiter=',')
i25hh=np.genfromtxt('/home/jovyan/3/HH/i25.csv', delimiter=',')
i35hh=np.genfromtxt('/home/jovyan/3/HH/i35.csv', delimiter=',')

ina5hh=np.genfromtxt('/home/jovyan/3/HH/ina5.csv', delimiter=',')
ina15hh=np.genfromtxt('/home/jovyan/3/HH/ina15.csv', delimiter=',')
ina25hh=np.genfromtxt('/home/jovyan/3/HH/ina25.csv', delimiter=',')
ina35hh=np.genfromtxt('/home/jovyan/3/HH/ina35.csv', delimiter=',')

ik5hh=np.genfromtxt('/home/jovyan/3/HH/ik5.csv', delimiter=',')
ik15hh=np.genfromtxt('/home/jovyan/3/HH/ik15.csv', delimiter=',')
ik25hh=np.genfromtxt('/home/jovyan/3/HH/ik25.csv', delimiter=',')
ik35hh=np.genfromtxt('/home/jovyan/3/HH/ik35.csv', delimiter=',')


ina=[]
ik=[]

i5_na = i5hh.copy()
i5_k=i5hh.copy()
i5_na[i5hh>0]=0
i5_k[i5hh<0]=0

plt.plot(i5_na)
plt.plot(i5_k)
#plt.show()


area_norm=np.trapz(i5_na,t,dx=0.001)
ina.append(i5_na)
ik.append(i5_k)

print (area_norm)

i15_na = i15hh.copy()
i15_k=i15hh.copy()
i15_na[i15_na>0]=0
i15_k[i15_na<0]=0
ina.append(i15_na)
ik.append(i15_k)

i25_na = i25hh.copy()
i25_k=i25hh.copy()
i25_na[i25_na>0]=0
i25_k[i25_na<0]=0
ina.append(i25_na)
ik.append(i25_k)

i35_na = i35hh.copy()
i35_k=i35hh.copy()
i35_na[i35_na>0]=0
i35_k[i35_na<0]=0
ina.append(i35_na)
ik.append(i35_k)



T = [5,15,25,35]

T=np.array(T)
T2=T+width


areas_na=[]
areas_k=[]
k=0
for temp in T:
	area_na=np.trapz(ina[k],t,dx=0.001)/area_norm
	area_k=np.trapz(ik[k],t,dx=0.001)/np.abs(area_norm)
	areas_na.append(area_na)
	areas_k.append(area_k)
	k=k+1


#############################################################
ax = plt.subplot(3,4,1)
plt.plot(t,ina5hh,label='$T_x = 5^o$C',lw=1)
#plt.text(56,0.8,'$T_x = 5^o$C')
#plt.legend()
#plt.xlabel('Time (ms)')
#plt.ylabel('Na Current (mA/cm$^2$)')
#plt.fill_between(t,hot5, where= hot5>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot5, where= hot5<0, facecolor='green', interpolate=True)
plt.xlim(50,60)
plt.ylim(-5,15)
#ax.set_xticks([50,60])
ax.set_yticks([-5,15])
ax.set_xticklabels([])
ax.set_xticks([])
#plt.ylabel ('nA/cm$^2$', fontsize = 16)
#ax.set_yticklabels([])
#ax.set_yticks([])


ax = plt.subplot(3,4,2)
plt.plot(t,ik5hh,label='$T_x = 5^o$C',lw=1)
#plt.legend()
plt.xlim(50,60)
plt.ylim(-5,15)
#plt.xlabel('Time (ms)')
#plt.ylabel('Na Current (mA/cm$^2$)')
#plt.fill_between(t,hot15, where= hot15>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot15, where= hot15<0, facecolor='green', interpolate=True)
#plt.text(56,0.8,'$T_x = 15^o$C')
ax.set_xticklabels([])
ax.set_xticks([])

ax.set_yticklabels([])
ax.set_yticks([])

ax = plt.subplot(3,4,3)
plt.plot(t,ina5hh + ik5hh,label='$T_x = 25^o$C',lw=1)
#plt.legend()
#plt.xlabel('Time (ms)')
#plt.ylabel('Na Current (mA/cm$^2$)')
#plt.text(56,0.8,'$T_x = 25^o$C')
plt.xlim(50,60)
#plt.fill_between(t,hot25, where= hot25>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot25, where= hot25<0, facecolor='green', interpolate=True)
plt.ylim(-5,15)
#plt.axhline(y=0, xmin=0, xmax=10, linewidth=1, color = 'k', linestyle = '--')
ax.set_xticklabels([])
ax.set_xticks([])
ax.set_yticklabels([])
ax.set_yticks([])

T = 5
ax = plt.subplot(3,4,4)
plt.bar(T-width,areas_na[0], width , color = 'green', label='Net inward current')
plt.bar(T,areas_k[0],  width, color='indianred', label='Net positive current')
ax.set_xticklabels([])
ax.set_xticks([])
plt.ylim(0,50)
ax.set_yticks([0,50])
#plt.text(T -width/2, 8,str(float(areas_k[0]/areas_na[0]))[:4])

print (areas_k[0]/areas_na[0])
#############################################################################

ax = plt.subplot(3,4,5)
plt.plot(t,ina15hh,label='$T_x = 5^o$C',lw=1)
#plt.text(56,0.8,'$T_x = 5^o$C')
#plt.legend()
#plt.xlabel('Time (ms)')
#plt.ylabel('K Current (mA/cm$^2$)')
#plt.fill_between(t,hot5, where= hot5>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot5, where= hot5<0, facecolor='green', interpolate=True)
plt.xlim(50,60)
plt.ylim(-5,15)
ax.set_yticks([-5,15])
ax.set_xticklabels([])
ax.set_xticks([])
#ax.set_yticklabels([])
#ax.set_yticks([])


ax = plt.subplot(3,4,6)
plt.plot(t,ik15hh,label='$T_x = 15^o$C',lw=1)
#plt.legend()
plt.xlim(50,60)
plt.ylim(-5,15)
#plt.xlabel('Time (ms)')
#plt.ylabel('K Current (mA/cm$^2$)')
#plt.fill_between(t,hot15, where= hot15>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot15, where= hot15<0, facecolor='green', interpolate=True)
#plt.text(56,0.8,'$T_x = 15^o$C')

ax.set_xticklabels([])
ax.set_xticks([])
ax.set_yticklabels([])
ax.set_yticks([])

ax = plt.subplot(3,4,7)
plt.plot(t,ik15hh + ina15hh ,label='$T_x = 25^o$C',lw=1)
#plt.axhline(y=0, xmin=0, xmax=10, linewidth=1, color = 'k', linestyle = '--')
#plt.legend()
#plt.xlabel('Time (ms)')
#plt.ylabel('K Current (mA/cm$^2$)')
#plt.text(56,0.8,'$T_x = 25^o$C')
plt.xlim(50,60)
#plt.fill_between(t,hot25, where= hot25>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot25, where= hot25<0, facecolor='green', interpolate=True)
plt.ylim(-5,15)
ax.set_xticklabels([])
ax.set_xticks([])
ax.set_yticklabels([])
ax.set_yticks([])

T = 15
ax = plt.subplot(3,4,8)
plt.bar(T-width,areas_na[1], width , color = 'green', label='Net inward current')
plt.bar(T,areas_k[1],  width, color='indianred', label='Net positive current')
ax.set_xticklabels([])
ax.set_xticks([])
plt.ylim(0,50)
ax.set_yticks([0,50])
print (areas_k[1]/areas_na[1])
#plt.text(T - width/2, 8,str(float(areas_k[1]/areas_na[1]))[:4])

#########################################################################################

ax = plt.subplot(3,4,9)
plt.plot(t,ina25hh,label='$T_x = 5^o$C',lw=1)
#plt.text(56,0.8,'$T_x = 5^o$C')
#plt.legend()
#plt.xlabel('Time (ms)')
#plt.ylabel('Total Current (mA/cm$^2$)')
#plt.fill_between(t,hot5, where= hot5>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot5, where= hot5<0, facecolor='green', interpolate=True)
plt.xlim(50,60)
plt.ylim(-5,15)
ax.set_xticklabels([])
ax.set_xticks([])
ax.set_yticks([-5,15])
#ax.set_yticklabels([])
#ax.set_yticks([])


ax = plt.subplot(3,4,10)
plt.plot(t,ik25hh,label='$T_x = 15^o$C',lw=1)
#plt.legend()
plt.xlim(50,60)
plt.ylim(-5,15)
#plt.xlabel('Time (ms)')
#plt.ylabel('Total Current (mA/cm$^2$)')
#plt.fill_between(t,hot15, where= hot15>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot15, where= hot15<0, facecolor='green', interpolate=True)
#plt.text(56,0.8,'$T_x = 15^o$C')
ax.set_xticklabels([])
ax.set_xticks([])
ax.set_yticklabels([])
ax.set_yticks([])

ax = plt.subplot(3,4,11)
plt.plot(t,i25hh,label='$T_x = 25^o$C',lw=1)
#plt.axhline(y=0, xmin=0, xmax=10, linewidth=1, color = 'k', linestyle = '--')
#plt.legend()
#plt.xlabel('Time (ms)')
#plt.ylabel('Total Current (mA/cm$^2$)')
#plt.text(56,0.8,'$T_x = 25^o$C')
plt.xlim(50,60)
#plt.fill_between(t,hot25, where= hot25>0, facecolor='red', interpolate=True)
#plt.fill_between(t,hot25, where= hot25<0, facecolor='green', interpolate=True)
plt.ylim(-5,15)
ax.set_xticklabels([])
ax.set_xticks([])
ax.set_yticklabels([])
ax.set_yticks([])

T = 25
ax = plt.subplot(3,4,12)
plt.bar(T-width,areas_na[2], width , color = 'green', label='Net inward current')
plt.bar(T,areas_k[2],  width, color='indianred', label='Net positive current')
ax.set_xticklabels([])
ax.set_xticks([])
plt.ylim(0,50)
ax.set_yticks([0,50])
print (areas_k[2]/areas_na[2])
#plt.text(T - width/2 , 8,str(float(areas_k[2]/areas_na[2]))[:4])
plt.subplots_adjust(wspace = 0.36, hspace = 0.29)
#plt.savefig('fig-3l.png',dpi=600, format = 'png', bbox_inches = 'tight')
plt.show()

#############################################################################################



    


















    






/srv/conda/envs/notebook/lib/python3.7/site-packages/ipykernel_launcher.py:12: DeprecationWarning: object of type <class 'float'> cannot be safely interpreted as an integer.
  if sys.path[0] == '':










    






-0.2711516889706277
1.3670119821438473
10.203735965640533
18.826004836430663










    
























In [8]:



    


pwd



    


















    
Out[8]:








'/home/jovyan/3'

















In [ ]:

Content source: mohitganguly/IRBlock_Vanderbilt

Similar notebooks:

notebook.community | gallery | about