print('{0}: {1:.0e}% at {2:.0f} GPa and {3:.0f} K'.format(variables[j], percentage_diff[i, j], PT[i][0], PT[i][1]))
print('')
print('Check excess entropy and landau in SLB2011')
file = '../../burnman/data/input_perplex/test_SLB_entropy_and_landau.dat'
f = open(file, 'r')
mins = {}
mins['Wus'] = SLB_2011.wuestite()
mins['Per'] = SLB_2011.periclase()
mins['Wad'] = SLB_2011.mg_wadsleyite()
mins['Ring'] = SLB_2011.mg_ringwoodite()
mins['Stv'] = SLB_2011.stishovite()
variables = ['H', 'S', 'V', 'C_p', 'alpha', 'K_T', 'rho']
percentage_diff = []
mineral_names = []
print(variables)
datalines = [line.strip().split()
for idx, line in enumerate(f.read().split('\n')) if line.strip() and idx > 0]
for line in datalines:
m = mins[line[0]]
mineral_names.append(m.name)
data = map(float, line[1:])
P, T, M, H, S, V, C_p, alpha, beta, C_p_over_C_v, rho = data
PT.append([P, T])
P = P * 1.e9
run_gt_solvus = True
run_fper_ol = True
run_fixed_ol_composition = True
run_upper_mantle = True
run_lower_mantle = True
run_olivine_polymorphs = True
gt = SLB_2011.garnet()
ol = SLB_2011.mg_fe_olivine()
wad = SLB_2011.mg_fe_wadsleyite()
rw = SLB_2011.mg_fe_ringwoodite()
bdg = SLB_2011.mg_fe_bridgmanite()
ppv = SLB_2011.post_perovskite()
fper = SLB_2011.ferropericlase()
opx = SLB_2011.orthopyroxene()
stv = SLB_2011.stishovite()
coe = SLB_2011.coesite()
cpv = SLB_2011.ca_perovskite()
if __name__ == "__main__" and run_aluminosilicates:
"""
Creates the classic aluminosilicate diagram involving
univariate reactions between andalusite, sillimanite and kyanite.
"""
sillimanite = HP_2011_ds62.sill()
andalusite = HP_2011_ds62.andalusite()
kyanite = HP_2011_ds62.ky()
# First, find the pressure and temperature of the invariant point
composition = sillimanite.formula
assemblage = burnman.Composite([sillimanite, andalusite, kyanite])
print('{0}: {1:.0e}% at {2:.0f} GPa and {3:.0f} K'.format(
variables[j], percentage_diff[i, j], PT[i][0], PT[i][1]))
print('')
print('Check excess entropy and landau in SLB2011')
file = '../../burnman/data/input_perplex/test_SLB_entropy_and_landau.dat'
f = open(file, 'r')
mins = {}
mins['Wus'] = SLB_2011.wuestite()
mins['Per'] = SLB_2011.periclase()
mins['Wad'] = SLB_2011.mg_wadsleyite()
mins['Ring'] = SLB_2011.mg_ringwoodite()
mins['Stv'] = SLB_2011.stishovite()
variables = ['H', 'S', 'V', 'C_p', 'alpha', 'K_T', 'rho']
percentage_diff = []
mineral_names = []
print(variables)
datalines = [
line.strip().split() for idx, line in enumerate(f.read().split('\n'))
if line.strip() and idx > 0
]
for line in datalines:
m = mins[line[0]]
mineral_names.append(m.name)
data = map(float, line[1:])
P, T, M, H, S, V, C_p, alpha, beta, C_p_over_C_v, rho = data