{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Two Feldspar thermometry\n", "- First, we consider analyses of touching Kspar-Plag pairs\n", "- Then we consider the situation common in lavas where you have a lot of individual Plag-Kspar crystals, and you want to assess all possible combinations from a single sample or deposit\n", "- You can download the excel spreadsheet this reads from here: https://github.com/PennyWieser/Thermobar/blob/main/docs/Examples/Feldspar_Thermobarometry/Two_Feldspar_input.xlsx" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### You need to install Thermobar once on your machine, if you haven't done this yet, uncomment the line below (remove the #)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "#!pip install Thermobar" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Import python things and set plotting parameters" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import Thermobar as pt\n", "pd.options.display.max_columns = None\n", "\n", "# This sets some plotting things\n", "plt.rcParams[\"font.family\"] = 'arial'\n", "plt.rcParams[\"font.size\"] =12\n", "plt.rcParams[\"mathtext.default\"] = \"regular\"\n", "plt.rcParams[\"mathtext.fontset\"] = \"dejavusans\"\n", "plt.rcParams['patch.linewidth'] = 1\n", "plt.rcParams['axes.linewidth'] = 1 \n", "plt.rcParams[\"xtick.direction\"] = \"in\"\n", "plt.rcParams[\"ytick.direction\"] = \"in\"\n", "plt.rcParams[\"ytick.direction\"] = \"in\"\n", "plt.rcParams[\"xtick.major.size\"] = 6 # Sets length of ticks\n", "plt.rcParams[\"ytick.major.size\"] = 4 # Sets length of ticks\n", "plt.rcParams[\"ytick.labelsize\"] = 12 # Sets size of numbers on tick marks\n", "plt.rcParams[\"xtick.labelsize\"] = 12 # Sets size of numbers on tick marks\n", "plt.rcParams[\"axes.titlesize\"] = 14 # Overall title\n", "plt.rcParams[\"axes.labelsize\"] = 14 # Axes labels" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Example 1 - Touching Pairs pressure calculations\n", "- Make an excel spreadsheet where you have rows for each touching pair, with the column headings SiO2_Kspar for SiO2 in the Kspar, SiO2_Plag for SiO2 in the Plag etc. The column order doesnt matter. \n", "- FeO is always as FeOt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Loading in excel spreadsheet" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "out=pt.import_excel('Two_Feldspar_input.xlsx', sheet_name=\"Paired_Plag_Kspar\")\n", "# This extracts a dataframe of all inputs\n", "my_input=out['my_input']\n", "# This extracts a dataframe of kspar compositions from the dictionary \"out\"\n", "Kspars=out['Kspars']\n", "# This extracts a dataframe of plag compositions from the dictionary \"out\"\n", "Plags=out['Plags']" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### check that your inputs have been read in properly, to do this, do \"name\".head() to look at the first few rows" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SiO2_PlagTiO2_PlagAl2O3_PlagFeOt_PlagMnO_PlagMgO_PlagCaO_PlagNa2O_PlagK2O_PlagCr2O3_PlagSample_ID_Plag
056.30.031.40.310.00.0514.602.700.710.00
160.50.024.80.070.00.005.927.361.420.01
257.80.026.70.080.00.008.066.081.130.02
346.20.033.50.310.00.0017.401.220.300.03
450.30.031.10.540.00.0614.502.760.790.04
\n", "
" ], "text/plain": [ " SiO2_Plag TiO2_Plag Al2O3_Plag FeOt_Plag MnO_Plag MgO_Plag CaO_Plag \\\n", "0 56.3 0.0 31.4 0.31 0.0 0.05 14.60 \n", "1 60.5 0.0 24.8 0.07 0.0 0.00 5.92 \n", "2 57.8 0.0 26.7 0.08 0.0 0.00 8.06 \n", "3 46.2 0.0 33.5 0.31 0.0 0.00 17.40 \n", "4 50.3 0.0 31.1 0.54 0.0 0.06 14.50 \n", "\n", " Na2O_Plag K2O_Plag Cr2O3_Plag Sample_ID_Plag \n", "0 2.70 0.71 0.0 0 \n", "1 7.36 1.42 0.0 1 \n", "2 6.08 1.13 0.0 2 \n", "3 1.22 0.30 0.0 3 \n", "4 2.76 0.79 0.0 4 " ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SiO2_KsparTiO2_KsparAl2O3_KsparFeOt_KsparMnO_KsparMgO_KsparCaO_KsparNa2O_KsparK2O_KsparCr2O3_KsparSample_ID_Kspar
065.00.019.10.110.00.030.512.0213.400.00
165.50.019.60.070.00.000.754.819.360.01
265.40.019.40.050.00.000.593.1311.500.02
364.60.018.80.090.00.000.391.1514.800.03
461.80.019.20.510.00.030.661.7112.900.04
\n", "
" ], "text/plain": [ " SiO2_Kspar TiO2_Kspar Al2O3_Kspar FeOt_Kspar MnO_Kspar MgO_Kspar \\\n", "0 65.0 0.0 19.1 0.11 0.0 0.03 \n", "1 65.5 0.0 19.6 0.07 0.0 0.00 \n", "2 65.4 0.0 19.4 0.05 0.0 0.00 \n", "3 64.6 0.0 18.8 0.09 0.0 0.00 \n", "4 61.8 0.0 19.2 0.51 0.0 0.03 \n", "\n", " CaO_Kspar Na2O_Kspar K2O_Kspar Cr2O3_Kspar Sample_ID_Kspar \n", "0 0.51 2.02 13.40 0.0 0 \n", "1 0.75 4.81 9.36 0.0 1 \n", "2 0.59 3.13 11.50 0.0 2 \n", "3 0.39 1.15 14.80 0.0 3 \n", "4 0.66 1.71 12.90 0.0 4 " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "display(Plags.head())\n", "display(Kspars.head())" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- If at any point you want to find out more about a function, and what equations it offers, you can use the help function" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Help on function calculate_plag_kspar_temp in module Thermobar.feldspar:\n", "\n", "calculate_plag_kspar_temp(*, plag_comps=None, kspar_comps=None, Two_Fspar_Match=None, equationT=None, P=None, eq_tests=False)\n", " Two feldspar thermometer (Kspar-Plag), returns temperature in Kelvin\n", " \n", " Parameters\n", " ----------\n", " \n", " plag_comps: pandas.DataFrame\n", " Plag compositions with column headings SiO2_Plag, MgO_Plag etc.\n", " \n", " kspar_comps: pandas.DataFrame\n", " Kspar compositions with column headings SiO2_Kspar, MgO_Kspar etc.\n", " \n", " EquationT: str\n", " choose from:\n", " \n", " | T_Put2008_eq27a (P-dependent, H2O-independent)\n", " | T_Put2008_eq27b (P-dependent, H2O-independent)\n", " | T_Put_Global_2Fspar (P-dependent, H2O-independent)\n", " \n", " P: float, int, pandas.Series, str\n", " Pressure in kbar to perform calculations at.\n", " Only needed for P-sensitive thermometers.\n", " If P=\"Solve\", returns a partial function,\n", " else, enter an integer, float, or panda series.\n", " \n", " Returns\n", " -------\n", " If eq_tests is False\n", " pandas.Series: Temperature in Kelvin\n", " If eq_tests is True\n", " pandas.DataFrame: Temperature in Kelvin+eq Tests + input compositions\n", "\n" ] } ], "source": [ "help(pt.calculate_plag_kspar_temp)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "- You can also do help on the equations themselves to get a bit more information, \n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Help on function T_Put2008_eq27a in module Thermobar.feldspar:\n", "\n", "T_Put2008_eq27a(P, *, K_Barth, Si_Kspar_cat_frac, Ca_Kspar_cat_frac, An_Kspar, An_Plag, Ab_Plag)\n", " Two feldspar thermometer: Equation 27a of Putirka (2008).\n", " :cite:`putirka2008thermometers`\n", " \n", " SEE±23°C for calibration\n", " \n", " SEE±44°C for test data\n", "\n" ] } ], "source": [ "help(pt.T_Put2008_eq27a)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 1a - Calculating T, fixed P\n", "- The function calculate_plag_kspar_temp calculates T in Kelvin using inputted Kspar, Plag compositions, a user-selected equation, and a pressure estimate (5 kbar here). " ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0 944.011837\n", "1 912.626406\n", "2 912.835560\n", "3 902.422454\n", "4 908.475430\n", "5 799.429554\n", "6 818.616248\n", "7 796.533005\n", "8 728.113848\n", "9 727.752376\n", "10 703.051858\n", "11 709.813651\n", "12 730.346965\n", "13 754.219939\n", "14 692.524231\n", "15 647.728144\n", "dtype: float64" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "T_27a_5kbar=pt.calculate_plag_kspar_temp(plag_comps=Plags, kspar_comps=Kspars, \n", " equationT=\"T_Put2008_eq27a\", P=5)-273.15\n", "T_27a_5kbar" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 1b - Can also specify a panda series of pressures for the original spreadsheet\n", "- here, presure stored in a column in GPa, so need to multiply by 10 to get kbar for the function" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "0 989.600535\n", "1 955.165833\n", "2 1053.643565\n", "3 943.994277\n", "4 950.619830\n", "5 831.883113\n", "6 930.082357\n", "7 988.688652\n", "8 969.169184\n", "9 754.549692\n", "10 728.029781\n", "11 735.283052\n", "12 757.339241\n", "13 783.040474\n", "14 716.746820\n", "15 668.870563\n", "dtype: float64" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "T_27a_5kbar=pt.calculate_plag_kspar_temp(plag_comps=Plags, kspar_comps=Kspars, \n", " equationT=\"T_Put2008_eq27a\", \n", " P=10*my_input['P_GPa'])-273.15\n", "T_27a_5kbar" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### 1c - Can also specify Eq_Test=True, this calculates activities based on Elkins and Grove (1990), following the protocol in the supplementary spreadsheets of Putirka. \n", "- Note, Putirka (2008) don't give any values of equilibrium tests unlike Cpx-Opx, or Cpx-Liq etc. Instead, they say \n", "\"As a test for equilibrium, one can compare the activities between those calculated for alkali and plagioclase feldspar, \n", "using the models of Elkins and Grove (1990) (columns AX – AZ). These should nominally be zero, but tests are required to \n", "determine whether this measure of equilibrium is effective at reducing error for T calculated using experimental data.\"\n", "- Thus, we don't apply any filters by default, but show users how to filter based on these values below" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
T_K_calcDelta_AnDelta_AbDelta_Ora_Ab_plga_An_plga_Or_plga_Ab_kspara_An_kspara_Or_ksparSiO2_PlagTiO2_PlagAl2O3_PlagFeOt_PlagMnO_PlagMgO_PlagCaO_PlagNa2O_PlagK2O_PlagCr2O3_PlagSample_ID_PlagSi_Plag_cat_propMg_Plag_cat_propFet_Plag_cat_propCa_Plag_cat_propAl_Plag_cat_propNa_Plag_cat_propK_Plag_cat_propMn_Plag_cat_propTi_Plag_cat_propCr_Plag_cat_propsumSi_Plag_cat_fracMg_Plag_cat_fracFet_Plag_cat_fracCa_Plag_cat_fracAl_Plag_cat_fracNa_Plag_cat_fracK_Plag_cat_fracMn_Plag_cat_fracTi_Plag_cat_fracCr_Plag_cat_fracAn_PlagAb_PlagOr_PlagSiO2_KsparTiO2_KsparAl2O3_KsparFeOt_KsparMnO_KsparMgO_KsparCaO_KsparNa2O_KsparK2O_KsparCr2O3_KsparSample_ID_KsparSi_Kspar_cat_propMg_Kspar_cat_propFet_Kspar_cat_propCa_Kspar_cat_propAl_Kspar_cat_propNa_Kspar_cat_propK_Kspar_cat_propMn_Kspar_cat_propTi_Kspar_cat_propCr_Kspar_cat_propsumSi_Kspar_cat_fracMg_Kspar_cat_fracFet_Kspar_cat_fracCa_Kspar_cat_fracAl_Kspar_cat_fracNa_Kspar_cat_fracK_Kspar_cat_fracMn_Kspar_cat_fracTi_Kspar_cat_fracCr_Kspar_cat_fracAn_KsparAb_KsparOr_KsparK_BarthTP
01262.7505350.0459870.1878720.2089130.5821540.9195221.0240460.3942820.8735350.81513356.30.031.40.310.00.0514.602.700.710.000.9370170.0012410.0043150.2603540.6159220.0871260.0150750.00.00.01.9210500.4877630.0006460.0022460.1355270.3206170.0453540.0078470.00.00.00.7181090.2403120.04158065.00.019.10.110.00.030.512.0213.400.001.0818130.0007440.0015310.0090950.3746530.0651830.2845130.00.00.01.8175330.5952100.0004100.0008420.0050040.2061330.0358640.1565380.00.00.00.0253480.1816750.7929770.7559981262.75053510
11228.3158330.248834-0.042074-0.0754900.6538090.8395070.6026590.6958830.5906720.67814960.50.024.80.070.00.005.927.361.420.011.0069190.0000000.0009740.1055680.4864610.2375000.0301500.00.00.01.8675720.5391590.0000000.0005220.0565270.2604780.1271710.0161440.00.00.00.2828590.6363570.08078465.50.019.60.070.00.000.754.819.360.011.0901350.0000000.0009740.0133740.3844610.1552140.1987350.00.00.01.8428930.5915350.0000000.0005290.0072570.2086180.0842230.1078380.00.00.00.0364100.4225550.5410350.6640221228.31583310
21326.7935650.412296-0.015576-0.2335320.6102830.8615560.5160500.6258590.4492610.74958257.80.026.70.080.00.008.066.081.130.020.9619820.0000000.0011130.1437300.5237300.1961960.0239930.00.00.01.8507430.5197810.0000000.0006020.0776610.2829830.1060090.0129640.00.00.00.3949510.5391210.06592865.40.019.40.050.00.000.593.1311.500.021.0884710.0000000.0006960.0105210.3805380.1010020.2441720.00.00.01.8253990.5962920.0000000.0003810.0057640.2084680.0553320.1337630.00.00.00.0295790.2839570.6864640.5267041326.79356520
31217.144277-0.1250970.3441290.3797720.6006150.9255941.2655070.2564871.0506910.88573546.20.033.50.310.00.0017.401.220.300.030.7689200.0000000.0043150.3102850.6571140.0393680.0063700.00.00.01.7863720.4304370.0000000.0024150.1736960.3678480.0220380.0035660.00.00.00.8715310.1105780.01789164.60.018.80.090.00.000.391.1514.800.031.0751560.0000000.0012530.0069550.3687680.0371090.3142380.00.00.01.8034800.5961560.0000000.0006950.0038560.2044760.0205770.1742400.00.00.00.0194100.1035700.8770200.9366271217.14427710
41223.769830-0.4429380.2232170.4456150.5890900.9322911.2751400.3658731.3752290.82952550.30.031.10.540.00.0614.502.760.790.040.8371570.0014890.0075160.2585710.6100370.0890630.0167740.00.00.01.8206060.4598230.0008180.0041280.1420250.3350740.0489190.0092130.00.00.00.7095660.2444040.04603061.80.019.20.510.00.030.661.7112.900.041.0285550.0007440.0070980.0117690.3766150.0551800.2738970.00.00.01.7538590.5864530.0004240.0040470.0067110.2147350.0314620.1561680.00.00.00.0345300.1618910.8035790.6623921223.76983010
\n", "
" ], "text/plain": [ " T_K_calc Delta_An Delta_Ab Delta_Or a_Ab_plg a_An_plg a_Or_plg \\\n", "0 1262.750535 0.045987 0.187872 0.208913 0.582154 0.919522 1.024046 \n", "1 1228.315833 0.248834 -0.042074 -0.075490 0.653809 0.839507 0.602659 \n", "2 1326.793565 0.412296 -0.015576 -0.233532 0.610283 0.861556 0.516050 \n", "3 1217.144277 -0.125097 0.344129 0.379772 0.600615 0.925594 1.265507 \n", "4 1223.769830 -0.442938 0.223217 0.445615 0.589090 0.932291 1.275140 \n", "\n", " a_Ab_kspar a_An_kspar a_Or_kspar SiO2_Plag TiO2_Plag Al2O3_Plag \\\n", "0 0.394282 0.873535 0.815133 56.3 0.0 31.4 \n", "1 0.695883 0.590672 0.678149 60.5 0.0 24.8 \n", "2 0.625859 0.449261 0.749582 57.8 0.0 26.7 \n", "3 0.256487 1.050691 0.885735 46.2 0.0 33.5 \n", "4 0.365873 1.375229 0.829525 50.3 0.0 31.1 \n", "\n", " FeOt_Plag MnO_Plag MgO_Plag CaO_Plag Na2O_Plag K2O_Plag Cr2O3_Plag \\\n", "0 0.31 0.0 0.05 14.60 2.70 0.71 0.0 \n", "1 0.07 0.0 0.00 5.92 7.36 1.42 0.0 \n", "2 0.08 0.0 0.00 8.06 6.08 1.13 0.0 \n", "3 0.31 0.0 0.00 17.40 1.22 0.30 0.0 \n", "4 0.54 0.0 0.06 14.50 2.76 0.79 0.0 \n", "\n", " Sample_ID_Plag Si_Plag_cat_prop Mg_Plag_cat_prop Fet_Plag_cat_prop \\\n", "0 0 0.937017 0.001241 0.004315 \n", "1 1 1.006919 0.000000 0.000974 \n", "2 2 0.961982 0.000000 0.001113 \n", "3 3 0.768920 0.000000 0.004315 \n", "4 4 0.837157 0.001489 0.007516 \n", "\n", " Ca_Plag_cat_prop Al_Plag_cat_prop Na_Plag_cat_prop K_Plag_cat_prop \\\n", "0 0.260354 0.615922 0.087126 0.015075 \n", "1 0.105568 0.486461 0.237500 0.030150 \n", "2 0.143730 0.523730 0.196196 0.023993 \n", "3 0.310285 0.657114 0.039368 0.006370 \n", "4 0.258571 0.610037 0.089063 0.016774 \n", "\n", " Mn_Plag_cat_prop Ti_Plag_cat_prop Cr_Plag_cat_prop sum \\\n", "0 0.0 0.0 0.0 1.921050 \n", "1 0.0 0.0 0.0 1.867572 \n", "2 0.0 0.0 0.0 1.850743 \n", "3 0.0 0.0 0.0 1.786372 \n", "4 0.0 0.0 0.0 1.820606 \n", "\n", " Si_Plag_cat_frac Mg_Plag_cat_frac Fet_Plag_cat_frac Ca_Plag_cat_frac \\\n", "0 0.487763 0.000646 0.002246 0.135527 \n", "1 0.539159 0.000000 0.000522 0.056527 \n", "2 0.519781 0.000000 0.000602 0.077661 \n", "3 0.430437 0.000000 0.002415 0.173696 \n", "4 0.459823 0.000818 0.004128 0.142025 \n", "\n", " Al_Plag_cat_frac Na_Plag_cat_frac K_Plag_cat_frac Mn_Plag_cat_frac \\\n", "0 0.320617 0.045354 0.007847 0.0 \n", "1 0.260478 0.127171 0.016144 0.0 \n", "2 0.282983 0.106009 0.012964 0.0 \n", "3 0.367848 0.022038 0.003566 0.0 \n", "4 0.335074 0.048919 0.009213 0.0 \n", "\n", " Ti_Plag_cat_frac Cr_Plag_cat_frac An_Plag Ab_Plag Or_Plag \\\n", "0 0.0 0.0 0.718109 0.240312 0.041580 \n", "1 0.0 0.0 0.282859 0.636357 0.080784 \n", "2 0.0 0.0 0.394951 0.539121 0.065928 \n", "3 0.0 0.0 0.871531 0.110578 0.017891 \n", "4 0.0 0.0 0.709566 0.244404 0.046030 \n", "\n", " SiO2_Kspar TiO2_Kspar Al2O3_Kspar FeOt_Kspar MnO_Kspar MgO_Kspar \\\n", "0 65.0 0.0 19.1 0.11 0.0 0.03 \n", "1 65.5 0.0 19.6 0.07 0.0 0.00 \n", "2 65.4 0.0 19.4 0.05 0.0 0.00 \n", "3 64.6 0.0 18.8 0.09 0.0 0.00 \n", "4 61.8 0.0 19.2 0.51 0.0 0.03 \n", "\n", " CaO_Kspar Na2O_Kspar K2O_Kspar Cr2O3_Kspar Sample_ID_Kspar \\\n", "0 0.51 2.02 13.40 0.0 0 \n", "1 0.75 4.81 9.36 0.0 1 \n", "2 0.59 3.13 11.50 0.0 2 \n", "3 0.39 1.15 14.80 0.0 3 \n", "4 0.66 1.71 12.90 0.0 4 \n", "\n", " Si_Kspar_cat_prop Mg_Kspar_cat_prop Fet_Kspar_cat_prop \\\n", "0 1.081813 0.000744 0.001531 \n", "1 1.090135 0.000000 0.000974 \n", "2 1.088471 0.000000 0.000696 \n", "3 1.075156 0.000000 0.001253 \n", "4 1.028555 0.000744 0.007098 \n", "\n", " Ca_Kspar_cat_prop Al_Kspar_cat_prop Na_Kspar_cat_prop K_Kspar_cat_prop \\\n", "0 0.009095 0.374653 0.065183 0.284513 \n", "1 0.013374 0.384461 0.155214 0.198735 \n", "2 0.010521 0.380538 0.101002 0.244172 \n", "3 0.006955 0.368768 0.037109 0.314238 \n", "4 0.011769 0.376615 0.055180 0.273897 \n", "\n", " Mn_Kspar_cat_prop Ti_Kspar_cat_prop Cr_Kspar_cat_prop sum \\\n", "0 0.0 0.0 0.0 1.817533 \n", "1 0.0 0.0 0.0 1.842893 \n", "2 0.0 0.0 0.0 1.825399 \n", "3 0.0 0.0 0.0 1.803480 \n", "4 0.0 0.0 0.0 1.753859 \n", "\n", " Si_Kspar_cat_frac Mg_Kspar_cat_frac Fet_Kspar_cat_frac \\\n", "0 0.595210 0.000410 0.000842 \n", "1 0.591535 0.000000 0.000529 \n", "2 0.596292 0.000000 0.000381 \n", "3 0.596156 0.000000 0.000695 \n", "4 0.586453 0.000424 0.004047 \n", "\n", " Ca_Kspar_cat_frac Al_Kspar_cat_frac Na_Kspar_cat_frac K_Kspar_cat_frac \\\n", "0 0.005004 0.206133 0.035864 0.156538 \n", "1 0.007257 0.208618 0.084223 0.107838 \n", "2 0.005764 0.208468 0.055332 0.133763 \n", "3 0.003856 0.204476 0.020577 0.174240 \n", "4 0.006711 0.214735 0.031462 0.156168 \n", "\n", " Mn_Kspar_cat_frac Ti_Kspar_cat_frac Cr_Kspar_cat_frac An_Kspar \\\n", "0 0.0 0.0 0.0 0.025348 \n", "1 0.0 0.0 0.0 0.036410 \n", "2 0.0 0.0 0.0 0.029579 \n", "3 0.0 0.0 0.0 0.019410 \n", "4 0.0 0.0 0.0 0.034530 \n", "\n", " Ab_Kspar Or_Kspar K_Barth T P \n", "0 0.181675 0.792977 0.755998 1262.750535 10 \n", "1 0.422555 0.541035 0.664022 1228.315833 10 \n", "2 0.283957 0.686464 0.526704 1326.793565 20 \n", "3 0.103570 0.877020 0.936627 1217.144277 10 \n", "4 0.161891 0.803579 0.662392 1223.769830 10 " ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "T_27a_5kbar_EqTests=pt.calculate_plag_kspar_temp(plag_comps=Plags, kspar_comps=Kspars, \n", " equationT=\"T_Put2008_eq27a\", P=10*my_input['P_GPa'], eq_tests=True)\n", "T_27a_5kbar_EqTests.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Example 2 - Combining all possible Plag-Kspar Pairs\n", "- Say you don't have touching Plag-Kspar pairs, but instead a list of Kspar compositions you have measured on the probe in one sheet of an excel spreadsheet, and a set of Plag compositions in another sheet\n", "- The functions ending in \"_matching\" in Thermobar consider all possible combinations of Plags and Kspars\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Import Kspars" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SiO2_KsparTiO2_KsparAl2O3_KsparFeOt_KsparMnO_KsparMgO_KsparCaO_KsparNa2O_KsparK2O_KsparCr2O3_KsparSample_ID_Kspar
065.50.019.60.070.00.000.754.819.3601
165.40.019.40.050.00.000.593.1311.5002
264.60.018.80.090.00.000.391.1514.8003
361.80.019.20.510.00.030.661.7112.9004
465.10.019.20.050.00.000.362.8712.6005
\n", "
" ], "text/plain": [ " SiO2_Kspar TiO2_Kspar Al2O3_Kspar FeOt_Kspar MnO_Kspar MgO_Kspar \\\n", "0 65.5 0.0 19.6 0.07 0.0 0.00 \n", "1 65.4 0.0 19.4 0.05 0.0 0.00 \n", "2 64.6 0.0 18.8 0.09 0.0 0.00 \n", "3 61.8 0.0 19.2 0.51 0.0 0.03 \n", "4 65.1 0.0 19.2 0.05 0.0 0.00 \n", "\n", " CaO_Kspar Na2O_Kspar K2O_Kspar Cr2O3_Kspar Sample_ID_Kspar \n", "0 0.75 4.81 9.36 0 1 \n", "1 0.59 3.13 11.50 0 2 \n", "2 0.39 1.15 14.80 0 3 \n", "3 0.66 1.71 12.90 0 4 \n", "4 0.36 2.87 12.60 0 5 " ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "out2=pt.import_excel('Two_Feldspar_input.xlsx', sheet_name=\"Kspar\")\n", "Kspars2=out2['Kspars']\n", "Kspars2.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Import plags" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SiO2_PlagTiO2_PlagAl2O3_PlagFeOt_PlagMnO_PlagMgO_PlagCaO_PlagNa2O_PlagK2O_PlagCr2O3_PlagSample_ID_Plag
057.30.0926.60.4300.038.336.110.4900
156.50.1226.90.4700.058.955.660.4701
257.60.1126.30.5000.078.506.270.4002
357.20.1627.00.6200.069.035.580.8403
456.70.1427.60.6900.119.465.580.4804
\n", "
" ], "text/plain": [ " SiO2_Plag TiO2_Plag Al2O3_Plag FeOt_Plag MnO_Plag MgO_Plag CaO_Plag \\\n", "0 57.3 0.09 26.6 0.43 0 0.03 8.33 \n", "1 56.5 0.12 26.9 0.47 0 0.05 8.95 \n", "2 57.6 0.11 26.3 0.50 0 0.07 8.50 \n", "3 57.2 0.16 27.0 0.62 0 0.06 9.03 \n", "4 56.7 0.14 27.6 0.69 0 0.11 9.46 \n", "\n", " Na2O_Plag K2O_Plag Cr2O3_Plag Sample_ID_Plag \n", "0 6.11 0.49 0 0 \n", "1 5.66 0.47 0 1 \n", "2 6.27 0.40 0 2 \n", "3 5.58 0.84 0 3 \n", "4 5.58 0.48 0 4 " ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" } ], "source": [ "out3=pt.import_excel('Two_Feldspar_input.xlsx', sheet_name=\"Plag\")\n", "Plags3=out3['Plags']\n", "Plags3.head()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Calculate T for all matches using equation 27a at 5 kbar" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Considering 77 Kspar-Plag pairs, be patient if this is >>1 million!\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
Delta_AnDelta_AbDelta_Ora_Ab_plga_An_plga_Or_plga_Ab_kspara_An_kspara_Or_ksparT_K_calcSiO2_KsparTiO2_KsparAl2O3_KsparFeOt_KsparMnO_KsparMgO_KsparCaO_KsparNa2O_KsparK2O_KsparCr2O3_KsparSample_ID_KsparSi_Kspar_cat_propMg_Kspar_cat_propFet_Kspar_cat_propCa_Kspar_cat_propAl_Kspar_cat_propNa_Kspar_cat_propK_Kspar_cat_propMn_Kspar_cat_propTi_Kspar_cat_propCr_Kspar_cat_propsumSi_Kspar_cat_fracMg_Kspar_cat_fracFet_Kspar_cat_fracCa_Kspar_cat_fracAl_Kspar_cat_fracNa_Kspar_cat_fracK_Kspar_cat_fracMn_Kspar_cat_fracTi_Kspar_cat_fracCr_Kspar_cat_fracAn_KsparAb_KsparOr_KsparID_KsparSiO2_PlagTiO2_PlagAl2O3_PlagFeOt_PlagMnO_PlagMgO_PlagCaO_PlagNa2O_PlagK2O_PlagCr2O3_PlagSample_ID_PlagSi_Plag_cat_propMg_Plag_cat_propFet_Plag_cat_propCa_Plag_cat_propAl_Plag_cat_propNa_Plag_cat_propK_Plag_cat_propMn_Plag_cat_propTi_Plag_cat_propCr_Plag_cat_propsumSi_Plag_cat_fracMg_Plag_cat_fracFet_Plag_cat_fracCa_Plag_cat_fracAl_Plag_cat_fracNa_Plag_cat_fracK_Plag_cat_fracMn_Plag_cat_fracTi_Plag_cat_fracCr_Plag_cat_fracAn_PlagAb_PlagOr_PlagID_PlagK_BarthTP
00.2955310.002892-0.3260510.6389160.9029590.3227020.6360230.6074270.6487541269.74882265.50.019.60.070.00.000.754.819.36011.0901350.0000000.0009740.0133740.3844610.1552140.1987350.00.00.01.8428930.5915350.0000000.0005290.0072570.2086180.0842230.1078380.00.00.00.0364100.4225550.5410350.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.7632081269.7488225
1-0.0926240.110226-0.3293780.6412230.9561650.4105940.5309971.0487880.7399721197.80028265.40.019.40.050.00.000.593.1311.50021.0884710.0000000.0006960.0105210.3805380.1010020.2441720.00.00.01.8253990.5962920.0000000.0003810.0057640.2084680.0553320.1337630.00.00.00.0295790.2839570.6864641.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.5128761197.8002825
2-2.3249930.3357390.0005170.6486921.1492210.8901090.3129543.4742140.8895921013.35492964.60.018.80.090.00.000.391.1514.80031.0751560.0000000.0012530.0069550.3687680.0371090.3142380.00.00.01.8034800.5961560.0000000.0006950.0038560.2044760.0205770.1742400.00.00.00.0194100.1035700.8770202.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.1870661013.3549295
3-1.9331090.255822-0.1877060.6456051.0653630.6471450.3897832.9984710.8348511081.99801861.80.019.20.510.00.030.661.7112.90041.0285550.0007440.0070980.0117690.3766150.0551800.2738970.00.00.01.7538590.5864530.0004240.0040470.0067110.2147350.0314620.1561680.00.00.00.0345300.1618910.8035793.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.2924041081.9980185
40.0644530.110505-0.2694220.6430951.0014630.4988680.5325900.9370100.7682901145.33298665.10.019.20.050.00.000.362.8712.60051.0834780.0000000.0006960.0064200.3766150.0926120.2675270.00.00.01.8273470.5929240.0000000.0003810.0035130.2060990.0506810.1464020.00.00.00.0175130.2526530.7298344.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.4563351145.3329865
..................................................................................................................................................................................................................................................................
72-0.0064070.168807-0.2632740.6423761.0291670.5422680.4735691.0355740.8055421120.94221464.90.019.30.060.00.030.292.3013.30071.0801490.0007440.0008350.0051710.3785760.0742190.2823900.00.00.01.8220850.5928090.0004090.0004580.0028380.2077710.0407330.1549820.00.00.00.0142940.2051490.7805576.056.10.2127.80.5600.099.945.530.38060.9336880.0022330.0077940.1772550.5453070.1784480.0080680.00.0026290.01.8554220.5032210.0012040.0042010.0955340.2938990.0961760.0043480.00.0014170.00.4872710.490550.0221860.4182021120.9422145
730.5270460.050488-0.2193970.6422281.0280530.5390650.5917410.5010070.7584621122.27711364.90.018.60.000.00.000.173.0212.10081.0801490.0000000.0000000.0030320.3648450.0974530.2569110.00.00.01.8023900.5992870.0000000.0000000.0016820.2024230.0540690.1425390.00.00.00.0084820.2726740.7188437.056.10.2127.80.5600.099.945.530.38060.9336880.0022330.0077940.1772550.5453070.1784480.0080680.00.0026290.01.8554220.5032210.0012040.0042010.0955340.2938990.0961760.0043480.00.0014170.00.4872710.490550.0221860.5558551122.2771135
740.643794-0.002624-0.2636000.6378590.9955270.4521400.6404820.3517330.7157391163.39376666.00.018.90.050.00.000.183.7311.30091.0984570.0000000.0006960.0032100.3707300.1203640.2399250.00.00.01.8333810.5991430.0000000.0003800.0017510.2022110.0656510.1308650.00.00.00.0088300.3311250.6600448.056.10.2127.80.5600.099.945.530.38060.9336880.0022330.0077940.1772550.5453070.1784480.0080680.00.0026290.01.8554220.5032210.0012040.0042010.0955340.2938990.0961760.0043480.00.0014170.00.4872710.490550.0221860.6750081163.3937665
750.713796-0.043818-0.2768880.6350870.9753230.4041720.6789050.2615270.6810601191.22821866.00.019.20.000.00.000.184.3110.400101.0984570.0000000.0000000.0032100.3766150.1390800.2208160.00.00.01.8381770.5975790.0000000.0000000.0017460.2048850.0756620.1201280.00.00.00.0088400.3830280.6081329.056.10.2127.80.5600.099.945.530.38060.9336880.0022330.0077940.1772550.5453070.1784480.0080680.00.0026290.01.8554220.5032210.0012040.0042010.0955340.2938990.0961760.0043480.00.0014170.00.4872710.490550.0221860.7808141191.2282185
760.4069450.094842-0.1951440.6444931.0452310.5902080.5496510.6382860.7853531102.19822165.50.018.90.110.00.040.172.5812.500111.0901350.0009920.0015310.0030320.3707300.0832540.2654040.00.00.01.8150780.6005990.0005470.0008440.0016700.2042500.0458680.1462220.00.00.00.0086200.2367260.75465410.056.10.2127.80.5600.099.945.530.38060.9336880.0022330.0077940.1772550.5453070.1784480.0080680.00.0026290.01.8554220.5032210.0012040.0042010.0955340.2938990.0961760.0043480.00.0014170.00.4872710.490550.0221860.4825731102.1982215
\n", "

77 rows × 85 columns

\n", "
" ], "text/plain": [ " Delta_An Delta_Ab Delta_Or a_Ab_plg a_An_plg a_Or_plg a_Ab_kspar \\\n", "0 0.295531 0.002892 -0.326051 0.638916 0.902959 0.322702 0.636023 \n", "1 -0.092624 0.110226 -0.329378 0.641223 0.956165 0.410594 0.530997 \n", "2 -2.324993 0.335739 0.000517 0.648692 1.149221 0.890109 0.312954 \n", "3 -1.933109 0.255822 -0.187706 0.645605 1.065363 0.647145 0.389783 \n", "4 0.064453 0.110505 -0.269422 0.643095 1.001463 0.498868 0.532590 \n", ".. ... ... ... ... ... ... ... \n", "72 -0.006407 0.168807 -0.263274 0.642376 1.029167 0.542268 0.473569 \n", "73 0.527046 0.050488 -0.219397 0.642228 1.028053 0.539065 0.591741 \n", "74 0.643794 -0.002624 -0.263600 0.637859 0.995527 0.452140 0.640482 \n", "75 0.713796 -0.043818 -0.276888 0.635087 0.975323 0.404172 0.678905 \n", "76 0.406945 0.094842 -0.195144 0.644493 1.045231 0.590208 0.549651 \n", "\n", " a_An_kspar a_Or_kspar T_K_calc SiO2_Kspar TiO2_Kspar Al2O3_Kspar \\\n", "0 0.607427 0.648754 1269.748822 65.5 0.0 19.6 \n", "1 1.048788 0.739972 1197.800282 65.4 0.0 19.4 \n", "2 3.474214 0.889592 1013.354929 64.6 0.0 18.8 \n", "3 2.998471 0.834851 1081.998018 61.8 0.0 19.2 \n", "4 0.937010 0.768290 1145.332986 65.1 0.0 19.2 \n", ".. ... ... ... ... ... ... \n", "72 1.035574 0.805542 1120.942214 64.9 0.0 19.3 \n", "73 0.501007 0.758462 1122.277113 64.9 0.0 18.6 \n", "74 0.351733 0.715739 1163.393766 66.0 0.0 18.9 \n", "75 0.261527 0.681060 1191.228218 66.0 0.0 19.2 \n", "76 0.638286 0.785353 1102.198221 65.5 0.0 18.9 \n", "\n", " FeOt_Kspar MnO_Kspar MgO_Kspar CaO_Kspar Na2O_Kspar K2O_Kspar \\\n", "0 0.07 0.0 0.00 0.75 4.81 9.36 \n", "1 0.05 0.0 0.00 0.59 3.13 11.50 \n", "2 0.09 0.0 0.00 0.39 1.15 14.80 \n", "3 0.51 0.0 0.03 0.66 1.71 12.90 \n", "4 0.05 0.0 0.00 0.36 2.87 12.60 \n", ".. ... ... ... ... ... ... \n", "72 0.06 0.0 0.03 0.29 2.30 13.30 \n", "73 0.00 0.0 0.00 0.17 3.02 12.10 \n", "74 0.05 0.0 0.00 0.18 3.73 11.30 \n", "75 0.00 0.0 0.00 0.18 4.31 10.40 \n", "76 0.11 0.0 0.04 0.17 2.58 12.50 \n", "\n", " Cr2O3_Kspar Sample_ID_Kspar Si_Kspar_cat_prop Mg_Kspar_cat_prop \\\n", "0 0 1 1.090135 0.000000 \n", "1 0 2 1.088471 0.000000 \n", "2 0 3 1.075156 0.000000 \n", "3 0 4 1.028555 0.000744 \n", "4 0 5 1.083478 0.000000 \n", ".. ... ... ... ... \n", "72 0 7 1.080149 0.000744 \n", "73 0 8 1.080149 0.000000 \n", "74 0 9 1.098457 0.000000 \n", "75 0 10 1.098457 0.000000 \n", "76 0 11 1.090135 0.000992 \n", "\n", " Fet_Kspar_cat_prop Ca_Kspar_cat_prop Al_Kspar_cat_prop \\\n", "0 0.000974 0.013374 0.384461 \n", "1 0.000696 0.010521 0.380538 \n", "2 0.001253 0.006955 0.368768 \n", "3 0.007098 0.011769 0.376615 \n", "4 0.000696 0.006420 0.376615 \n", ".. ... ... ... \n", "72 0.000835 0.005171 0.378576 \n", "73 0.000000 0.003032 0.364845 \n", "74 0.000696 0.003210 0.370730 \n", "75 0.000000 0.003210 0.376615 \n", "76 0.001531 0.003032 0.370730 \n", "\n", " Na_Kspar_cat_prop K_Kspar_cat_prop Mn_Kspar_cat_prop Ti_Kspar_cat_prop \\\n", "0 0.155214 0.198735 0.0 0.0 \n", "1 0.101002 0.244172 0.0 0.0 \n", "2 0.037109 0.314238 0.0 0.0 \n", "3 0.055180 0.273897 0.0 0.0 \n", "4 0.092612 0.267527 0.0 0.0 \n", ".. ... ... ... ... \n", "72 0.074219 0.282390 0.0 0.0 \n", "73 0.097453 0.256911 0.0 0.0 \n", "74 0.120364 0.239925 0.0 0.0 \n", "75 0.139080 0.220816 0.0 0.0 \n", "76 0.083254 0.265404 0.0 0.0 \n", "\n", " Cr_Kspar_cat_prop sum Si_Kspar_cat_frac Mg_Kspar_cat_frac \\\n", "0 0.0 1.842893 0.591535 0.000000 \n", "1 0.0 1.825399 0.596292 0.000000 \n", "2 0.0 1.803480 0.596156 0.000000 \n", "3 0.0 1.753859 0.586453 0.000424 \n", "4 0.0 1.827347 0.592924 0.000000 \n", ".. ... ... ... ... \n", "72 0.0 1.822085 0.592809 0.000409 \n", "73 0.0 1.802390 0.599287 0.000000 \n", "74 0.0 1.833381 0.599143 0.000000 \n", "75 0.0 1.838177 0.597579 0.000000 \n", "76 0.0 1.815078 0.600599 0.000547 \n", "\n", " Fet_Kspar_cat_frac Ca_Kspar_cat_frac Al_Kspar_cat_frac \\\n", "0 0.000529 0.007257 0.208618 \n", "1 0.000381 0.005764 0.208468 \n", "2 0.000695 0.003856 0.204476 \n", "3 0.004047 0.006711 0.214735 \n", "4 0.000381 0.003513 0.206099 \n", ".. ... ... ... \n", "72 0.000458 0.002838 0.207771 \n", "73 0.000000 0.001682 0.202423 \n", "74 0.000380 0.001751 0.202211 \n", "75 0.000000 0.001746 0.204885 \n", "76 0.000844 0.001670 0.204250 \n", "\n", " Na_Kspar_cat_frac K_Kspar_cat_frac Mn_Kspar_cat_frac Ti_Kspar_cat_frac \\\n", "0 0.084223 0.107838 0.0 0.0 \n", "1 0.055332 0.133763 0.0 0.0 \n", "2 0.020577 0.174240 0.0 0.0 \n", "3 0.031462 0.156168 0.0 0.0 \n", "4 0.050681 0.146402 0.0 0.0 \n", ".. ... ... ... ... \n", "72 0.040733 0.154982 0.0 0.0 \n", "73 0.054069 0.142539 0.0 0.0 \n", "74 0.065651 0.130865 0.0 0.0 \n", "75 0.075662 0.120128 0.0 0.0 \n", "76 0.045868 0.146222 0.0 0.0 \n", "\n", " Cr_Kspar_cat_frac An_Kspar Ab_Kspar Or_Kspar ID_Kspar SiO2_Plag \\\n", "0 0.0 0.036410 0.422555 0.541035 0.0 57.3 \n", "1 0.0 0.029579 0.283957 0.686464 1.0 57.3 \n", "2 0.0 0.019410 0.103570 0.877020 2.0 57.3 \n", "3 0.0 0.034530 0.161891 0.803579 3.0 57.3 \n", "4 0.0 0.017513 0.252653 0.729834 4.0 57.3 \n", ".. ... ... ... ... ... ... \n", "72 0.0 0.014294 0.205149 0.780557 6.0 56.1 \n", "73 0.0 0.008482 0.272674 0.718843 7.0 56.1 \n", "74 0.0 0.008830 0.331125 0.660044 8.0 56.1 \n", "75 0.0 0.008840 0.383028 0.608132 9.0 56.1 \n", "76 0.0 0.008620 0.236726 0.754654 10.0 56.1 \n", "\n", " TiO2_Plag Al2O3_Plag FeOt_Plag MnO_Plag MgO_Plag CaO_Plag Na2O_Plag \\\n", "0 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", "1 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", "2 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", "3 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", "4 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", ".. ... ... ... ... ... ... ... \n", "72 0.21 27.8 0.56 0 0.09 9.94 5.53 \n", "73 0.21 27.8 0.56 0 0.09 9.94 5.53 \n", "74 0.21 27.8 0.56 0 0.09 9.94 5.53 \n", "75 0.21 27.8 0.56 0 0.09 9.94 5.53 \n", "76 0.21 27.8 0.56 0 0.09 9.94 5.53 \n", "\n", " K2O_Plag Cr2O3_Plag Sample_ID_Plag Si_Plag_cat_prop Mg_Plag_cat_prop \\\n", "0 0.49 0 0 0.95366 0.000744 \n", "1 0.49 0 0 0.95366 0.000744 \n", "2 0.49 0 0 0.95366 0.000744 \n", "3 0.49 0 0 0.95366 0.000744 \n", "4 0.49 0 0 0.95366 0.000744 \n", ".. ... ... ... ... ... \n", "72 0.38 0 6 0.933688 0.002233 \n", "73 0.38 0 6 0.933688 0.002233 \n", "74 0.38 0 6 0.933688 0.002233 \n", "75 0.38 0 6 0.933688 0.002233 \n", "76 0.38 0 6 0.933688 0.002233 \n", "\n", " Fet_Plag_cat_prop Ca_Plag_cat_prop Al_Plag_cat_prop Na_Plag_cat_prop \\\n", "0 0.005985 0.148545 0.521768 0.197164 \n", "1 0.005985 0.148545 0.521768 0.197164 \n", "2 0.005985 0.148545 0.521768 0.197164 \n", "3 0.005985 0.148545 0.521768 0.197164 \n", "4 0.005985 0.148545 0.521768 0.197164 \n", ".. ... ... ... ... \n", "72 0.007794 0.177255 0.545307 0.178448 \n", "73 0.007794 0.177255 0.545307 0.178448 \n", "74 0.007794 0.177255 0.545307 0.178448 \n", "75 0.007794 0.177255 0.545307 0.178448 \n", "76 0.007794 0.177255 0.545307 0.178448 \n", "\n", " K_Plag_cat_prop Mn_Plag_cat_prop Ti_Plag_cat_prop Cr_Plag_cat_prop \\\n", "0 0.010404 0.0 0.001127 0.0 \n", "1 0.010404 0.0 0.001127 0.0 \n", "2 0.010404 0.0 0.001127 0.0 \n", "3 0.010404 0.0 0.001127 0.0 \n", "4 0.010404 0.0 0.001127 0.0 \n", ".. ... ... ... ... \n", "72 0.008068 0.0 0.002629 0.0 \n", "73 0.008068 0.0 0.002629 0.0 \n", "74 0.008068 0.0 0.002629 0.0 \n", "75 0.008068 0.0 0.002629 0.0 \n", "76 0.008068 0.0 0.002629 0.0 \n", "\n", " sum Si_Plag_cat_frac Mg_Plag_cat_frac Fet_Plag_cat_frac \\\n", "0 1.839397 0.518464 0.000405 0.003254 \n", "1 1.839397 0.518464 0.000405 0.003254 \n", "2 1.839397 0.518464 0.000405 0.003254 \n", "3 1.839397 0.518464 0.000405 0.003254 \n", "4 1.839397 0.518464 0.000405 0.003254 \n", ".. ... ... ... ... \n", "72 1.855422 0.503221 0.001204 0.004201 \n", "73 1.855422 0.503221 0.001204 0.004201 \n", "74 1.855422 0.503221 0.001204 0.004201 \n", "75 1.855422 0.503221 0.001204 0.004201 \n", "76 1.855422 0.503221 0.001204 0.004201 \n", "\n", " Ca_Plag_cat_frac Al_Plag_cat_frac Na_Plag_cat_frac K_Plag_cat_frac \\\n", "0 0.080757 0.283663 0.107189 0.005656 \n", "1 0.080757 0.283663 0.107189 0.005656 \n", "2 0.080757 0.283663 0.107189 0.005656 \n", "3 0.080757 0.283663 0.107189 0.005656 \n", "4 0.080757 0.283663 0.107189 0.005656 \n", ".. ... ... ... ... \n", "72 0.095534 0.293899 0.096176 0.004348 \n", "73 0.095534 0.293899 0.096176 0.004348 \n", "74 0.095534 0.293899 0.096176 0.004348 \n", "75 0.095534 0.293899 0.096176 0.004348 \n", "76 0.095534 0.293899 0.096176 0.004348 \n", "\n", " Mn_Plag_cat_frac Ti_Plag_cat_frac Cr_Plag_cat_frac An_Plag Ab_Plag \\\n", "0 0.0 0.000613 0.0 0.417129 0.553656 \n", "1 0.0 0.000613 0.0 0.417129 0.553656 \n", "2 0.0 0.000613 0.0 0.417129 0.553656 \n", "3 0.0 0.000613 0.0 0.417129 0.553656 \n", "4 0.0 0.000613 0.0 0.417129 0.553656 \n", ".. ... ... ... ... ... \n", "72 0.0 0.001417 0.0 0.487271 0.49055 \n", "73 0.0 0.001417 0.0 0.487271 0.49055 \n", "74 0.0 0.001417 0.0 0.487271 0.49055 \n", "75 0.0 0.001417 0.0 0.487271 0.49055 \n", "76 0.0 0.001417 0.0 0.487271 0.49055 \n", "\n", " Or_Plag ID_Plag K_Barth T P \n", "0 0.029215 0 0.763208 1269.748822 5 \n", "1 0.029215 0 0.512876 1197.800282 5 \n", "2 0.029215 0 0.187066 1013.354929 5 \n", "3 0.029215 0 0.292404 1081.998018 5 \n", "4 0.029215 0 0.456335 1145.332986 5 \n", ".. ... ... ... ... .. \n", "72 0.02218 6 0.418202 1120.942214 5 \n", "73 0.02218 6 0.555855 1122.277113 5 \n", "74 0.02218 6 0.675008 1163.393766 5 \n", "75 0.02218 6 0.780814 1191.228218 5 \n", "76 0.02218 6 0.482573 1102.198221 5 \n", "\n", "[77 rows x 85 columns]" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "T_27a_5kbar_All=pt.calculate_plag_kspar_temp_matching(plag_comps=Plags3, kspar_comps=Kspars2, \n", " equationT=\"T_Put2008_eq27a\", P=5)\n", "T_27a_5kbar_All" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Examine equilibrium tests\n", "- As it is currently unclear what delta values are acceptable, we leave it up to the users to investigate the experimental data, but provide an example here of how you can filter based on these values\n", "- First, plot up the delta values to get an idea of the degree of disequilibrium" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, '# of measurements')" ] }, "execution_count": 21, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3YAAAE9CAYAAABHrfALAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAkg0lEQVR4nO3deZxkZX3v8c/XGVBU3ECDE5XBBTQioLaSuGuuS8CouGJYRMSW6KgQic69JkZIDNwrEuJ1AxHB3asSAcddgxu4dBAQRS8YRyCgXhEXFJWB3/2jarBoejk93dWnTvXn/XrVa+qcOn3qe6aqnu5fPc95TqoKSZIkSVJ33aLtAJIkSZKkxbGwkyRJkqSOs7CTJEmSpI6zsJMkSZKkjrOwkyRJkqSOs7CTJEmSpI5b3XaAppJ4XQZpDFVV2s6wGLZN0viyfZI0imZrmzpT2AF4zT1pvCSd/pvpRrZN0vixfZI0iuZqmxyKKUmSJEkdZ2EnSZIkSR1nYSdJkiRJHWdhJ0mSJEkdZ2EnSZIkSR1nYSdJkiRJHWdhJ0mSJEkdZ2EnSZIkSR1nYSdJkiRJHWdhJ0mSJEkdt7rtANI4WLt+w1D2u/GYvYeyX2mc+XmUVgY/69JN2WMnSZIkSR1nYSdJkiRJHWdhJ0mSJEkdZ2EnSZIkSR1nYSdJkiRJHWdhJ0mSJEkdZ2EnSZIkSR1nYSdJkiRJHWdhJ0mSJEkdZ2EnSZIkSR1nYSdJkiRJHWdhJ0mSJEkdZ2EnSZIkSR1nYSdJQHpOTXLEDI+dluRNbeSSJElqwsJO0oqX5H7A54BnzvDYK4FHLnsoSZKkBVjddgBJGgEvAU4CLh1cmeQxwJOAtwF3XPZUkiRJDdljJ2nFq6p1VfW+wXVJ1gD/CuwHXN9KMEmSpIbssZOkaZJsBbwfOLyqrkwy5/YTExM33p+cnGRycnK4ASVJkqaxsJOkm5sA7gkc1y/qdgBWJblVVR0yfeOpqalljidJknRTFnaSNE1VnQPcffNyktcC21fVutZCSZIkzcFz7CRJkkZQkv2TnJ/kvCRnJ5lIsirJ8Um+m+SSJIe2nVPSaLDHTpL6quqgWda/dnmTSFrpkuwCvB54UP9c372A04BjgJ2BXYFtgXOSnFtVX28vraRRYGEnrVBr128Yyn43HrP3UPYrSSvM74BDqurK/vIUvfN9nwW8pao2AVcn+QCwP2BhJ61wFnaSJEkjpqo2AhsB0pvF6TjgDHo9dZcNbHo5sNts+3HWXmnlsLCTJEkaUUluA5xCb0KnJ9HrmavBTZjjWpvO2iutHE6eIkmSNIKS3AM4m17h9tiq+jlwKbBmYLM19HrtJK1wFnaSJEkjJsm2wFnAaVW1b1Vd23/odODgJKuT3AHYF/hoKyEljRSHYkqSJI2edcCOwD5J9hlY/0TgXsD5wNbACVX1hRbySRoxFnaSJEkjpqqOBo6e5eHDljGKpI5wKKYkSZIkddyyFnZJ9k9yfpLzkpydZCLJqiTHJ/lukkuSHLqcmSRJkiSp65ZtKGaSXYDXAw+qqiuT7AWcBhwD7EzvuizbAuckObeqvNCmJEmSJDWwnD12vwMOqaor+8tTwA7As4B3VtWmqroa+ACw/zLmkiRJkqROW7Yeu6raCGwESBLgOOAMej11lw1sejmw20z7mJiYuPH+5OQkk5OTwwkrSZIkSR2y7LNiJrkNcApwd+BJwNeBGtyE3oU4b2ZqamrY8SRJkiSpc5a1sEtyD+BM4CLgsVV1bZJLgTUDm62h12snSZKkDlu7fkPbEaQVY9nOsUuyLXAWcFpV7VtV1/YfOh04OMnqJHcA9gU+uly5JEmSJKnrlrPHbh2wI7BPkn0G1j8RuBdwPrA1cEJVfWEZc0mSJElSpy3n5ClHA0fP8vBhy5VDkiRJksbNsl6gXJIkSZK09CzsJEmSJKnjLOwkSZIkqeMs7CRJkiSp4yzsJEmSJKnjlvUC5ZIWxgu7SpIkqQl77CRJkiSp4yzsJEmSJKnjLOwkSZIkqeMs7CQJSM+pSY7oL2+T5OQkFyb5dv/+Nm3nlCRJmomFnaQVL8n9gM8BzxxY/Wp6E0zt1r9tA/z35U8nSZI0P2fFlCR4CXAScOnAui8CG6vqBoAk3wTu30I2SZKkeVnYSVrxqmodQJInDKz79Ob7SXYEDgMmZ/r5iYmJG+9PTk4yOTnjZkM1rEtjbDxm76Hst4uGefkR/58lSYtlYSdJc0jyYODfgDdV1cdm2mZqamp5Q0mSJE1jYSdJs0iyL/AWYF1Vva/tPJIkSbOxsJOkGST5S+CNwBOqyi45SZI00izsJGlmxwIBTkqyed1Xquol7UWSJEmamYWdJPVV1UED93dpMYokSdKCeB07SZIkSeo4CztJkiRJ6jgLO0mSJEnqOAs7SZIkSeo4CztJkiRJ6jgLO0mSJEnqOAs7SZIkSeo4CztJkiRJ6rjGhV2Sxye5S//+85JsSHJkkq2GF0+SJEmSNJ9GhV2S9cBHgXsmeThwEnAF8Bzg2KGlkyRJkiTNq2mP3aHAs6vqq8ABwNlV9ULgQGDfYYWTJEmSJM2vaWF3F+CC/v0nA2f2718FbLPUoSRJkiRJza1uuN13gIOS/ARYA5yeZGvgb4HzhpRNkiRJktRA08LuCODDwJ2Ao6vq4iRvBp5CrwdPkiRJHbV2/Ya2I0hapEaFXVWd1Z8R8/ZVdXV/9VHAYcDth5RNkiRJktRA01kxrwe2GyjqqKofA38MbBxONEmSJElSE7P22CV5HvCCzYvAGUmum7bZXeld9kCSJEmS1JK5hmJ+GNiRXlH3CODLwDUDj1d/+SNDSydJkiRJmteshV1V/ZreeXQk2Qh8sKp+u0y5JEmSJEkNNZ085dQk900yAWxFrxdv8PGThxFOkiRJkjS/RoVdkvXAPwM/A3417eECLOwkSZIkqSVNr2N3KPDqqjp6mGEkSZIkSQvX6HIHwB3oTaYiSZKkZZKeU5McMbDup0nOG7jt12ZGSaOhaY/du4FDkxxRVTXMQJIkSYIk9wPeDOwJfKu/bhfgZ1W1R4vRJI2gpoXd9sDTgf36M2T+fvDBqnrUEueSJEla6V4CnARcOrDuYcD1Sb4E3J7eiKrXVdX1LeSTNEKaFnbfpTd5iiSNpSQBTgG+VVXHJlkFvAF4Er228tiqeluLESWtMFW1DiDJEwZWrwY+C6ynN1P5BuCXwPEz7WNiYuLG+5OTk0xOTg4praS2Nb3cwZGb7ydZDVzvkExJ42Km4U7Ai4CdgV2BbYFzkpxbVV9vJ6UkQVW9fXA5yXHAy5ilsJuamlqGVJJGQdPJU0jykiSXAL8BdkpyYpJ/TtJ4H5I0ojYPd/rQwLp9gHdW1aaquhr4ALB/G+EkabMkByTZbXAVcF1beSSNjkZFWZK/AV4JHMUfGo9PAwcD/zScaJK0PKpqXVW9b9rquwOXDSxfDtxt+VJJ0ox2BY5KsirJNsA64IMtZ5I0Apr2th0KTFbVu4AbAKrqw8ABwIFNn8wpeyV1yC2AwSHnAWacnGBiYuLG24knnrgs4SStWEcCP6M3bPwC4Gx6Iw4krXBNJ0+5G3DxDOsvA+7YZAdO2SupYy4F1gwsr6HXa3cznsMiaZiq6qCB+7+hN2JKkm6iaY/dOcDzB5arf27dq4CvNdzHTOew3Dhlb5ILkrymPxOdJLXtdODgJKuT3AHYF/hoq4kkSZJm0bTH7mXAJ5I8GbgVcAK92eJuSW8q8Hk5Za+kjnkrcC/gfGBr4ISq+kK7kSRJkmbW9HIH306yM7AfcN/+z50GvKeqrtnSJ3fKXkmjZNpwp03AYa2FkSRJWoCmPXZU1W+Bdyzlkyc5ADi/qi7YvAqn7JUkSZKkBWl6uYOHJflGkl8l+f302yKe3yl7JUmSJGmRmvbYvZPerJivA65dwuc/EngTvVkyt6I3sYpT9kqSJEnSAjQt7NYAT6mq7y32CZ2yV5IkSZKWVtPLHbwfeNYwg0iSJEmStkzTHrtjgKkkLwB+CNww+GBVPW6pg0mSJEmSmmla2L0X+ClwBkt7jp0kSZIkaZGaFna7Aw+qqu8OM4wkSZIkaeGanmP3ZXqXJpAkSZIkjZimPXafBd6R5JnAfwKbBh+sqtcsdTBJkiRJUjNNC7u/AM4F/qh/G1RLmkiSJEmStCCNCruqeuywg0iSJEmStkzTHjuS7Aa8FLgPsB+wD3BJVX1ySNkkSWNs7foNbUeQJGlsNJo8JckTgXOAVcCewC2BOwNnJPmr4cWTJEmSJM2n6ayYrwMOr6qD6U+cUlX/ALwM+PshZZMkSZIkNdC0sLsf8JkZ1n8GWLtkaSRJkiRJC9a0sPsB8GczrP9Lepc/kCRJkiS1pOnkKX8HvCvJQ/o/8/wk9wSeRW8iFUmSJElSSxr12FXVR4FHANsBFwJPpjeRyiOr6kNDSydJkiRJmlejHrskJwHHVNWBQ84jSZIkSVqgpufYPR24fphBJEmSJElbpuk5dscBb0tyPPBD4LeDD1aVE6hIkiRJUkuaFnZH9f99fP/f6v+b/v1VSxlKkiRJktRc08Jup6GmkCRJkiRtsaaF3WVDTSFJkiRJ2mJNC7tN/GH45UwciilJkiRJLWla2D12hp+7F3A48OolTSRJIyTJPsCRwA3Az4AXVtX3200lSZJ0U40Ku6r6wgyrP5fke8C/AqctaSpJGgFJtgHeA+xeVZckORx4I7B3u8kkSZJuqul17GbzU2CXpQgiSSNoFb3Zf2/fX74t0y73IkmSNAoa9dglOXiG1bcDDgXOWdJEkjQiquqaJIcCZye5il6h9/Dp201MTNx4f3JyksnJyVn3uXb9hiEklSRJK13Tc+z+ftpyAb8HvgH83ZImkqQRkeQBwGuAP6mq7yd5GfCRJHtU1Y0TSk1NTbWWUZIkCZqfY+d17CStRE8EvjIwWcqbgX8BtqM3FF2SJGkkND7HLsnjk9ylf/95STYkOTLJVsOLJ0mtOhd4dJI/6i8/DfhBVVnUSZKkkdKosEuyHvgocM8kDwdOAq4AngMcO7R0ktSiqvo88HrgrCTnA+uAp7abSpIk6eaanmN3KPDsqvpqkrcBZ1fVC5M8FDgTePnQEkpSi6rqzfSGYEqSJI2spkMx7wJc0L//ZHrFHMBVwDZLHUqSJEmS1FzTHrvvAAcl+QmwBjg9ydbA3wLnDSmbJEmSJKmBpoXdEcCHgTsBR1fVxUneDDyFXg+eJEmSJKklTS93cFZ/RszbV9XV/dVHAYdV1XVDSydJkiRJmlfTHjvonWe3S5JV/eUAt0zywKp63dJHkyRJkiQ10aiwS3Io8L+BVUDRK+ro3/8qYGEnSZIkSS1pOivmeuCf6M2A+WNgR2BX4JvA6cOJJkmSJElqomlhtwY4tap+B5wL/FlVfYfe9eteOKxwkiRJkqT5NS3sfgzcuX//u8AD+/evoFf0SZIkSZJa0nTylA8A70pyCPBJ4L1JzqN3qYOLh5RNkiRJktRA08LuvwM/B7arqjOSvB14E3AV8PwhZZMkSZIkNdBoKGZVbaqq11XVGf3lv6uqO1fVfavqnOFGlCRJWpnSc2qSI/rLq5Icn+S7SS7pz1wuSc2vY5dkb+Bw4N7Ao4FDgMuq6sQhZZMktWzt+g1tR5BWrCT3A94M7Al8q7/6RcDO9GYn3xY4J8m5VfX1dlJKGhWNeuyS7A+8G/givQuVrwJ+BByX5PDhxZMkSVqxXgKcBHxoYN0+wDv7o6mupjcPwv5thJM0WprOivkq4EVVdRRwPUBVvRk4CHjZcKJJkiStXFW1rqreN2313YHLBpYvB+62fKkkjaqmQzHvBUzNsP48YIclSyNJkqS53AKogeXQ/9J9JhMTEzfen5ycZHJycnjJJLWqaWH3LWBvejNhwh8alIOBC5o+WZIApwDfqqpjk6wC3gA8qZ/l2Kp6W9P9SZIkrTCXctNrCK+h12s3o6mpmb6XlzSOmhZ2rwA2JPlzYGvgNUl2Bvagdy27eXkCsCRJ0qKdDhyc5EzgtsC+gDNjSmpW2FXVl5PsArwY2ATcAfgS8NyqurThc20+AXhw+32AE6tqE3B1ks0nAFvYSZIk3dxb6Z0icz69L9tPqKovtBtJ0ihofLmDqvoR8JotfaKqWgeQ5AkDq2c6AXi3LX0OSZKkcVNVBw3c3wQc1loYSSOrUWHX7637J+C+wC2nP15VO2/h83sCsCRJkiQtUtMeu/fTK7jeAVy7hM/vCcCSJEmStEhNC7tdgIdU1XeW+Pk9AViSJEmSFqlpYfdJ4M+ApS7sPAFYkiRJkhapaWF3OPDNJPsDPwRuGHywqg5u+oSeACxJkiRJS+sWDbc7gd4kJ1cBq4Ctpt0kaSwleUCSs5J8M8lUkge3nUmSJGm6pj12jwYeUVXnDjOMJI2SJLcGPg28oKo+nuSpwHvpzRAsSZI0MpoWdhcCdxxmEEkaQU8Avl9VH+8vnwH8oMU8kiRJM2pa2L0deHeSU+n9UbNp8MGqOnmpg0nSCNgZ+FGSdwC7Az8HXtlqIkmSpBk0Lez+B/A7epcjmK4ACztJ42grYC/gsVX1tf5QzI8n2bGqfrd5o4mJiRt/YHJyksnJyeVPKs1g7foNQ9v3xmP2Hsp+h5V5WHklaVQ0KuyqaqdhB5GkEXQFcFFVfQ2gqk5PchJwT+CizRtNTU21FE+SJKmn6ayYkrQSfQLYafNMmEkeRW+UgufZSZKkkdJ0KKYkrThV9aMkTwPekuQ29IakP72qfttuMkmSpJuysJOkOVTVF4E9284hSZI0l1mHYiY5Ncl2/fuPSmIRKEmSJEkjaK5z7J4N3K1//9/xOnaSJEmSNJLm6oX7BPC1JP8PCPDNJDfMtGFV3WMY4SRJkiRJ85ursHsW8GTgDsA7gf8F/GIZMkmSJEmSFmDWwq6qrgdOB0iyFjipqn6zTLkkSZIkSQ01uo5dVR0J/EmSDya5IMm3k5yW5HFDzidJkiRJmkejwi7JM4GzgRuAt/dv1wGfTPLU4cWTJEmSJM2n6SUMjgReWVXHD6w7Pslh/cdOX+JckiRJkqSGGvXYATsBH5th/ceAXZYujiRJkiRpoZoWdhcBfzHD+r2BHyxdHEmSJEnSQjUdivkPwEeSPAz4Wn/dnsDTgb8aRjBJkiRJUjNNZ8X8GL0eu62BSeBAehOpPKyqPjK8eJIkSZKk+TTtsaOqPg98fohZJEmSJElboOk5dpIkSZKkEWVhJ0mSJEkdZ2EnSZIkSR3X6By7JCcDL6+qX/WX3wK8pqp+OsxwkiRJ0nJau37DUPa78Zi9h7JfabNZe+yS/DTJp5IcDTwPuH+SzdvvD9xuOQJKkiRJkuY211DMBwBvBK4FApwGXJPkXOBWwAFJ9kyy9fBjSpIkSZJmM2thV1VXVtWGqjqqv+qBwA7Ay4BNwK7AW4Grhp5SkiRJkjSrWc+xS/IpYAo4Fyjghqr6JfDlJNcDr6qq/1xpPXaOu5YkSZI0auaaPOVNwIPonV8X4PtJLgQuBLYCdk9yZVVdO/yYkiRJkqTZzFrYVdWZwJkASW4A9gS2A+4PHAj8I/CeJBdX1R7DjypJkiRJmkmjyx0APwR+UlUX0RuK+VLgL4AfAbsPK5wkSZIkaX6NCruq2mna8q4Di1NLmkiSJEmStCBzXe5AkiRJktQBFnaSNI8kT0vyq7ZzSJIkzcbCTpLmkOQ+wLH0ZgeWJEkaSRZ2kjSLJLcG3gP8TdtZJEmS5tJ0VkxJWolO6N8umGujiYmJG+9PTk4yOTk55FiSJEk3ZWEnSTNI8mJgU1WdnGTtXNtOTTk5sCRJapeFnSTN7CDg1knOA7YGtunf36uqrmgxlyRJ0s1Y2EnSDKrqoZvv93vsLqyqPVoLJEmSNAcnT5EkSZKkjrPHTpLmUVUbgdu2nUOSJGk29thJkiRJUsdZ2EmSJHVMkjckuTTJef3bB9vOJKldDsWUJEnqnocB+1bV2W0HkTQa7LGTJEnqkCS3BB4IvDLJt5J8JMk92s4lqV0j0WOX5A3As4Cf9Vd9r6qe02IkSZKkUbUG+Dzwd8C3gSOA05M8qKpqcMOJiYkb709OTjI5ObmcOTVg7foNQ9v3xmP2Htq+1R0jUdjhcAJJkqRGquoHwF6bl5McC/w9sBb4weC2U1NTy5pNUntaH4rpcAJJkqTmkuyW5IDpq4Hr2sgjaTSMQo/dkg8nGGZXtyRJUstuAN6Y5Mv93ru/Bi6oqstbziWpRa0Xdg4nkCRJaq6qLkzyUuDMJKuAy4HnthxLUstaL+yS7AbsXlXvHlyNwwkkSZJmVFXvAd7Tdg5Jo6P1c+z4w3CCnfrLDieQJEmSpAVovcfO4QSSJEmStDitF3bgcAJJkiRJWoxRGIopSZIkSVoECztJkiRJ6jgLO0mSJEnqOAs7SZIkSeo4CztJkiRJ6jgLO0mSJEnqOAs7SZIkSeo4CztJkiRJ6jgLO0mSJEnqOAs7SZIkSeo4CztJkiRJ6jgLO0mSJEnqOAs7SZpDkv2TnJ/kvCRnJ5loO5MkSdJ0q9sOIEmjKskuwOuBB1XVlUn2Ak4D7tFuMkmSpJuyx06SZvc74JCqurK/PAXskGTrFjNJkiTdjD12kjSLqtoIbARIEuA44Iyq+n2LsSRJuom16ze0HWHBNh6zd9sRxo6FnSTNI8ltgFOAuwNPmv74xMQfTrubnJxkcnJy2bJJbeniH5KSNM4s7CRpDknuAZwJXAQ8tqqunb7N1NTUsueSJEka5Dl2kjSLJNsCZwGnVdW+MxV1kiRJo8AeO0ma3TpgR2CfJPsMrP/zqrqqpUySJEk3Y2EnSbOoqqOBo9vOIUmSNB+HYkqSJElSx1nYSZIkSVLHWdhJkiRJUsdZ2EmSJElSx1nYSZIkSVLHWdhJkiRJUsdZ2EmSJElSx3kduxGxdv2Goex34zF7D2W/kiRJkkaHPXaSJEmS1HEWdpIkSZLUcRZ2kiRJktRxFnaSJEmS1HEWdpIkSZLUcRZ2kiRJktRxFnaSJEmS1HEWdpIkSZLUcRZ2kiRJktRxFnaSJEmS1HGr2w4gSZIkaWVZu37DUPa78Zi9h7LfYeWFpctsj50kSZIkdZyFnSRJkiR1nIWdJEmSJHWchZ0kSZIkdZyF3Rb41XmfbDvCyDvxxBPbjtAJvpfGy7i978fp/TlOxzJO77Nxel26blxeC49jtIzLcXSl3bWw2wLXnD8eb9Jh6soHoG2+l8bLuL3vx+n9OU7HMk7vs3F6XbpuXF4Lj2O0jMtxdKXdtbCTJEmSpI4bicIuyd5JLkjyvSQfSnK7tjNJEtg+SRpNtk2Spmu9sEtyZ+CdwDOqahfgP4Fj2k0lSbZPkkaTbZOkmbRe2AFPAL5RVRf3l98K7JckLWaSJLB9kjSabJsk3Uyqqt0AyXpgbVUd2l9eDVwH3L6qfjmwXbtBJQ1FVY3sHyJN2ifbJml8jWr75N9O0so2W9u0ermDzOAWwEwNz/WDC6PauEoaa/O2T7ZNklrg306SbmYUhmJeCqwZWP5j4Oqq+nVLeSRpM9snSaPItknSzYxCYfdp4E+T3Ke/fChweot5JGkz2ydJo8i2SdLNtF7YVdVPgOcDH05yEfAA4BXtpppfkv2TnJ/kvCRnJ5loO9MoSs+pSY5oO8socZrqbpivfVrI6ziqn4Umx9CV92vTnKP6Wgxq+Lp04vdQw2NZl+TbSS5McnqSu7SRtSuW4m+nhX6ukzwtya+2PPVwLOBzP3Kfl3Fpf8elvVrg7/SR/DxQVd4WeAN2Aa4E7tpf3gu4tO1co3YD7gd8Hvg1cETbeUblBtwZ+Alwn/7y/wTe0nYub8N7HUf1s9DkGLryfm2ac1Rfiy14XTrxe6jhsTwY2Ehv4g+AY4ET2s4+zreFfq6B+wCXANe0nX1LjmMUPy/j0v6OS3u1wN/pI/l5qKr2e+w66nfAIVV1ZX95CtghydYtZhpFLwFOAj7UdpAR4zTV42Ehr+OofhaaHENX3q9Nc47qazGoybF05ffQvMdSVf9B74+pXyS5Fb3zxa5a/qgrSuPPdZJbA+8B/mYZ8zXV9DhG8fMyLu3vuLRXjf6vR/zzYGE3lyR7Jdk0/QY8qqo29LcJcBxwRlX9vtXALZjt/yjJgVW1rqre13bGEXR34LKB5cuB2wHbthNHc5mjHbgXDV/HEf4sNHkvduX92ijnCL8Wg+Y9lqra2JHfQ01fl+uSPK3/+KPoXXxbi7QU7RdwQv92wfATz2yxxzGin5dxaX/Hpb1q+n/d+udhLhZ2c6iqj1fV6hlu7wJIchvg/wD3Bg5pNWxL5vs/0owaTVOt0TDbexzYRPdfxybvxa68X7uSs4nGx9KB30ONj6WqPlpV2wOvBT6VxL9RFmmx7VeSFwObqurk5cg7m6Vqh0fs8zIu7e+4tFfzHseofB7mYqO5hZLcAzib3gv+2Kr6ebuJ1CFOUz0exuF1bHIMXTnOruRsotGxdOT30LzHkuTeSR4xsM3JwI7AHZcn4orU9PNyEPCQJOcBHwe26U9+sYbR0PhzP4Kfl3Fpf8elvWpyHAcx2p8HC7stkWRb4CzgtKrat6qubTmSusVpqsfDOLyOTY6hK8fZlZxNzHssHfo91OR1uSvwgSTb95f3Ay6sKs+zG55Gn5eqemhV7VpVe9Cb8OLaqtqjqq5YvqhzanQcI/p5GZf2d1zaq3mPowOfBwu7LbSO3reJ+/Qr9c237doOptFXHb3Eh25qrtcxyUT/G72RNtsxDObvyvu1ybF0RcNj6cTvoYbvsS8BrwPO6q/bF3haK4FXiHFov2BBxzFyn5dxaX/Hpb0al98hqZppOKkkSZIkqSvssZMkSZKkjrOwkyRJkqSOs7CTJEmSpI6zsJMkSZKkjrOwkyRJkqSOs7CTJEmSpI6zsNOCJNmYpPq3G5Jck+QrSZ64gH08pv/zq/vLeyR55CJzndDf54MXsx9J3dTFtqm//r8tZv+SRt8It0/PSHJ2kl8n+WmSDUkevph9ql0WdtoSrwDuCtwN+FPgK8CGRfyB8m/ALlsaJslWwDOAS4Dnbel+JHWebZOkUTVq7dOrgXcBHwH2AB4DXAh8Psm+W7pftWt12wHUSb+sqh/1718BvDLJXYF/AR6wBfvLIvM8EdgaOA44Kskrquq6Re5TUvfYNkkaVSPTPiV5IHAUsFdVfWrgoVcluQp4W5LPV9VPtvQ51A577LRUTgR2TXJvgCS3T3Jqkl8k+VGSE5NsO/2HkpwF7Ai8Pckp/XVPTnJukt/2f/6DSW43x3M/F/gycDqwHfDkac9RSQ5MckF/+MO/J7nXUhy0pJE3sm1T3yOT/N8kv0ny3iS3XdTRSuqSttqnFwDnTivqNvuX/r/P2fxcSd6U5JIk/5Vk+0UdsYbKwk5L5Tv9f/+k/+/JwPbAI4G96Q0XOGWGn3s6cDm9IQovT7ITvWEBbwPuCzwLeBxw6ExPmuTWwFOBM6vqCmCKmYc8/QNwOL2hBjsA/7yQg5PUWaPeNv018PJ+nl37+5e0MrTSPgEPAb4+0wP9UQVfpTdcdLPnAwcBT62qn85/WGqLQzG1VH7R/3fbfm/YPsD2VfUzgCQHAhuT3H3wh6rqZ0mupzdE4RdJ7gK8vKpO7G+yMclngfvP8rxPAW5N7xtxgNPoDXnaflrjc3xVfa6f5a3AYYs5WEmdMept0z9W1Sf6WV4GfC7Juqr6+WIOWlIntNU+3Qm4eo5cV9MrMDf7RFV9eUFHplbYY6elsrm7/5fA/eiN/b60P/TxGuDb/cd3nmsnVXUx8PEkr07y/iQXAM8GVs3yI38FfKP/jTj0/njaqr9+0PcH7v+yv42k8TfqbdM3Bu6f29/ffeY/LEljoK326Wf0Ri/NZg1w1cDyxjmPQiPDHjstld36/14I7A78mt4sS9NdSW8IwIyS7E5vpqgzgS/Rm3TgsFm2vSO9yQm2SrJp2sMHAW8cWP799B+fLYOksTLqbdMNA/c3f9k6vb2SNJ6WvX3q+yrw6Fn2tQ29nr7TBlb/do59aYRY2GmpHAz8R1X9IMmtgNsAq6rqewD9E4OPA140w8/WwP0DgK9U1XM3r0hyH+DiGX7umfTew4+h9+3TZs8AXpvkAVX1rS0/JEljYNTbpgfwh3NdHkqvqBscYSBpfLXRPgG8HXhxkmdW1YenPfZyerP5fmBLDkjtsrDTlrhdkh3o9XptT292pX2BxwNU1UVJPgm8O8lL6X3T81Z6jdWVSaZfd+Ua4L5J7kSv63/XJHvS+4PoUHrfUl06Q47nAp+pqi8OrkxyGfBKehMVHLEUByypE7rYNv1jkh8Cv6LXk3diVV2zmP8ESSNpVNonqurCJOuBU/rn753JH4aK/y3w/Kr68RIeu5aJ59hpS7yB3rCAK4DP0pu16XFV9YWBbQ6g903Rp4EvAP9Fb4a4mbyJ3rdRb6f3h81XgM8AZwNrgSOZNjQhvWu/PBp4x/SdVdUvgPcD+yXxywtp5ehi2/T6/v4/B5xDr/CTNH5ab58GVdUb6H0J9TTgP+hdmuWBwJ9X1fsXfngaBamq+beSJEmSJI0se+wkSZIkqeMs7CRJkiSp4yzsJEmSJKnjLOwkSZIkqeMs7CRJkiSp4yzsJEmSJKnjLOwkSZIkqeMs7CRJkiSp4/4/J1mWxU3nkEEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ((ax0, ax1, ax2)) = plt.subplots(1,3, figsize = (15,5)) # adjust dimensions of figure here\n", "ax0.hist(T_27a_5kbar_All['Delta_An'])\n", "ax1.hist(T_27a_5kbar_All['Delta_Ab'])\n", "ax2.hist(T_27a_5kbar_All['Delta_Or'])\n", "ax0.set_xlabel('Delta An')\n", "ax1.set_xlabel('Delta Ab')\n", "ax2.set_xlabel('Delta Or')\n", "ax0.set_ylabel('# of measurements')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Filtering outputs\n", "- Here we make filters, the \"filt_a_...\" considers all avlues of Delta An <1, and is a bool, so a series of true and false statements" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [], "source": [ "# This says accept pairs within Delta An between +-1\n", "filt_a_An=abs(T_27a_5kbar_All['Delta_An'])<1 \n", "\n", "# This says accept pairs with Delta Ab between +-0.15\n", "filt_a_Ab=abs(T_27a_5kbar_All['Delta_Ab'])<0.15 \n", "\n", "# This says accept pairs with Delta Or between +-0.2\n", "filt_a_Or=abs(T_27a_5kbar_All['Delta_Or'])<0.2 " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Apply filters\n", "- Here, we apply these filters to only consider pairs meeting all 3 criteria above" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
Delta_AnDelta_AbDelta_Ora_Ab_plga_An_plga_Or_plga_Ab_kspara_An_kspara_Or_ksparT_K_calcSiO2_KsparTiO2_KsparAl2O3_KsparFeOt_KsparMnO_KsparMgO_KsparCaO_KsparNa2O_KsparK2O_KsparCr2O3_KsparSample_ID_KsparSi_Kspar_cat_propMg_Kspar_cat_propFet_Kspar_cat_propCa_Kspar_cat_propAl_Kspar_cat_propNa_Kspar_cat_propK_Kspar_cat_propMn_Kspar_cat_propTi_Kspar_cat_propCr_Kspar_cat_propsumSi_Kspar_cat_fracMg_Kspar_cat_fracFet_Kspar_cat_fracCa_Kspar_cat_fracAl_Kspar_cat_fracNa_Kspar_cat_fracK_Kspar_cat_fracMn_Kspar_cat_fracTi_Kspar_cat_fracCr_Kspar_cat_fracAn_KsparAb_KsparOr_KsparID_KsparSiO2_PlagTiO2_PlagAl2O3_PlagFeOt_PlagMnO_PlagMgO_PlagCaO_PlagNa2O_PlagK2O_PlagCr2O3_PlagSample_ID_PlagSi_Plag_cat_propMg_Plag_cat_propFet_Plag_cat_propCa_Plag_cat_propAl_Plag_cat_propNa_Plag_cat_propK_Plag_cat_propMn_Plag_cat_propTi_Plag_cat_propCr_Plag_cat_propsumSi_Plag_cat_fracMg_Plag_cat_fracFet_Plag_cat_fracCa_Plag_cat_fracAl_Plag_cat_fracNa_Plag_cat_fracK_Plag_cat_fracMn_Plag_cat_fracTi_Plag_cat_fracCr_Plag_cat_fracAn_PlagAb_PlagOr_PlagID_PlagK_BarthTP
70.4831750.034127-0.1498590.6450621.0512350.6118000.6109350.5680600.7616581095.06777764.90.018.60.000.00.000.173.0212.10081.0801490.0000000.0000000.0030320.3648450.0974530.2569110.00.00.01.8023900.5992870.0000000.0000000.0016820.2024230.0540690.1425390.00.00.00.0084820.2726740.7188437.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.4924981095.0677775
80.615621-0.015126-0.1978610.6435411.0125410.5225000.6586660.3969200.7203611133.53189266.00.018.90.050.00.000.183.7311.30091.0984570.0000000.0006960.0032100.3707300.1203640.2399250.00.00.01.8333810.5991430.0000000.0003800.0017510.2022110.0656510.1308650.00.00.00.0088300.3311250.6600448.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.598071133.5318925
100.3459560.076966-0.1241670.6458481.0717570.6636430.5688820.7258010.7878101076.24081465.50.018.90.110.00.040.172.5812.500111.0901350.0009920.0015310.0030320.3707300.0832540.2654040.00.00.01.8150780.6005990.0005470.0008440.0016700.2042500.0458680.1462220.00.00.00.0086200.2367260.75465410.057.30.0926.60.4300.038.336.110.49000.953660.0007440.0059850.1485450.5217680.1971640.0104040.00.0011270.01.8393970.5184640.0004050.0032540.0807570.2836630.1071890.0056560.00.0006130.00.4171290.5536560.02921500.4275691076.2408145
180.5125590.037062-0.1322060.6377681.0441040.6277600.6007060.5315450.7599661109.29333864.90.018.60.000.00.000.173.0212.10081.0801490.0000000.0000000.0030320.3648450.0974530.2569110.00.00.01.8023900.5992870.0000000.0000000.0016820.2024230.0540690.1425390.00.00.00.0084820.2726740.7188437.056.50.1226.90.4700.058.955.660.47010.9403450.0012410.0065420.1596010.5276530.1826430.0099790.00.0015020.01.8295060.5139890.0006780.0035760.0872370.2884130.0998320.0054550.00.0008210.00.4531250.5185430.02833210.5258471109.2933385
190.635801-0.013936-0.1861240.6350491.0081420.5317890.6489850.3723400.7179131149.13269966.00.018.90.050.00.000.183.7311.30091.0984570.0000000.0006960.0032100.3707300.1203640.2399250.00.00.01.8333810.5991430.0000000.0003800.0017510.2022110.0656510.1308650.00.00.00.0088300.3311250.6600448.056.50.1226.90.4700.058.955.660.47010.9403450.0012410.0065420.1596010.5276530.1826430.0099790.00.0015020.01.8295060.5139890.0006780.0035760.0872370.2884130.0998320.0054550.00.0008210.00.4531250.5185430.02833210.6385681149.1326995
210.3850290.080548-0.1026970.6391751.0631410.6838120.5586270.6781110.7865091089.81675165.50.018.90.110.00.040.172.5812.500111.0901350.0009920.0015310.0030320.3707300.0832540.2654040.00.00.01.8150780.6005990.0005470.0008440.0016700.2042500.0458680.1462220.00.00.00.0086200.2367260.75465410.056.50.1226.90.4700.058.955.660.47010.9403450.0012410.0065420.1596010.5276530.1826430.0099790.00.0015020.01.8295060.5139890.0006780.0035760.0872370.2884130.0998320.0054550.00.0008210.00.4531250.5185430.02833210.4565211089.8167515
330.342311-0.030924-0.1412470.5995130.9202770.5034330.6304370.5779660.6446801287.49642065.50.019.60.070.00.000.754.819.36011.0901350.0000000.0009740.0133740.3844610.1552140.1987350.00.00.01.8428930.5915350.0000000.0005290.0072570.2086180.0842230.1078380.00.00.00.0364100.4225550.5410350.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.8422971287.4964205
34-0.0166920.078517-0.0976870.6023350.9738610.6403150.5238170.9905530.7380021213.05822165.40.019.40.050.00.000.593.1311.50021.0884710.0000000.0006960.0105210.3805380.1010020.2441720.00.00.01.8253990.5962920.0000000.0003810.0057640.2084680.0553320.1337630.00.00.00.0295790.2839570.6864641.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.5660241213.0582215
370.1355720.0799690.0108470.6046261.0194480.7777280.5246570.8838760.7668821158.89579465.10.019.20.050.00.000.362.8712.60051.0834780.0000000.0006960.0064200.3766150.0926120.2675270.00.00.01.8273470.5929240.0000000.0003810.0035130.2060990.0506810.1464020.00.00.00.0175130.2526530.7298344.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.5036231158.8957945
38-0.0066810.1130730.0441220.6054231.0357640.8320180.4923501.0424450.7878961141.21107765.00.019.20.000.00.000.352.5112.80061.0818130.0000000.0000000.0062410.3766150.0809950.2717740.00.00.01.8174380.5952410.0000000.0000000.0034340.2072230.0445660.1495370.00.00.00.0173850.2256070.7570085.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.4497121141.2110775
39-0.0433940.1238340.1518490.6070991.0708420.9585240.4832651.1142350.8066761105.81391164.90.019.30.060.00.030.292.3013.30071.0801490.0007440.0008350.0051710.3785760.0742190.2823900.00.00.01.8220850.5928090.0004090.0004580.0028380.2077710.0407330.1549820.00.00.00.0142940.2051490.7805576.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.4089321105.8139115
400.5325500.0047820.1932270.6070361.0695050.9534510.6022540.5369550.7602241107.10155964.90.018.60.000.00.000.173.0212.10081.0801490.0000000.0000000.0030320.3648450.0974530.2569110.00.00.01.8023900.5992870.0000000.0000000.0016820.2024230.0540690.1425390.00.00.00.0084820.2726740.7188437.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.5435341107.1015595
410.654607-0.0452790.0962220.6051721.0305920.8145070.6504510.3759860.7182851146.72738366.00.018.90.050.00.000.183.7311.30091.0984570.0000000.0006960.0032100.3707300.1203640.2399250.00.00.01.8333810.5991430.0000000.0003800.0017510.2022110.0656510.1308650.00.00.00.0088300.3311250.6600448.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.6600461146.7273835
420.727768-0.0841430.0522770.6039861.0065190.7366680.6881290.2787520.6843911173.51372666.00.019.20.000.00.000.184.3110.400101.0984570.0000000.0000000.0032100.3766150.1390800.2208160.00.00.01.8381770.5975790.0000000.0000000.0017460.2048850.0756620.1201280.00.00.00.0088400.3830280.6081329.057.20.1627.00.6200.069.035.580.84030.9519960.0014890.008630.1610270.5296140.1800610.0178350.00.0020030.01.8526550.5138550.0008040.0046580.0869170.2858680.0971910.0096270.00.0010810.00.448640.501670.04969130.7635061173.5137265
510.5231650.038410-0.1141090.6346691.0393920.6451130.5962580.5162270.7592221115.67233464.90.018.60.000.00.000.173.0212.10081.0801490.0000000.0000000.0030320.3648450.0974530.2569110.00.00.01.8023900.5992870.0000000.0000000.0016820.2024230.0540690.1425390.00.00.00.0084820.2726740.7188437.056.70.1427.60.6900.119.465.580.48040.9436740.0027290.0096040.1686950.5413830.1800610.0101920.00.0017530.01.8580920.5078730.0014690.0051690.090790.2913650.0969070.0054850.00.0009430.00.4699710.5016360.02839340.543571115.6723345
520.642779-0.013464-0.1723840.6313041.0047890.5444540.6447690.3620090.7168381156.13659766.00.018.90.050.00.000.183.7311.30091.0984570.0000000.0006960.0032100.3707300.1203640.2399250.00.00.01.8333810.5991430.0000000.0003800.0017510.2022110.0656510.1308650.00.00.00.0088300.3311250.6600448.056.70.1427.60.6900.119.465.580.48040.9436740.0027290.0096040.1686950.5413830.1800610.0101920.00.0017530.01.8580920.5078730.0014690.0051690.090790.2913650.0969070.0054850.00.0009430.00.4699710.5016360.02839340.6600911156.1365975
530.714493-0.053705-0.1938720.6291680.9833250.4886240.6828730.2688320.6824971183.51211566.00.019.20.000.00.000.184.3110.400101.0984570.0000000.0000000.0032100.3766150.1390800.2208160.00.00.01.8381770.5975790.0000000.0000000.0017460.2048850.0756620.1201280.00.00.00.0088400.3830280.6081329.056.70.1427.60.6900.119.465.580.48040.9436740.0027290.0096040.1686950.5413830.1800610.0101920.00.0017530.01.8580920.5078730.0014690.0051690.090790.2913650.0969070.0054850.00.0009430.00.4699710.5016360.02839340.7635581183.5121155
540.3995640.082237-0.0818720.6364101.0576920.7040650.5541730.6581280.7859371095.90095065.50.018.90.110.00.040.172.5812.500111.0901350.0009920.0015310.0030320.3707300.0832540.2654040.00.00.01.8150780.6005990.0005470.0008440.0016700.2042500.0458680.1462220.00.00.00.0086200.2367260.75465410.056.70.1427.60.6900.119.465.580.48040.9436740.0027290.0096040.1686950.5413830.1800610.0101920.00.0017530.01.8580920.5078730.0014690.0051690.090790.2913650.0969070.0054850.00.0009430.00.4699710.5016360.02839340.4719081095.9009505
620.5112630.039141-0.1564120.6398211.0427180.6035500.6006800.5314540.7599621109.33042964.90.018.60.000.00.000.173.0212.10081.0801490.0000000.0000000.0030320.3648450.0974530.2569110.00.00.01.8023900.5992870.0000000.0000000.0016820.2024230.0540690.1425390.00.00.00.0084820.2726740.7188437.056.00.2126.60.5600.18.995.70.45050.9320240.0024810.0077940.1603140.5217680.1839340.0095550.00.0026290.01.8204990.5119610.0013630.0042810.0880610.2866070.1010350.0052480.00.0014440.00.4531180.5198770.02700550.5244981109.3304295
650.3836810.082652-0.1289060.6412531.0616740.6576000.5586010.6779930.7865061089.85213465.50.018.90.110.00.040.172.5812.500111.0901350.0009920.0015310.0030320.3707300.0832540.2654040.00.00.01.8150780.6005990.0005470.0008440.0016700.2042500.0458680.1462220.00.00.00.0086200.2367260.75465410.056.00.2126.60.5600.18.995.70.45050.9320240.0024810.0077940.1603140.5217680.1839340.0095550.00.0026290.01.8204990.5119610.0013630.0042810.0880610.2866070.1010350.0052480.00.0014440.00.4531180.5198770.02700550.455351089.8521345
760.4069450.094842-0.1951440.6444931.0452310.5902080.5496510.6382860.7853531102.19822165.50.018.90.110.00.040.172.5812.500111.0901350.0009920.0015310.0030320.3707300.0832540.2654040.00.00.01.8150780.6005990.0005470.0008440.0016700.2042500.0458680.1462220.00.00.00.0086200.2367260.75465410.056.10.2127.80.5600.099.945.530.38060.9336880.0022330.0077940.1772550.5453070.1784480.0080680.00.0026290.01.8554220.5032210.0012040.0042010.0955340.2938990.0961760.0043480.00.0014170.00.4872710.490550.0221860.4825731102.1982215
\n", "
" ], "text/plain": [ " Delta_An Delta_Ab Delta_Or a_Ab_plg a_An_plg a_Or_plg a_Ab_kspar \\\n", "7 0.483175 0.034127 -0.149859 0.645062 1.051235 0.611800 0.610935 \n", "8 0.615621 -0.015126 -0.197861 0.643541 1.012541 0.522500 0.658666 \n", "10 0.345956 0.076966 -0.124167 0.645848 1.071757 0.663643 0.568882 \n", "18 0.512559 0.037062 -0.132206 0.637768 1.044104 0.627760 0.600706 \n", "19 0.635801 -0.013936 -0.186124 0.635049 1.008142 0.531789 0.648985 \n", "21 0.385029 0.080548 -0.102697 0.639175 1.063141 0.683812 0.558627 \n", "33 0.342311 -0.030924 -0.141247 0.599513 0.920277 0.503433 0.630437 \n", "34 -0.016692 0.078517 -0.097687 0.602335 0.973861 0.640315 0.523817 \n", "37 0.135572 0.079969 0.010847 0.604626 1.019448 0.777728 0.524657 \n", "38 -0.006681 0.113073 0.044122 0.605423 1.035764 0.832018 0.492350 \n", "39 -0.043394 0.123834 0.151849 0.607099 1.070842 0.958524 0.483265 \n", "40 0.532550 0.004782 0.193227 0.607036 1.069505 0.953451 0.602254 \n", "41 0.654607 -0.045279 0.096222 0.605172 1.030592 0.814507 0.650451 \n", "42 0.727768 -0.084143 0.052277 0.603986 1.006519 0.736668 0.688129 \n", "51 0.523165 0.038410 -0.114109 0.634669 1.039392 0.645113 0.596258 \n", "52 0.642779 -0.013464 -0.172384 0.631304 1.004789 0.544454 0.644769 \n", "53 0.714493 -0.053705 -0.193872 0.629168 0.983325 0.488624 0.682873 \n", "54 0.399564 0.082237 -0.081872 0.636410 1.057692 0.704065 0.554173 \n", "62 0.511263 0.039141 -0.156412 0.639821 1.042718 0.603550 0.600680 \n", "65 0.383681 0.082652 -0.128906 0.641253 1.061674 0.657600 0.558601 \n", "76 0.406945 0.094842 -0.195144 0.644493 1.045231 0.590208 0.549651 \n", "\n", " a_An_kspar a_Or_kspar T_K_calc SiO2_Kspar TiO2_Kspar Al2O3_Kspar \\\n", "7 0.568060 0.761658 1095.067777 64.9 0.0 18.6 \n", "8 0.396920 0.720361 1133.531892 66.0 0.0 18.9 \n", "10 0.725801 0.787810 1076.240814 65.5 0.0 18.9 \n", "18 0.531545 0.759966 1109.293338 64.9 0.0 18.6 \n", "19 0.372340 0.717913 1149.132699 66.0 0.0 18.9 \n", "21 0.678111 0.786509 1089.816751 65.5 0.0 18.9 \n", "33 0.577966 0.644680 1287.496420 65.5 0.0 19.6 \n", "34 0.990553 0.738002 1213.058221 65.4 0.0 19.4 \n", "37 0.883876 0.766882 1158.895794 65.1 0.0 19.2 \n", "38 1.042445 0.787896 1141.211077 65.0 0.0 19.2 \n", "39 1.114235 0.806676 1105.813911 64.9 0.0 19.3 \n", "40 0.536955 0.760224 1107.101559 64.9 0.0 18.6 \n", "41 0.375986 0.718285 1146.727383 66.0 0.0 18.9 \n", "42 0.278752 0.684391 1173.513726 66.0 0.0 19.2 \n", "51 0.516227 0.759222 1115.672334 64.9 0.0 18.6 \n", "52 0.362009 0.716838 1156.136597 66.0 0.0 18.9 \n", "53 0.268832 0.682497 1183.512115 66.0 0.0 19.2 \n", "54 0.658128 0.785937 1095.900950 65.5 0.0 18.9 \n", "62 0.531454 0.759962 1109.330429 64.9 0.0 18.6 \n", "65 0.677993 0.786506 1089.852134 65.5 0.0 18.9 \n", "76 0.638286 0.785353 1102.198221 65.5 0.0 18.9 \n", "\n", " FeOt_Kspar MnO_Kspar MgO_Kspar CaO_Kspar Na2O_Kspar K2O_Kspar \\\n", "7 0.00 0.0 0.00 0.17 3.02 12.10 \n", "8 0.05 0.0 0.00 0.18 3.73 11.30 \n", "10 0.11 0.0 0.04 0.17 2.58 12.50 \n", "18 0.00 0.0 0.00 0.17 3.02 12.10 \n", "19 0.05 0.0 0.00 0.18 3.73 11.30 \n", "21 0.11 0.0 0.04 0.17 2.58 12.50 \n", "33 0.07 0.0 0.00 0.75 4.81 9.36 \n", "34 0.05 0.0 0.00 0.59 3.13 11.50 \n", "37 0.05 0.0 0.00 0.36 2.87 12.60 \n", "38 0.00 0.0 0.00 0.35 2.51 12.80 \n", "39 0.06 0.0 0.03 0.29 2.30 13.30 \n", "40 0.00 0.0 0.00 0.17 3.02 12.10 \n", "41 0.05 0.0 0.00 0.18 3.73 11.30 \n", "42 0.00 0.0 0.00 0.18 4.31 10.40 \n", "51 0.00 0.0 0.00 0.17 3.02 12.10 \n", "52 0.05 0.0 0.00 0.18 3.73 11.30 \n", "53 0.00 0.0 0.00 0.18 4.31 10.40 \n", "54 0.11 0.0 0.04 0.17 2.58 12.50 \n", "62 0.00 0.0 0.00 0.17 3.02 12.10 \n", "65 0.11 0.0 0.04 0.17 2.58 12.50 \n", "76 0.11 0.0 0.04 0.17 2.58 12.50 \n", "\n", " Cr2O3_Kspar Sample_ID_Kspar Si_Kspar_cat_prop Mg_Kspar_cat_prop \\\n", "7 0 8 1.080149 0.000000 \n", "8 0 9 1.098457 0.000000 \n", "10 0 11 1.090135 0.000992 \n", "18 0 8 1.080149 0.000000 \n", "19 0 9 1.098457 0.000000 \n", "21 0 11 1.090135 0.000992 \n", "33 0 1 1.090135 0.000000 \n", "34 0 2 1.088471 0.000000 \n", "37 0 5 1.083478 0.000000 \n", "38 0 6 1.081813 0.000000 \n", "39 0 7 1.080149 0.000744 \n", "40 0 8 1.080149 0.000000 \n", "41 0 9 1.098457 0.000000 \n", "42 0 10 1.098457 0.000000 \n", "51 0 8 1.080149 0.000000 \n", "52 0 9 1.098457 0.000000 \n", "53 0 10 1.098457 0.000000 \n", "54 0 11 1.090135 0.000992 \n", "62 0 8 1.080149 0.000000 \n", "65 0 11 1.090135 0.000992 \n", "76 0 11 1.090135 0.000992 \n", "\n", " Fet_Kspar_cat_prop Ca_Kspar_cat_prop Al_Kspar_cat_prop \\\n", "7 0.000000 0.003032 0.364845 \n", "8 0.000696 0.003210 0.370730 \n", "10 0.001531 0.003032 0.370730 \n", "18 0.000000 0.003032 0.364845 \n", "19 0.000696 0.003210 0.370730 \n", "21 0.001531 0.003032 0.370730 \n", "33 0.000974 0.013374 0.384461 \n", "34 0.000696 0.010521 0.380538 \n", "37 0.000696 0.006420 0.376615 \n", "38 0.000000 0.006241 0.376615 \n", "39 0.000835 0.005171 0.378576 \n", "40 0.000000 0.003032 0.364845 \n", "41 0.000696 0.003210 0.370730 \n", "42 0.000000 0.003210 0.376615 \n", "51 0.000000 0.003032 0.364845 \n", "52 0.000696 0.003210 0.370730 \n", "53 0.000000 0.003210 0.376615 \n", "54 0.001531 0.003032 0.370730 \n", "62 0.000000 0.003032 0.364845 \n", "65 0.001531 0.003032 0.370730 \n", "76 0.001531 0.003032 0.370730 \n", "\n", " Na_Kspar_cat_prop K_Kspar_cat_prop Mn_Kspar_cat_prop Ti_Kspar_cat_prop \\\n", "7 0.097453 0.256911 0.0 0.0 \n", "8 0.120364 0.239925 0.0 0.0 \n", "10 0.083254 0.265404 0.0 0.0 \n", "18 0.097453 0.256911 0.0 0.0 \n", "19 0.120364 0.239925 0.0 0.0 \n", "21 0.083254 0.265404 0.0 0.0 \n", "33 0.155214 0.198735 0.0 0.0 \n", "34 0.101002 0.244172 0.0 0.0 \n", "37 0.092612 0.267527 0.0 0.0 \n", "38 0.080995 0.271774 0.0 0.0 \n", "39 0.074219 0.282390 0.0 0.0 \n", "40 0.097453 0.256911 0.0 0.0 \n", "41 0.120364 0.239925 0.0 0.0 \n", "42 0.139080 0.220816 0.0 0.0 \n", "51 0.097453 0.256911 0.0 0.0 \n", "52 0.120364 0.239925 0.0 0.0 \n", "53 0.139080 0.220816 0.0 0.0 \n", "54 0.083254 0.265404 0.0 0.0 \n", "62 0.097453 0.256911 0.0 0.0 \n", "65 0.083254 0.265404 0.0 0.0 \n", "76 0.083254 0.265404 0.0 0.0 \n", "\n", " Cr_Kspar_cat_prop sum Si_Kspar_cat_frac Mg_Kspar_cat_frac \\\n", "7 0.0 1.802390 0.599287 0.000000 \n", "8 0.0 1.833381 0.599143 0.000000 \n", "10 0.0 1.815078 0.600599 0.000547 \n", "18 0.0 1.802390 0.599287 0.000000 \n", "19 0.0 1.833381 0.599143 0.000000 \n", "21 0.0 1.815078 0.600599 0.000547 \n", "33 0.0 1.842893 0.591535 0.000000 \n", "34 0.0 1.825399 0.596292 0.000000 \n", "37 0.0 1.827347 0.592924 0.000000 \n", "38 0.0 1.817438 0.595241 0.000000 \n", "39 0.0 1.822085 0.592809 0.000409 \n", "40 0.0 1.802390 0.599287 0.000000 \n", "41 0.0 1.833381 0.599143 0.000000 \n", "42 0.0 1.838177 0.597579 0.000000 \n", "51 0.0 1.802390 0.599287 0.000000 \n", "52 0.0 1.833381 0.599143 0.000000 \n", "53 0.0 1.838177 0.597579 0.000000 \n", "54 0.0 1.815078 0.600599 0.000547 \n", "62 0.0 1.802390 0.599287 0.000000 \n", "65 0.0 1.815078 0.600599 0.000547 \n", "76 0.0 1.815078 0.600599 0.000547 \n", "\n", " Fet_Kspar_cat_frac Ca_Kspar_cat_frac Al_Kspar_cat_frac \\\n", "7 0.000000 0.001682 0.202423 \n", "8 0.000380 0.001751 0.202211 \n", "10 0.000844 0.001670 0.204250 \n", "18 0.000000 0.001682 0.202423 \n", "19 0.000380 0.001751 0.202211 \n", "21 0.000844 0.001670 0.204250 \n", "33 0.000529 0.007257 0.208618 \n", "34 0.000381 0.005764 0.208468 \n", "37 0.000381 0.003513 0.206099 \n", "38 0.000000 0.003434 0.207223 \n", "39 0.000458 0.002838 0.207771 \n", "40 0.000000 0.001682 0.202423 \n", "41 0.000380 0.001751 0.202211 \n", "42 0.000000 0.001746 0.204885 \n", "51 0.000000 0.001682 0.202423 \n", "52 0.000380 0.001751 0.202211 \n", "53 0.000000 0.001746 0.204885 \n", "54 0.000844 0.001670 0.204250 \n", "62 0.000000 0.001682 0.202423 \n", "65 0.000844 0.001670 0.204250 \n", "76 0.000844 0.001670 0.204250 \n", "\n", " Na_Kspar_cat_frac K_Kspar_cat_frac Mn_Kspar_cat_frac Ti_Kspar_cat_frac \\\n", "7 0.054069 0.142539 0.0 0.0 \n", "8 0.065651 0.130865 0.0 0.0 \n", "10 0.045868 0.146222 0.0 0.0 \n", "18 0.054069 0.142539 0.0 0.0 \n", "19 0.065651 0.130865 0.0 0.0 \n", "21 0.045868 0.146222 0.0 0.0 \n", "33 0.084223 0.107838 0.0 0.0 \n", "34 0.055332 0.133763 0.0 0.0 \n", "37 0.050681 0.146402 0.0 0.0 \n", "38 0.044566 0.149537 0.0 0.0 \n", "39 0.040733 0.154982 0.0 0.0 \n", "40 0.054069 0.142539 0.0 0.0 \n", "41 0.065651 0.130865 0.0 0.0 \n", "42 0.075662 0.120128 0.0 0.0 \n", "51 0.054069 0.142539 0.0 0.0 \n", "52 0.065651 0.130865 0.0 0.0 \n", "53 0.075662 0.120128 0.0 0.0 \n", "54 0.045868 0.146222 0.0 0.0 \n", "62 0.054069 0.142539 0.0 0.0 \n", "65 0.045868 0.146222 0.0 0.0 \n", "76 0.045868 0.146222 0.0 0.0 \n", "\n", " Cr_Kspar_cat_frac An_Kspar Ab_Kspar Or_Kspar ID_Kspar SiO2_Plag \\\n", "7 0.0 0.008482 0.272674 0.718843 7.0 57.3 \n", "8 0.0 0.008830 0.331125 0.660044 8.0 57.3 \n", "10 0.0 0.008620 0.236726 0.754654 10.0 57.3 \n", "18 0.0 0.008482 0.272674 0.718843 7.0 56.5 \n", "19 0.0 0.008830 0.331125 0.660044 8.0 56.5 \n", "21 0.0 0.008620 0.236726 0.754654 10.0 56.5 \n", "33 0.0 0.036410 0.422555 0.541035 0.0 57.2 \n", "34 0.0 0.029579 0.283957 0.686464 1.0 57.2 \n", "37 0.0 0.017513 0.252653 0.729834 4.0 57.2 \n", "38 0.0 0.017385 0.225607 0.757008 5.0 57.2 \n", "39 0.0 0.014294 0.205149 0.780557 6.0 57.2 \n", "40 0.0 0.008482 0.272674 0.718843 7.0 57.2 \n", "41 0.0 0.008830 0.331125 0.660044 8.0 57.2 \n", "42 0.0 0.008840 0.383028 0.608132 9.0 57.2 \n", "51 0.0 0.008482 0.272674 0.718843 7.0 56.7 \n", "52 0.0 0.008830 0.331125 0.660044 8.0 56.7 \n", "53 0.0 0.008840 0.383028 0.608132 9.0 56.7 \n", "54 0.0 0.008620 0.236726 0.754654 10.0 56.7 \n", "62 0.0 0.008482 0.272674 0.718843 7.0 56.0 \n", "65 0.0 0.008620 0.236726 0.754654 10.0 56.0 \n", "76 0.0 0.008620 0.236726 0.754654 10.0 56.1 \n", "\n", " TiO2_Plag Al2O3_Plag FeOt_Plag MnO_Plag MgO_Plag CaO_Plag Na2O_Plag \\\n", "7 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", "8 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", "10 0.09 26.6 0.43 0 0.03 8.33 6.11 \n", "18 0.12 26.9 0.47 0 0.05 8.95 5.66 \n", "19 0.12 26.9 0.47 0 0.05 8.95 5.66 \n", "21 0.12 26.9 0.47 0 0.05 8.95 5.66 \n", "33 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "34 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "37 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "38 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "39 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "40 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "41 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "42 0.16 27.0 0.62 0 0.06 9.03 5.58 \n", "51 0.14 27.6 0.69 0 0.11 9.46 5.58 \n", "52 0.14 27.6 0.69 0 0.11 9.46 5.58 \n", "53 0.14 27.6 0.69 0 0.11 9.46 5.58 \n", "54 0.14 27.6 0.69 0 0.11 9.46 5.58 \n", "62 0.21 26.6 0.56 0 0.1 8.99 5.7 \n", "65 0.21 26.6 0.56 0 0.1 8.99 5.7 \n", "76 0.21 27.8 0.56 0 0.09 9.94 5.53 \n", "\n", " K2O_Plag Cr2O3_Plag Sample_ID_Plag Si_Plag_cat_prop Mg_Plag_cat_prop \\\n", "7 0.49 0 0 0.95366 0.000744 \n", "8 0.49 0 0 0.95366 0.000744 \n", "10 0.49 0 0 0.95366 0.000744 \n", "18 0.47 0 1 0.940345 0.001241 \n", "19 0.47 0 1 0.940345 0.001241 \n", "21 0.47 0 1 0.940345 0.001241 \n", "33 0.84 0 3 0.951996 0.001489 \n", "34 0.84 0 3 0.951996 0.001489 \n", "37 0.84 0 3 0.951996 0.001489 \n", "38 0.84 0 3 0.951996 0.001489 \n", "39 0.84 0 3 0.951996 0.001489 \n", "40 0.84 0 3 0.951996 0.001489 \n", "41 0.84 0 3 0.951996 0.001489 \n", "42 0.84 0 3 0.951996 0.001489 \n", "51 0.48 0 4 0.943674 0.002729 \n", "52 0.48 0 4 0.943674 0.002729 \n", "53 0.48 0 4 0.943674 0.002729 \n", "54 0.48 0 4 0.943674 0.002729 \n", "62 0.45 0 5 0.932024 0.002481 \n", "65 0.45 0 5 0.932024 0.002481 \n", "76 0.38 0 6 0.933688 0.002233 \n", "\n", " Fet_Plag_cat_prop Ca_Plag_cat_prop Al_Plag_cat_prop Na_Plag_cat_prop \\\n", "7 0.005985 0.148545 0.521768 0.197164 \n", "8 0.005985 0.148545 0.521768 0.197164 \n", "10 0.005985 0.148545 0.521768 0.197164 \n", "18 0.006542 0.159601 0.527653 0.182643 \n", "19 0.006542 0.159601 0.527653 0.182643 \n", "21 0.006542 0.159601 0.527653 0.182643 \n", "33 0.00863 0.161027 0.529614 0.180061 \n", "34 0.00863 0.161027 0.529614 0.180061 \n", "37 0.00863 0.161027 0.529614 0.180061 \n", "38 0.00863 0.161027 0.529614 0.180061 \n", "39 0.00863 0.161027 0.529614 0.180061 \n", "40 0.00863 0.161027 0.529614 0.180061 \n", "41 0.00863 0.161027 0.529614 0.180061 \n", "42 0.00863 0.161027 0.529614 0.180061 \n", "51 0.009604 0.168695 0.541383 0.180061 \n", "52 0.009604 0.168695 0.541383 0.180061 \n", "53 0.009604 0.168695 0.541383 0.180061 \n", "54 0.009604 0.168695 0.541383 0.180061 \n", "62 0.007794 0.160314 0.521768 0.183934 \n", "65 0.007794 0.160314 0.521768 0.183934 \n", "76 0.007794 0.177255 0.545307 0.178448 \n", "\n", " K_Plag_cat_prop Mn_Plag_cat_prop Ti_Plag_cat_prop Cr_Plag_cat_prop \\\n", "7 0.010404 0.0 0.001127 0.0 \n", "8 0.010404 0.0 0.001127 0.0 \n", "10 0.010404 0.0 0.001127 0.0 \n", "18 0.009979 0.0 0.001502 0.0 \n", "19 0.009979 0.0 0.001502 0.0 \n", "21 0.009979 0.0 0.001502 0.0 \n", "33 0.017835 0.0 0.002003 0.0 \n", "34 0.017835 0.0 0.002003 0.0 \n", "37 0.017835 0.0 0.002003 0.0 \n", "38 0.017835 0.0 0.002003 0.0 \n", "39 0.017835 0.0 0.002003 0.0 \n", "40 0.017835 0.0 0.002003 0.0 \n", "41 0.017835 0.0 0.002003 0.0 \n", "42 0.017835 0.0 0.002003 0.0 \n", "51 0.010192 0.0 0.001753 0.0 \n", "52 0.010192 0.0 0.001753 0.0 \n", "53 0.010192 0.0 0.001753 0.0 \n", "54 0.010192 0.0 0.001753 0.0 \n", "62 0.009555 0.0 0.002629 0.0 \n", "65 0.009555 0.0 0.002629 0.0 \n", "76 0.008068 0.0 0.002629 0.0 \n", "\n", " sum Si_Plag_cat_frac Mg_Plag_cat_frac Fet_Plag_cat_frac \\\n", "7 1.839397 0.518464 0.000405 0.003254 \n", "8 1.839397 0.518464 0.000405 0.003254 \n", "10 1.839397 0.518464 0.000405 0.003254 \n", "18 1.829506 0.513989 0.000678 0.003576 \n", "19 1.829506 0.513989 0.000678 0.003576 \n", "21 1.829506 0.513989 0.000678 0.003576 \n", "33 1.852655 0.513855 0.000804 0.004658 \n", "34 1.852655 0.513855 0.000804 0.004658 \n", "37 1.852655 0.513855 0.000804 0.004658 \n", "38 1.852655 0.513855 0.000804 0.004658 \n", "39 1.852655 0.513855 0.000804 0.004658 \n", "40 1.852655 0.513855 0.000804 0.004658 \n", "41 1.852655 0.513855 0.000804 0.004658 \n", "42 1.852655 0.513855 0.000804 0.004658 \n", "51 1.858092 0.507873 0.001469 0.005169 \n", "52 1.858092 0.507873 0.001469 0.005169 \n", "53 1.858092 0.507873 0.001469 0.005169 \n", "54 1.858092 0.507873 0.001469 0.005169 \n", "62 1.820499 0.511961 0.001363 0.004281 \n", "65 1.820499 0.511961 0.001363 0.004281 \n", "76 1.855422 0.503221 0.001204 0.004201 \n", "\n", " Ca_Plag_cat_frac Al_Plag_cat_frac Na_Plag_cat_frac K_Plag_cat_frac \\\n", "7 0.080757 0.283663 0.107189 0.005656 \n", "8 0.080757 0.283663 0.107189 0.005656 \n", "10 0.080757 0.283663 0.107189 0.005656 \n", "18 0.087237 0.288413 0.099832 0.005455 \n", "19 0.087237 0.288413 0.099832 0.005455 \n", "21 0.087237 0.288413 0.099832 0.005455 \n", "33 0.086917 0.285868 0.097191 0.009627 \n", "34 0.086917 0.285868 0.097191 0.009627 \n", "37 0.086917 0.285868 0.097191 0.009627 \n", "38 0.086917 0.285868 0.097191 0.009627 \n", "39 0.086917 0.285868 0.097191 0.009627 \n", "40 0.086917 0.285868 0.097191 0.009627 \n", "41 0.086917 0.285868 0.097191 0.009627 \n", "42 0.086917 0.285868 0.097191 0.009627 \n", "51 0.09079 0.291365 0.096907 0.005485 \n", "52 0.09079 0.291365 0.096907 0.005485 \n", "53 0.09079 0.291365 0.096907 0.005485 \n", "54 0.09079 0.291365 0.096907 0.005485 \n", "62 0.088061 0.286607 0.101035 0.005248 \n", "65 0.088061 0.286607 0.101035 0.005248 \n", "76 0.095534 0.293899 0.096176 0.004348 \n", "\n", " Mn_Plag_cat_frac Ti_Plag_cat_frac Cr_Plag_cat_frac An_Plag Ab_Plag \\\n", "7 0.0 0.000613 0.0 0.417129 0.553656 \n", "8 0.0 0.000613 0.0 0.417129 0.553656 \n", "10 0.0 0.000613 0.0 0.417129 0.553656 \n", "18 0.0 0.000821 0.0 0.453125 0.518543 \n", "19 0.0 0.000821 0.0 0.453125 0.518543 \n", "21 0.0 0.000821 0.0 0.453125 0.518543 \n", "33 0.0 0.001081 0.0 0.44864 0.50167 \n", "34 0.0 0.001081 0.0 0.44864 0.50167 \n", "37 0.0 0.001081 0.0 0.44864 0.50167 \n", "38 0.0 0.001081 0.0 0.44864 0.50167 \n", "39 0.0 0.001081 0.0 0.44864 0.50167 \n", "40 0.0 0.001081 0.0 0.44864 0.50167 \n", "41 0.0 0.001081 0.0 0.44864 0.50167 \n", "42 0.0 0.001081 0.0 0.44864 0.50167 \n", "51 0.0 0.000943 0.0 0.469971 0.501636 \n", "52 0.0 0.000943 0.0 0.469971 0.501636 \n", "53 0.0 0.000943 0.0 0.469971 0.501636 \n", "54 0.0 0.000943 0.0 0.469971 0.501636 \n", "62 0.0 0.001444 0.0 0.453118 0.519877 \n", "65 0.0 0.001444 0.0 0.453118 0.519877 \n", "76 0.0 0.001417 0.0 0.487271 0.49055 \n", "\n", " Or_Plag ID_Plag K_Barth T P \n", "7 0.029215 0 0.492498 1095.067777 5 \n", "8 0.029215 0 0.59807 1133.531892 5 \n", "10 0.029215 0 0.427569 1076.240814 5 \n", "18 0.028332 1 0.525847 1109.293338 5 \n", "19 0.028332 1 0.638568 1149.132699 5 \n", "21 0.028332 1 0.456521 1089.816751 5 \n", "33 0.049691 3 0.842297 1287.496420 5 \n", "34 0.049691 3 0.566024 1213.058221 5 \n", "37 0.049691 3 0.503623 1158.895794 5 \n", "38 0.049691 3 0.449712 1141.211077 5 \n", "39 0.049691 3 0.408932 1105.813911 5 \n", "40 0.049691 3 0.543534 1107.101559 5 \n", "41 0.049691 3 0.660046 1146.727383 5 \n", "42 0.049691 3 0.763506 1173.513726 5 \n", "51 0.028393 4 0.54357 1115.672334 5 \n", "52 0.028393 4 0.660091 1156.136597 5 \n", "53 0.028393 4 0.763558 1183.512115 5 \n", "54 0.028393 4 0.471908 1095.900950 5 \n", "62 0.027005 5 0.524498 1109.330429 5 \n", "65 0.027005 5 0.45535 1089.852134 5 \n", "76 0.02218 6 0.482573 1102.198221 5 " ] }, "execution_count": 23, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# This filters pairs to only those meeting these three criteria defined above\n", "Filt_pairs=T_27a_5kbar_All.loc[filt_a_An&filt_a_Ab&filt_a_Or] \n", "Filt_pairs" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.12" } }, "nbformat": 4, "nbformat_minor": 4 }