%====================================================== % AX file for NCFMASCrO involving ol cpx opx g sp pl liq %====================================================== % ================== di-jd-ct-fo-fa-q melt ======================== liq 8 1 % approx model works well from 1bar to 60 kbar ct(L) 0.3 jd(L) 0.15 di(L) 0.2 fa(L) 0.05 fo(L) 0.1 hm(L) 0.001 ek(L) 0.001 %------------------------------------------------------------------- p(q) 1 1 1 7 -1 ct -1 di -1 jd -1 fa -1 fo -1 hm -1 ek p(di) 1 1 0 1 1 di p(jd) 1 1 0 1 1 jd p(ct) 1 1 0 1 1 ct p(fo) 1 1 0 1 1 fo p(fa) 1 1 0 1 1 fa p(hm) 1 1 0 1 1 hm p(ek) 1 1 0 1 1 ek % -------------------------------------------------- sf W(q,di) 26 0 -0.4 W(q,jd) -10 0 0 W(q,ct) -10 0 0 W(q,fo) -25 0 -0.1 W(q,fa) -7 0 0 W(q,hm) 0 0 0 W(q,ek) 0 0 0 W(di,jd) 0 0 0 W(di,ct) -2 0 0 W(di,fo) 24 0 0.2 W(di,fa) 17 0 0 W(di,hm) 0 0 0 W(di,ek) 0 0 0 W(jd,ct) 0 0 0 W(jd,fo) 0 0 0 W(jd,fa) 0 0 0 W(jd,hm) 0 0 0 W(jd,ek) 0 0 0 W(ct,fo) -1 0 0.1 W(ct,fa) 0 0 0 W(ct,hm) 0 0 0 W(ct,ek) 0 0 0 W(fo,fa) 9 0 0 W(fo,hm) 0 0 0 W(fo,ek) 0 0 0 W(fa,hm) 0 0 0 W(fa,ek) 0 0 0 W(hm,ek) 0 0 0 % -------------------------------------------------- 12 pq 1 1 1 7 -1 ct -1 di -1 jd -1 fa -1 fo -1 hm -1 ek pdi 1 1 0 1 1 di pjd 1 1 0 1 1 jd pct 1 1 0 1 1 ct pfa 1 1 0 1 1 fa pfo 1 1 0 1 1 fo phm 1 1 0 1 1 hm pek 1 1 0 1 1 ek feM 1 1 0 1 2 fa mgM 1 1 0 2 2 fo 1 di ol 1 1 0 2 1 fo 1 fa sumM 1 1 0 3 2 fa 2 di 2 fo % -------------------------------------------------- % ideal mixing activities q 1 1 pq 1 make 1 qL 2 DQF 4.0 0 -0.43 di 4 3 pdi 2 mgM 1 sumM -2 make 1 diL 1 DQF -6.7 0 0.17 jd 1 1 pjd 1 make 2 abL 1 qL -1 DQF 18.7 0 -0.16 ct 1 1 pct 1 make 2 anL 1 qL -1 DQF 2.5 0 0.29 fo 1 3 mgM 2 ol 1 sumM -2 make 1 foL 1 DQF 2.0 0 -0.03 fa 1 3 feM 2 ol 1 sumM -2 make 1 faL 1 DQF -2.5 0 0.18 hm 1 1 phm 1 make 1 hemL 1/2 DQF 0.0 0 0.20 ek 1 1 pek 1 make 1 eskL 1/2 DQF 21.0 0 0 % __________________________________________________________________ % ================== plag 2T ASF ======================== pl 2 1 Ca(pl) 0.8 %------------------------------------------------------------------- p(an) 1 1 0 1 1 Ca p(ab) 1 1 1 1 -1 Ca % -------------------------------------------------- asf W(an,ab) 22.4 0 0 an 0.39 0 0 ab 1.00 0 0 % -------------------------------------------------- 4 xNaA 1 1 1 1 -1 Ca xCaA 1 1 0 1 1 Ca xAlT 1 1 1/2 1 1/2 Ca xSiT 1 1 1/2 1 -1/2 Ca % -------------------------------------------------- % Al-avoidance mixing activities an 1 2 xCaA 1 xAlT 2 ab 4 3 xNaA 1 xAlT 1 xSiT 1 make 1 abh 1 %__________________________________________________________________ % ================================================= % ------------------------------------------------- % cpx: NCFMASO % clinopyroxene % % M1 M2 T % Mg Fe Al Fe3 Cr Mg Fe Ca Na Si Al % di 1 0 0 0 0 0 0 1 0 2 0 % fs 0 1 0 0 0 0 1 0 0 2 0 % cats 0 0 1 0 0 0 0 1 0 1 1 % crdi 0 0 0 0 1 0 0 1 0 1 1 % cess 0 0 0 1 0 0 0 1 0 1 1 % jd 0 0 1 0 0 0 0 0 1 2 0 % cen 1 0 0 0 0 1 0 0 0 2 0 % cfm 1 0 0 0 0 0 1 0 0 2 0 % % xFeM1 + xFeM2 % x -> ----------------------------- % xFeM1 + xFeM2 + xMgM1 + xMgM2 % % y -> 2 xAlT % % o -> xFeM2 + xMgM2 % % n -> xNaM2 % % xFeM1 % Q -> -x + ------------- % xFeM1 + xMgM1 % % f -> xFe3M1 % % cr -> xCrM1 % ------------------------------------------------- cpx 8 1 x(cpx) 0.075 y(cpx) 0.112 o(cpx) 0.05 n(cpx) 0.11 Q(cpx) -0.0005 range -1 1 f(cpx) 0.001 cr(cpx) 0.001 % ------------------------------------------------- p(di) 1 1 1 3 -1 n -1 o -1 y p(fs) 5 1 0 2 1 Q 1 x 2 0 1 -1 n 0 1 1 Q 2 0 1 -1 n 0 1 1 x 2 0 1 -1 Q 0 1 1 y 2 0 1 -1 x 0 1 1 y p(cats) 1 1 0 3 -1 cr -1 f 1 y p(crdi) 1 1 0 1 1 cr p(cess) 1 1 0 1 1 f p(jd) 1 1 0 1 1 n p(cen) 4 1 0 2 1 o 1 Q 2 0 1 -1 n 0 1 1 Q 2 0 1 -1 o 0 1 1 x 2 0 1 -1 Q 0 1 1 y p(cfm) 6 1 0 2 -1 x -2 Q 2 0 1 2 n 0 1 1 Q 2 0 1 1 n 0 1 1 x 2 0 1 1 o 0 1 1 x 2 0 1 2 Q 0 1 1 y 2 0 1 1 x 0 1 1 y % ------------------------------------------------- asf W(di,fs) 20 0 0 W(di,cats) 12.3 0 -0.1 W(di,crdi) 8 0 0 W(di,cess) 8 0 0 W(di,jd) 26 0 0 W(di,cen) 29.8 0 -0.03 W(di,cfm) 18 0 0 W(fs,cats) 25 0 -0.1 W(fs,crdi) 34 0 0 % Wfscrdi = 24 - x W(fs,cess) 34 0 0 % Wfscess = 24 - x W(fs,jd) 24 0 0 W(fs,cen) 7 0 0 W(fs,cfm) 4 0 0 W(cats,crdi) 2 0 0 W(cats,cess) 2 0 0 W(cats,jd) 6 0 0 W(cats,cen) 45.7 0 -0.29 W(cats,cfm) 27 0 -0.1 W(crdi,cess) 2 0 0 W(crdi,jd) 3 0 0 W(crdi,cen) 48 0 0 % Wcrdicen = 38 - x W(crdi,cfm) 36 0 0 % Wcrdicfm = 26 - x W(cess,jd) 3 0 0 W(cess,cen) 58 0 0 % Wcesscen = 48 - x W(cess,cfm) 36 0 0 % Wcesscfm = 26 - x W(jd,cen) 40 0 0 W(jd,cfm) 40 0 0 W(cen,cfm) 4 0 0 di 1.2 0 0 fs 1.0 0 0 cats 1.9 0 0 crdi 1.9 0 0 cess 1.9 0 0 jd 1.2 0 0 cen 1.0 0 0 cfm 1.0 0 0 % ------------------------------------------------- 11 xMgM1 5 1 1 4 -1 n -1 Q -1 x -1 y 2 0 1 1 n 0 1 1 Q 2 0 1 1 n 0 1 1 x 2 0 1 1 Q 0 1 1 y 2 0 1 1 x 0 1 1 y xFeM1 5 1 0 2 1 Q 1 x 2 0 1 -1 n 0 1 1 Q 2 0 1 -1 n 0 1 1 x 2 0 1 -1 Q 0 1 1 y 2 0 1 -1 x 0 1 1 y xAlM1 1 1 0 4 -1 cr -1 f 1 n 1 y xFe3M1 1 1 0 1 1 f xCrM1 1 1 0 1 1 cr xMgM2 4 1 0 2 1 o 1 Q 2 0 1 -1 n 0 1 1 Q 2 0 1 -1 o 0 1 1 x 2 0 1 -1 Q 0 1 1 y xFeM2 4 1 0 1 -1 Q 2 0 1 1 n 0 1 1 Q 2 0 1 1 o 0 1 1 x 2 0 1 1 Q 0 1 1 y xCaM2 1 1 1 2 -1 n -1 o xNaM2 1 1 0 1 1 n xSiT 1 1 1 1 -1/2 y xAlT 1 1 0 1 1/2 y % ------------------------------------------------- di 1 3 xMgM1 1 xCaM2 1 xSiT 1/2 check 0 0 0 0 0 0 0 fs 1 3 xFeM1 1 xFeM2 1 xSiT 1/2 check 1 0 1 0 0 0 0 make 1 fs 1 DQF 3.8 -0.003 0.03 % Lindsley 1965 cats 1.4142 4 xAlM1 1 xCaM2 1 xSiT 1/4 xAlT 1/4 check 0 1 0 0 0 0 0 crdi 1.4142 4 xCrM1 1 xCaM2 1 xSiT 1/4 xAlT 1/4 check 0 1 0 0 0 0 1 make 3 cats 1 kos 1 jd -1 DQF -3 0 0 cess 1.4142 4 xFe3M1 1 xCaM2 1 xSiT 1/4 xAlT 1/4 check 0 1 0 0 0 1 0 make 3 cats 1 acm 1 jd -1 DQF -6 0 0 jd 1 3 xAlM1 1 xNaM2 1 xSiT 1/2 check 0 0 0 1 0 0 0 cen 1 3 xMgM1 1 xMgM2 1 xSiT 1/2 check 0 0 1 0 0 0 0 make 1 en 1 DQF 3.5 -0.0020 0.048 cfm 1 3 xMgM1 1 xFeM2 1 xSiT 1/2 check 1/2 0 1 0 -1/2 0 0 make 2 en 1/2 fs 1/2 DQF -3.0 0 0 % ================================================= % ================== fo-fa olivine ======================== ol 2 1 x(ol) 0.25 p(fo) 1 1 1 1 -1 x p(fa) 1 1 0 1 1 x sf W(fo,fa) 9.0 0 0 2 pfo 1 1 1 1 -1 x pfa 1 1 0 1 1 x % ideal mixing activities fo 1 1 pfo 2 fa 1 1 pfa 2 %__________________________________________________________________ % ================================================= % ------------------------------------------------- % spinel: FMACr % % % M % Fe Fe3 Al Mg Cr % sp 0 0 2 1 0 % hc 1 0 2 0 0 % mt 1 2 0 0 0 % pcr 0 0 0 1 2 % % xFeM % x -> ----------- % xFeM + xMgM % % 3 xFe3M % f -> ------- % 2 % % 3 xCrM % cr -> ------ % 2 % ------------------------------------------------- sp 4 1 % aluminous spinel x(sp) 0.3 f(sp) 0.003 cr(sp) 0.05 % ------------------------------------------------- p(sp) 1 1 1 2 -1 cr -1 x p(hc) 1 1 0 2 -1 f 1 x p(mt) 1 1 0 1 1 f p(pcr) 1 1 0 1 1 cr % ------------------------------------------------- sf % from Bryndzia & Wood 90 W(sp,hc) 4 0 0 % changed from zero to 4 W(sp,mt) 56 0 0 W(sp,pcr) 39 0 0 W(hc,mt) 32 0 0 W(hc,pcr) 27 0 0 W(mt,pcr) 36 0 0 % ------------------------------------------------- 5 xFeM 1 1 0 1 1/3 x xFe3M 1 1 0 1 2/3 f xAlM 1 1 2/3 2 -2/3 cr -2/3 f xMgM 1 1 1/3 1 -1/3 x xCrM 1 1 0 1 2/3 cr % ------------------------------------------------- sp 27/4 2 xAlM 2 xMgM 1 check 0 0 0 hc 27/4 2 xFeM 1 xAlM 2 make 1 herc 1 check 1 0 0 mt 27/4 2 xFeM 1 xFe3M 2 check 1 1 0 pcr 27/4 2 xMgM 1 xCrM 2 make 1 picr 1 check 0 0 1 % ================================================= % ------------------------------------------------- % chromite: FMACr % % % M % Fe Fe3 Al Mg Cr % sp 0 0 2 1 0 % hc 1 0 2 0 0 % mt 1 2 0 0 0 % pcr 0 0 0 1 2 % % xFeM % x -> ----------- % xFeM + xMgM % % 3 xFe3M % f -> ------- % 2 % % 3 xCrM % cr -> ------ % 2 % ------------------------------------------------- cm 4 1 % Cr-rich spinel x(cm) 0.4 f(cm) 0.01 cr(cm) 0.7 % ------------------------------------------------- p(sp) 1 1 1 2 -1 cr -1 x p(hc) 1 1 0 2 -1 f 1 x p(mt) 1 1 0 1 1 f p(pcr) 1 1 0 1 1 cr % ------------------------------------------------- sf % from Bryndzia & Wood 90 W(sp,hc) 4 0 0 % changed from zero to 4 W(sp,mt) 56 0 0 W(sp,pcr) 39 0 0 W(hc,mt) 32 0 0 W(hc,pcr) 27 0 0 W(mt,pcr) 36 0 0 % ------------------------------------------------- 5 xFeM 1 1 0 1 1/3 x xFe3M 1 1 0 1 2/3 f xAlM 1 1 2/3 2 -2/3 cr -2/3 f xMgM 1 1 1/3 1 -1/3 x xCrM 1 1 0 1 2/3 cr % ------------------------------------------------- sp 27/4 2 xAlM 2 xMgM 1 check 0 0 0 hc 27/4 2 xFeM 1 xAlM 2 make 1 herc 1 check 1 0 0 mt 27/4 2 xFeM 1 xFe3M 2 check 1 1 0 pcr 27/4 2 xMgM 1 xCrM 2 make 1 picr 1 check 0 0 1 % ================================================= % ------------------------------------------------- % gt: CFMASCr % garnet % % M1 M2 % Mg Fe Ca Al Cr Fe3 % py 3 0 0 2 0 0 % alm 0 3 0 2 0 0 % gr 0 0 3 2 0 0 % andr 0 0 3 0 0 2 % knor 3 0 0 0 2 0 % % xFeM1 % x -> ------------- % xFeM1 + xMgM1 % % c -> xCaM1 % % f -> xFe3M2 % % cr -> xCrM2 % ------------------------------------------------- g 5 1 x(g) 0.3 c(g) 0.2 f(g) 0.01 cr(g) 0.01 % ------------------------------------------------- p(py) 2 1 1 3 -1 c -1 cr -1 x 2 0 1 1 c 0 1 1 x p(alm) 2 1 0 1 1 x 2 0 1 -1 c 0 1 1 x p(gr) 1 1 0 2 1 c -1 f p(andr) 1 1 0 1 1 f p(knor) 1 1 0 1 1 cr % ------------------------------------------------- sf W(py,alm) 4 0 0.1 W(py,gr) 35 0 0.1 W(py,andr) 91 0.0017 0.032 W(py,knor) 2 0 0 W(alm,gr) 4 0 0.1 W(alm,andr) 60 0.0017 0.032 W(alm,knor) 6 0 0.01 W(gr,andr) 2 0 0 W(gr,knor) 47 -0.0338 0.221 W(andr,knor) 101 -0.0321 0.153 % ------------------------------------------------- 6 xMgM1 2 1 1 2 -1 c -1 x 2 0 1 1 c 0 1 1 x xFeM1 2 1 0 1 1 x 2 0 1 -1 c 0 1 1 x xCaM1 1 1 0 1 1 c xAlM2 1 1 1 2 -1 cr -1 f xCrM2 1 1 0 1 1 cr xFe3M2 1 1 0 1 1 f % ------------------------------------------------- py 1 2 xMgM1 3 xAlM2 2 check 0 0 0 0 alm 1 2 xFeM1 3 xAlM2 2 check 1 0 0 0 gr 1 2 xCaM1 3 xAlM2 2 check 0 1 0 0 andr 1 2 xCaM1 3 xFe3M2 2 check 0 1 1 0 knor 1 2 xMgM1 3 xCrM2 2 check 0 0 0 1 % ================================================= % ------------------------------------------------- % opx: NCFMASTO % orthopyroxene % % M1 M2 T % Mg Fe Al Fe3 Cr Mg Fe Ca Si Al % en 1 0 0 0 0 1 0 0 2 0 % fs 0 1 0 0 0 0 1 0 2 0 % fm 1 0 0 0 0 0 1 0 2 0 % odi 1 0 0 0 0 0 0 1 2 0 % mgts 0 0 1 0 0 1 0 0 1 1 % cren 0 0 0 0 1 1 0 0 1 1 % mess 0 0 0 1 0 1 0 0 1 1 % % xFeM1 + xFeM2 % x -> ----------------------------- % xFeM1 + xFeM2 + xMgM1 + xMgM2 % % y -> 2 xAlT % % c -> xCaM2 % % xFeM1 % Q -> -x + ------------- % xFeM1 + xMgM1 % % f -> xFe3M1 % % cr -> xCrM1 % ------------------------------------------------- opx 7 1 x(opx) 0.05 y(opx) 0.006 c(opx) 0.025 Q(opx) 0.032 range -1 1 f(opx) 0.001 cr(opx) 0.001 % ------------------------------------------------- p(en) 3 1 1 4 -1 c 1 Q -1 x -1 y 2 0 1 1 c 0 1 1 x 2 0 1 -1 Q 0 1 1 y p(fs) 3 1 0 2 1 Q 1 x 2 0 1 -1 Q 0 1 1 y 2 0 1 -1 x 0 1 1 y p(fm) 4 1 0 1 -2 Q 2 0 1 -1 c 0 1 1 x 2 0 1 2 Q 0 1 1 y 2 0 1 1 x 0 1 1 y p(odi) 1 1 0 1 1 c p(mgts) 1 1 0 3 -1 cr -1 f 1 y p(cren) 1 1 0 1 1 cr p(mess) 1 1 0 1 1 f % ------------------------------------------------- asf W(en,fs) 5.2 0 0 W(en,fm) 4 0 0 W(en,odi) 32.2 0 0.12 W(en,mgts) 13 0 -0.15 W(en,cren) 8 0 0 W(en,mess) 8 0 0 W(fs,fm) 4 0 0 W(fs,odi) 24 0 0 W(fs,mgts) 7 0 -0.15 W(fs,cren) 10 0 0 W(fs,mess) 10 0 0 W(fm,odi) 18 0 0 W(fm,mgts) 2 0 -0.15 W(fm,cren) 12 0 0 W(fm,mess) 12 0 0 W(odi,mgts) 75.4 0 -0.94 W(odi,cren) 30 0 0 % Wodicren = 40 + x W(odi,mess) 30 0 0 % Wodimess = 40 + x W(mgts,cren) 2 0 0 W(mgts,mess) 2 0 0 W(cren,mess) 2 0 0 en 1.0 0 0 fs 1.0 0 0 fm 1.0 0 0 odi 1.2 0 0 mgts 1.0 0 0 cren 1.0 0 0 mess 1.0 0 0 % ------------------------------------------------- 10 xMgM1 3 1 1 3 -1 Q -1 x -1 y 2 0 1 1 Q 0 1 1 y 2 0 1 1 x 0 1 1 y xFeM1 3 1 0 2 1 Q 1 x 2 0 1 -1 Q 0 1 1 y 2 0 1 -1 x 0 1 1 y xAlM1 1 1 0 3 -1 cr -1 f 1 y xFe3M1 1 1 0 1 1 f xCrM1 1 1 0 1 1 cr xMgM2 3 1 1 3 -1 c 1 Q -1 x 2 0 1 1 c 0 1 1 x 2 0 1 -1 Q 0 1 1 y xFeM2 3 1 0 2 -1 Q 1 x 2 0 1 -1 c 0 1 1 x 2 0 1 1 Q 0 1 1 y xCaM2 1 1 0 1 1 c xSiT 1 1 1 1 -1/2 y xAlT 1 1 0 1 1/2 y % ------------------------------------------------- en 1 3 xMgM1 1 xMgM2 1 xSiT 1/2 check 0 0 0 0 0 0 fs 1 3 xFeM1 1 xFeM2 1 xSiT 1/2 check 1 0 0 0 0 0 fm 1 3 xMgM1 1 xFeM2 1 xSiT 1/2 check 1/2 0 0 -1/2 0 0 make 2 en 1/2 fs 1/2 DQF -6.0 0 0 odi 1 3 xMgM1 1 xCaM2 1 xSiT 1/2 check 0 0 1 0 0 0 make 1 di 1 DQF -0.1 0.000211 0.005 mgts 1.4142 4 xAlM1 1 xMgM2 1 xSiT 1/4 xAlT 1/4 check 0 1 0 0 0 0 cren 1.4142 4 xCrM1 1 xMgM2 1 xSiT 1/4 xAlT 1/4 check 0 1 0 0 0 1 make 3 mgts 1 kos 1 jd -1 DQF 3 0 0 mess 1.4142 4 xFe3M1 1 xMgM2 1 xSiT 1/4 xAlT 1/4 check 0 1 0 0 1 0 make 3 mgts 1 acm 1 jd -1 DQF -15 0 0.15 *