In [1]:
%cat 0Source_Citation.txt
In [2]:
%matplotlib inline
# %matplotlib notebook # for interactive
For high dpi displays.
In [3]:
%config InlineBackend.figure_format = 'retina'
This example compares pressure calculated from pytheos
and original publication for the gold scale by Shim 2001.
In [4]:
import matplotlib.pyplot as plt
import numpy as np
from uncertainties import unumpy as unp
import pytheos as eos
In [5]:
eta = np.linspace(0., 0.34, 18)
print(eta)
In [6]:
shim_au = eos.gold.Shim2002()
In [7]:
shim_au.print_equations()
In [8]:
shim_au.print_equations()
In [9]:
shim_au.print_parameters()
In [10]:
v0 = 67.84742110765599
In [11]:
shim_au.three_r
Out[11]:
In [12]:
v = v0 * (1.-eta)
temp = 3000.
In [13]:
p = shim_au.cal_p(v, temp * np.ones_like(v))
In [14]:
print('for T = ', temp)
for eta_i, p_i in zip(eta, p):
print("{0: .3f} {1: .2f}".format(eta_i, p_i))
In [15]:
v = shim_au.cal_v(p, temp * np.ones_like(p), min_strain=0.6)
print(1.-(v/v0))
In [ ]: