plt.rc('text', usetex=True)
plt.rcParams['text.latex.preamble'] = '\usepackage{relsize}'
plt.rc('font', family='sanserif')
figure=plt.figure(dpi=100,figsize=figsize)
#plot v_s
plt.plot(pressure/1.e9,seis_vs/1.e3,linestyle="-",color='k',linewidth=2.0,label='PREM')
#plot v_phi
plt.plot(pressure/1.e9,seis_vphi/1.e3,linestyle="-",color='k',linewidth=2.0)
#plot density
plt.plot(pressure/1.e9,seis_rho/1.e3,linestyle="-",color='k',linewidth=2.0)
plt.plot(pressure/1.e9,vphi/1.e3,linestyle="-",color=colors.color(3),linewidth=1.0,marker='s',markerfacecolor=colors.color(3), label="vphi")
plt.plot(pressure/1.e9,vs/1.e3,linestyle="-",color=colors.color(4),linewidth=1.0,marker='v',markerfacecolor=colors.color(4), label="vs")
plt.plot(pressure/1.e9,rho/1.e3,linestyle="-",color=colors.color(2),linewidth=1.0,marker='o',markerfacecolor=colors.color(2), label="rho")
rock, anchor_t = array_to_rock(lit, names)
temperature = burnman.geotherm.adiabatic(pressure, anchor_t, rock)
rho, vp, vs, vphi, K, G = \
burnman.velocities_from_rock(rock, pressure, temperature, burnman.averaging_schemes.HashinShtrikmanAverage())
plt.plot(pressure/1.e9,vs/1.e3,dashes=dashstyle2,color=colors.color(4),linewidth=1.0)
plt.plot(pressure/1.e9,vphi/1.e3,dashes=dashstyle2,color=colors.color(3),linewidth=1.0, label="literature")
plt.plot(pressure/1.e9,rho/1.e3,dashes=dashstyle2,color=colors.color(2),linewidth=1.0)
plt.xlabel("Pressure (GPa)")
plt.ylabel("Velocities (km/s) and Density (kg/m$^3$)")
plt.legend(bbox_to_anchor=(1.0, 0.9),prop={'size':12})
seis_vphi / 1.e3,
linestyle="-",
color='k',
linewidth=2.0)
# plot density
plt.plot(pressure / 1.e9,
seis_rho / 1.e3,
linestyle="-",
color='k',
linewidth=2.0)
plt.plot(pressure / 1.e9,
vphi / 1.e3,
linestyle="-",
color=colors.color(3),
linewidth=1.0,
marker='s',
markerfacecolor=colors.color(3),
label="vphi")
plt.plot(pressure / 1.e9,
vs / 1.e3,
linestyle="-",
color=colors.color(4),
linewidth=1.0,
marker='v',
markerfacecolor=colors.color(4),
label="vs")
plt.plot(pressure / 1.e9,
rho / 1.e3,
linestyle="-",
plt.rcParams['text.latex.preamble'] = '\usepackage{relsize}'
plt.rc('font', family='sanserif')
figure=plt.figure(dpi=100,figsize=figsize)
#plot v_s
plt.plot(pressure/1.e9,seis_vs/1.e3,linestyle="-",color='k',linewidth=2.0,label='PREM')
#plot v_phi
plt.plot(pressure/1.e9,seis_vphi/1.e3,linestyle="-",color='k',linewidth=2.0)
#plot density
plt.plot(pressure/1.e9,seis_rho/1.e3,linestyle="-",color='k',linewidth=2.0)
plt.plot(pressure/1.e9,vphi/1.e3,linestyle="-",color=colors.color(3),linewidth=1.0,marker='s',markerfacecolor=colors.color(3), label="vphi")
plt.plot(pressure/1.e9,vs/1.e3,linestyle="-",color=colors.color(4),linewidth=1.0,marker='v',markerfacecolor=colors.color(4), label="vs")
plt.plot(pressure/1.e9,rho/1.e3,linestyle="-",color=colors.color(2),linewidth=1.0,marker='o',markerfacecolor=colors.color(2), label="rho")
rock, anchor_t = array_to_rock(lit, names)
temperature = burnman.geotherm.adiabatic(pressure, anchor_t, rock)
rho, vp, vs, vphi, K, G = \
burnman.velocities_from_rock(rock, pressure, temperature, burnman.averaging_schemes.hashin_shtrikman_average())
plt.plot(pressure/1.e9,vs/1.e3,dashes=dashstyle2,color=colors.color(4),linewidth=1.0)
plt.plot(pressure/1.e9,vphi/1.e3,dashes=dashstyle2,color=colors.color(3),linewidth=1.0, label="literature")
plt.plot(pressure/1.e9,rho/1.e3,dashes=dashstyle2,color=colors.color(2),linewidth=1.0)
plt.xlabel("Pressure (GPa)")
plt.ylabel("Velocities (km/s) and Density (kg/m$^3$)")