import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
from scipy import interpolate
import math
import sys
# import plot_func

# reload (plot_func)
import os


data1 = pd.read_csv('CO2(pO2).csv')
data2 = pd.read_csv('CO2(DPG).csv')
data3 = pd.read_csv('CO2(pH).csv')
data4 = pd.read_csv('CO2(Temp).csv')
data5 = pd.read_csv('CO2(Hct).csv')

Y1 = data1[data1.keys()[0]]
Y2 = data1[data1.keys()[1]]
Y3 = data1[data1.keys()[2]]
Y4 = data1[data1.keys()[3]]
Y5 = data1[data1.keys()[4]]

Y6 = data2[data2.keys()[0]]
Y7 = data2[data2.keys()[1]]
Y8 = data2[data2.keys()[2]]
Y9 = data2[data2.keys()[3]]
Y10 = data2[data2.keys()[4]]

Y11 = data3[data3.keys()[0]]
Y12 = data3[data3.keys()[1]]
Y13 = data3[data3.keys()[2]]
Y14 = data3[data3.keys()[3]]
Y15 = data3[data3.keys()[4]]

Y16 = data4[data4.keys()[0]]
Y17 = data4[data4.keys()[1]]
Y18 = data4[data4.keys()[2]]
Y19 = data4[data4.keys()[3]]
Y20 = data4[data4.keys()[4]]

Y21 = data5[data5.keys()[0]]
Y22 = data5[data5.keys()[1]]
Y23 = data5[data5.keys()[2]]
Y24 = data5[data5.keys()[3]]
Y25 = data5[data5.keys()[4]]
# print(Y1)

PCO2 = list(range(1,102))
# print(PO2)


# Y_values = []
data = [data1, data2, data3, data4, data5]
plt.figure(figsize=(14,10))
plt.subplot(2,2,1)

plt.plot(PCO2, Y1, 'navy' ,linestyle='-',  label = 'pH=6.92',linewidth=3)
plt.plot(PCO2, Y2, 'green' ,linestyle='-',  label = 'pH=7.08',linewidth=3)
plt.plot(PCO2, Y3, 'red' ,linestyle='-',  label = 'pH=7.24',linewidth=3)
plt.plot(PCO2, Y4, 'black' ,linestyle='-',  label = 'pH=7.40',linewidth=3)
plt.plot(PCO2, Y5, 'purple' ,linestyle='-',  label = 'pH=7.56',linewidth=3)
plt.plot(PCO2, Y6, 'navy' ,linestyle='-',  label = 'DPG$_{rbc}$=0.003',linewidth=3)
plt.plot(PCO2, Y7, 'green' ,linestyle='-',  label = 'DPG$_{rbc}$=0.0045',linewidth=3)
plt.plot(PCO2, Y8, 'red' ,linestyle='-',  label = 'DPG$_{rbc}$=0.006',linewidth=3)
plt.plot(PCO2, Y9, 'black' ,linestyle='-',  label = 'DPG$_{rbc}$=0.0075',linewidth=3)
plt.plot(PCO2, Y10, 'purple' ,linestyle='-',  label = 'DPG$_{rbc}$=0.009',linewidth=3)

plt.xlim(0, 100)
plt.ylim(0.00, 150)
# plt.yticks(np.arange(0,0.8,0.15))
plt.ylabel ('CO$_2$',fontsize=12)
plt.xlabel ('P$_{CO_2}$(mmHg)',fontsize=12)
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.title('A')
plt.legend(loc='upper left')

#######

plt.subplot(2,2,2)

plt.plot(PCO2, Y11, 'navy' ,linestyle='-',  label = 'pH=6.92',linewidth=3)
plt.plot(PCO2, Y12, 'green' ,linestyle='-',  label = 'pH=7.08',linewidth=3)
plt.plot(PCO2, Y13, 'red' ,linestyle='-',  label = 'pH=7.24',linewidth=3)
plt.plot(PCO2, Y14, 'black' ,linestyle='-',  label = 'pH=7.40',linewidth=3)
plt.plot(PCO2, Y15, 'purple' ,linestyle='-',  label = 'pH=7.56',linewidth=3)

plt.xlim(0, 100)
plt.ylim(0.00, 250)
# plt.yticks(np.arange(0,0.8,0.15))
plt.ylabel ('CO$_2$',fontsize=12)
plt.xlabel ('P$_{CO_2}$(mmHg)',fontsize=12)
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.title('B')
plt.legend(loc='upper left')

#######
plt.subplot(2,2,3)

plt.plot(PCO2, Y16, 'navy' ,linestyle='-',  label = 'T=27',linewidth=3)
plt.plot(PCO2, Y17, 'green' ,linestyle='-',  label = 'T=32',linewidth=3)
plt.plot(PCO2, Y18, 'red' ,linestyle='-',  label = 'T=37',linewidth=3)
plt.plot(PCO2, Y19, 'black' ,linestyle='-',  label = 'T=42',linewidth=3)
plt.plot(PCO2, Y20, 'purple' ,linestyle='-',  label = 'T=47',linewidth=3)

plt.xlim(0, 100)
plt.ylim(0.00, 150)
# plt.yticks(np.arange(0,0.8,0.15))
plt.ylabel ('CO$_2$',fontsize=12)
plt.xlabel ('P$_{CO_2}$(mmHg)',fontsize=12)
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.title('C')
plt.legend(loc='upper left')
#######
plt.subplot(2,2,4)

plt.plot(PCO2, Y21, 'navy' ,linestyle='-',  label = 'Hct=0.1',linewidth=3)
plt.plot(PCO2, Y22, 'green' ,linestyle='-',  label = 'Hct=0.2',linewidth=3)
plt.plot(PCO2, Y23, 'red' ,linestyle='-',  label = 'Hct=0.3',linewidth=3)
plt.plot(PCO2, Y24, 'black' ,linestyle='-',  label = 'Hct=0.4',linewidth=3)
plt.plot(PCO2, Y25, 'purple' ,linestyle='-',  label = 'Hct=0.5',linewidth=3)

plt.xlim(0, 100)
plt.ylim(0.00, 150)
# plt.yticks(np.arange(0,0.8,0.15))
plt.ylabel ('CO$_2$',fontsize=12)
plt.xlabel ('P$_{CO_2}$(mmHg)',fontsize=12)
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.title('D')
plt.legend(loc='upper left')
#######

plt.savefig('Figure06.png')
plt.show()