This directory contains the source code for the Fe-Ti oxide geothermometer
described in:

Ghiorso MS, Sack RO (1991) Fe-Ti oxide Geothermometry: Thermodynamic 
   Formulation and the Estimation of Intensive Variables in Silicic Magmas.
   Contrib. Mineral. Petrol. (in press)

To get the program working:

(1a) Copy the source file GEOTHERMOMETER.FOR to your machine.
(1b) Copy the BinHex file geothermometer-sit.hqx to your Macintosh.
     Use Stuffit (vers 1.5 or latter) to de-BinHex the file and
                 unstuff the various files. You will get all the
                 source files, a THINK C project library and two versions
                 of the executable (one for co-processor machines and
                 one for non-co-processor machines)
     Proceed to test the code with the executables provided, or...

(2) Compile and link the code.

(3) Run the test case indicated below.

Example (Vax 3200 running VMS 5.4):

$ fort geothermometer
$ link geothermometer
$ run geothermometer
***** CUBIC OXIDE PHASE*****
Input Wt % SiO2  in spinel [default:   0.00]: .03
Input Wt % TiO2  in spinel [default:   0.00]: 9.2
Input Wt % Al2O3 in spinel [default:   0.00]: 1.13
Input Wt % V2O3  in spinel [default:   0.00]: .12
Input Wt % Cr2O3 in spinel [default:   0.00]: .01
Input Wt % FeO   in spinel [default:   0.00]: 83.29
Input Wt % MnO   in spinel [default:   0.00]: .56
Input Wt % MgO   in spinel [default:   0.00]: .60
Input Wt % CaO   in spinel [default:   0.00]: .01
Input Wt % ZnO   in spinel [default:   0.00]: .17
Input Wt % NiO   in spinel [default:   0.00]: 0

***** RHOMBOHEDRAL OXIDE PHASE*****
Input Wt % SiO2  in rhom ox [default:   0.00]:
Input Wt % TiO2  in rhom ox [default:   0.00]: 47.26
Input Wt % Al2O3 in rhom ox [default:   0.00]: .11
Input Wt % V2O3  in rhom ox [default:   0.00]:
Input Wt % Cr2O3 in rhom ox [default:   0.00]:
Input Wt % FeO   in rhom ox [default:   0.00]: 49.51
Input Wt % MnO   in rhom ox [default:   0.00]: 1.05
Input Wt % MgO   in rhom ox [default:   0.00]: 1.20
Input Wt % CaO   in rhom ox [default:   0.00]:
Input Wt % ZnO   in rhom ox [default:   0.00]: .05
Input Wt % NiO   in rhom ox [default:   0.00]:

*****CALCULATING THE TEMPERATURE*****

*****CALCULATING THE FO2*****

*******************************
T (C):  758.55 log f O2: -14.57
*******************************

Do another calculation (Y or N)? n
FORTRAN STOP
$

If you get a different answer, or you have a problem, contact:

Mark S. Ghiorso
Department of Geological Sciences, AJ-20
University of Washington
Seattle, WA 98l95

Telephone: (206)-685-2482
Fax:       (206)-543-3836
e-mail:    ghiorso@fondue.geology.washington.edu  (internet)

The code has been tested on a VAX 3200 using the DIGITAL FORTRAN compiler
and on a Zenith 386-SX running Microsoft FORTRAN. The C version has been
tested on a Macintosh SE (no co-processor, i.e. very, very, very slow) 
and on a Macintosh SE/30 and Mac IIx (with co-processor, acceptably fast).

*END-OF-FILE README*