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

Plag-Liq and Kspar-Liq Matching¶

This notebook shows how to assess all possible matches for Plag-Liq and Kspar-Liq pairs where there aren’t obvious pairwise analyses
You can download the excel spreadsheet here: https://github.com/PennyWieser/Thermobar/blob/main/docs/Examples/Feldspar_Thermobarometry/Feldspar_Liquid.xlsx

[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	MgO_Plag	CaO_Plag	Na2O_Plag	K2O_Plag	Sample_ID_Plag
0	57.3	0.09	26.6	0.43	0.03	8.33	6.11	0.49	Plag1
1	56.5	0.12	26.9	0.47	0.05	8.95	5.66	0.47	Plag2
2	57.6	0.11	26.3	0.50	0.07	8.50	6.27	0.40	Plag3_core
3	57.2	0.16	27.0	0.62	0.06	9.03	5.58	0.84	Plag4
4	56.7	0.14	27.6	0.69	0.11	9.46	5.58	0.48	Plag5

	SiO2_Liq	TiO2_Liq	Al2O3_Liq	FeOt_Liq	MnO_Liq	MgO_Liq	CaO_Liq	Na2O_Liq	K2O_Liq	P2O5_Liq	Sample_ID_Liq
0	49.1	3.22	14.4	14.8	0.14	3.20	6.72	3.34	1.70	1.13	Glass1
1	49.2	3.89	15.3	13.7	0.12	3.88	6.76	3.44	1.22	0.83	Glass2
2	49.6	3.79	15.8	13.0	0.14	4.26	6.59	3.65	1.04	0.63	Glass3
3	47.1	4.21	12.0	17.8	0.18	3.40	7.28	2.93	2.02	2.32	Glass4
4	48.1	3.88	13.2	16.4	0.16	4.02	6.51	3.36	1.36	1.59	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_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_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_P2O5_Liq	Mean_H2O_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_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_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_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_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_P2O5_Liq	Std_H2O_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_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_Mg_Number_Liq_NoFe3	Std_Mg_Number_Liq_Fe3	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	1.839397	0.518464	0.000405	0.003254	0.080757	0.283663	0.107189	0.005656	0.000000	0.000613	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1417.824824	1417.824834	24	Plag1	0.000035	0.575867	9.742417e-07	0.092707	0.091742	0.000277	0.092707	0.091742	0.000277	0.280843	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.829506	0.513989	0.000678	0.003576	0.087237	0.288413	0.099832	0.005455	0.000000	0.000821	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1427.141564	1427.141578	24	Plag2	0.000035	0.575619	9.396090e-07	0.093896	0.089875	0.000290	0.093896	0.089875	0.000290	0.242137	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.847006	0.519031	0.000940	0.003768	0.082066	0.279308	0.109543	0.004598	0.000000	0.000746	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1414.323909	1414.323916	24	Plag3_core	0.000036	0.575946	9.857607e-07	0.092260	0.092460	0.000272	0.092260	0.092460	0.000272	0.282434	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.852655	0.513855	0.000804	0.004658	0.086917	0.285868	0.097191	0.009627	0.000000	0.001081	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1439.081728	1439.081751	24	Plag4	0.000034	0.575317	8.979065e-07	0.095419	0.087572	0.000308	0.095419	0.087572	0.000308	0.236600	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.858092	0.507873	0.001469	0.005169	0.090790	0.291365	0.096907	0.005485	0.000000	0.000943	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1431.625552	1431.625568	24	Plag5	0.000034	0.575482	9.210654e-07	0.093678	0.087349	0.000297	0.094467	0.088998	0.000297	0.225845	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.820499	0.511961	0.001363	0.004281	0.088061	0.286607	0.101035	0.005248	0.000000	0.001444	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1426.316444	1426.316458	24	Plag6	0.000035	0.575639	9.424861e-07	0.093791	0.090038	0.000289	0.093791	0.090038	0.000289	0.242763	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.855422	0.503221	0.001204	0.004201	0.095534	0.293899	0.096176	0.004348	0.000000	0.001417	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1432.417065	1432.417081	24	Plag7	0.000034	0.575432	9.147264e-07	0.089663	0.083709	0.000298	0.094567	0.088844	0.000298	0.213013	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.841894	0.509628	0.001078	0.004005	0.093040	0.290733	0.096007	0.004150	0.000000	0.001359	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1430.527996	1430.528010	24	Plag8	0.000034	0.575491	9.227664e-07	0.091093	0.085495	0.000295	0.094327	0.089211	0.000295	0.218336	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.827502	0.509999	0.001222	0.004874	0.094456	0.289802	0.092878	0.005809	0.000000	0.000959	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1437.451418	1437.451437	24	Plag9	0.000034	0.575293	8.958051e-07	0.090884	0.081445	0.000306	0.095209	0.087879	0.000306	0.208055	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.821202	0.498055	0.003270	0.005732	0.108687	0.296190	0.083986	0.004080	0.000000	0.000000	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1447.055772	1447.055793	24	Plag10	0.000032	0.574843	8.389790e-07	0.055535	0.062015	0.000320	0.096428	0.086080	0.000320	0.163502	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.826772	0.490159	0.002716	0.005105	0.116165	0.303877	0.078607	0.003371	0.000000	0.000000	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1452.074689	1452.074710	24	Plag11	0.000031	0.574554	8.039550e-07	0.035455	0.045289	0.000328	0.097062	0.085162	0.000328	0.143180	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.838783	0.489672	0.002969	0.005677	0.115406	0.307226	0.076163	0.002887	0.000000	0.000000	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1452.426711	1452.426732	24	Plag12	0.000031	0.574514	7.992196e-07	0.033126	0.042248	0.000329	0.097106	0.085098	0.000329	0.139640	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.820575	0.475189	0.002862	0.004664	0.125768	0.315686	0.071430	0.003499	0.000077	0.000825	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1459.095295	1459.095318	24	Plag13	0.000029	0.574108	7.517030e-07	0.038506	0.031519	0.000340	0.097948	0.083904	0.000340	0.120173	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.830065	0.481001	0.011253	0.009051	0.116930	0.304188	0.071236	0.005221	0.000231	0.000889	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1459.708251	1459.708275	24	Plag14	0.000029	0.574187	7.601939e-07	0.032584	0.034128	0.000341	0.098029	0.083799	0.000341	0.128904	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.821194	0.456750	0.002997	0.001529	0.141000	0.340351	0.055123	0.001632	0.000000	0.000619	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1465.114311	1465.114329	24	Plag15	0.000026	0.573293	6.632173e-07	0.083145	0.068832	0.000350	0.098694	0.082837	0.000350	0.082719	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.834463	0.467963	0.003381	0.006146	0.133175	0.325913	0.060335	0.002199	0.000000	0.000887	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1463.450180	1463.450200	24	Plag16	0.000027	0.573609	6.966147e-07	0.071049	0.053764	0.000347	0.098491	0.083132	0.000347	0.095861	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.822984	0.474379	0.005580	0.005650	0.126286	0.317313	0.066734	0.003028	0.000000	0.001030	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1460.773734	1460.773756	24	Plag17	0.000028	0.573947	7.335450e-07	0.047745	0.034627	0.000343	0.098159	0.083606	0.000343	0.111811	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.825786	0.476294	0.007610	0.006099	0.123748	0.311132	0.069813	0.004070	0.000000	0.001234	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1460.303106	1460.303128	24	Plag18	0.000029	0.574059	7.459361e-07	0.037773	0.031492	0.000342	0.098102	0.083691	0.000342	0.119368	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.815544	0.483840	0.002460	0.005673	0.114428	0.305648	0.080515	0.006315	0.000155	0.000965	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1455.479984	1455.480008	24	Plag19	0.000031	0.574500	7.966720e-07	0.048180	0.048440	0.000334	0.097497	0.084553	0.000334	0.148881	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.799816	0.450617	0.003033	0.002707	0.152682	0.342868	0.046616	0.001062	0.000000	0.000417	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1465.883918	1465.883932	24	Plag20	0.000024	0.572803	6.138485e-07	0.097877	0.081522	0.000351	0.098778	0.082692	0.000351	0.064601	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.809442	0.461189	0.003428	0.002692	0.144872	0.331504	0.054393	0.001291	0.000078	0.000553	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1465.135556	1465.135573	24	Plag21	0.000026	0.573220	6.557926e-07	0.085731	0.071653	0.000350	0.098695	0.082831	0.000350	0.079442	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.824063	0.469445	0.003128	0.003892	0.135206	0.323147	0.062448	0.001979	0.000000	0.000755	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1462.681697	1462.681717	24	Plag22	0.000027	0.573669	7.031835e-07	0.069589	0.051140	0.000346	0.098394	0.083267	0.000346	0.097728	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.828654	0.474819	0.003121	0.005024	0.129015	0.317294	0.066703	0.002903	0.000231	0.000890	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1461.160790	1461.160812	24	Plag23	0.000028	0.573900	7.284139e-07	0.051121	0.036659	0.000343	0.098207	0.083537	0.000343	0.109396	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.841300	0.482224	0.004312	0.006274	0.122706	0.307126	0.072028	0.003690	0.000077	0.001564	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1458.319397	1458.319420	24	Plag24	0.000029	0.574180	7.598295e-07	0.034909	0.032264	0.000338	0.097851	0.084042	0.000338	0.124203	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.829599	0.494860	0.000814	0.001217	0.107701	0.308768	0.085364	0.001277	0.000000	0.000000	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1440.189524	1440.189540	24	Plag25	0.000033	0.575052	8.661768e-07	0.049493	0.067070	0.000310	0.095553	0.087357	0.000310	0.167707	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.838583	0.479768	0.000405	0.002120	0.122887	0.315794	0.078102	0.000924	0.000000	0.000000	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1450.575320	1450.575339	24	Plag26	0.000031	0.574522	8.008386e-07	0.033397	0.038950	0.000326	0.096869	0.085431	0.000326	0.134478	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.834272	0.490877	0.000271	0.001821	0.111412	0.306912	0.087434	0.001273	0.000000	0.000000	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1440.667623	1440.667639	24	Plag27	0.000033	0.575030	8.633897e-07	0.047826	0.065796	0.000311	0.095614	0.087267	0.000311	0.166050	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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	1.837634	0.493601	0.000405	0.001969	0.110432	0.305283	0.086923	0.001386	0.000000	0.000000	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.782917	0.0925	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.472184	0.068771	0.001474	0.10486	0.08669	0.159639	0.059923	0.015287	0.024757	0.006414	0.395638	0.395638	5.0	1440.739581	1440.739598	24	Plag28	0.000033	0.575031	8.634728e-07	0.048698	0.066159	0.000311	0.095623	0.087253	0.000311	0.166545	3.680864	1.171897	1.065274	2.169584	0.036584	1.213167	1.941711	0.394098	0.443442	0.480756	0.124035	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.031852	0.016971	0.000299	0.01849	0.019822	0.011072	0.007131	0.005473	0.008687	0.003992	0.079518	0.079518	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