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 Dorogokupets 2007.
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(1., 0.65, 8)
print(eta)
In [6]:
dorogokupets2007_au = eos.gold.Dorogokupets2007()
In [7]:
help(dorogokupets2007_au)
In [8]:
dorogokupets2007_au.print_equations()
In [9]:
dorogokupets2007_au.print_equations()
In [10]:
dorogokupets2007_au.print_parameters()
In [11]:
v0 = 67.84742110765599
In [12]:
dorogokupets2007_au.three_r
Out[12]:
In [13]:
v = v0 * (eta)
temp = 2500.
In [14]:
p = dorogokupets2007_au.cal_p(v, temp * np.ones_like(v))
In [15]:
print('for T = ', temp)
for eta_i, p_i in zip(eta, p):
print("{0: .3f} {1: .2f} ".format(eta_i, p_i))
In [16]:
v = dorogokupets2007_au.cal_v(p, temp * np.ones_like(p), min_strain=0.6)
print(1.-(v/v0))