!     This session illustrates the component transformations necessary
!     to compute the "classical Thompson AFM diagram".

!     In the classical AFM projection, compositions are reduced through
!     the projection heirarchy H2O, SIO2 (q), KAl3O5 (mu) into the 
!     subcomposition MGO-FEO-AL2O3. The difficulty with this projection
!     is that biotite has a negative composition and would normally be
!     excluded by VERTEX. This problem is circumvented here by defining
!     two components, phl (=3 MgO - Al2O3) and ann (=3 Feo - Al2O3), to 
!     replace the normal FeO and MgO components so that Biotite plots
!     inside the composition space considered by Vertex. 

!     Note that the user must also define a muscovite component 
!     (mu = 0.5 K2O - 1.5 Al2O3, after projection through H2O and SiO2) so 
!     that the system can be saturated with respect to this component.
!     This has the result that K-feldpar phases, which have negative 
!     negative compositions along the KAl3O5-SiO2 join, will be rejected.
!     Consequently, vertex will not recognize any thermal limit for 
!     muscovite stability and the calculations will only be valid within
!     the stability field for Mu+Q.

!     run build to create input file in5.dat for vertex:

ostrich{jamie}84: build


 NO is the default answer to all Y/N prompts


 Enter the name of the input file for Vertex (i.e. the file
 OUTPUT by this program) < 15 characters and left justified:

in5.dat

 The file in5.dat        exists, do you really want to overwrite it (y/n)? 
y

 Enter the thermo data file name (e.g. hp90ver.dat),
 < 15 characters, left justified: 
hp90ver.dat
 Do you want to generate a   print  file (Y/N)? 
y
 Do you want to generate a graphics file (Y/N)? 
y
 Enter a name for the   print  file, <15 characters, left justified: 
print5.out
 Enter a name for the graphics file, <15 characters, left justified: 
plot5.out

 Specify the kind of phase diagram calculation:

  0 - for a Composition diagram
  1 - for a Schreinemakers-type diagram
  3 - for a Mixed-variable diagram

0
 Print dependent potentials for chemographies?
 Answer no if you do not know what this means.

n

 The data base components are:

  NA2O  MGO   AL2O3 SIO2  K2O   CAO   TIO2  MNO   FEO   O     H2O   CO2  

 Do you want to redefine them (Y/N)? 
y

 Enter the new component name (< 5 characters left justified, blank to finish):

mu

 Enter the component you wish to replace with mu   :
K2O
 Enter the stoichiometric coefficient of NA2O  in component mu   :
0
 Enter the stoichiometric coefficient of MGO   in component mu   :
0
 Enter the stoichiometric coefficient of AL2O3 in component mu   :
1.5
 Enter the stoichiometric coefficient of SIO2  in component mu   :
0
 Enter the stoichiometric coefficient of K2O   in component mu   :
0.5
 Enter the stoichiometric coefficient of CAO   in component mu   :
0
 Enter the stoichiometric coefficient of TIO2  in component mu   :
0
 Enter the stoichiometric coefficient of MNO   in component mu   :
0
 Enter the stoichiometric coefficient of FEO   in component mu   :
0
 Enter the stoichiometric coefficient of O     in component mu   :
0
 Enter the stoichiometric coefficient of H2O   in component mu   :
0
 Enter the stoichiometric coefficient of CO2   in component mu   :
0

 mu    = +  0.00 NA2O +  0.00 MGO  +  1.50 AL2O3+  0.00 SIO2 +  0.50 K2O  +  0.00 CAO  
         +  0.00 TIO2 +  0.00 MNO  +  0.00 FEO  +  0.00 O    +  0.00 H2O  +  0.00 CO2  

 Is this correct (Y/N)? 
y

 Enter the new component name (< 5 characters left justified, blank to finish):

phl

 Enter the component you wish to replace with phl  :
MGO
 Enter the stoichiometric coefficient of NA2O  in component phl  :
0
 Enter the stoichiometric coefficient of MGO   in component phl  :
3
 Enter the stoichiometric coefficient of AL2O3 in component phl  :
-1
 Enter the stoichiometric coefficient of SIO2  in component phl  :
0
 Enter the stoichiometric coefficient of mu    in component phl  :
0
 Enter the stoichiometric coefficient of CAO   in component phl  :
0
 Enter the stoichiometric coefficient of TIO2  in component phl  :
0
 Enter the stoichiometric coefficient of MNO   in component phl  :
0
 Enter the stoichiometric coefficient of FEO   in component phl  :
0
 Enter the stoichiometric coefficient of O     in component phl  :
