This page was generated from docs/Examples/Feldspar_Thermobarometry/Fspar_Liq_Matching.ipynb. Interactive online version: Binder badge.

Python Notebook Download

Plag-Liq and Kspar-Liq Matching

[1]:
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import Thermobar as pt
pd.options.display.max_columns = None

Example 1 - Plag and Liquid.

[2]:
## Load in the liquids
out_Liq=pt.import_excel('Feldspar_Liquid.xlsx', sheet_name="Liq_only")
Liqs=out_Liq['Liqs']
# Load in the Plags
out_Plag=pt.import_excel('Feldspar_Liquid.xlsx', sheet_name="Plag_only")
Plags=out_Plag['Plags']
[3]:
## Inspect them to check they read in right
display(Plags.head())
display(Liqs.head())
SiO2_Plag TiO2_Plag Al2O3_Plag FeOt_Plag MnO_Plag MgO_Plag CaO_Plag Na2O_Plag K2O_Plag Cr2O3_Plag Sample_ID_Plag
0 57.3 0.09 26.6 0.43 0.0 0.03 8.33 6.11 0.49 0.0 Plag1
1 56.5 0.12 26.9 0.47 0.0 0.05 8.95 5.66 0.47 0.0 Plag2
2 57.6 0.11 26.3 0.50 0.0 0.07 8.50 6.27 0.40 0.0 Plag3_core
3 57.2 0.16 27.0 0.62 0.0 0.06 9.03 5.58 0.84 0.0 Plag4
4 56.7 0.14 27.6 0.69 0.0 0.11 9.46 5.58 0.48 0.0 Plag5
SiO2_Liq TiO2_Liq Al2O3_Liq FeOt_Liq MnO_Liq MgO_Liq CaO_Liq Na2O_Liq K2O_Liq Cr2O3_Liq P2O5_Liq H2O_Liq Fe3Fet_Liq NiO_Liq CoO_Liq CO2_Liq Sample_ID_Liq
0 49.1 3.22 14.4 14.8 0.14 3.20 6.72 3.34 1.70 0.0 1.13 0.0 0.0 0.0 0.0 0.0 Glass1
1 49.2 3.89 15.3 13.7 0.12 3.88 6.76 3.44 1.22 0.0 0.83 0.0 0.0 0.0 0.0 0.0 Glass2
2 49.6 3.79 15.8 13.0 0.14 4.26 6.59 3.65 1.04 0.0 0.63 0.0 0.0 0.0 0.0 0.0 Glass3
3 47.1 4.21 12.0 17.8 0.18 3.40 7.28 2.93 2.02 0.0 2.32 0.0 0.0 0.0 0.0 0.0 Glass4
4 48.1 3.88 13.2 16.4 0.16 4.02 6.51 3.36 1.36 0.0 1.59 0.0 0.0 0.0 0.0 0.0 Glass5

Example 1a - All possible Liq-Plag matches, An-Ab equilibrium

  • You can specify to use the An-Ab equilibrium test of Putirka (2008) as in this example

[4]:
MM_dict=pt.calculate_fspar_liq_temp_matching(liq_comps=Liqs, plag_comps=Plags,
                                             equationT="T_Put2008_eq24a", P=5,
                                                Ab_An_P2008=True, H2O_Liq=2)
Av_Matches=MM_dict['Av_PTs']
All_Matches=MM_dict['All_PTs']
Considering N=28 Fspar & N=24 Liqs, which is a total of N=672 Liq-Fspar pairs, be patient if this is >>1 million!
Done!!! I found a total of N=336 Fspar-Liq matches using the specified filter. N=28 Fspar out of the N=28 Fspar that you input matched to 1 or more liquids

Example 1b - All possible Liq-Plag matches, no filter

  • This could be used if you want to develop your own equilibrium filters

[5]:
MM_dict=pt.calculate_fspar_liq_temp_matching(liq_comps=Liqs, plag_comps=Plags,
                                             equationT="T_Put2008_eq24a", P=5,
                                                Ab_An_P2008=False, H2O_Liq=2)
Av_Matches=MM_dict['Av_PTs']
All_Matches=MM_dict['All_PTs']
Considering N=28 Fspar & N=24 Liqs, which is a total of N=672 Liq-Fspar pairs, be patient if this is >>1 million!
We are returning all pairs, if you want to use the Ab-An equilibrium test of Putirka (2008), enter Ab_An_P2008=True
Done!!! I found a total of N=672 Fspar-Liq matches using the specified filter. N=28 Fspar out of the N=28 Fspar that you input matched to 1 or more liquids

