def realize_pyrolite():
# approximate four component pyrolite model
x_pv = 0.67
x_fp = 0.33
pv_fe_num = 0.07
fp_fe_num = 0.2
mg_perovskite = minerals.SLB_2011_ZSB_2013.mg_perovskite()
realize_mineral(mg_perovskite)
fe_perovskite = minerals.SLB_2011_ZSB_2013.fe_perovskite()
realize_mineral(fe_perovskite)
wuestite = minerals.SLB_2011_ZSB_2013.wuestite()
realize_mineral(wuestite)
periclase = minerals.SLB_2011_ZSB_2013.periclase()
realize_mineral(periclase)
perovskite = HelperSolidSolution(
[mg_perovskite, fe_perovskite], [1.0 - pv_fe_num, pv_fe_num])
ferropericlase = HelperSolidSolution(
[periclase, wuestite], [1.0 - fp_fe_num, fp_fe_num])
pyrolite = burnman.Composite([perovskite, ferropericlase], [x_pv, x_fp])
pyrolite.set_method('slb3')
anchor_temperature = normal(loc=1935.0, scale=200.0)
return pyrolite, anchor_temperature
def make_rock():
# approximate four component pyrolite model
x_pv = 0.67
x_fp = 0.33
pv_fe_num = 0.07
fp_fe_num = 0.2
mg_perovskite = minerals.SLB_2011_ZSB_2013.mg_perovskite()
fe_perovskite = minerals.SLB_2011_ZSB_2013.fe_perovskite()
wuestite = minerals.SLB_2011_ZSB_2013.wuestite()
periclase = minerals.SLB_2011_ZSB_2013.periclase()
perovskite = HelperSolidSolution([mg_perovskite, fe_perovskite],
[1.0 - pv_fe_num, pv_fe_num])
ferropericlase = HelperSolidSolution([periclase, wuestite],
[1.0 - fp_fe_num, fp_fe_num])
pyrolite = burnman.Composite([perovskite, ferropericlase], [x_pv, x_fp])
pyrolite.set_method('slb3')
anchor_temperature = 1935.0
return pyrolite, anchor_temperature
def __init__(self, fe_num):
endmembers = [mg_perovskite, fe_perovskite]
molar_fractions = [1. - fe_num, 0.0 + fe_num]
HelperSolidSolution.__init__(self, endmembers, molar_fractions)
def __init__(self, fe_num):
endmembers = [periclase, wustite]
molar_fractions = [1. - fe_num, 0.0 + fe_num]
HelperSolidSolution.__init__(self, endmembers, molar_fractions)