0
 Enter the stoichiometric coefficient of H2O   in component phl  :
0
 Enter the stoichiometric coefficient of CO2   in component phl  :
0

 phl   = +  0.00 NA2O +  3.00 MGO  + -1.00 AL2O3+  0.00 SIO2 +  0.00 mu   +  0.00 CAO  
         +  0.00 TIO2 +  0.00 MNO  +  0.00 FEO  +  0.00 O    +  0.00 H2O  +  0.00 CO2  

 Is this correct (Y/N)? 
y

 Enter the new component name (< 5 characters left justified, blank to finish):

ann

 Enter the component you wish to replace with ann  :
FEO
 Enter the stoichiometric coefficient of NA2O  in component ann  :
0
 Enter the stoichiometric coefficient of phl   in component ann  :
0
 Enter the stoichiometric coefficient of AL2O3 in component ann  :
-1
 Enter the stoichiometric coefficient of SIO2  in component ann  :
0
 Enter the stoichiometric coefficient of mu    in component ann  :
0
 Enter the stoichiometric coefficient of CAO   in component ann  :
0
 Enter the stoichiometric coefficient of TIO2  in component ann  :
0
 Enter the stoichiometric coefficient of MNO   in component ann  :
0
 Enter the stoichiometric coefficient of FEO   in component ann  :
3
 Enter the stoichiometric coefficient of O     in component ann  :
0
 Enter the stoichiometric coefficient of H2O   in component ann  :
0
 Enter the stoichiometric coefficient of CO2   in component ann  :
0

 ann   = +  0.00 NA2O +  0.00 phl  + -1.00 AL2O3+  0.00 SIO2 +  0.00 mu   +  0.00 CAO  
         +  0.00 TIO2 +  0.00 MNO  +  3.00 FEO  +  0.00 O    +  0.00 H2O  +  0.00 CO2  

 Is this correct (Y/N)? 
y

 Enter the new component name (< 5 characters left justified, blank to finish):



 Do you want to do calculations with a saturated phase (Y/N)?
 The phase is:   FLUID 
 Its components can be: H2O   CO2  
 Its compositional variable is called: X(CO2)  

y

 Enter the number of components in the   FLUID 
 (1 or 2 for COH buffered fluids): 

1
 Enter the name of the component, left justified: 
H2O

 Do you want to do calculations with saturated components (Y/N)? 
y

 Select saturated components from the following set:

 NA2O  phl   AL2O3 SIO2  mu    CAO   TIO2  MNO   ann   O     CO2  

 How many saturated components (maximum 3)? 
2

 **warning ver015** the order you enter the components determines the
 saturation heirarchy and will effect your results (see Connolly 1990).

 Enter component names, left justified, one per line:

SIO2
mu

 Do you want to treat the chemical potential of a
 component as an INDEPENDENT variable (Y/N)? 
n

 Select remaining components from the following set:

 NA2O  phl   AL2O3 CAO   TIO2  MNO   ann   O     CO2  

 How many thermodynamic components (minimum 2, maximum 7)? 
3
 Enter component names, left justified, one per line:

ann 
AL2O3
phl

 working...


 Do you want to exclude phases from your calculations (Y/N)? 
n

 Select the equation of state for the   FLUID  phase

  1 - MRK (DeSantis et al 1974)
  2 - HSMRK (Kerrick and Jacobs 1981)
  3 - Hybrid MRK-HSMRK
  4 - Saxena and Fei 1987, pseudo-virial expansion
  5 - Bottinga and Richet 1982, RK
  6 - Holland and Powell 1990, CORK
  7 - Hybrid Haar/HSMRK 
  8 - Graphite buffered COH-MRK fluid

6

 Do you want to treat solution phases (Y/N)? 
y
 Enter the name of the file which contains the solution models you want to
 use (e.g.solut.dat), left justified, < 15 characters: 
