Julia-MELTS interface demo

This notebook runs the StatGeochem.jl package, which implements an interface for interacting with MELTS from the Julia programming language, including from Jupyter notebooks such as this.

StatGeochem.jl also includes some of the codes and utilities used in Keller & Schoene 2012, Keller et al. 2015 and Keller & Schoene 2018.

Launch Binder notebook

If running this notebook as an online Binder notebook and the webpage times out, click the badge at left to relaunch (refreshing will not work). Note that any changes will be lost!

Hint: shift-enter to run a single cell, or from the Cell menu select Run All to run the whole file. Any code from this notebook can be copied and pasted into the Julia REPL or a .jl script.


Load required Julia packages

In [1]:
## --- Load (and install if neccesary) the StatGeochem package
try
    using StatGeochem
catch
    using Pkg
    Pkg.add(PackageSpec(url="https://github.com/brenhinkeller/StatGeochem.jl"))
    using StatGeochem
end

using Statistics, StatsBase, DelimitedFiles
using Plots; gr();

AlphaMELTS calculations

Try to download and install alphaMELTS

In [2]:
# Download precompiled executable from caltech
alphameltsversion = "linux_alphamelts_1-9"
download("https://magmasource.caltech.edu/alphamelts/zipfiles/$alphameltsversion.zip","./$alphameltsversion.zip")
run(`unzip -o $alphameltsversion.zip`);
run(`mv $alphameltsversion/alphamelts_linux64 $alphameltsversion/alphamelts`);

# The absolute path to the alphamelts perl script must be input here:
  #  meltspath = "/usr/local/bin/run_alphamelts.command"
meltspath = pwd()*"/$alphameltsversion/run_alphamelts.command";

# The absolute or relative path to a direcory where you want to run MELTS -
# A number of configuration and output files will be written here.
  #  scratchdir = "scratch/"
scratchdir = pwd()*"/$alphameltsversion/";

Run alphamelts

In [3]:
## --- # # # # # # # # # # # pMelts equil. batch melting # # # # # # # # # # # #

# Conditions
P_range = [20000,20000]
T_range = [1700,800]
# Starting composition
elements = ["SiO2",  "TiO2","Al2O3","Fe2O3","Cr2O3",  "FeO",  "MnO",  "MgO",   "NiO",  "CoO",  "CaO",  "Na2O", "K2O", "P2O5", "H2O",]
composition=[44.8030, 0.1991, 4.4305, 0.9778, 0.3823, 7.1350, 0.1344, 37.6345, 0.2489, 0.0129, 3.5345, 0.3584, 0.0289, 0.0209, 0.15,] #mcdbse (McDonough Pyrolite)
# Run simulation
melts_configure(meltspath, scratchdir, composition, elements, T_range, P_range,
    batchstring="1\nsc.melts\n10\n1\n3\n1\nliquid\n1\n1.0\n0\n10\n0\n4\n0\n",
    dT=-10, dP=0, index=1, version="pMELTS",mode="isobaric",fo2path="FMQ")

# Read results
melt_comp = melts_query_liquid(scratchdir, index=1)
solid_comp = melts_query_solid(scratchdir, index=1)
modes = melts_query_modes(scratchdir, index=1)
ALPHAMELTS_VERSION pMELTS
ALPHAMELTS_DELTAT -10.0
ALPHAMELTS_MINP 2.0
ALPHAMELTS_MODE isobaric
ALPHAMELTS_MINT 800.0
ALPHAMELTS_CELSIUS_OUTPUT true
ALPHAMELTS_DELTAP 0.0
ALPHAMELTS_MAXP 90000.0
ALPHAMELTS_MAXT 3000.0
Checking for updates...


*** alphaMELTS 1.2 -- pMELTS isobar w/ or w/o liquid ***

This front end is the work of Paul Asimow and Paula Antoshechkina
(nee Smith) and it uses the MELTS and pMELTS algorithms developed by
Mark Ghiorso & co-workers. You are welcome to use and distribute this
program, under the condition that you acknowledge all the contributors
by citing the appropriate references with any results:

Smith & Asimow (2005) -- front end
Ghiorso et al (2002) -- pMELTS thermodynamic model
Asimow & Ghiorso (1998) -- subsolidus calculations (if used)

There may be others, depending on the exact calculation performed
and the source and / or partition coefficients (e.g. D(P,T,X)) used.
See Smith & Asimow (2005), documentation and the forum for details.

Unable to automatically check for updates; check for updates at
http://magmasource.caltech.edu/alphamelts/ and
http://magmasource.caltech.edu/forum/

Choose:
 1. Read MELTS file to set composition of system
 2. Twiddle starting or continuation parameters
 3. Single (batch) calculation
 4. Execute (follow path, mineral isograd or melt contour)
 5. Set fO2 buffer
 6. Set H2O (ppm) or aH2O
 7. Impose initial entropy, enthalpy or volume
 8. Adjust solid phase setting(s)
 9. Turn liquid on / off
10. Turn phase diagram mode on / off
11. Update state using restart file
12. Update composition using MELTS file
13. Write out restart file
14. Write out MELTS file
15. Write thermodynamic output for all phases
16. Calculate integrated melt and output file(s)
17. Fit parental melt composition (amoeba)
18. Cumulate Invertor (not yet implemented)
-1. Turn off menu display for options 1-18
 0. QUIT
Your choice: MELTS filename: input file open
Successfully read 'initial composition: sio2 44.7801'
Successfully read 'initial composition: tio2 0.1989'
Successfully read 'initial composition: al2o3 4.4282'
Successfully read 'initial composition: fe2o3 0.9773'
Successfully read 'initial composition: cr2o3 0.3821'
Successfully read 'initial composition: feo 7.1313'
Successfully read 'initial composition: mno 0.1343'
Successfully read 'initial composition: mgo 37.6152'
Successfully read 'initial composition: nio 0.2487'
Successfully read 'initial composition: coo 0.0128'
Successfully read 'initial composition: cao 3.5326'
Successfully read 'initial composition: na2o 0.3582'
Successfully read 'initial composition: k2o 0.0288'
Successfully read 'initial composition: p2o5 0.0208'
Successfully read 'initial composition: h2o 0.1499'
Successfully read 'initial temperature: 1700.0'
Successfully read 'initial pressure: 20000.0'
Successfully read 'log fo2 path: fmq'

Choose:
 1. Read MELTS file to set composition of system
 2. Twiddle starting or continuation parameters
 3. Single (batch) calculation
 4. Execute (follow path, mineral isograd or melt contour)
 5. Set fO2 buffer
 6. Set H2O (ppm) or aH2O
 7. Impose initial entropy, enthalpy or volume
 8. Adjust solid phase setting(s)
 9. Turn liquid on / off
10. Turn phase diagram mode on / off
11. Update state using restart file
12. Update composition using MELTS file
13. Write out restart file
14. Write out MELTS file
15. Write thermodynamic output for all phases
16. Calculate integrated melt and output file(s)
17. Fit parental melt composition (amoeba)
18. Cumulate Invertor (not yet implemented)
-1. Turn off menu display for options 1-18
 0. QUIT
Your choice: Turn phase diagram mode on (1) or off (0)? 
Choose:
 1. Read MELTS file to set composition of system
 2. Twiddle starting or continuation parameters
 3. Single (batch) calculation
 4. Execute (follow path, mineral isograd or melt contour)
 5. Set fO2 buffer
 6. Set H2O (ppm) or aH2O
 7. Impose initial entropy, enthalpy or volume
 8. Adjust solid phase setting(s)
 9. Turn liquid on / off
10. Turn phase diagram mode on / off
11. Update state using restart file
12. Update composition using MELTS file
13. Write out restart file
14. Write out MELTS file
15. Write thermodynamic output for all phases
16. Calculate integrated melt and output file(s)
17. Fit parental melt composition (amoeba)
18. Cumulate Invertor (not yet implemented)
-1. Turn off menu display for options 1-18
 0. QUIT
Your choice: Superliquidus (1) or subsolidus (0) initial guess ? Phase to track boundary of (by name, lower case): Type of melt contour to track:
0. Phi (melt fraction by volume)
1. F (melt fraction by mass)
2. aH2O (activity of water in the melt)
Choose: Type the F value to set (or < 0.0 for default): Use special monotonic search for troublesome phases (1) or quick search (0)? Initial Guess (not an isograd solution): P 20000.000000, T 1700.000000
liquid: 99.959 g 44.80 0.20 4.43 0.57 0.38 7.50 0.13 37.63 0.25 0.01 3.53 0.36 0.03 0.02 0.15 
Activity of H2O = 0.00111967  Melt fraction = 1
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
Isograd solution at: P 20000.000000, T 1624.070312
liquid: 99.957 g 44.80 0.20 4.43 0.56 0.38 7.51 0.13 37.63 0.25 0.01 3.53 0.36 0.03 0.02 0.15 
Activity of H2O = 0.00119344  Melt fraction = 1

Choose:
 1. Read MELTS file to set composition of system
 2. Twiddle starting or continuation parameters
 3. Single (batch) calculation
 4. Execute (follow path, mineral isograd or melt contour)
 5. Set fO2 buffer
 6. Set H2O (ppm) or aH2O
 7. Impose initial entropy, enthalpy or volume
 8. Adjust solid phase setting(s)
 9. Turn liquid on / off
10. Turn phase diagram mode on / off
11. Update state using restart file
12. Update composition using MELTS file
13. Write out restart file
14. Write out MELTS file
15. Write thermodynamic output for all phases
16. Calculate integrated melt and output file(s)
17. Fit parental melt composition (amoeba)
18. Cumulate Invertor (not yet implemented)
-1. Turn off menu display for options 1-18
 0. QUIT
Your choice: Turn phase diagram mode on (1) or off (0)? 
Choose:
 1. Read MELTS file to set composition of system
 2. Twiddle starting or continuation parameters
 3. Single (batch) calculation
 4. Execute (follow path, mineral isograd or melt contour)
 5. Set fO2 buffer
 6. Set H2O (ppm) or aH2O
 7. Impose initial entropy, enthalpy or volume
 8. Adjust solid phase setting(s)
 9. Turn liquid on / off
10. Turn phase diagram mode on / off
11. Update state using restart file
12. Update composition using MELTS file
13. Write out restart file
14. Write out MELTS file
15. Write thermodynamic output for all phases
16. Calculate integrated melt and output file(s)
17. Fit parental melt composition (amoeba)
18. Cumulate Invertor (not yet implemented)
-1. Turn off menu display for options 1-18
 0. QUIT
Your choice: ...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
...Adding the solid phase olivine to the assemblage.
...Dropping phase olivine from the assemblage.
Initial alphaMELTS calculation at: P 20000.000000, T 1624.070312
liquid: 99.957 g 44.80 0.20 4.43 0.56 0.38 7.51 0.13 37.63 0.25 0.01 3.53 0.36 0.03 0.02 0.15 
Activity of H2O = 0.00119344  Melt fraction = 1
...Adding the solid phase olivine to the assemblage.
alphaMELTS at: P 20000.000000, T 1614.070312
liquid: 91.072 g 45.12 0.22 4.86 0.58 0.42 7.80 0.12 36.23 0.17 0.01 3.85 0.39 0.03 0.02 0.16 
Activity of H2O = 0.00143887  Melt fraction = 0.911137
olivine: 8.882234 g, composition (Ca0.00Mg0.94Fe''0.05Mn0.00Co0.00Ni0.01)2SiO4
alphaMELTS at: P 20000.000000, T 1604.070312
liquid: 83.344 g 45.46 0.24 5.31 0.60 0.46 8.05 0.11 34.78 0.13 0.01 4.19 0.43 0.03 0.02 0.18 
Activity of H2O = 0.00171279  Melt fraction = 0.833847
olivine: 16.607219 g, composition (Ca0.00Mg0.94Fe''0.05Mn0.00Co0.00Ni0.01)2SiO4
...Adding the solid phase spinel to the assemblage.
alphaMELTS at: P 20000.000000, T 1594.070312
liquid: 76.418 g 45.90 0.26 5.76 0.62 0.41 8.25 0.10 33.32 0.10 0.01 4.54 0.47 0.04 0.03 0.20 
Activity of H2O = 0.00202505  Melt fraction = 0.764573
olivine: 23.396165 g, composition (Ca0.00Mg0.94Fe''0.05Mn0.00Co0.00Ni0.01)2SiO4
spinel: 0.134547 g, composition Fe''0.15Mg0.85Fe'''0.07Al0.75Cr1.17Ti0.00O4
alphaMELTS at: P 20000.000000, T 1584.070312
liquid: 70.600 g 46.34 0.28 6.20 0.63 0.38 8.41 0.09 31.90 0.08 0.01 4.89 0.51 0.04 0.03 0.21 
Activity of H2O = 0.00235381  Melt fraction = 0.706376
olivine: 29.111411 g, composition (Ca0.00Mg0.93Fe''0.06Mn0.00Co0.00Ni0.01)2SiO4
spinel: 0.235384 g, composition Fe''0.16Mg0.84Fe'''0.07Al0.77Cr1.15Ti0.00O4
...Adding the solid phase orthopyroxene to the assemblage.
...Dropping phase orthopyroxene from the assemblage.
alphaMELTS at: P 20000.000000, T 1574.070312
liquid: 65.643 g 46.79 0.30 6.63 0.64 0.35 8.51 0.08 30.53 0.07 0.01 5.24 0.55 0.04 0.03 0.23 
Activity of H2O = 0.00269699  Melt fraction = 0.656799
olivine: 33.979596 g, composition (Ca0.00Mg0.93Fe''0.06Mn0.00Co0.00Ni0.01)2SiO4
spinel: 0.321473 g, composition Fe''0.17Mg0.84Fe'''0.08Al0.79Cr1.12Ti0.00O4
...Adding the solid phase orthopyroxene to the assemblage.
alphaMELTS at: P 20000.000000, T 1564.070312
liquid: 55.506 g 45.68 0.35 7.43 0.69 0.32 9.11 0.08 29.32 0.06 0.01 5.96 0.64 0.05 0.04 0.27 
Activity of H2O = 0.00378705  Melt fraction = 0.555371
olivine: 36.012180 g, composition (Ca0.00Mg0.93Fe''0.06Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 7.960639 g, composition opx Na0.00Ca0.05Fe''0.10Mg1.79Fe'''0.01Ti0.00Al0.09Si1.95O6
spinel: 0.464989 g, composition Fe''0.17Mg0.84Fe'''0.08Al0.89Cr1.02Ti0.01O4
alphaMELTS at: P 20000.000000, T 1554.070312
liquid: 49.596 g 45.27 0.39 8.03 0.72 0.29 9.40 0.07 28.13 0.05 0.01 6.52 0.72 0.06 0.04 0.30 
Activity of H2O = 0.00472055  Melt fraction = 0.496244
olivine: 38.364714 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 11.410543 g, composition opx Na0.00Ca0.06Fe''0.10Mg1.78Fe'''0.01Ti0.00Al0.10Si1.95O6
spinel: 0.571432 g, composition Fe''0.17Mg0.84Fe'''0.08Al0.95Cr0.96Ti0.01O4
alphaMELTS at: P 20000.000000, T 1544.070312
liquid: 44.678 g 44.90 0.42 8.60 0.74 0.26 9.63 0.07 26.99 0.05 0.01 7.09 0.80 0.06 0.05 0.34 
Activity of H2O = 0.00577711  Melt fraction = 0.447037
olivine: 40.307067 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 14.282365 g, composition opx Na0.00Ca0.06Fe''0.11Mg1.77Fe'''0.01Ti0.00Al0.11Si1.94O6
spinel: 0.674624 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.01Cr0.89Ti0.01O4
alphaMELTS at: P 20000.000000, T 1534.070312
liquid: 40.507 g 44.57 0.46 9.14 0.76 0.24 9.80 0.06 25.90 0.04 0.01 7.66 0.87 0.07 0.05 0.37 
Activity of H2O = 0.00696605  Melt fraction = 0.405313
olivine: 41.933870 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 16.722170 g, composition opx Na0.00Ca0.07Fe''0.11Mg1.76Fe'''0.01Ti0.00Al0.11Si1.94O6
spinel: 0.777506 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.08Cr0.83Ti0.01O4
alphaMELTS at: P 20000.000000, T 1524.070312
liquid: 36.917 g 44.27 0.49 9.64 0.78 0.21 9.93 0.06 24.85 0.04 0.00 8.22 0.96 0.08 0.06 0.41 
Activity of H2O = 0.00829785  Melt fraction = 0.369393
olivine: 43.311147 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 18.829911 g, composition opx Na0.00Ca0.07Fe''0.12Mg1.75Fe'''0.01Ti0.00Al0.12Si1.93O6
spinel: 0.882041 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.14Cr0.77Ti0.01O4
alphaMELTS at: P 20000.000000, T 1514.070312
liquid: 33.788 g 44.00 0.52 10.11 0.80 0.19 10.01 0.06 23.84 0.04 0.00 8.80 1.04 0.09 0.06 0.44 
Activity of H2O = 0.00978407  Melt fraction = 0.338081
olivine: 44.486449 g, composition (Ca0.00Mg0.92Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 20.676180 g, composition opx Na0.00Ca0.07Fe''0.12Mg1.73Fe'''0.01Ti0.00Al0.13Si1.93O6
spinel: 0.989456 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.19Cr0.71Ti0.01O4
...Adding the solid phase clinopyroxene to the assemblage.
...Dropping phase clinopyroxene from the assemblage.
alphaMELTS at: P 20000.000000, T 1504.070312
liquid: 31.032 g 43.76 0.56 10.55 0.82 0.17 10.05 0.05 22.87 0.03 0.00 9.37 1.13 0.09 0.07 0.48 
Activity of H2O = 0.0114373  Melt fraction = 0.310505
olivine: 45.495020 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 22.312332 g, composition opx Na0.00Ca0.08Fe''0.12Mg1.72Fe'''0.01Ti0.00Al0.13Si1.93O6
spinel: 1.100424 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.25Cr0.66Ti0.01O4
...Adding the solid phase clinopyroxene to the assemblage.
...Dropping phase clinopyroxene from the assemblage.
alphaMELTS at: P 20000.000000, T 1494.070312
liquid: 28.583 g 43.54 0.59 10.96 0.83 0.15 10.05 0.05 21.93 0.03 0.00 9.94 1.22 0.10 0.07 0.52 
Activity of H2O = 0.0132713  Melt fraction = 0.286008
olivine: 46.363698 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 23.776873 g, composition opx Na0.00Ca0.08Fe''0.13Mg1.71Fe'''0.01Ti0.00Al0.14Si1.92O6
spinel: 1.215221 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.30Cr0.61Ti0.01O4
...Adding the solid phase clinopyroxene to the assemblage.
...Dropping phase clinopyroxene from the assemblage.
alphaMELTS at: P 20000.000000, T 1484.070312
liquid: 26.392 g 43.35 0.62 11.34 0.84 0.13 10.01 0.05 21.03 0.03 0.00 10.51 1.32 0.11 0.08 0.57 
Activity of H2O = 0.0153017  Melt fraction = 0.26408
olivine: 47.113450 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4
orthopyroxene: 25.099600 g, composition opx Na0.00Ca0.09Fe''0.13Mg1.70Fe'''0.01Ti0.00Al0.14Si1.92O6
spinel: 1.333875 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.34Cr0.56Ti0.01O4
...Adding the solid phase clinopyroxene to the assemblage.
...Dropping phase clinopyroxene from the assemblage.
alphaMELTS at: P 20000.000000, T 1474.070312
liquid: 24.417 g 43.18 0.66 11.70 0.85 0.12 9.95 0.04 20.16 0.03 0.00 11.06 1.42 0.12 0.09 0.61 
Activity of H2O = 0.0175466  Melt fraction = 0.244318
olivine: 47.761057 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 26.304323 g, composition opx Na0.00Ca0.09Fe''0.13Mg1.69Fe'''0.01Ti0.00Al0.15Si1.92O6
spinel: 1.456306 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.39Cr0.52Ti0.01O4
...Adding the solid phase clinopyroxene to the assemblage.
alphaMELTS at: P 20000.000000, T 1464.070312
liquid: 21.132 g 43.00 0.70 12.12 0.86 0.11 10.08 0.04 19.68 0.03 0.00 10.83 1.60 0.14 0.10 0.71 
Activity of H2O = 0.022457  Melt fraction = 0.211448
olivine: 48.588784 g, composition (Ca0.00Mg0.91Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 25.514705 g, composition opx Na0.00Ca0.09Fe''0.14Mg1.69Fe'''0.01Ti0.00Al0.15Si1.92O6
clinopyroxene: 3.049418 g, composition cpx Na0.02Ca0.57Fe''0.11Mg1.18Fe'''0.02Ti0.01Al0.21Si1.89O6
spinel: 1.653782 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.45Cr0.46Ti0.01O4
alphaMELTS at: P 20000.000000, T 1454.070312
liquid: 18.290 g 42.87 0.74 12.48 0.86 0.10 10.16 0.04 19.23 0.03 0.00 10.59 1.80 0.16 0.11 0.82 
Activity of H2O = 0.0285482  Melt fraction = 0.183011
olivine: 49.290361 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.942667 g, composition opx Na0.00Ca0.09Fe''0.14Mg1.68Fe'''0.01Ti0.00Al0.16Si1.91O6
clinopyroxene: 5.569546 g, composition cpx Na0.02Ca0.58Fe''0.11Mg1.16Fe'''0.02Ti0.01Al0.22Si1.89O6
spinel: 1.846146 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.49Cr0.42Ti0.01O4
alphaMELTS at: P 20000.000000, T 1444.070312
liquid: 15.881 g 42.80 0.78 12.81 0.86 0.09 10.19 0.04 18.78 0.03 0.00 10.34 2.01 0.18 0.13 0.94 
Activity of H2O = 0.0358609  Melt fraction = 0.158908
olivine: 49.880834 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.566195 g, composition opx Na0.00Ca0.08Fe''0.14Mg1.68Fe'''0.01Ti0.00Al0.17Si1.91O6
clinopyroxene: 7.586043 g, composition cpx Na0.02Ca0.59Fe''0.11Mg1.14Fe'''0.02Ti0.01Al0.22Si1.88O6
spinel: 2.024578 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.53Cr0.38Ti0.01O4
alphaMELTS at: P 20000.000000, T 1434.070312
liquid: 13.853 g 42.78 0.81 13.10 0.86 0.09 10.17 0.04 18.33 0.03 0.00 10.10 2.23 0.21 0.15 1.08 
Activity of H2O = 0.0444289  Melt fraction = 0.138614
olivine: 50.378410 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.332580 g, composition opx Na0.00Ca0.08Fe''0.15Mg1.68Fe'''0.01Ti0.01Al0.17Si1.91O6
clinopyroxene: 9.189705 g, composition cpx Na0.02Ca0.60Fe''0.12Mg1.12Fe'''0.02Ti0.01Al0.23Si1.88O6
spinel: 2.185259 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.55Cr0.35Ti0.01O4
alphaMELTS at: P 20000.000000, T 1424.070312
liquid: 12.149 g 42.80 0.83 13.38 0.86 0.09 10.11 0.04 17.89 0.03 0.00 9.86 2.46 0.24 0.17 1.23 
Activity of H2O = 0.0542477  Melt fraction = 0.121566
olivine: 50.798768 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.200332 g, composition opx Na0.00Ca0.08Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.17Si1.90O6
clinopyroxene: 10.463339 g, composition cpx Na0.03Ca0.61Fe''0.11Mg1.11Fe'''0.02Ti0.01Al0.24Si1.88O6
spinel: 2.327302 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.57Cr0.33Ti0.01O4
...Adding the solid phase garnet to the assemblage.
...Dropping phase garnet from the assemblage.
alphaMELTS at: P 20000.000000, T 1414.070312
liquid: 10.717 g 42.86 0.85 13.64 0.85 0.09 10.00 0.05 17.44 0.03 0.00 9.63 2.70 0.27 0.19 1.40 
Activity of H2O = 0.0652847  Melt fraction = 0.107232
olivine: 51.155082 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.138395 g, composition opx Na0.01Ca0.08Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6
clinopyroxene: 11.477210 g, composition cpx Na0.03Ca0.62Fe''0.11Mg1.09Fe'''0.02Ti0.01Al0.24Si1.87O6
spinel: 2.451633 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.59Cr0.32Ti0.01O4
...Adding the solid phase garnet to the assemblage.
...Dropping phase garnet from the assemblage.
alphaMELTS at: P 20000.000000, T 1404.070312
liquid: 9.508 g 42.94 0.86 13.90 0.85 0.08 9.86 0.05 16.99 0.03 0.00 9.40 2.95 0.30 0.22 1.58 
Activity of H2O = 0.0774864  Melt fraction = 0.0951415
olivine: 51.458212 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.124315 g, composition opx Na0.01Ca0.08Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6
clinopyroxene: 12.288111 g, composition cpx Na0.03Ca0.63Fe''0.11Mg1.08Fe'''0.02Ti0.01Al0.25Si1.87O6
spinel: 2.560042 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.60Cr0.30Ti0.01O4
...Adding the solid phase garnet to the assemblage.
...Dropping phase garnet from the assemblage.
alphaMELTS at: P 20000.000000, T 1394.070312
liquid: 8.485 g 43.05 0.86 14.15 0.84 0.08 9.68 0.05 16.54 0.03 0.00 9.16 3.20 0.34 0.25 1.77 
Activity of H2O = 0.0907855  Melt fraction = 0.0848974
olivine: 51.717009 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.142219 g, composition opx Na0.01Ca0.07Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6
clinopyroxene: 12.940628 g, composition cpx Na0.04Ca0.63Fe''0.11Mg1.07Fe'''0.02Ti0.01Al0.25Si1.87O6
spinel: 2.654583 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.62Cr0.29Ti0.01O4
...Adding the solid phase garnet to the assemblage.
...Dropping phase garnet from the assemblage.
alphaMELTS at: P 20000.000000, T 1384.070312
liquid: 7.613 g 43.17 0.85 14.41 0.82 0.09 9.48 0.05 16.08 0.03 0.00 8.94 3.46 0.38 0.27 1.97 
Activity of H2O = 0.105106  Melt fraction = 0.076174
olivine: 51.938661 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
orthopyroxene: 24.181036 g, composition opx Na0.01Ca0.07Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6
clinopyroxene: 13.469233 g, composition cpx Na0.04Ca0.64Fe''0.11Mg1.05Fe'''0.02Ti0.01Al0.26Si1.87O6
spinel: 2.737252 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.62Cr0.28Ti0.01O4
...Adding the solid phase garnet to the assemblage.
alphaMELTS at: P 20000.000000, T 1374.070312
liquid: 6.863 g 43.31 0.85 14.64 0.81 0.09 9.26 0.05 15.62 0.03 0.00 8.71 3.72 0.42 0.30 2.18 
Activity of H2O = 0.120429  Melt fraction = 0.0686709
olivine: 52.150643 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 0.089132 g, composition (Ca0.18Fe''0.14Mg0.68)3Al2Si3O12
orthopyroxene: 24.175149 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6
clinopyroxene: 13.870519 g, composition cpx Na0.04Ca0.65Fe''0.11Mg1.04Fe'''0.02Ti0.01Al0.26Si1.87O6
spinel: 2.790367 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.63Cr0.28Ti0.01O4
alphaMELTS at: P 20000.000000, T 1364.070312
liquid: 6.194 g 43.52 0.86 14.78 0.80 0.10 8.99 0.05 15.16 0.03 0.00 8.46 4.03 0.46 0.34 2.42 
Activity of H2O = 0.137078  Melt fraction = 0.0619795
olivine: 52.510356 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 0.903057 g, composition (Ca0.18Fe''0.14Mg0.68)3Al2Si3O12
orthopyroxene: 23.701151 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6
clinopyroxene: 13.953296 g, composition cpx Na0.04Ca0.65Fe''0.11Mg1.03Fe'''0.02Ti0.01Al0.26Si1.87O6
spinel: 2.675436 g, composition Fe''0.17Mg0.83Fe'''0.07Al1.62Cr0.29Ti0.01O4
alphaMELTS at: P 20000.000000, T 1354.070312
liquid: 5.622 g 43.74 0.87 14.92 0.78 0.11 8.70 0.06 14.70 0.03 0.00 8.22 4.33 0.51 0.37 2.67 
Activity of H2O = 0.154397  Melt fraction = 0.0562602
olivine: 52.826652 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 1.641180 g, composition (Ca0.18Fe''0.15Mg0.68)3Al2Si3O12
orthopyroxene: 23.273605 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6
clinopyroxene: 14.003803 g, composition cpx Na0.04Ca0.66Fe''0.11Mg1.02Fe'''0.02Ti0.01Al0.26Si1.87O6
spinel: 2.568454 g, composition Fe''0.18Mg0.83Fe'''0.07Al1.60Cr0.30Ti0.01O4
alphaMELTS at: P 20000.000000, T 1344.070312
liquid: 5.130 g 43.95 0.87 15.05 0.76 0.12 8.41 0.06 14.24 0.03 0.00 7.97 4.64 0.56 0.41 2.92 
Activity of H2O = 0.172301  Melt fraction = 0.051337
olivine: 53.106072 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 2.314755 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12
orthopyroxene: 22.884334 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.90O6
clinopyroxene: 14.030766 g, composition cpx Na0.05Ca0.67Fe''0.11Mg1.02Fe'''0.02Ti0.01Al0.26Si1.87O6
spinel: 2.468664 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.59Cr0.32Ti0.01O4
alphaMELTS at: P 20000.000000, T 1334.070312
liquid: 4.704 g 44.16 0.87 15.18 0.75 0.14 8.11 0.06 13.78 0.04 0.01 7.73 4.95 0.61 0.44 3.19 
Activity of H2O = 0.190714  Melt fraction = 0.0470717
olivine: 53.353971 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 2.932822 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12
orthopyroxene: 22.527017 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.90O6
clinopyroxene: 14.040591 g, composition cpx Na0.05Ca0.67Fe''0.11Mg1.01Fe'''0.02Ti0.01Al0.26Si1.87O6
spinel: 2.375397 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.57Cr0.33Ti0.01O4
alphaMELTS at: P 20000.000000, T 1324.070312
liquid: 4.333 g 44.35 0.86 15.30 0.73 0.15 7.81 0.06 13.33 0.04 0.01 7.49 5.25 0.66 0.48 3.46 
Activity of H2O = 0.209569  Melt fraction = 0.0433543
olivine: 53.574796 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 3.502746 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12
orthopyroxene: 22.196672 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.91O6
clinopyroxene: 14.038032 g, composition cpx Na0.05Ca0.68Fe''0.11Mg1.00Fe'''0.02Ti0.01Al0.26Si1.87O6
spinel: 2.288062 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.56Cr0.34Ti0.01O4
alphaMELTS at: P 20000.000000, T 1314.070312
liquid: 4.007 g 44.54 0.86 15.42 0.71 0.17 7.52 0.07 12.88 0.04 0.01 7.26 5.55 0.72 0.52 3.74 
Activity of H2O = 0.22881  Melt fraction = 0.0400961
olivine: 53.772272 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 4.030605 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12
orthopyroxene: 21.889309 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.91O6
clinopyroxene: 14.026655 g, composition cpx Na0.05Ca0.68Fe''0.10Mg0.99Fe'''0.02Ti0.02Al0.26Si1.87O6
spinel: 2.206131 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.54Cr0.36Ti0.01O4
alphaMELTS at: P 20000.000000, T 1304.070312
liquid: 3.720 g 44.71 0.85 15.54 0.69 0.19 7.23 0.07 12.45 0.04 0.01 7.03 5.85 0.77 0.56 4.03 
Activity of H2O = 0.248388  Melt fraction = 0.037225
olivine: 53.949552 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 4.521471 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12
orthopyroxene: 21.601681 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.17Si1.91O6
clinopyroxene: 14.009152 g, composition cpx Na0.06Ca0.69Fe''0.10Mg0.99Fe'''0.02Ti0.02Al0.26Si1.87O6
spinel: 2.129135 g, composition Fe''0.19Mg0.83Fe'''0.08Al1.53Cr0.37Ti0.01O4
alphaMELTS at: P 20000.000000, T 1294.070312
liquid: 3.466 g 44.87 0.84 15.64 0.67 0.21 6.94 0.07 12.02 0.04 0.01 6.80 6.13 0.83 0.60 4.33 
Activity of H2O = 0.268264  Melt fraction = 0.0346821
olivine: 54.109313 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 4.979616 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12
orthopyroxene: 21.331113 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6
clinopyroxene: 13.987567 g, composition cpx Na0.06Ca0.69Fe''0.10Mg0.98Fe'''0.02Ti0.02Al0.26Si1.88O6
spinel: 2.056649 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.51Cr0.38Ti0.01O4
alphaMELTS at: P 20000.000000, T 1284.070312
liquid: 3.240 g 45.01 0.83 15.74 0.65 0.23 6.66 0.08 11.61 0.04 0.01 6.58 6.41 0.89 0.64 4.63 
Activity of H2O = 0.288408  Melt fraction = 0.0324187
olivine: 54.253844 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 5.408668 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12
orthopyroxene: 21.075380 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6
clinopyroxene: 13.963453 g, composition cpx Na0.06Ca0.70Fe''0.10Mg0.97Fe'''0.02Ti0.02Al0.26Si1.88O6
spinel: 1.988296 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.50Cr0.40Ti0.01O4
alphaMELTS at: P 20000.000000, T 1274.070312
liquid: 3.037 g 45.13 0.81 15.83 0.63 0.25 6.39 0.08 11.20 0.04 0.01 6.37 6.69 0.95 0.68 4.94 
Activity of H2O = 0.308796  Melt fraction = 0.0303942
olivine: 54.385101 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 5.811732 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12
orthopyroxene: 20.832617 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6
clinopyroxene: 13.937990 g, composition cpx Na0.06Ca0.70Fe''0.10Mg0.97Fe'''0.02Ti0.02Al0.26Si1.88O6
spinel: 1.923733 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.48Cr0.41Ti0.01O4
alphaMELTS at: P 20000.000000, T 1264.070312
liquid: 2.855 g 45.22 0.80 15.91 0.61 0.27 6.13 0.08 10.81 0.05 0.01 6.17 6.95 1.01 0.73 5.25 
Activity of H2O = 0.32941  Melt fraction = 0.0285752
olivine: 54.504761 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.191479 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12
orthopyroxene: 20.601247 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6
clinopyroxene: 13.912073 g, composition cpx Na0.07Ca0.70Fe''0.10Mg0.96Fe'''0.02Ti0.02Al0.26Si1.88O6
spinel: 1.862651 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.47Cr0.43Ti0.01O4
alphaMELTS at: P 20000.000000, T 1254.070312
liquid: 2.691 g 45.30 0.78 15.98 0.60 0.30 5.88 0.09 10.43 0.05 0.01 5.98 7.20 1.07 0.77 5.57 
Activity of H2O = 0.350238  Melt fraction = 0.0269334
olivine: 54.614264 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.550223 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12
orthopyroxene: 20.379933 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.15Si1.91O6
clinopyroxene: 13.886375 g, composition cpx Na0.07Ca0.71Fe''0.10Mg0.96Fe'''0.02Ti0.02Al0.26Si1.88O6
spinel: 1.804774 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.45Cr0.44Ti0.01O4
alphaMELTS at: P 20000.000000, T 1244.070312
liquid: 2.543 g 45.35 0.77 16.05 0.58 0.32 5.64 0.09 10.07 0.05 0.01 5.79 7.44 1.13 0.82 5.90 
Activity of H2O = 0.37127  Melt fraction = 0.0254453
olivine: 54.714847 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.889977 g, composition (Ca0.16Fe''0.16Mg0.67)3Al2Si3O12
orthopyroxene: 20.167534 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.15Si1.91O6
clinopyroxene: 13.861396 g, composition cpx Na0.07Ca0.71Fe''0.10Mg0.95Fe'''0.02Ti0.02Al0.25Si1.88O6
spinel: 1.749849 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.44Cr0.46Ti0.01O4
alphaMELTS at: P 20000.000000, T 1234.070312
liquid: 2.407 g 45.39 0.76 16.10 0.56 0.35 5.41 0.09 9.72 0.05 0.01 5.62 7.67 1.20 0.86 6.23 
Activity of H2O = 0.392501  Melt fraction = 0.0240909
olivine: 54.807573 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.212501 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12
orthopyroxene: 19.963075 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.15Si1.92O6
clinopyroxene: 13.837503 g, composition cpx Na0.07Ca0.71Fe''0.09Mg0.95Fe'''0.02Ti0.02Al0.25Si1.88O6
spinel: 1.697649 g, composition Fe''0.20Mg0.82Fe'''0.08Al1.42Cr0.47Ti0.01O4
alphaMELTS at: P 20000.000000, T 1224.070312
liquid: 2.284 g 45.40 0.74 16.15 0.54 0.38 5.18 0.10 9.38 0.05 0.01 5.45 7.89 1.26 0.91 6.56 
Activity of H2O = 0.413926  Melt fraction = 0.0228534
olivine: 54.893359 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.519342 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12
orthopyroxene: 19.765716 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.15Si1.92O6
clinopyroxene: 13.814956 g, composition cpx Na0.08Ca0.72Fe''0.09Mg0.94Fe'''0.02Ti0.02Al0.25Si1.88O6
spinel: 1.647966 g, composition Fe''0.20Mg0.82Fe'''0.08Al1.41Cr0.49Ti0.01O4
alphaMELTS at: P 20000.000000, T 1214.070312
liquid: 2.170 g 45.38 0.73 16.20 0.53 0.41 4.97 0.10 9.05 0.06 0.01 5.29 8.09 1.33 0.96 6.91 
Activity of H2O = 0.435546  Melt fraction = 0.0217186
olivine: 54.972993 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.811862 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12
orthopyroxene: 19.574738 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.14Si1.92O6
clinopyroxene: 13.793936 g, composition cpx Na0.08Ca0.72Fe''0.09Mg0.94Fe'''0.02Ti0.02Al0.25Si1.88O6
spinel: 1.600614 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.39Cr0.50Ti0.01O4
...Adding the solid phase hornblende to the assemblage.
...Dropping phase hornblende from the assemblage.
alphaMELTS at: P 20000.000000, T 1204.070312
liquid: 2.066 g 45.35 0.71 16.24 0.51 0.44 4.77 0.11 8.74 0.06 0.01 5.14 8.29 1.39 1.01 7.26 
Activity of H2O = 0.457358  Melt fraction = 0.0206743
olivine: 55.047157 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.091268 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12
orthopyroxene: 19.389518 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.14Si1.92O6
clinopyroxene: 13.774559 g, composition cpx Na0.08Ca0.72Fe''0.09Mg0.93Fe'''0.02Ti0.02Al0.25Si1.89O6
spinel: 1.555422 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.37Cr0.52Ti0.01O4
...Adding the solid phase hornblende to the assemblage.
...Dropping phase hornblende from the assemblage.
...Adding the solid phase biotite to the assemblage.
alphaMELTS at: P 20000.000000, T 1194.070312
liquid: 1.951 g 45.19 0.69 16.24 0.50 0.48 4.61 0.11 8.49 0.06 0.01 5.04 8.50 1.36 1.07 7.64 
Activity of H2O = 0.48032  Melt fraction = 0.0195234
olivine: 55.101579 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.353334 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12
orthopyroxene: 19.217252 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.70Fe'''0.01Ti0.01Al0.14Si1.92O6
clinopyroxene: 13.768639 g, composition cpx Na0.08Ca0.73Fe''0.09Mg0.93Fe'''0.02Ti0.02Al0.25Si1.89O6
biotite: 0.019341 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.512253 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.36Cr0.53Ti0.01O4
...Adding the solid phase hornblende to the assemblage.
...Dropping phase hornblende from the assemblage.
alphaMELTS at: P 20000.000000, T 1184.070312
liquid: 1.826 g 44.91 0.67 16.21 0.49 0.55 4.50 0.11 8.33 0.06 0.01 4.98 8.75 1.22 1.14 8.07 
Activity of H2O = 0.504807  Melt fraction = 0.0182712
olivine: 55.138126 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.599245 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12
orthopyroxene: 19.055767 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.70Fe'''0.01Ti0.01Al0.14Si1.92O6
clinopyroxene: 13.775658 g, composition cpx Na0.09Ca0.73Fe''0.09Mg0.92Fe'''0.02Ti0.02Al0.25Si1.89O6
biotite: 0.057551 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.470694 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.34Cr0.55Ti0.01O4
...Adding the solid phase hornblende to the assemblage.
...Adding the solid phase apatite to the assemblage.
alphaMELTS at: P 20000.000000, T 1174.070312
liquid: 0.688 g 44.57 0.70 16.01 0.50 0.62 4.72 0.13 8.34 0.07 0.01 5.05 8.26 1.02 1.06 8.95 
Activity of H2O = 0.561124  Melt fraction = 0.00688762
olivine: 54.216269 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.034694 g, composition (Ca0.16Fe''0.18Mg0.66)3Al2Si3O12
orthopyroxene: 20.990514 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.14Si1.92O6
clinopyroxene: 11.608657 g, composition cpx Na0.08Ca0.74Fe''0.09Mg0.92Fe'''0.02Ti0.02Al0.24Si1.89O6
hornblende: 3.724794 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.192739 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.433585 g, composition Fe''0.21Mg0.80Fe'''0.09Al1.31Cr0.58Ti0.01O4
apatite: 0.032 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1164.070312
liquid: 0.269 g 44.34 0.69 15.88 0.49 0.68 4.83 0.14 8.33 0.07 0.01 5.08 7.76 0.87 0.94 9.88 
Activity of H2O = 0.618219  Melt fraction = 0.00269047
olivine: 53.845070 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.554240 g, composition (Ca0.16Fe''0.18Mg0.66)3Al2Si3O12
orthopyroxene: 21.717746 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.13Si1.92O6
clinopyroxene: 10.592297 g, composition cpx Na0.08Ca0.75Fe''0.09Mg0.92Fe'''0.02Ti0.02Al0.24Si1.89O6
hornblende: 5.275140 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.234543 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.389359 g, composition Fe''0.22Mg0.80Fe'''0.09Al1.28Cr0.60Ti0.02O4
apatite: 0.043 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1154.070312
liquid: 0.107 g 44.14 0.68 15.79 0.48 0.74 4.84 0.15 8.26 0.08 0.01 5.08 7.41 0.75 0.86 10.74 
Activity of H2O = 0.667748  Melt fraction = 0.00107007
olivine: 53.698077 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.479341 g, composition (Ca0.16Fe''0.18Mg0.66)3Al2Si3O12
orthopyroxene: 21.952849 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.13Si1.92O6
clinopyroxene: 10.087276 g, composition cpx Na0.07Ca0.76Fe''0.09Mg0.92Fe'''0.02Ti0.03Al0.23Si1.89O6
hornblende: 5.952685 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.248069 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.347398 g, composition Fe''0.22Mg0.79Fe'''0.09Al1.25Cr0.63Ti0.02O4
apatite: 0.047 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1144.070312
liquid: 0.027 g 43.93 0.66 15.73 0.47 0.81 4.81 0.17 8.16 0.09 0.01 5.03 7.15 0.66 0.79 11.54 
Activity of H2O = 0.712375  Melt fraction = 0.000273933
olivine: 53.628726 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.530873 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12
orthopyroxene: 22.018594 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.13Si1.93O6
clinopyroxene: 9.772084 g, composition cpx Na0.07Ca0.76Fe''0.09Mg0.92Fe'''0.01Ti0.03Al0.23Si1.89O6
hornblende: 6.330544 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.253575 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.308516 g, composition Fe''0.23Mg0.79Fe'''0.09Al1.22Cr0.65Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
...Dropping a liquid from the assemblage.
alphaMELTS at: P 20000.000000, T 1134.070312
olivine: 53.623437 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.672821 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12
orthopyroxene: 21.948711 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6
clinopyroxene: 9.620627 g, composition cpx Na0.07Ca0.77Fe''0.09Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6
hornblende: 6.475084 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255172 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.273360 g, composition Fe''0.23Mg0.79Fe'''0.10Al1.20Cr0.67Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1124.070312
olivine: 53.658848 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 6.867438 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12
orthopyroxene: 21.794215 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6
clinopyroxene: 9.576887 g, composition cpx Na0.07Ca0.77Fe''0.09Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6
hornblende: 6.474978 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255171 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.241208 g, composition Fe''0.23Mg0.79Fe'''0.10Al1.18Cr0.69Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1114.070312
olivine: 53.692537 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.056088 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12
orthopyroxene: 21.643001 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6
clinopyroxene: 9.536212 g, composition cpx Na0.07Ca0.78Fe''0.08Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6
hornblende: 6.474870 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255170 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.210412 g, composition Fe''0.23Mg0.78Fe'''0.10Al1.16Cr0.71Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1104.070312
olivine: 53.724590 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.239090 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12
orthopyroxene: 21.494946 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6
clinopyroxene: 9.498400 g, composition cpx Na0.07Ca0.78Fe''0.08Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6
hornblende: 6.474759 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255169 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.180899 g, composition Fe''0.24Mg0.78Fe'''0.10Al1.14Cr0.73Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1094.070312
olivine: 53.755083 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.416727 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12
orthopyroxene: 21.349934 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.11Si1.93O6
clinopyroxene: 9.463270 g, composition cpx Na0.07Ca0.78Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.89O6
hornblende: 6.474645 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255168 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.152600 g, composition Fe''0.24Mg0.78Fe'''0.10Al1.12Cr0.75Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1084.070312
olivine: 53.784087 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.589261 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12
orthopyroxene: 21.207863 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.93O6
clinopyroxene: 9.430657 g, composition cpx Na0.07Ca0.79Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.89O6
hornblende: 6.474529 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255167 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.125454 g, composition Fe''0.24Mg0.77Fe'''0.10Al1.09Cr0.77Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1074.070312
olivine: 53.811662 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.756924 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12
orthopyroxene: 21.068640 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.93O6
clinopyroxene: 9.400412 g, composition cpx Na0.07Ca0.79Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.90O6
hornblende: 6.474411 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255167 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.099406 g, composition Fe''0.25Mg0.77Fe'''0.10Al1.07Cr0.79Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1064.070312
olivine: 53.837863 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 7.919928 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12
orthopyroxene: 20.932182 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.94O6
clinopyroxene: 9.372399 g, composition cpx Na0.07Ca0.79Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.90O6
hornblende: 6.474290 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255166 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.074409 g, composition Fe''0.25Mg0.77Fe'''0.10Al1.05Cr0.81Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1054.070312
olivine: 53.862741 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.078465 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12
orthopyroxene: 20.798414 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.94O6
clinopyroxene: 9.346495 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.90O6
hornblende: 6.474167 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255165 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.050419 g, composition Fe''0.25Mg0.77Fe'''0.10Al1.03Cr0.84Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1044.070312
olivine: 53.886339 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.232706 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12
orthopyroxene: 20.667270 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.10Si1.94O6
clinopyroxene: 9.322588 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.20Si1.90O6
hornblende: 6.474041 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255164 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.027398 g, composition Fe''0.26Mg0.76Fe'''0.11Al1.00Cr0.86Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
alphaMELTS at: P 20000.000000, T 1034.070312
olivine: 53.908697 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.382809 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12
orthopyroxene: 20.538689 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.73Fe'''0.01Ti0.01Al0.10Si1.94O6
clinopyroxene: 9.300575 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6
hornblende: 6.473912 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255163 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 1.005312 g, composition Fe''0.26Mg0.76Fe'''0.11Al0.98Cr0.88Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
...Adding the solid phase water to the assemblage.
alphaMELTS at: P 20000.000000, T 1024.070312
olivine: 53.930719 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.530180 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12
orthopyroxene: 20.410844 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.10Si1.94O6
clinopyroxene: 9.282727 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6
hornblende: 6.470963 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255163 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.984166 g, composition Fe''0.26Mg0.75Fe'''0.11Al0.95Cr0.90Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.000 g, composition H2O
alphaMELTS at: P 20000.000000, T 1014.070313
olivine: 53.992584 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.733804 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12
orthopyroxene: 20.201784 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.10Si1.94O6
clinopyroxene: 9.377289 g, composition cpx Na0.08Ca0.81Fe''0.07Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6
hornblende: 6.335583 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255162 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.965428 g, composition Fe''0.27Mg0.75Fe'''0.11Al0.93Cr0.92Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.003 g, composition H2O
alphaMELTS at: P 20000.000000, T 1004.070313
olivine: 54.056671 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 8.939458 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12
orthopyroxene: 19.988425 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.10Si1.94O6
clinopyroxene: 9.481035 g, composition cpx Na0.08Ca0.81Fe''0.07Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6
hornblende: 6.190129 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255161 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.947424 g, composition Fe''0.27Mg0.75Fe'''0.11Al0.91Cr0.94Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.006 g, composition H2O
alphaMELTS at: P 20000.000000, T 994.070313
olivine: 54.123351 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 9.147760 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12
orthopyroxene: 19.770058 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.09Si1.94O6
clinopyroxene: 9.594762 g, composition cpx Na0.08Ca0.81Fe''0.07Mg0.89Fe'''0.01Ti0.03Al0.20Si1.90O6
hornblende: 6.033525 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255160 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.930133 g, composition Fe''0.27Mg0.74Fe'''0.11Al0.89Cr0.97Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.009 g, composition H2O
alphaMELTS at: P 20000.000000, T 984.070313
olivine: 54.193047 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 9.359398 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12
orthopyroxene: 19.545873 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6
clinopyroxene: 9.719410 g, composition cpx Na0.09Ca0.81Fe''0.07Mg0.89Fe'''0.01Ti0.03Al0.20Si1.91O6
hornblende: 5.864521 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255159 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.913535 g, composition Fe''0.28Mg0.74Fe'''0.11Al0.87Cr0.99Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.013 g, composition H2O
alphaMELTS at: P 20000.000000, T 974.070313
olivine: 54.266249 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 9.575158 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12
orthopyroxene: 19.314933 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6
clinopyroxene: 9.856092 g, composition cpx Na0.09Ca0.80Fe''0.07Mg0.89Fe'''0.01Ti0.03Al0.20Si1.91O6
hornblende: 5.681651 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255159 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.897614 g, composition Fe''0.28Mg0.74Fe'''0.11Al0.85Cr1.01Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.017 g, composition H2O
alphaMELTS at: P 20000.000000, T 964.070313
olivine: 54.343534 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 9.795942 g, composition (Ca0.14Fe''0.23Mg0.64)3Al2Si3O12
orthopyroxene: 19.076147 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6
clinopyroxene: 10.006133 g, composition cpx Na0.09Ca0.80Fe''0.07Mg0.88Fe'''0.01Ti0.03Al0.20Si1.91O6
hornblende: 5.483185 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255158 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.882352 g, composition Fe''0.28Mg0.73Fe'''0.11Al0.82Cr1.03Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.021 g, composition H2O
alphaMELTS at: P 20000.000000, T 954.070313
olivine: 54.425581 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 10.022804 g, composition (Ca0.14Fe''0.23Mg0.64)3Al2Si3O12
orthopyroxene: 18.828229 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6
clinopyroxene: 10.171117 g, composition cpx Na0.10Ca0.80Fe''0.07Mg0.88Fe'''0.01Ti0.03Al0.20Si1.91O6
hornblende: 5.267059 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255157 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.867734 g, composition Fe''0.29Mg0.73Fe'''0.12Al0.80Cr1.05Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.026 g, composition H2O
alphaMELTS at: P 20000.000000, T 944.070313
olivine: 54.513207 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 10.256990 g, composition (Ca0.13Fe''0.23Mg0.64)3Al2Si3O12
orthopyroxene: 18.569643 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.08Si1.95O6
clinopyroxene: 10.352961 g, composition cpx Na0.10Ca0.80Fe''0.07Mg0.87Fe'''0.01Ti0.03Al0.20Si1.91O6
hornblende: 5.030797 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255156 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.853744 g, composition Fe''0.29Mg0.72Fe'''0.12Al0.78Cr1.07Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.031 g, composition H2O
alphaMELTS at: P 20000.000000, T 934.070313
olivine: 54.607395 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 10.499991 g, composition (Ca0.13Fe''0.23Mg0.64)3Al2Si3O12
orthopyroxene: 18.298542 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.08Si1.95O6
clinopyroxene: 10.553993 g, composition cpx Na0.10Ca0.80Fe''0.07Mg0.87Fe'''0.01Ti0.03Al0.20Si1.91O6
hornblende: 4.771393 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255156 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.840369 g, composition Fe''0.30Mg0.72Fe'''0.12Al0.76Cr1.09Ti0.02O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.036 g, composition H2O
alphaMELTS at: P 20000.000000, T 924.070313
olivine: 54.709348 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 10.753619 g, composition (Ca0.13Fe''0.24Mg0.64)3Al2Si3O12
orthopyroxene: 18.012671 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.95O6
clinopyroxene: 10.777075 g, composition cpx Na0.11Ca0.80Fe''0.07Mg0.87Fe'''0.01Ti0.03Al0.20Si1.91O6
hornblende: 4.485163 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255155 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.827593 g, composition Fe''0.30Mg0.72Fe'''0.12Al0.74Cr1.11Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.042 g, composition H2O
alphaMELTS at: P 20000.000000, T 914.070313
olivine: 54.820555 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 11.020101 g, composition (Ca0.13Fe''0.24Mg0.64)3Al2Si3O12
orthopyroxene: 17.709249 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.96O6
clinopyroxene: 11.025755 g, composition cpx Na0.11Ca0.80Fe''0.07Mg0.86Fe'''0.01Ti0.03Al0.21Si1.92O6
hornblende: 4.167540 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255154 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.815404 g, composition Fe''0.30Mg0.71Fe'''0.12Al0.72Cr1.13Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.049 g, composition H2O
alphaMELTS at: P 20000.000000, T 904.070313
olivine: 54.942876 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4
garnet: 11.302216 g, composition (Ca0.12Fe''0.24Mg0.64)3Al2Si3O12
orthopyroxene: 17.384800 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.96O6
clinopyroxene: 11.304489 g, composition cpx Na0.12Ca0.79Fe''0.06Mg0.86Fe'''0.01Ti0.03Al0.21Si1.92O6
hornblende: 3.812798 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255153 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.803787 g, composition Fe''0.31Mg0.71Fe'''0.12Al0.70Cr1.15Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.057 g, composition H2O
alphaMELTS at: P 20000.000000, T 894.070313
olivine: 55.078676 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 11.603480 g, composition (Ca0.12Fe''0.24Mg0.64)3Al2Si3O12
orthopyroxene: 17.034920 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.96O6
clinopyroxene: 11.618940 g, composition cpx Na0.12Ca0.79Fe''0.06Mg0.85Fe'''0.01Ti0.03Al0.21Si1.92O6
hornblende: 3.413659 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255153 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.792730 g, composition Fe''0.31Mg0.70Fe'''0.12Al0.68Cr1.17Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.065 g, composition H2O
alphaMELTS at: P 20000.000000, T 884.070313
olivine: 55.230995 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 11.928412 g, composition (Ca0.12Fe''0.24Mg0.64)3Al2Si3O12
orthopyroxene: 16.653946 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.07Si1.96O6
clinopyroxene: 11.976403 g, composition cpx Na0.13Ca0.79Fe''0.06Mg0.84Fe'''0.01Ti0.03Al0.21Si1.92O6
hornblende: 2.960746 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255152 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.782220 g, composition Fe''0.31Mg0.70Fe'''0.12Al0.66Cr1.19Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.075 g, composition H2O
alphaMELTS at: P 20000.000000, T 874.070313
olivine: 55.403811 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 12.282912 g, composition (Ca0.12Fe''0.25Mg0.64)3Al2Si3O12
orthopyroxene: 16.234481 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.07Si1.96O6
clinopyroxene: 12.386425 g, composition cpx Na0.14Ca0.78Fe''0.06Mg0.84Fe'''0.01Ti0.03Al0.22Si1.92O6
hornblende: 2.441783 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255151 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.772247 g, composition Fe''0.32Mg0.70Fe'''0.12Al0.64Cr1.21Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.086 g, composition H2O
alphaMELTS at: P 20000.000000, T 864.070313
olivine: 55.602415 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 12.674821 g, composition (Ca0.11Fe''0.25Mg0.64)3Al2Si3O12
orthopyroxene: 15.766697 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6
clinopyroxene: 12.861706 g, composition cpx Na0.14Ca0.78Fe''0.06Mg0.83Fe'''0.01Ti0.03Al0.22Si1.93O6
hornblende: 1.840431 g, composition NaCa2Mg3.73Fe2+0.27Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255150 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.762800 g, composition Fe''0.32Mg0.69Fe'''0.12Al0.62Cr1.23Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.099 g, composition H2O
alphaMELTS at: P 20000.000000, T 854.070313
olivine: 55.833981 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 13.114759 g, composition (Ca0.11Fe''0.25Mg0.64)3Al2Si3O12
orthopyroxene: 15.237290 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6
clinopyroxene: 13.419453 g, composition cpx Na0.15Ca0.77Fe''0.06Mg0.83Fe'''0.01Ti0.03Al0.23Si1.93O6
hornblende: 1.134537 g, composition NaCa2Mg3.73Fe2+0.27Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255150 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.753873 g, composition Fe''0.32Mg0.69Fe'''0.12Al0.61Cr1.25Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.114 g, composition H2O
alphaMELTS at: P 20000.000000, T 844.070313
olivine: 56.108425 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 13.617393 g, composition (Ca0.11Fe''0.25Mg0.64)3Al2Si3O12
orthopyroxene: 14.627896 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6
clinopyroxene: 14.083439 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.82Fe'''0.01Ti0.03Al0.23Si1.93O6
hornblende: 0.293464 g, composition NaCa2Mg3.73Fe2+0.27Al1.00Fe3+0.00Al2Si6O22(OH)2
biotite: 0.255149 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.745463 g, composition Fe''0.32Mg0.69Fe'''0.12Al0.59Cr1.27Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.132 g, composition H2O
...Dropping phase hornblende from the assemblage.
alphaMELTS at: P 20000.000000, T 834.070313
olivine: 56.206913 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 13.859489 g, composition (Ca0.11Fe''0.26Mg0.64)3Al2Si3O12
orthopyroxene: 14.359445 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6
clinopyroxene: 14.307738 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6
biotite: 0.255148 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.736116 g, composition Fe''0.33Mg0.68Fe'''0.12Al0.57Cr1.29Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.138 g, composition H2O
alphaMELTS at: P 20000.000000, T 824.070313
olivine: 56.210501 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 13.960153 g, composition (Ca0.11Fe''0.26Mg0.63)3Al2Si3O12
orthopyroxene: 14.273795 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.06Si1.97O6
clinopyroxene: 14.298327 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6
biotite: 0.255148 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.726683 g, composition Fe''0.34Mg0.67Fe'''0.12Al0.55Cr1.31Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.138 g, composition H2O
alphaMELTS at: P 20000.000000, T 814.070313
olivine: 56.213232 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 14.058245 g, composition (Ca0.11Fe''0.26Mg0.63)3Al2Si3O12
orthopyroxene: 14.189909 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.77Fe'''0.00Ti0.01Al0.06Si1.97O6
clinopyroxene: 14.290043 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6
biotite: 0.255147 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.717796 g, composition Fe''0.34Mg0.67Fe'''0.12Al0.53Cr1.33Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.138 g, composition H2O
alphaMELTS at: P 20000.000000, T 804.070313
olivine: 56.215135 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4
garnet: 14.153870 g, composition (Ca0.10Fe''0.27Mg0.63)3Al2Si3O12
orthopyroxene: 14.107733 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.77Fe'''0.00Ti0.01Al0.06Si1.97O6
clinopyroxene: 14.282833 g, composition cpx Na0.16Ca0.77Fe''0.05Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6
biotite: 0.255147 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2
spinel: 0.709425 g, composition Fe''0.35Mg0.66Fe'''0.12Al0.50Cr1.35Ti0.01O4
apatite: 0.049 g, composition Ca5(PO4)3OH
water: 0.138 g, composition H2O
Minimum Temperature reached
Successful return from alphamelts...

Choose:
 1. Read MELTS file to set composition of system
 2. Twiddle starting or continuation parameters
 3. Single (batch) calculation
 4. Execute (follow path, mineral isograd or melt contour)
 5. Set fO2 buffer
 6. Set H2O (ppm) or aH2O
 7. Impose initial entropy, enthalpy or volume
 8. Adjust solid phase setting(s)
 9. Turn liquid on / off
10. Turn phase diagram mode on / off
11. Update state using restart file
12. Update composition using MELTS file
13. Write out restart file
14. Write out MELTS file
15. Write thermodynamic output for all phases
16. Calculate integrated melt and output file(s)
17. Fit parental melt composition (amoeba)
18. Cumulate Invertor (not yet implemented)
-1. Turn off menu display for options 1-18
 0. QUIT
Your choice: 

Out[3]:
Dict{Any,Any} with 15 entries:
  "Temperature"     => [1624.07, 1614.07, 1604.07, 1594.07, 1584.07, 1574.07, 1…
  "water"           => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  …  0.…
  "orthopyroxene_0" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 7.96064, 11.4105, 14.2824…
  "elements"        => Any["Pressure", "Temperature", "mass", "liquid_0", "oliv…
  "clinopyroxene_0" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  …  11…
  "mass"            => [99.9575, 99.9544, 99.9515, 99.949, 99.9467, 99.9445, 99…
  "olivine_0"       => [0.0, 8.88223, 16.6072, 23.3962, 29.1114, 33.9796, 36.01…
  "hornblende_0"    => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  …  3.…
  "garnet_0"        => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  …  11…
  "spinel_0"        => [0.0, 0.0, 0.0, 0.134547, 0.235384, 0.321473, 0.464989, …
  "liquid_0"        => [99.9575, 91.0722, 83.3443, 76.4183, 70.5999, 65.6434, 5…
  "apatite"         => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  …  0.…
  "Pressure"        => [20000.0, 20000.0, 20000.0, 20000.0, 20000.0, 20000.0, 2…
  ""                => [NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN  …  Na…
  "biotite_0"       => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  …  0.…

Plot results

In [4]:
## --- Plot melt composition
melt_comp = melts_query_liquid(scratchdir, index=1)

h = plot(xlabel="Percent melt",ylabel="Abundance (wt. %) in melt")
for e in ["SiO2","Al2O3","CaO","MgO","FeO","Na2O","K2O"]
    plot!(h,melt_comp["mass"],melt_comp[e],label=e)
end
plot!(h,xlims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:right)
Out[4]:
0 20 40 60 80 100 0 10 20 30 40 Percent melt Abudance (wt. %) in melt SiO2 Al2O3 CaO MgO FeO Na2O K2O
In [5]:
## --- Plot solid composition
solid_comp = melts_query_solid(scratchdir, index=1)

h = plot(xlabel="Percent melt",ylabel="Abundance (wt. %) in solid")
for e in ["SiO2","Al2O3","CaO","MgO","FeO","Na2O","K2O"]
    plot!(h,100 .- solid_comp["mass"],solid_comp[e],label=e)
end
plot!(h,xlims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:right)
Out[5]:
0 20 40 60 80 100 0 10 20 30 40 50 Percent melt Abudance (wt. %) in solid SiO2 Al2O3 CaO MgO FeO Na2O K2O
In [6]:
## --- Plot phase modes
modes = melts_query_modes(scratchdir, index=1)

h = plot(xlabel="Temperature (C)",ylabel="Abundance (wt. %)")
for m in modes["elements"][4:end]
    plot!(h,modes["Temperature"],modes[m],label=m)
end
plot!(h,ylims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:topleft)
Out[6]:
800 1000 1200 1400 1600 0 20 40 60 80 100 Temperature (C) Abudance (wt. %) liquid_0 olivine_0 garnet_0 orthopyroxene_0 clinopyroxene_0 hornblende_0 biotite_0 spinel_0 apatite water
In [ ]: