Worked examples
This page summarizes the different examples available. If you want a specific example, we can always include more.
Liquids and Olivines
This example shows various options for calculating liquid-only temperatures, as well as sensitivity to temperature and pressure.
This example shows:
how to perform olivine-liquid thermometry,
both for paired olivines and liquids, and showing how to consider all possible pairings. - How to calculate equilibrium olivine-liquid contents for different liquid compositions using a variety of Kd models (Roeder, Matzen, Toplis, Putirka) - How to plot olivines and liquids on Rhodes diagrams - How to calculate Fe3/FeT from a buffer position.
This example shows how to perform olivine-liquid hygrometry calculations.
This example shows how to plot a Rhodes diagram for a single sample with an olivine population.
This example shows how to plot a Rhodes diagram when you have different matrix glass/whole rock samples, and an olivine population from each sample.
This example shows how to perform CHOMPI calculations following Blundy (2022).
Clinopyroxene (Cpx)
This example shows how to perform Cpx-Liq and Cpx-only calculations when you have either just Cpx, or paired Cpx-Liqs. Melt matching is covered below. It also shows how to plot Cpx Liq comps on a Rhodes diagram.
This example shows how to perform Cpx-Liq and Cpx-only calculations when you have using the Petrelli and Jorgenson machine learning models. There are two options, using onnx and using the voting method of Jorgenson. Note, if you use these models, you have to install the .pkl and .onnx files from a zip from github (instructions in notebook) and check your sklearn installation matches the setup.py (should check automatically but may need manual upgrading).
This example shows how to plot a pyroxene ternary diagram for pyroxene data from Kilauea
This example is a simple introduction to melt matching. It shows how to combine all possible liquid and cpx compositions using various equilibrium tests.
This example is a more complex example following the method of Scruggs and Putirka (2018) - if you dont have many liquids (perhaps just mafic and silicic end members) it shows how to make synthetic liquids.
Orthopyroxene (Opx)
This example shows how to plot pyroxene compositions on a ternary diagram
This example shows how to consider matches between all possible liq-opx pairs.
Two Pyroxene (Opx-Cpx)
This example shows how to perform Opx-Cpx calcs on pre-matched pairs.
This example shows how to consider all possible Cpx-Opx pairs using equilibrium tests to calculate P and T.
Amphiboles
This example shows how to perform Amp-only and Amp-Liq thermobarometry and chemometry calculations.
This example shows how to do Amp-Liq melt matching.
This example shows how to plot amphibole classification diagrams.
Feldspars
This example shows how to perform feldspar-liquid thermobarometry.
This example shows how to perform plagioclase-liquid hygrometry calculations.
This :doc:`example <Examples/Feldspar_Thermobarometry/Two_Feldspar_All_Possible_Pairs>’ shows how to perform two feldspar thermobarometry.
This :doc:`example <Examples/Feldspar_Thermobarometry/Fspar_Liq_Matching>’ shows how to perform Kspar-liq and plag-liq melt matching.
This :doc:`example <Examples/Feldspar_Thermobarometry/Fspar_Ternary_Plot>’ shows how to plot a feldspar ternary diagram.
This :doc:`example <Examples/General_Plotting/Plagioclase_Classification_Kilauea>’ shows how to plot a plagioclase ternary diagram segmented by sample.
Garnet
This :doc:`example <Examples/Garnet_Geotherms/Garnet_Functions>’ shows how to perform garnet thermobarometry and plot a garnet compositional section with a geotherm etc.
This :doc:`example <Examples/Garnet_Geotherms/Geotherm_functions>’ shows how to calculate a garnet geotherm
Error Propagation
We do not include examples for all phases, but each example below shows some different data vis strategies, so you should be able to adapt them for your problem of interest (e.g. look at them all, even if you are only interested in Cpx!)
This :doc:`example <Examples/Error_propagation/Liquid_Thermometry_Error_prop>’ shows how to propagate analytical error when performing liquid thermometry. It also shows how to plot error ellipses.
This :doc:`example <Examples/Error_propagation/Cpx_only_contour_plot>’ shows how to propagate analytical error when performing Cpx-Liq thermobarometry.
Melt inclusions
This :doc:`example <Examples/Melt_Inclusion_FeMg_Equilibrium/Ol_Melt_Inclusion_Mg_Fe_Eq_SingleSamples>’ shows how to assess Fe-Mg disequilibrium between melt inclusion glasses and the host olivine for a single sample, this :doc:`example <Examples/Melt_Inclusion_FeMg_Equilibrium/Ol_Melt_Inclusion_Mg_Fe_Eq_MultipleSamples> shows the same but for multiple different samples
This :doc:`example <Examples/Integration_with_VESIcal/Combining_VESIcal_Thermobar_SatPs>’ shows how to integrate Thermobar with VESical to get a saturation pressure for each melt inclusion.
Plotting Mineral classification diagrams.
This example shows how to plot pyroxene compositions on a ternary diagram
This example shows how to plot a pyroxene ternary diagram for pyroxene data from Kilauea
This example shows how to plot amphibole classification diagrams.
This :doc:`example <Examples/Feldspar_Thermobarometry/Fspar_Ternary_Plot>’ shows how to plot a feldspar ternary diagram.
This :doc:`example <Examples/General_Plotting/Plagioclase_Classification_Kilauea>’ shows how to plot a plagioclase ternary diagram segmented by sample.
Other Useful examples
Calculating equilibrium mineral contents from a liquid line of descent
This example shows how to calculate equilibrium mineral contents from a Petrolog3 liquid line of descent (it could be an output from any other model, or experiments)
Assesing the calibration range of models
This example shows how to get and plot the data used to calibrate various models in Thermobar
Calculating viscosity
This example shows how to calculate viscosity using the model of Giordano et al. (2008).
Converting pressures to depths
This example shows how to convert pressures you have calculated into depths using a variety of crustal density models.