- Author:
- nima <nafs080@aucklanduni.ac.nz>
- Date:
- 2021-03-26 09:25:46+13:00
- Desc:
- figure 06 and 07 are updated
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/648/rawfile/506de3694dd9b28f81a1199be8c55d57648493b7/Fig06-plot.py
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('Fig06.csv')
X1 = data[data.keys()[4]]*1000
Y1 = data[data.keys()[0]]
Y2 = data[data.keys()[1]]
Y3 = data[data.keys()[2]]
Y4 = data[data.keys()[3]]
Y5 = data[data.keys()[5]]*1e6
Y6 = data[data.keys()[6]]*1e6
Y7 = data[data.keys()[7]]*1e6
Y8 = data[data.keys()[8]]*1e6
Y9 = data[data.keys()[9]]*1000/57
Y10 = data[data.keys()[10]]*1000/57
Y11 = data[data.keys()[11]]*1000/57
Y12 = data[data.keys()[12]]*1000/57
plt.figure(figsize=(16,16))
plt.subplot(2,2,1)
plt.plot(X1, Y1, 'navy',linestyle='-', label = 'Control', linewidth= 3)
plt.plot(X1, Y2, 'red',linestyle='-', label = 'Block 50%', linewidth= 3)
plt.plot(X1, Y3, 'green',linestyle='-', label = 'Block 75%', linewidth= 3)
plt.plot(X1, Y4, 'purple',linestyle='-', label = 'Block 90%', linewidth= 3)
plt.grid()
plt.xlim(0, 1900)
plt.ylim(-70,30)
plt.xlabel ('time(ms)',fontsize=16)
plt.tick_params(axis='both', labelsize=14)
plt.ylabel ('V$_m$(mV)', fontsize=16)
plt.title('A', fontsize='14')
plt.legend(bbox_to_anchor=(0., 1.02, 1.02, 0.1), loc='best',fontsize=14,
ncol=10, mode="expand")
plt.subplot(2,2,2)
plt.plot(X1, Y1, 'navy',linestyle='-', label = 'Control', linewidth= 3)
plt.plot(X1, Y2, 'red',linestyle='-', label = 'Block 50%', linewidth= 3)
plt.plot(X1, Y3, 'green',linestyle='-', label = 'Block 75%', linewidth= 3)
plt.plot(X1, Y4, 'purple',linestyle='-', label = 'Block 90%', linewidth= 3)
plt.grid()
plt.xlim(400, 1150)
plt.ylim(-70,-35)
plt.xlabel ('time(ms)',fontsize=16)
plt.tick_params(axis='both', labelsize=14)
plt.ylabel ('V$_m$(mV)', fontsize=16)
plt.title('B', fontsize='14')
plt.legend(bbox_to_anchor=(0., 1.02, 1.02, 0.1), loc='best',fontsize=14,
ncol=10, mode="expand")
plt.subplot(2,2,3)
plt.plot(X1, Y5, 'navy',linestyle='-', label = 'Control', linewidth= 3)
plt.plot(X1, Y6, 'red',linestyle='-', label = 'Block 50%', linewidth= 3)
plt.plot(X1, Y7, 'green',linestyle='-', label = 'Block 75%', linewidth= 3)
plt.plot(X1, Y8, 'purple',linestyle='-', label = 'Block 90%', linewidth= 3)
plt.grid()
plt.xlim(0, 1900)
plt.ylim(50,450)
plt.xlabel ('time(ms)',fontsize=16)
plt.tick_params(axis='both', labelsize=14)
plt.ylabel ('[Ca$^{2+}]_i$ (nM)', fontsize=16)
plt.title('A', fontsize='14')
plt.legend(bbox_to_anchor=(0., 1.02, 1.02, 0.1), loc='best',fontsize=14,
ncol=10, mode="expand")
plt.subplot(2,2,4)
plt.plot(X1, Y9, 'navy',linestyle='-', label = 'Control', linewidth= 3)
plt.plot(X1, Y10, 'red',linestyle='-', label = 'Block 50%', linewidth= 3)
plt.plot(X1, Y11, 'green',linestyle='-', label = 'Block 75%', linewidth= 3)
plt.plot(X1, Y12, 'purple',linestyle='-', label = 'Block 90%', linewidth= 3)
plt.grid()
plt.xlim(400, 1150)
plt.ylim(-0.16, 0)
plt.xlabel ('time(ms)',fontsize=16)
plt.tick_params(axis='both', labelsize=14)
plt.ylabel ('I$_{NaCa}$ (pA/pF)', fontsize=16)
plt.title('D', fontsize='14')
plt.legend(bbox_to_anchor=(0., 1.02, 1.02, 0.1), loc='best',fontsize=14,
ncol=10, mode="expand")
plt.tight_layout(pad=0.4, w_pad=0.5, h_pad=1.0)
plt.savefig('C:/Nima/ABI/Physiome Journal/sinus/Python_codes/figure06.png')
plt.show()
#~ plt.legend(loc=0)