Example 1c - Iterating T and H2O for all possible Plag-Liq matches

[9]:
T_H_Iter_dict=pt.calculate_fspar_liq_temp_hygr_matching(liq_comps=Liqs, plag_comps=Plags,
        equationT="T_Put2008_eq24a", equationH="H_Waters2015", P=5)
T_H_Calc_Av=T_H_Iter_dict.get('Av_HTs')
T_H_Calc_All=T_H_Iter_dict.get('All_HTs')
T_H_Evol=T_H_Iter_dict.get('T_H_Evolution')
Considering N=28 Fspar & N=24 Liqs, which is a total of N=672 Liq- Fspar pairs, be patient if this is >>1 million!
20
Done!!! I found a total of N=672 Fspar-Liq matches using the specified filter. N=28 Fspar out of the N=28 Fspar that you input matched to 1 or more liquids
[10]:
T_H_Calc_Av
[10]:
ID_Fspar Mean_Sample_ID_Plag Mean_# of iterations Mean_Delta T (last 2 iters) Mean_H2O_calc Mean_Delta H (last 2 iters) Mean_Delta_An Mean_Delta_Ab Mean_Delta_Or Mean_Pred_An_EqE Mean_Pred_Ab_EqF Mean_Pred_Or_EqG Mean_Obs_Kd_Ab_An Mean_Si_Plag_cat_prop Mean_Mg_Plag_cat_prop Mean_Fet_Plag_cat_prop Mean_Ca_Plag_cat_prop Mean_Al_Plag_cat_prop Mean_Na_Plag_cat_prop Mean_K_Plag_cat_prop Mean_Mn_Plag_cat_prop Mean_Ti_Plag_cat_prop Mean_Cr_Plag_cat_prop Mean_sum Mean_Si_Plag_cat_frac Mean_Mg_Plag_cat_frac Mean_Fet_Plag_cat_frac Mean_Ca_Plag_cat_frac Mean_Al_Plag_cat_frac Mean_Na_Plag_cat_frac Mean_K_Plag_cat_frac Mean_Mn_Plag_cat_frac Mean_Ti_Plag_cat_frac Mean_Cr_Plag_cat_frac Mean_An_Plag Mean_Ab_Plag Mean_Or_Plag Mean_SiO2_Liq Mean_TiO2_Liq Mean_Al2O3_Liq Mean_FeOt_Liq Mean_MnO_Liq Mean_MgO_Liq Mean_CaO_Liq Mean_Na2O_Liq Mean_K2O_Liq Mean_Cr2O3_Liq Mean_P2O5_Liq Mean_H2O_Liq Mean_Fe3Fet_Liq Mean_NiO_Liq Mean_CoO_Liq Mean_CO2_Liq Mean_ID_liq Mean_Sample_ID_liq Mean_SiO2_Liq_mol_frac Mean_MgO_Liq_mol_frac Mean_MnO_Liq_mol_frac Mean_FeOt_Liq_mol_frac Mean_CaO_Liq_mol_frac Mean_Al2O3_Liq_mol_frac Mean_Na2O_Liq_mol_frac Mean_K2O_Liq_mol_frac Mean_TiO2_Liq_mol_frac Mean_P2O5_Liq_mol_frac Mean_Cr2O3_Liq_mol_frac Mean_Si_Liq_cat_frac Mean_Mg_Liq_cat_frac Mean_Mn_Liq_cat_frac Mean_Fet_Liq_cat_frac Mean_Ca_Liq_cat_frac Mean_Al_Liq_cat_frac Mean_Na_Liq_cat_frac Mean_K_Liq_cat_frac Mean_Ti_Liq_cat_frac Mean_P_Liq_cat_frac Mean_Cr_Liq_cat_frac Mean_Mg_Number_Liq_NoFe3 Mean_Mg_Number_Liq_Fe3 Mean_P Mean_T Mean_T_K_calc # of Liqs Averaged Std_Sample_ID_Plag Std_# of iterations Std_Delta T (last 2 iters) Std_H2O_calc Std_Delta H (last 2 iters) Std_Delta_An Std_Delta_Ab Std_Delta_Or Std_Pred_An_EqE Std_Pred_Ab_EqF Std_Pred_Or_EqG Std_Obs_Kd_Ab_An Std_Si_Plag_cat_prop Std_Mg_Plag_cat_prop Std_Fet_Plag_cat_prop Std_Ca_Plag_cat_prop Std_Al_Plag_cat_prop Std_Na_Plag_cat_prop Std_K_Plag_cat_prop Std_Mn_Plag_cat_prop Std_Ti_Plag_cat_prop Std_Cr_Plag_cat_prop Std_sum Std_Si_Plag_cat_frac Std_Mg_Plag_cat_frac Std_Fet_Plag_cat_frac Std_Ca_Plag_cat_frac Std_Al_Plag_cat_frac Std_Na_Plag_cat_frac Std_K_Plag_cat_frac Std_Mn_Plag_cat_frac Std_Ti_Plag_cat_frac Std_Cr_Plag_cat_frac Std_An_Plag Std_Ab_Plag Std_Or_Plag Std_SiO2_Liq Std_TiO2_Liq Std_Al2O3_Liq Std_FeOt_Liq Std_MnO_Liq Std_MgO_Liq Std_CaO_Liq Std_Na2O_Liq Std_K2O_Liq Std_Cr2O3_Liq Std_P2O5_Liq Std_H2O_Liq Std_Fe3Fet_Liq Std_NiO_Liq Std_CoO_Liq Std_CO2_Liq Std_ID_liq Std_Sample_ID_liq Std_SiO2_Liq_mol_frac Std_MgO_Liq_mol_frac Std_MnO_Liq_mol_frac Std_FeOt_Liq_mol_frac Std_CaO_Liq_mol_frac Std_Al2O3_Liq_mol_frac Std_Na2O_Liq_mol_frac Std_K2O_Liq_mol_frac Std_TiO2_Liq_mol_frac Std_P2O5_Liq_mol_frac Std_Cr2O3_Liq_mol_frac Std_Si_Liq_cat_frac Std_Mg_Liq_cat_frac Std_Mn_Liq_cat_frac Std_Fet_Liq_cat_frac Std_Ca_Liq_cat_frac Std_Al_Liq_cat_frac Std_Na_Liq_cat_frac Std_K_Liq_cat_frac Std_Ti_Liq_cat_frac Std_P_Liq_cat_frac Std_Cr_Liq_cat_frac Std_Mg_Number_Liq_NoFe3 Std_Mg_Number_Liq_Fe3 Std_P Std_T Std_T_K_calc Sample_ID_Plag
0 0 Plag1 20.0 0.000009 -0.152865 -2.516342e-07 0.156529 0.164705 0.028989 0.573657 0.388952 0.000226 0.687211 0.953660 0.000744 0.005985 0.148545 0.521768 0.197164 0.010404 0.000000 0.001127 0.0 1.839397 0.518464 0.000405 0.003254 0.080757 0.283663 0.107189 0.005656 0.000000 0.000613 0.0 0.417129 0.553656 0.029215 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1417.824824 1417.824834 24 Plag1 0.0 0.000035 0.575867 9.742417e-07 0.092707 0.091742 0.000277 0.092707 0.091742 0.000277 0.280843 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.126137 11.126151 Plag1
1 1 Plag2 20.0 0.000014 -0.231324 -3.722000e-07 0.127930 0.137478 0.028095 0.581055 0.381065 0.000237 0.592499 0.940345 0.001241 0.006542 0.159601 0.527653 0.182643 0.009979 0.000000 0.001502 0.0 1.829506 0.513989 0.000678 0.003576 0.087237 0.288413 0.099832 0.005455 0.000000 0.000821 0.0 0.453125 0.518543 0.028332 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1427.141564 1427.141578 24 Plag2 0.0 0.000035 0.575619 9.396090e-07 0.093896 0.089875 0.000290 0.093896 0.089875 0.000290 0.242137 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.249813 11.249827 Plag2
2 2 Plag3_core 20.0 0.000006 -0.104766 -1.716115e-07 0.152616 0.166319 0.023213 0.570877 0.391984 0.000222 0.691103 0.958653 0.001737 0.006959 0.151576 0.515884 0.202327 0.008493 0.000000 0.001377 0.0 1.847006 0.519031 0.000940 0.003768 0.082066 0.279308 0.109543 0.004598 0.000000 0.000746 0.0 0.418261 0.558303 0.023435 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1414.323909 1414.323916 24 Plag3_core 0.0 0.000036 0.575946 9.857607e-07 0.092260 0.092460 0.000272 0.092260 0.092460 0.000272 0.282434 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.078454 11.078469 Plag3_core
3 3 Plag4 20.0 0.000023 -0.389840 -6.073598e-07 0.141893 0.130335 0.049439 0.590532 0.371335 0.000252 0.578949 0.951996 0.001489 0.008630 0.161027 0.529614 0.180061 0.017835 0.000000 0.002003 0.0 1.852655 0.513855 0.000804 0.004658 0.086917 0.285868 0.097191 0.009627 0.000000 0.001081 0.0 0.448640 0.501670 0.049691 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1439.081728 1439.081751 24 Plag4 0.0 0.000034 0.575317 8.979065e-07 0.095419 0.087572 0.000308 0.095419 0.087572 0.000308 0.236600 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.410293 11.410306 Plag4
4 4 Plag5 20.0 0.000016 -0.269530 -4.269059e-07 0.115262 0.125390 0.028150 0.584614 0.377362 0.000243 0.552634 0.943674 0.002729 0.009604 0.168695 0.541383 0.180061 0.010192 0.000000 0.001753 0.0 1.858092 0.507873 0.001469 0.005169 0.090790 0.291365 0.096907 0.005485 0.000000 0.000943 0.0 0.469971 0.501636 0.028393 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1431.625552 1431.625568 24 Plag5 0.0 0.000034 0.575482 9.210654e-07 0.093678 0.087349 0.000297 0.094467 0.088998 0.000297 0.225845 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.307907 11.307921 Plag5
5 5 Plag6 20.0 0.000013 -0.220612 -3.555738e-07 0.127282 0.138124 0.026769 0.580400 0.381753 0.000236 0.594031 0.932024 0.002481 0.007794 0.160314 0.521768 0.183934 0.009555 0.000000 0.002629 0.0 1.820499 0.511961 0.001363 0.004281 0.088061 0.286607 0.101035 0.005248 0.000000 0.001444 0.0 0.453118 0.519877 0.027005 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1426.316444 1426.316458 24 Plag6 0.0 0.000035 0.575639 9.424861e-07 0.093791 0.090038 0.000289 0.093791 0.090038 0.000289 0.242763 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.238705 11.238719 Plag6
6 6 Plag7 20.0 0.000015 -0.260039 -4.086991e-07 0.102296 0.117505 0.021936 0.585243 0.376715 0.000244 0.521234 0.933688 0.002233 0.007794 0.177255 0.545307 0.178448 0.008068 0.000000 0.002629 0.0 1.855422 0.503221 0.001204 0.004201 0.095534 0.293899 0.096176 0.004348 0.000000 0.001417 0.0 0.487271 0.490550 0.022180 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1432.417065 1432.417081 24 Plag7 0.0 0.000034 0.575432 9.147264e-07 0.089663 0.083709 0.000298 0.094567 0.088844 0.000298 0.213013 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.315809 11.315823 Plag7
7 7 Plag8 20.0 0.000014 -0.242243 -3.833287e-07 0.104935 0.121268 0.021239 0.583743 0.378263 0.000241 0.534258 0.938681 0.001985 0.007377 0.171370 0.535499 0.176834 0.007644 0.000000 0.002504 0.0 1.841894 0.509628 0.001078 0.004005 0.093040 0.290733 0.096007 0.004150 0.000000 0.001359 0.0 0.481582 0.496938 0.021480 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1430.527996 1430.528010 24 Plag8 0.0 0.000034 0.575491 9.227664e-07 0.091093 0.085495 0.000295 0.094327 0.089211 0.000295 0.218336 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.291525 11.291539 Plag8
8 8 Plag9 20.0 0.000019 -0.321922 -4.982051e-07 0.103970 0.112959 0.029827 0.589238 0.372639 0.000250 0.509102 0.932024 0.002233 0.008908 0.172619 0.529614 0.169735 0.010616 0.000000 0.001753 0.0 1.827502 0.509999 0.001222 0.004874 0.094456 0.289802 0.092878 0.005809 0.000000 0.000959 0.0 0.489046 0.480877 0.030077 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1437.451418 1437.451437 24 Plag9 0.0 0.000034 0.575293 8.958051e-07 0.090884 0.081445 0.000306 0.095209 0.087879 0.000306 0.208055 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.382352 11.382366 Plag9
9 9 Plag10 20.0 0.000021 -0.371962 -5.414099e-07 0.089059 0.085056 0.020477 0.596858 0.365061 0.000262 0.400083 0.907059 0.005955 0.010439 0.197941 0.539422 0.152955 0.007431 0.000000 0.000000 0.0 1.821202 0.498055 0.003270 0.005732 0.108687 0.296190 0.083986 0.004080 0.000000 0.000000 0.0 0.552402 0.426859 0.020739 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1447.055772 1447.055793 24 Plag10 0.0 0.000032 0.574843 8.389790e-07 0.055535 0.062015 0.000320 0.096428 0.086080 0.000320 0.163502 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.492334 11.492347 Plag10
10 10 Plag11 20.0 0.000021 -0.396379 -5.540146e-07 0.089645 0.079034 0.016743 0.600838 0.361201 0.000268 0.350354 0.895409 0.004962 0.009325 0.212207 0.555114 0.143597 0.006157 0.000000 0.000000 0.0 1.826772 0.490159 0.002716 0.005105 0.116165 0.303877 0.078607 0.003371 0.000000 0.000000 0.0 0.586269 0.396720 0.017011 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1452.074689 1452.074710 24 Plag11 0.0 0.000031 0.574554 8.039550e-07 0.035455 0.045289 0.000328 0.097062 0.085162 0.000328 0.143180 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.544579 11.544592 Plag11
11 11 Plag12 20.0 0.000021 -0.392990 -5.459771e-07 0.089685 0.078577 0.014576 0.601117 0.360932 0.000269 0.341694 0.900402 0.005458 0.010439 0.212207 0.564922 0.140048 0.005308 0.000000 0.000000 0.0 1.838783 0.489672 0.002969 0.005677 0.115406 0.307226 0.076163 0.002887 0.000000 0.000000 0.0 0.593482 0.391673 0.014845 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1452.426711 1452.426732 24 Plag12 0.0 0.000031 0.574514 7.992196e-07 0.033126 0.042248 0.000329 0.097106 0.085098 0.000329 0.139640 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.546313 11.546325 Plag12
12 12 Plag13 20.0 0.000022 -0.434295 -5.693605e-07 0.090462 0.076121 0.017155 0.606403 0.355913 0.000278 0.294059 0.865118 0.005210 0.008490 0.228969 0.574730 0.130044 0.006370 0.000141 0.001502 0.0 1.820575 0.475189 0.002862 0.004664 0.125768 0.315686 0.071430 0.003499 0.000077 0.000825 0.0 0.626655 0.355912 0.017433 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1459.095295 1459.095318 24 Plag13 0.0 0.000029 0.574108 7.517030e-07 0.038506 0.031519 0.000340 0.097948 0.083904 0.000340 0.120173 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.613181 11.613193 Plag13
13 13 Plag14 20.0 0.000024 -0.464925 -6.173124e-07 0.090537 0.076026 0.026719 0.606889 0.355458 0.000278 0.315424 0.880263 0.020593 0.016563 0.213990 0.556683 0.130367 0.009555 0.000423 0.001627 0.0 1.830065 0.481001 0.011253 0.009051 0.116930 0.304188 0.071236 0.005221 0.000231 0.000889 0.0 0.604643 0.368360 0.026997 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1459.708251 1459.708275 24 Plag14 0.0 0.000029 0.574187 7.601939e-07 0.032584 0.034128 0.000341 0.098029 0.083799 0.000341 0.128904 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.630428 11.630439 Plag14
14 14 Plag15 20.0 0.000018 -0.397618 -4.595993e-07 0.114366 0.085592 0.007968 0.611171 0.351477 0.000286 0.202410 0.831831 0.005458 0.002784 0.256788 0.619845 0.100389 0.002973 0.000000 0.001127 0.0 1.821194 0.456750 0.002997 0.001529 0.141000 0.340351 0.055123 0.001632 0.000000 0.000619 0.0 0.713004 0.278743 0.008254 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1465.114311 1465.114329 24 Plag15 0.0 0.000026 0.573293 6.632173e-07 0.083145 0.068832 0.000350 0.098694 0.082837 0.000350 0.082719 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.629039 11.629049 Plag15
15 15 Plag16 20.0 0.000020 -0.421155 -5.118339e-07 0.097190 0.076320 0.010953 0.609853 0.352695 0.000283 0.234567 0.858461 0.006203 0.011274 0.244305 0.597876 0.110683 0.004034 0.000000 0.001627 0.0 1.834463 0.467963 0.003381 0.006146 0.133175 0.325913 0.060335 0.002199 0.000000 0.000887 0.0 0.680474 0.308290 0.011236 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1463.450180 1463.450200 24 Plag16 0.0 0.000027 0.573609 6.966147e-07 0.071049 0.053764 0.000347 0.098491 0.083132 0.000347 0.095861 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.633754 11.633765 Plag16
16 16 Plag17 20.0 0.000022 -0.433805 -5.551707e-07 0.091520 0.075851 0.015166 0.607733 0.354667 0.000280 0.273596 0.864785 0.010173 0.010300 0.230218 0.578456 0.121654 0.005520 0.000000 0.001878 0.0 1.822984 0.474379 0.005580 0.005650 0.126286 0.317313 0.066734 0.003028 0.000000 0.001030 0.0 0.644159 0.340394 0.015446 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1460.773734 1460.773756 24 Plag17 0.0 0.000028 0.573947 7.335450e-07 0.047745 0.034627 0.000343 0.098159 0.083606 0.000343 0.111811 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.624195 11.624206 Plag17
17 17 Plag18 20.0 0.000023 -0.448259 -5.836822e-07 0.090603 0.075928 0.020316 0.607360 0.355016 0.000279 0.292089 0.869612 0.013894 0.011135 0.225938 0.568060 0.127463 0.007431 0.000000 0.002253 0.0 1.825786 0.476294 0.007610 0.006099 0.123748 0.311132 0.069813 0.004070 0.000000 0.001234 0.0 0.626158 0.353247 0.020595 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1460.303106 1460.303128 24 Plag18 0.0 0.000029 0.574059 7.459361e-07 0.037773 0.031492 0.000342 0.098102 0.083691 0.000342 0.119368 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.627570 11.627581 Plag18
18 18 Plag19 20.0 0.000024 -0.455126 -6.324904e-07 0.090046 0.079497 0.031106 0.603538 0.358621 0.000273 0.364306 0.878432 0.004466 0.010300 0.207749 0.554918 0.146179 0.011465 0.000282 0.001753 0.0 1.815544 0.483840 0.002460 0.005673 0.114428 0.305648 0.080515 0.006315 0.000155 0.000965 0.0 0.568562 0.400059 0.031378 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1455.479984 1455.480008 24 Plag19 0.0 0.000031 0.574500 7.966720e-07 0.048180 0.048440 0.000334 0.097497 0.084553 0.000334 0.148881 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.593366 11.593378 Plag19
19 19 Plag20 20.0 0.000014 -0.332381 -3.537067e-07 0.150825 0.119029 0.005012 0.611779 0.350914 0.000287 0.158076 0.811027 0.005458 0.004872 0.274799 0.617099 0.083900 0.001911 0.000000 0.000751 0.0 1.799816 0.450617 0.003033 0.002707 0.152682 0.342868 0.046616 0.001062 0.000000 0.000417 0.0 0.762040 0.232661 0.005299 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1465.883918 1465.883932 24 Plag20 0.0 0.000024 0.572803 6.138485e-07 0.097877 0.081522 0.000351 0.098778 0.082692 0.000351 0.064601 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.592705 11.592714 Plag20
20 20 Plag21 20.0 0.000017 -0.385527 -4.402840e-07 0.121034 0.089973 0.006150 0.611187 0.351461 0.000286 0.194391 0.834494 0.006203 0.004872 0.262138 0.599837 0.098421 0.002336 0.000141 0.001002 0.0 1.809442 0.461189 0.003428 0.002692 0.144872 0.331504 0.054393 0.001291 0.000078 0.000553 0.0 0.722354 0.271210 0.006436 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1465.135556 1465.135573 24 Plag21 0.0 0.000026 0.573220 6.557926e-07 0.085731 0.071653 0.000350 0.098695 0.082831 0.000350 0.079442 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.622237 11.622247 Plag21
21 21 Plag22 20.0 0.000020 -0.416992 -5.113539e-07 0.096254 0.075985 0.009630 0.609244 0.353259 0.000282 0.239137 0.856297 0.005707 0.007098 0.246623 0.589441 0.113910 0.003609 0.000000 0.001377 0.0 1.824063 0.469445 0.003128 0.003892 0.135206 0.323147 0.062448 0.001979 0.000000 0.000755 0.0 0.677271 0.312816 0.009912 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1462.681697 1462.681717 24 Plag22 0.0 0.000027 0.573669 7.031835e-07 0.069589 0.051140 0.000346 0.098394 0.083267 0.000346 0.097728 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.627457 11.627467 Plag22
22 22 Plag23 20.0 0.000021 -0.432072 -5.490876e-07 0.091988 0.075789 0.014334 0.608039 0.354381 0.000280 0.267686 0.868280 0.005707 0.009186 0.235924 0.580222 0.121977 0.005308 0.000423 0.001627 0.0 1.828654 0.474819 0.003121 0.005024 0.129015 0.317294 0.066703 0.002903 0.000231 0.000890 0.0 0.649554 0.335831 0.014614 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1461.160790 1461.160812 24 Plag23 0.0 0.000028 0.573900 7.284139e-07 0.051121 0.036659 0.000343 0.098207 0.083537 0.000343 0.109396 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.625848 11.625859 Plag23
23 23 Plag24 20.0 0.000022 -0.434452 -5.756347e-07 0.090372 0.076247 0.018320 0.605788 0.356491 0.000277 0.303920 0.887919 0.007940 0.011552 0.225938 0.565510 0.132626 0.006794 0.000141 0.002879 0.0 1.841300 0.482224 0.004312 0.006274 0.122706 0.307126 0.072028 0.003690 0.000077 0.001564 0.0 0.618401 0.363002 0.018596 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1458.319397 1458.319420 24 Plag24 0.0 0.000029 0.574180 7.598295e-07 0.034909 0.032264 0.000338 0.097851 0.084042 0.000338 0.124203 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.607756 11.607767 Plag24
24 24 Plag25 20.0 0.000016 -0.286227 -4.272864e-07 0.088252 0.087951 0.006315 0.591410 0.370451 0.000253 0.410372 0.905395 0.001489 0.002227 0.197049 0.564922 0.156182 0.002336 0.000000 0.000000 0.0 1.829599 0.494860 0.000814 0.001217 0.107701 0.308768 0.085364 0.001277 0.000000 0.000000 0.0 0.554183 0.439248 0.006569 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1440.189524 1440.189540 24 Plag25 0.0 0.000033 0.575052 8.661768e-07 0.049493 0.067070 0.000310 0.095553 0.087357 0.000310 0.167707 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.403394 11.403408 Plag25
25 25 Plag26 20.0 0.000019 -0.355752 -4.940624e-07 0.089466 0.078352 0.004309 0.599648 0.362346 0.000266 0.329062 0.882094 0.000744 0.003897 0.225938 0.580614 0.143597 0.001699 0.000000 0.000000 0.0 1.838583 0.479768 0.000405 0.002120 0.122887 0.315794 0.078102 0.000924 0.000000 0.000000 0.0 0.608613 0.386811 0.004576 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1450.575320 1450.575339 24 Plag26 0.0 0.000031 0.574522 8.008386e-07 0.033397 0.038950 0.000326 0.096869 0.085431 0.000326 0.134478 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.517846 11.517858 Plag26
26 26 Plag27 20.0 0.000016 -0.289504 -4.309242e-07 0.088308 0.087273 0.006109 0.591790 0.370072 0.000254 0.406318 0.900402 0.000496 0.003340 0.204360 0.562960 0.160377 0.002336 0.000000 0.000000 0.0 1.834272 0.490877 0.000271 0.001821 0.111412 0.306912 0.087434 0.001273 0.000000 0.000000 0.0 0.556729 0.436908 0.006363 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1440.667623 1440.667639 24 Plag27 0.0 0.000033 0.575030 8.633897e-07 0.047826 0.065796 0.000311 0.095614 0.087267 0.000311 0.166050 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.408898 11.408911 Plag27
27 27 Plag28 20.0 0.000016 -0.291488 -4.339911e-07 0.088317 0.087392 0.006722 0.591847 0.370015 0.000254 0.407528 0.907059 0.000744 0.003619 0.202934 0.560999 0.159732 0.002548 0.000000 0.000000 0.0 1.837634 0.493601 0.000405 0.001969 0.110432 0.305283 0.086923 0.001386 0.000000 0.000000 0.0 0.555658 0.437366 0.006976 49.082083 3.40875 14.079167 13.005 0.180833 4.797917 8.408333 3.212083 1.243333 0.0 0.782917 0.0925 0.0 0.0 0.0 0.0 11.5 11.5 0.53705 0.078143 0.001675 0.119259 0.098482 0.090789 0.034093 0.008698 0.028164 0.003648 0.0 0.472184 0.068771 0.001474 0.10486 0.08669 0.159639 0.059923 0.015287 0.024757 0.006414 0.0 0.395638 0.395638 5.0 1440.739581 1440.739598 24 Plag28 0.0 0.000033 0.575031 8.634728e-07 0.048698 0.066159 0.000311 0.095623 0.087253 0.000311 0.166545 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 3.680864 1.171897 1.065274 2.169584 0.036584 1.213167 1.941711 0.394098 0.443442 0.0 0.480756 0.124035 0.0 0.0 0.0 0.0 7.071068 7.071068 0.037479 0.018976 0.000334 0.021092 0.022046 0.006564 0.004228 0.003124 0.009896 0.002268 0.0 0.031852 0.016971 0.000299 0.01849 0.019822 0.011072 0.007131 0.005473 0.008687 0.003992 0.0 0.079518 0.079518 0.0 11.410121 11.410134 Plag28
[ ]:
    FsparNumbers = Combo_fspar_liqs2['ID_Fspar'].unique()

    Combo_fspar_liqs3=Combo_fspar_liqs2.drop(['Pass An-Ab Eq Test Put2008?', 'T_K_calc'], axis=1)
    Combo_fspar_liqs3['T_K_calc']=Combo_fspar_liqs2['T_K_calc'].astype(float)
    if len(FsparNumbers) > 0:
        if plag_comps is not None:
            df1_Mean_nopref=Combo_fspar_liqs3.groupby(['ID_Fspar', 'Sample_ID_Plag'], as_index=False).mean()
            df1_Std_nopref=Combo_fspar_liqs3.groupby(['ID_Fspar', 'Sample_ID_Plag'], as_index=False).std()
            count=Combo_fspar_liqs2.groupby('ID_Fspar',as_index=False).count().iloc[:, 1]
            df1_Mean_nopref['# of Liqs Averaged']=count

            Sample_ID_Fspar_Mean=df1_Mean_nopref['Sample_ID_Plag']
            Sample_ID_Fspar_Std=df1_Std_nopref['Sample_ID_Plag']
            df1_Mean=df1_Mean_nopref.add_prefix('Mean_')
            df1_Std=df1_Std_nopref.add_prefix('Std_')

            df1_Mean.rename(columns={"Mean_ID_Fspar": "ID_Fspar"}, inplace=True)
            df1_Mean.rename(columns={"Mean_# of Liqs Averaged": "# of Liqs Averaged"}, inplace=True)
            df1_Std.rename(columns={"Std_ID_Fspar": "ID_Fspar"}, inplace=True)



            df1_M=pd.merge(df1_Mean, df1_Std, on=['ID_Fspar'])
            df1_M['Sample_ID_Plag']=Sample_ID_Fspar_Mean

    else:
        raise Exception(
            'No Matches - to set less strict filters, change our Kd filter')


    print('Done!!! I found a total of N='+str(len(Combo_fspar_liqs3)) + ' Fspar-Liq matches using the specified filter. N=' + str(len(df1_M)) + ' Fspar out of the N='+str(LenFspar)+' Fspar that you input matched to 1 or more liquids')

    return {'Av_PTs': df1_M, 'All_PTs': Combo_fspar_liqs2}
    # Now we do the averaging step for each feldspar crystal



