Author: Eric Nuzzi
Project Java Webmaster: Glenn A. Richard
Center for High Pressure Research
Stony Brook University
The tr660 program was orignally written in Fortran 77 by Don Weidner of the Center for High Pressure Research (CHiPR). The tr660 Java code was created by Eric Nuzzi, a computer science undergraduate student at SUNY Stony Brook, as part of Project Java under the supervision of Glenn Richard. The current version of the program was developed using Java 2 version 1.2.
This region of the Earth has been identified seismically as demonstrating changes in seismic velocities that cannot be explained by compression due to burial alone. Laboratory studies have demonstrated changes in mineralogy (phase) that will occur in candidate chemical compositions. The program utilizes published phase diagrams for these reactions along with published properties of the relevant phases. Some estimates, described by WW are necessary to complete the data set. This program is designed to be used for research purposed to test the effects of varying chemical composition, temperature, and physical properties on the inferred physical properties of the Earth through these phase transformations.
Input of Data
The chemical system is modeled in terms of 5 oxides, MgO, FeO, CaO, SiO2, and Al2O3. For most models of the Earth, these oxides constitute over 95% of the bulk composition. Variables XMG, XFE, XCA, XSI, XAL represent the number of the respective cation in the calculation. The appropriate number of oxygens are added automatically for charge neutrality. The absolute value of these variables is not important, they must simply be in the appropriate ratio that the user wishes. The default value is appropriate for the pyrolite model as discussed in WW. All of these variables can be changed by the user to test alternate models.
The phases that are included in this calculation are:
The properties of the phases are given in a table. The variables correspond to the following:
All values are the room pressure, temperature values
These parameters are set at default values as defined by the paper WW. Users can change any of the values.
Iterations indicates the total number of iterations or pressure steps for the calculation
T is the foot temperature corresponding to depth of the starting pressure P. 17 GPa is set as the default which corresponds to the depth of 500 km.
Output of Results
The program calculates the stable phases at each P,T condition along an adiabatic temperature gradient, starting at P, T for "Iterations steps of DTAP in pressure. At each condition, it then calculates the composition of each phase using the data of WW and it calculates the volume per cents of each phase. Then it calculates the properties of the phases and finally calculates the properties of the bulk mineralogy. Output is graphical and tabular. The graphing output is versatile allowing the user to view any of the calculated variables as a function of depth. The output can be imported to a file in the users platform for further analysis.
The opening screen of the program (data screen) contains the starting parameters from which the results will be generated. Clicking the Run tr660 button begins the calculations, and loads the results screen (results 1), containing the graph. At this point, the user may wish to resize the graph.
The Java-related images and names on this site are copyrights, trademarks, or registered trademarks of Sun Microsystems, Incorporated/JavaSoft. Please read JavaSoft's Copyright information.