import numpy as np
import pandas as pd
import sys
sys.path.append("..") # This allows you to put the python file in the folder above. So you can have lots of sub folders
import matplotlib.pyplot as plt
import Thermobar as pt
out=pt.import_excel('Waters2021.xlsx', sheet_name="Sheet1")
my_input=out['my_input']
myLiquids1=out['Liqs']
myCpxs1=out['Cpxs']
import os
import sys
os.path.dirname(sys.executable)
'C:\\Users\\penny\\anaconda3'
T34=pt.calculate_liq_only_temp(liq_comps=myLiquids1, equationT="T_Put2008_eq34_cpx_sat", P=10*my_input['P'])-273.15
plt.plot(T34, my_input['temp_out'], 'ok')
[<matplotlib.lines.Line2D at 0x1f69edc0dc0>]
T33=pt.calculate_cpx_liq_temp(liq_comps=myLiquids1, cpx_comps=myCpxs1, equationT="T_Put2008_eq33", P=my_input['P'])-273.15
plt.plot(T, my_input['temp_out'], 'ok')
C:\Users\penny\anaconda3\lib\site-packages\pandas\core\arraylike.py:358: RuntimeWarning: divide by zero encountered in log result = getattr(ufunc, method)(*inputs, **kwargs)
[<matplotlib.lines.Line2D at 0x1f69ed030d0>]
plt.plot(T33, T34, 'ok')
plt.plot([950, 1000], [950, 1000], '-r')
plt.xlabel('equation 33 (cpx-liq temp)')
plt.ylabel('equation 34 (liq-only sat temp)')
Text(0, 0.5, 'equation 34 (liq-only sat temp)')