Again, lets filter using the equilibrium test of Putirka (2008)

[ ]:
T_H_Calc_Filter=T_H_Calc.loc[T_H_Calc['Pass An-Ab Eq Test Put2008?'].str.contains('Yes')]
[ ]:
# Maybe you want to plot a histogram of H2O contents
plt.hist(T_H_Calc_Filter['H2O_calc'])
plt.xlabel('H$_2$O Calc (wt%)')

Example 2 - Kspar-Liquid calculations

[ ]:
# Load in the Kspars
out_Kspar=pt.import_excel('Feldspar_Liquid.xlsx', sheet_name="Kspar_only")
Kspars=out_Kspar['Kspars']

# Load in the liquids you want to match with
out_Liq2=pt.import_excel('Feldspar_Liquid.xlsx', sheet_name="Evolved_Liq_only")
Liqs=out_Liq2['Liqs']

2a - Lets calculate temperature for all liq-kspar matches

  • Currently no equilibrium tests exist for Kspar-Liquid that we are aware off. If you know of ones I can implement, please reach out!

[ ]:
melt_match_2=pt.calculate_fspar_liq_temp_matching(liq_comps=Liqs, kspar_comps=Kspars,
                                             equationT="T_Put2008_eq24b", P=5,
                                                eq_tests=True, H2O_Liq=2)
melt_match_2