import GuggenWater as gw
import numpy as np
T, P, Feature, Src = np.loadtxt('Features.dat',
dtype={
'names': ('T', 'P', 'F', 'S'),
'formats': (np.float, np.float, 'S9', 'S9')
},
usecols=(0, 1, 2, 3),
unpack=True,
skiprows=1)
Property = np.array([str(s, 'utf-8') for s in Feature])
Source = np.array([str(s, 'utf-8') for s in Src])
model = gw.GuggenheimWater()
for t, p, ppty, src in zip(T, P, Property, Source):
model.addPoint(t, p, ppty, src)
model.plot(savefig=True, skip=['Widom'], style='ggplot', wait=True)
model.writetoFile('All.dat')
model2 = gw.GuggenheimWater(gw.ST2I_mean_field())
T, P, Feature, Src = np.loadtxt('FeaturesST2.dat',
dtype={
'names': ('T', 'P', 'F', 'S'),
'formats': (np.float, np.float, 'S9', 'S9')
},
usecols=(0, 1, 2, 3),
unpack=True,
skiprows=1)