%matplotlib inline
# %matplotlib notebook # for interactive
For high dpi displays.
%config InlineBackend.figure_format = 'retina'
This example compares pressure calculated from pytheos
and original publication for the gold scale by Jamieson 1983.
import matplotlib.pyplot as plt
import numpy as np
from uncertainties import unumpy as unp
import pytheos as eos
eta = np.linspace(0., 0.225, 46)
print(eta)
jamieson_aul = eos.gold.Jamieson1982L()
jamieson_auh = eos.gold.Jamieson1982H()
jamieson_aul.print_equations()
jamieson_auh.print_equations()
jamieson_aul.print_parameters()
jamieson_auh.print_parameters()
v0 = 67.84747902176544
jamieson_aul.three_r
jamieson_auh.three_r
v = v0 * (1.-eta)
temp = 1500.
p = jamieson_auh.cal_p(v, temp * np.ones_like(v))
print('for T = ', temp)
for eta_i, p_i in zip(eta, p):
print("{0: .3f} {1: .2f}".format(eta_i, p_i))
v = jamieson_auh.cal_v(p, temp * np.ones_like(p), min_strain=0.6)
print(1.-(v/v0))