Wed 13 Feb, 2002 datafile. Changes added since 1998 paper.


¥ Added syv, hlt, syvL and hltL endmembers,
	(sylvite, halite and their molten equivalents).
	Solids data from Hemingway & Robie 95. 
	Melting curves from Clarke 1959.
	
¥ cc q wo CO2, mag en fo CO2, ta en q H2O 
    brackets from Aranovich & Newton (Am Min)
    
¥ new melt endmember data to agree with enthalpy
    of melting (Holland & Powell 2001). 
       
¥ updated chum data using experiments of Pawley,
    pers comm. New alpha, K, V, S, H
    
¥ updated phl data using experiments of Aranovich
    & Newton (Am Min 1998).
    
¥ updated chl data using experiments of Pawley (99
    pers com) and Fockenberg 1995, with new K. 
    S raised by 20 J for all chlorite endmembers.
    
¥ updated Prl for new K (Pawley 99 pers com). Had
    to raise alpha to satisfy experiments.
    
¥ The entropies of Mg-chloritoid and Fe-chloritoid have been updated so that
   that the thermodynamic dataset reproduces the important and careful study
   of FeÐMg partitioning between chlorite and chloritoid by Vidal et al. (1999)
   at temperatures down to 300¡C.
   
¥ The properties of the eastonite end-member of biotite were also modified
   slightly from those in Holland & Powell (1998), to allow agreement between
   calculated and natural alumina compositions of biotites in regional metamorphic
   pelites, such as those from the Al2SiO5 triple point at Mt Moosilauke (Hodges
   & Spear 1982). This was accomplished by raising the entropy of eastonite
   slightly and deriving its enthalpy from the same constraints used in Holland
   & Powell (1998).



Example datafile for KFMASH phengite

% ============= Phengite in KFMASH (with vacancies) ========================
mu 4  
      x(mu) 0.35
      y(mu) 0.75    % 2 XAlT1
      z(mu) 0.02    % vacancy on A site

   p(mu)    1 1      0  1  1  y
   
   p(cel)   1 2      1  1 -1  x       1  2  -1  y  -1  z
   
   p(fcel)  1 2      0  1  1  x       1  2  -1  y  -1  z

   p(prl)   1 1      0  1  1  z

 
  asf       % asymmetric formalism  (van Laar)

  W(mu,cel)      0  0  0.2 
  W(mu,fcel)     0  0  0.2
  W(mu,prl)     20  0  0.2
  W(cel,fcel)    0  0  0    
  W(cel,prl)    25  0  0.2    
  W(fcel,prl)   25  0  0.2
  
  vmu   0.55  0  0       % asymmetry parameters
  vcel  0.55  0  0    
  vfcel 0.55  0  0    
  vprl  0.45  0  0    
 
7  x(Al,M2A)   1 1      0  2   1  y  1  z

   x(Mg,M2A)   1 2      1  1  -1  x       1  2  -1  y  -1  z
    
   x(Fe,M2A)   1 2      0  1   1  x       1  2  -1  y  -1  z
    
   x(Al,T1)    1 1      0  1  1/2 y
   
   x(Si,T1)    1 1      1  1 -1/2 y

   x(V,A)      1 1      0  1  1  z

   x(K,A)      1 1      1  1 -1  z

      mu     4  4      x(K,A) 1   x(Al,M2A) 1    x(Al,T1) 1    x(Si,T1)  1
             
      cel    1  3      x(K,A) 1   x(Mg,M2A) 1    x(Si,T1) 2
      
      fcel   1  3      x(K,A) 1   x(Fe,M2A) 1    x(Si,T1) 2

      prl    1  3      x(V,A) 1   x(Al,M2A) 1    x(Si,T1) 2
              
%__________________________________________________________________