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.
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 the StatGeochem package
using StatGeochem
using Plots
# 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 = abspath("$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/";
## --- # # # # # # # # # # # 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_CELSIUS_OUTPUT true ALPHAMELTS_MINT 800.0 ALPHAMELTS_MINP 2.0 ALPHAMELTS_DELTAP 0.0 ALPHAMELTS_DELTAT -10.0 ALPHAMELTS_MAXP 90000.0 ALPHAMELTS_MODE isobaric 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:
Dict{String, Union{Vector{Float64}, Vector{String}}} with 14 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" => ["Pressure", "Temperature", "mass", "liquid_0", "olivine… "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… "biotite_0" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.…
## --- 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)
## --- 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)
## --- 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)