%======================================================
% AX file for NCFMASCrO involving ol cpx opx g sp pl liq
% Jennings & Holland (2015)
%======================================================



% ==================   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
  DQF equilibrium 0 0 0
  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
  DQF equilibrium 0 0 0
  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
  DQF equilibrium 0 0 0
  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
  DQF equilibrium 0 0 0
  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  
  
  
*