solut.dat
 
 **warning ver025** missing all endmembers for GlTr(i)   
 **warning ver025** missing all endmembers for Tschermaki
 **warning ver025** missing all endmembers for Pl(i)     
 **warning ver025** missing all endmembers for MaPa(i)   
 **warning ver025** missing all endmembers for Pl(hB)    
 **warning ver025** missing all endmembers for Pl(h)     
 **warning ver025** missing all endmembers for DiCats(i) 
 **warning ver026** only one endmember for Ab(h)     
 **warning ver026** only one endmember for Kf(h)     
 **warning ver026** only one endmember for Ab        
 **warning ver026** only one endmember for Kf        
 **warning ver026** only one endmember for Kspar(h)  
 **warning ver026** only one endmember for Kspar     
 **warning ver025** missing all endmembers for Sanidine  
 **warning ver025** missing all endmembers for Omph(i)   
 Fluid      is not a possible solution phase because of endmember H2O     
 **warning ver025** missing all endmembers for Tr(i)     
 **warning ver026** only one endmember for Mica      
 **warning ver026** only one endmember for Pa        
 **warning ver026** only one endmember for Mu        
 **warning ver026** only one endmember for gralad    
 **warning ver025** missing all endmembers for GrAd(E&W) 
 **warning ver025** missing all endmembers for GrAd(i)   
 **warning ver025** missing all endmembers for EpCz(i)   
 **warning ver025** missing all endmembers for h-EpCz(i) 
 **warning ver025** missing all endmembers for HeDi(i)   
 **warning ver025** missing all endmembers for Dol(i)    
 **warning ver026** only one endmember for Neph(F&B) 

 ** warning ver111 ** the following endmembers are missing for solution Gt(i-c)   

  gr      

 it will be treated as a simpler solution between the remaining endmembers


 ** warning ver111 ** the following endmembers are missing for solution Gt(i-b)   

  gr      

 it will be treated as a simpler solution between the remaining endmembers

 **warning ver025** missing all endmembers for tr-ed     
 **warning ver025** missing all endmembers for mgal-amph 
 **warning ver025** missing all endmembers for mg-fe-hb  

 ** warning ver111 ** the following endmembers are missing for solution GrPyAlSp(B

  gr      

 it will be treated as a simpler solution between the remaining endmembers


 In solution GrPyAlSp(B the following endmembers have fixed compositions:

  spess    X = 0.100000    

 To relax or change these constraints you must modify solut.dat     .


 ** warning ver111 ** the following endmembers are missing for solution GrPyAl(B) 

  gr      

 it will be treated as a simpler solution between the remaining endmembers

 **warning ver025** missing all endmembers for JdDi(G1)  
 **warning ver025** missing all endmembers for JdDi(G2)  
 **warning ver025** missing all endmembers for JdDi(W?)  

 Select phases from the following list, enter the
 names one at a time and left justified, ENTER A BLANK WHEN YOU ARE doNE.

  Chl(i)      Bio(i)      St(i)       Carp(i)     Crd(i)      Cumm(i)   
  Anth(i)     EnFs(i)     Talc(i)     Ol(i)       Sp(J&R)     Sp(G&S)   
  Sp(i)       Gt(i-c)     Gt(i-b)     Ctd(i)      GrPyAlSp(B  GrPyAl(B) 

Chl(i)
Bio(i)
St(i)
Carp(i)
Crd(i)
Cumm(i)
Anth(i)
EnFs(i)
Talc(i)
Sp(J&R)
Ol(i)
Ctd(i)
GrPyAl(i)
GrPyAl(i)  is not a valid name, try again
GrPyAl(B)

 Calculate high variance phase fields (Y/N)?
n
 Enter a one-line title for your calculation:
test 5
 how many composition diagrams do you want to calculate?
6

 Specify the values for :
          P(bars)   T(K)    
 For calculation number 1
4000 600

 Specify the values for :
          P(bars)   T(K)    
 For calculation number 2
4000 700

 Specify the values for :
          P(bars)   T(K)    
 For calculation number 3
4000 800

 Specify the values for :
          P(bars)   T(K)    
 For calculation number 4
4000 920

 Specify the values for :
          P(bars)   T(K)    
 For calculation number 5
4000 1000

 Specify the values for :
          P(bars)   T(K)    
 For calculation number 6
4000 1100
 
ostrich{jamie}85: vertex
 Enter the name of the computational option file (i.e. the
 file created with BUILD) < 15 characters, left justified: 
in5.dat

 Reading thermodynamic data from file: hp90ver.dat   

 Writing print output to file: print5.out    

 Writing plot output to file: plot5.out     

 Reading solution models from file: solut.dat     


 **warning ver117** the subdivision scheme for the high order solution GrPyAl(B) 
 will be interpreted as a binary model with the parameters:
           imod= 1
           xmax=   1.000
           xmin=  10.000
           xinc=   9.000
 If you wish to change this or an  error (probably from routine CART)
 follows then you must modify the solution model file  (documentation section 4.1)
 
ostrich{jamie}86: psvdraw 

 Enter plot file name: 

plot5.out

 PostScript will be written to file: plot5.out.ps  

ostrich{jamie}87: pview plot5.out.ps
[1] 502
ostrich{jamie}88: