PHEMGP


PHEMGP has moved to: bitbucket.org/inverseproblem/phemgp


A quick guide on how to generate tables for PHEMGP from Perple_X output

PHEMGP was originally intended to permit users to extract physical properties from Perple_X phase equilibrium calculations without learning how to run Perple_X; however, while it is relatively inexpensive to save the raw Perple_X output for a whole mantle calculation, summarizing the entire calculation in a table with the resolution required to model the upper mantle generates huge tables. For this reason, it is preferable for users to generate their own tables from Perple_X output in order to best meet the demands of their specific applications. To generate tables users must run the Perple_X program WERAMI. The essential steps in this process are:

1) Copy WERAMI by downloading Perple_X as described at WINDOWS, LINUX, or OSX. These guides describe how to download the entire package, you only need the program WERAMI (i.e., you can delete the other files/programs).

2) Obtain the i/o data for a Perple_X calculation. For example from raw_data_for_mantle_tables; this example is a series of whole mantle calculations made with the Stixrude and Lithgow-Bertelloni (2011) data base as a function of composition. Composition being defined in terms of a binary mixture of a basaltic and harzburgitic components. Each calculation is designated by a project name of the form hzmb_###, where ### is the weight fraction of the basaltic component multiplied by 1000. Save these files in the same directory as WERAMI.

3) Start WERAMI, enter the project name (e.g., hzmb_###), select computational mode #2, and specify property choice #38; an example of a similar calculation is illustrated here. The desired table will be written to the file hzmb_###_X.phm, where X is an integer that indexes the WERAMI calculation, once you delete the "_X" part of the file name, the file can be used as PHEMGP input. Experimenting first with computational mode #1, which permits you to interactively sample computational results, may give you a better idea of the programs capabilities.

Currently the following raw data for phemgp tables is available from this site:

1) raw_data_for_mantle_tables: a series of calculations made for mantle mineralogy with the Stixrude & Lithgow-Bertelloni (2011) data base as a function of mantle composition. Composition being defined in terms of a binary mixture of a basaltic and harzburgitic components. The calculations are designated by the project name hzmb_###, where ### is the weight fraction of the basaltic component multiplied by 1000. 

2) raw_data_for_crust_tables: a series of calculations made for crustal mineralogy with the Holland & Powell (1998, 2002 revision) data base as a function of bulk water content. The calculations are designated by the project name taylor#, where # is the weight % of water in the bulk composition.