Add Ti-in-biotite thermometer of Wu & Chen (2015)

Author: Pierre Lanari

Changelog.md

@@ -18,6 +18,7 @@
 - External functions:
   - Add thermobarometry functions: Phengite\_P\_All
   - Add multi-equilibrium functions: GrtOpx\_T\_All; GrtPh\_T\_All
+  - Add Ti-in-biotite thermometer of Wu & Chen (2015)
 - Other:
   - Fix an issue in the EPMA calibration assistant resulting the program to crash when no standard is available for the first phase in the list
   - Fix an issue preventing a function description to be displayed for external function with the standalone versions  

Program/External/Core_Call_TB.m

@@ -8,10 +8,10 @@
         InputData = GenerateInputData(Data,OxList,InputVariables);
         [OutputData,OutputVariables] = Amphibole_T_All(InputData,InputVariables,AddParameters,ElOxDataDef);
 
-    case 'Biotite_T_H05'
+    case 'Biotite_T_All'
         InputVariables = {'SiO2','TiO2','Al2O3','FeO','Fe2O3','MnO','MgO','CaO','Na2O','K2O'};
         InputData = GenerateInputData(Data,OxList,InputVariables);
-        [OutputData,OutputVariables] = Biotite_T_H05(InputData,InputVariables,AddParameters,ElOxDataDef);
+        [OutputData,OutputVariables] = Biotite_T_All(InputData,InputVariables,AddParameters,ElOxDataDef);
 
     case 'Chlorite_T_All'
         InputVariables = {'SiO2','TiO2','Al2O3','FeO','Fe2O3','MnO','MgO','CaO','Na2O','K2O'};

Program/External/Core_Function_Indexing.m

@@ -106,10 +106,10 @@
 Functions.Min(Min).TB = [];
 
 Fct = length(Functions.Min(Min).TB)+1;
-Functions.Min(Min).TB.Name{Fct} = 'T.Bt (Henry et al. 2005)';
-Functions.Min(Min).TB.FileName{Fct} = 'Biotite_T_H05';
-Functions.Min(Min).TB.Details(Fct).AddVar = {};
-Functions.Min(Min).TB.Details(Fct).Def = [];
+Functions.Min(Min).TB.Name{Fct} = 'T.Bt (all calibrations)';
+Functions.Min(Min).TB.FileName{Fct} = 'Biotite_T_All';
+Functions.Min(Min).TB.Details(Fct).AddVar = {'P_kbar'};
+Functions.Min(Min).TB.Details(Fct).Def = [7];
 
 % -------------------------------------------------------------------------
 % Multi-equilibrium thermobarometry

Program/External/TB/Biotite_T_H05.m renamed to Program/External/TB/Biotite_T_All.m

@@ -3,17 +3,21 @@
 % XMapTools External Function: biotite thermobarometry 
 %  
 % List of output variables:
-%       - T_H05      Henry et al. (2005)
+%       - T_W20     Wu & Chen (2015)
+%       - T_H05     Henry et al. (2005)
 %
 % 22 Oxygen-basis 
 %
-% P. Lanari - Last update 06.03.2021
+% P. Lanari - Last update 10.02.2022
 % Find out more at https://xmaptools.ch
 
-OutputVariables = {'T_H05','Ti','XMg'};
+OutputVariables = {'T_H05','T_W15','Ti','XMg'};
 
 OutputData = zeros(size(InputData,1),length(OutputVariables));
 
+P_kbar = AddParameters.Values(1);
+P_GPa = P_kbar/10;
+
 % General structural formula function for 22 oxygen
 WhereMin = find(sum(InputData,2) > 50);
 [MatrixSF,ElementsList] = SF_OxNorm(InputData(WhereMin,:),InputVariables,22,ElOxDataDef);
@@ -30,6 +34,11 @@
 
 XMg = Mg./(Mg+Fe);
 
+Alvi = Al-(6-Si);
+X_Ti = Ti./(Fe+Mg+Alvi+Ti);
+X_Fe = Fe./(Fe+Mg+Alvi+Ti);
+X_Mg = Mg./(Fe+Mg+Alvi+Ti);
+
 A = -2.3594;
 B = 4.6482e-9;
 C = -1.7283;
@@ -42,7 +51,11 @@
     end
 end
 
-OutputData(WhereMin,:) = [T_H05,Ti,XMg];
+T_W15 = exp(6.313+0.224.*log(X_Ti)-0.288.*log(X_Fe)-0.449.*log(X_Mg)+0.15.*P_GPa);
+Filter = find(T_W15 < 100 | T_W15 > 1000);
+T_W15(Filter) = zeros(size(Filter));
+
+OutputData(WhereMin,:) = [T_H05,T_W15,Ti,XMg];
 
 end