from pylab import *
import sys
sys.path.append('../equilibre') 
from thermochimie import *          
th=Thermochimie("../equilibre/thermodata.txt") 
G=EnthalpieLibre(th,"../equilibre/liste_reactions.txt")
figure(figsize=(10,12))
xlabel('T (K)')  
ylabel(r'$\Delta_rG^0\ (kJ/mol)$')
G.plot(10,'r')
G.plot(11,'b') 
G.plot(8,'g')
G.plot(9,'brown')   
G.plot(4,'black')  
G.plot(14,'black')
G.plot(15,'navy')
G.plot(16,'brown')
G.plot(17,'r')
axis([300,3000,-1200,300])
grid(True)
            
T=300
R=8.31
G0=G.gstandard(10,300)*1e3
            
figure(figsize=(7,7))
xlabel('T (K)') 
ylabel(r'$RT\ ln(P_{O_2}\ )$')  
G.plot(10,'r')
text(1000,-50,r'CuO')
text(1000,-200,r'Cu')
axis([300,2200,-300,0]) 
            
figure(figsize=(10,10))
xlabel('T (K)') 
ylabel(r'$\Delta_rG^0\ \rm(kJ/mol)$')
G.plot(1,'r')
G.plot(2,'b') 
G.plot(3,'g') 
#G.plot(7,'black')
#axis([300,3000,-600,0])
axis([300,3000,-1200,300]) 
grid(True) 
            
[H,S]=th.reaction("C(s)+CO_2(g)=2*CO(g)")
figure(figsize=(10,10))
xlabel('T (K)')  
ylabel(r'$\Delta_rG^0\ \rm (kJ/mol)$') 
G.plot(1,'r')   
G.plot(2,'b')  
G.plot(3,'g')  
G.plot(4,'black')    
axis([300,3000,-600,0])
grid(True) 
            
figure(figsize=(10,10))
xlabel('T (K)')  
ylabel(r'$\Delta_rG^0\ \rm(kJ/mol)$')    
G.plot(2,'b')  
G.plot(11,'black')
G.plot(12,'black')
G.plot(13,'black')
axis([300,3000,-600,0])
grid(True) 
            
[DH,DS]=th.reaction("FeO(s)+CO(g)=Fe(s)+CO_2(g)")
T=1300
DG=DH*1e3-T*DS
            
R=8.31
K=exp(-DG/(R*T))