from constants0 import *
from scipy.integrate import quad

def Cp_n2(T):
    return cpg[2]

def Cp_o2(T):
    return cpg[3]

def Cp_co2g(T):
    return cpg[0]

def Cp_mea(T): #kJ/kg
    return Aa + (Ba * T) + (Ca * (T**2))

def Cp_h2og(T):
    return cpg[1]

def Cp_h2o(T):  #kJ/kg
    return Aw + (Bw * T) + (Cw * (T**2))

def Cp_co2(T):  #kJ/kg for co2 absorbert i MEA
    return Ac + (Bc*T)

def Cp_mea_aq(T):   #kJ/kg
    w_mea = waMEA
    cp_h2o = Cp_h2o(T)
    cp_mea = Cp_mea(T)

    return (1-w_mea)*cp_h2o + (w_mea*cp_mea) + ((w_mea*(1-w_mea))*(As + Bs*T + Cs * w_mea * ((T - 273.15)**(-1.5859))))

def mean_Cp(T1,T2):
    I = quad(Cp_mea_aq,T1,T2)
    return I[0]/(T2-T1)