{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Julia-MELTS interface demo\n", "\n", "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.\n", "\n", "StatGeochem.jl also includes some of the codes and utilities used in Keller & Schoene 2012, Keller et al. 2015 and Keller & Schoene 2018.\n", "\n", "\"Launch \n", "

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!

\n", "\n", "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.\n", "***\n", "### Load required Julia packages" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "## --- Load (and install if neccesary) the StatGeochem package\n", "try\n", " using StatGeochem\n", "catch\n", " using Pkg\n", " Pkg.add(\"StatGeochem\")\n", " using StatGeochem\n", "end\n", "\n", "using Statistics, StatsBase, DelimitedFiles\n", "using Plots; gr();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### AlphaMELTS calculations\n", "#### Try to download and install alphaMELTS" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Download precompiled executable from caltech\n", "alphameltsversion = \"linux_alphamelts_1-9\"\n", "download(\"https://magmasource.caltech.edu/alphamelts/zipfiles/$alphameltsversion.zip\",\"./$alphameltsversion.zip\")\n", "run(`unzip -o $alphameltsversion.zip`);\n", "run(`mv $alphameltsversion/alphamelts_linux64 $alphameltsversion/alphamelts`);\n", "\n", "# The absolute path to the alphamelts perl script must be input here:\n", " # meltspath = \"/usr/local/bin/run_alphamelts.command\"\n", "meltspath = pwd()*\"/$alphameltsversion/run_alphamelts.command\";\n", "\n", "# The absolute or relative path to a direcory where you want to run MELTS -\n", "# A number of configuration and output files will be written here.\n", " # scratchdir = \"scratch/\"\n", "scratchdir = pwd()*\"/$alphameltsversion/\";" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Run alphamelts" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ALPHAMELTS_VERSION pMELTS\n", "ALPHAMELTS_DELTAT -10.0\n", "ALPHAMELTS_MINP 2.0\n", "ALPHAMELTS_MODE isobaric\n", "ALPHAMELTS_MINT 800.0\n", "ALPHAMELTS_CELSIUS_OUTPUT true\n", "ALPHAMELTS_DELTAP 0.0\n", "ALPHAMELTS_MAXP 90000.0\n", "ALPHAMELTS_MAXT 3000.0\n", "Checking for updates...\n", "\n", "\n", "*** alphaMELTS 1.2 -- pMELTS isobar w/ or w/o liquid ***\n", "\n", "This front end is the work of Paul Asimow and Paula Antoshechkina\n", "(nee Smith) and it uses the MELTS and pMELTS algorithms developed by\n", "Mark Ghiorso & co-workers. You are welcome to use and distribute this\n", "program, under the condition that you acknowledge all the contributors\n", "by citing the appropriate references with any results:\n", "\n", "Smith & Asimow (2005) -- front end\n", "Ghiorso et al (2002) -- pMELTS thermodynamic model\n", "Asimow & Ghiorso (1998) -- subsolidus calculations (if used)\n", "\n", "There may be others, depending on the exact calculation performed\n", "and the source and / or partition coefficients (e.g. D(P,T,X)) used.\n", "See Smith & Asimow (2005), documentation and the forum for details.\n", "\n", "Unable to automatically check for updates; check for updates at\n", "http://magmasource.caltech.edu/alphamelts/ and\n", "http://magmasource.caltech.edu/forum/\n", "\n", "Choose:\n", " 1. Read MELTS file to set composition of system\n", " 2. Twiddle starting or continuation parameters\n", " 3. Single (batch) calculation\n", " 4. Execute (follow path, mineral isograd or melt contour)\n", " 5. Set fO2 buffer\n", " 6. Set H2O (ppm) or aH2O\n", " 7. Impose initial entropy, enthalpy or volume\n", " 8. Adjust solid phase setting(s)\n", " 9. Turn liquid on / off\n", "10. Turn phase diagram mode on / off\n", "11. Update state using restart file\n", "12. Update composition using MELTS file\n", "13. Write out restart file\n", "14. Write out MELTS file\n", "15. Write thermodynamic output for all phases\n", "16. Calculate integrated melt and output file(s)\n", "17. Fit parental melt composition (amoeba)\n", "18. Cumulate Invertor (not yet implemented)\n", "-1. Turn off menu display for options 1-18\n", " 0. QUIT\n", "Your choice: MELTS filename: input file open\n", "Successfully read 'initial composition: sio2 44.7801'\n", "Successfully read 'initial composition: tio2 0.1989'\n", "Successfully read 'initial composition: al2o3 4.4282'\n", "Successfully read 'initial composition: fe2o3 0.9773'\n", "Successfully read 'initial composition: cr2o3 0.3821'\n", "Successfully read 'initial composition: feo 7.1313'\n", "Successfully read 'initial composition: mno 0.1343'\n", "Successfully read 'initial composition: mgo 37.6152'\n", "Successfully read 'initial composition: nio 0.2487'\n", "Successfully read 'initial composition: coo 0.0128'\n", "Successfully read 'initial composition: cao 3.5326'\n", "Successfully read 'initial composition: na2o 0.3582'\n", "Successfully read 'initial composition: k2o 0.0288'\n", "Successfully read 'initial composition: p2o5 0.0208'\n", "Successfully read 'initial composition: h2o 0.1499'\n", "Successfully read 'initial temperature: 1700.0'\n", "Successfully read 'initial pressure: 20000.0'\n", "Successfully read 'log fo2 path: fmq'\n", "\n", "Choose:\n", " 1. Read MELTS file to set composition of system\n", " 2. Twiddle starting or continuation parameters\n", " 3. Single (batch) calculation\n", " 4. Execute (follow path, mineral isograd or melt contour)\n", " 5. Set fO2 buffer\n", " 6. Set H2O (ppm) or aH2O\n", " 7. Impose initial entropy, enthalpy or volume\n", " 8. Adjust solid phase setting(s)\n", " 9. Turn liquid on / off\n", "10. Turn phase diagram mode on / off\n", "11. Update state using restart file\n", "12. Update composition using MELTS file\n", "13. Write out restart file\n", "14. Write out MELTS file\n", "15. Write thermodynamic output for all phases\n", "16. Calculate integrated melt and output file(s)\n", "17. Fit parental melt composition (amoeba)\n", "18. Cumulate Invertor (not yet implemented)\n", "-1. Turn off menu display for options 1-18\n", " 0. QUIT\n", "Your choice: Turn phase diagram mode on (1) or off (0)? \n", "Choose:\n", " 1. Read MELTS file to set composition of system\n", " 2. Twiddle starting or continuation parameters\n", " 3. Single (batch) calculation\n", " 4. Execute (follow path, mineral isograd or melt contour)\n", " 5. Set fO2 buffer\n", " 6. Set H2O (ppm) or aH2O\n", " 7. Impose initial entropy, enthalpy or volume\n", " 8. Adjust solid phase setting(s)\n", " 9. Turn liquid on / off\n", "10. Turn phase diagram mode on / off\n", "11. Update state using restart file\n", "12. Update composition using MELTS file\n", "13. Write out restart file\n", "14. Write out MELTS file\n", "15. Write thermodynamic output for all phases\n", "16. Calculate integrated melt and output file(s)\n", "17. Fit parental melt composition (amoeba)\n", "18. Cumulate Invertor (not yet implemented)\n", "-1. Turn off menu display for options 1-18\n", " 0. QUIT\n", "Your choice: Superliquidus (1) or subsolidus (0) initial guess ? Phase to track boundary of (by name, lower case): Type of melt contour to track:\n", "0. Phi (melt fraction by volume)\n", "1. F (melt fraction by mass)\n", "2. aH2O (activity of water in the melt)\n", "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\n", "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 \n", "Activity of H2O = 0.00111967 Melt fraction = 1\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "Isograd solution at: P 20000.000000, T 1624.070312\n", "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 \n", "Activity of H2O = 0.00119344 Melt fraction = 1\n", "\n", "Choose:\n", " 1. Read MELTS file to set composition of system\n", " 2. Twiddle starting or continuation parameters\n", " 3. Single (batch) calculation\n", " 4. Execute (follow path, mineral isograd or melt contour)\n", " 5. Set fO2 buffer\n", " 6. Set H2O (ppm) or aH2O\n", " 7. Impose initial entropy, enthalpy or volume\n", " 8. Adjust solid phase setting(s)\n", " 9. Turn liquid on / off\n", "10. Turn phase diagram mode on / off\n", "11. Update state using restart file\n", "12. Update composition using MELTS file\n", "13. Write out restart file\n", "14. Write out MELTS file\n", "15. Write thermodynamic output for all phases\n", "16. Calculate integrated melt and output file(s)\n", "17. Fit parental melt composition (amoeba)\n", "18. Cumulate Invertor (not yet implemented)\n", "-1. Turn off menu display for options 1-18\n", " 0. QUIT\n", "Your choice: Turn phase diagram mode on (1) or off (0)? \n", "Choose:\n", " 1. Read MELTS file to set composition of system\n", " 2. Twiddle starting or continuation parameters\n", " 3. Single (batch) calculation\n", " 4. Execute (follow path, mineral isograd or melt contour)\n", " 5. Set fO2 buffer\n", " 6. Set H2O (ppm) or aH2O\n", " 7. Impose initial entropy, enthalpy or volume\n", " 8. Adjust solid phase setting(s)\n", " 9. Turn liquid on / off\n", "10. Turn phase diagram mode on / off\n", "11. Update state using restart file\n", "12. Update composition using MELTS file\n", "13. Write out restart file\n", "14. Write out MELTS file\n", "15. Write thermodynamic output for all phases\n", "16. Calculate integrated melt and output file(s)\n", "17. Fit parental melt composition (amoeba)\n", "18. Cumulate Invertor (not yet implemented)\n", "-1. Turn off menu display for options 1-18\n", " 0. QUIT\n", "Your choice: ...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "...Adding the solid phase olivine to the assemblage.\n", "...Dropping phase olivine from the assemblage.\n", "Initial alphaMELTS calculation at: P 20000.000000, T 1624.070312\n", "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 \n", "Activity of H2O = 0.00119344 Melt fraction = 1\n", "...Adding the solid phase olivine to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1614.070312\n", "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 \n", "Activity of H2O = 0.00143887 Melt fraction = 0.911137\n", "olivine: 8.882234 g, composition (Ca0.00Mg0.94Fe''0.05Mn0.00Co0.00Ni0.01)2SiO4\n", "alphaMELTS at: P 20000.000000, T 1604.070312\n", "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 \n", "Activity of H2O = 0.00171279 Melt fraction = 0.833847\n", "olivine: 16.607219 g, composition (Ca0.00Mg0.94Fe''0.05Mn0.00Co0.00Ni0.01)2SiO4\n", "...Adding the solid phase spinel to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1594.070312\n", "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 \n", "Activity of H2O = 0.00202505 Melt fraction = 0.764573\n", "olivine: 23.396165 g, composition (Ca0.00Mg0.94Fe''0.05Mn0.00Co0.00Ni0.01)2SiO4\n", "spinel: 0.134547 g, composition Fe''0.15Mg0.85Fe'''0.07Al0.75Cr1.17Ti0.00O4\n", "alphaMELTS at: P 20000.000000, T 1584.070312\n", "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 \n", "Activity of H2O = 0.00235381 Melt fraction = 0.706376\n", "olivine: 29.111411 g, composition (Ca0.00Mg0.93Fe''0.06Mn0.00Co0.00Ni0.01)2SiO4\n", "spinel: 0.235384 g, composition Fe''0.16Mg0.84Fe'''0.07Al0.77Cr1.15Ti0.00O4\n", "...Adding the solid phase orthopyroxene to the assemblage.\n", "...Dropping phase orthopyroxene from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1574.070312\n", "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 \n", "Activity of H2O = 0.00269699 Melt fraction = 0.656799\n", "olivine: 33.979596 g, composition (Ca0.00Mg0.93Fe''0.06Mn0.00Co0.00Ni0.01)2SiO4\n", "spinel: 0.321473 g, composition Fe''0.17Mg0.84Fe'''0.08Al0.79Cr1.12Ti0.00O4\n", "...Adding the solid phase orthopyroxene to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1564.070312\n", "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 \n", "Activity of H2O = 0.00378705 Melt fraction = 0.555371\n", "olivine: 36.012180 g, composition (Ca0.00Mg0.93Fe''0.06Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 7.960639 g, composition opx Na0.00Ca0.05Fe''0.10Mg1.79Fe'''0.01Ti0.00Al0.09Si1.95O6\n", "spinel: 0.464989 g, composition Fe''0.17Mg0.84Fe'''0.08Al0.89Cr1.02Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1554.070312\n", "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 \n", "Activity of H2O = 0.00472055 Melt fraction = 0.496244\n", "olivine: 38.364714 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 11.410543 g, composition opx Na0.00Ca0.06Fe''0.10Mg1.78Fe'''0.01Ti0.00Al0.10Si1.95O6\n", "spinel: 0.571432 g, composition Fe''0.17Mg0.84Fe'''0.08Al0.95Cr0.96Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1544.070312\n", "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 \n", "Activity of H2O = 0.00577711 Melt fraction = 0.447037\n", "olivine: 40.307067 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 14.282365 g, composition opx Na0.00Ca0.06Fe''0.11Mg1.77Fe'''0.01Ti0.00Al0.11Si1.94O6\n", "spinel: 0.674624 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.01Cr0.89Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1534.070312\n", "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 \n", "Activity of H2O = 0.00696605 Melt fraction = 0.405313\n", "olivine: 41.933870 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 16.722170 g, composition opx Na0.00Ca0.07Fe''0.11Mg1.76Fe'''0.01Ti0.00Al0.11Si1.94O6\n", "spinel: 0.777506 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.08Cr0.83Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1524.070312\n", "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 \n", "Activity of H2O = 0.00829785 Melt fraction = 0.369393\n", "olivine: 43.311147 g, composition (Ca0.00Mg0.92Fe''0.07Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 18.829911 g, composition opx Na0.00Ca0.07Fe''0.12Mg1.75Fe'''0.01Ti0.00Al0.12Si1.93O6\n", "spinel: 0.882041 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.14Cr0.77Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1514.070312\n", "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 \n", "Activity of H2O = 0.00978407 Melt fraction = 0.338081\n", "olivine: 44.486449 g, composition (Ca0.00Mg0.92Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 20.676180 g, composition opx Na0.00Ca0.07Fe''0.12Mg1.73Fe'''0.01Ti0.00Al0.13Si1.93O6\n", "spinel: 0.989456 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.19Cr0.71Ti0.01O4\n", "...Adding the solid phase clinopyroxene to the assemblage.\n", "...Dropping phase clinopyroxene from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1504.070312\n", "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 \n", "Activity of H2O = 0.0114373 Melt fraction = 0.310505\n", "olivine: 45.495020 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 22.312332 g, composition opx Na0.00Ca0.08Fe''0.12Mg1.72Fe'''0.01Ti0.00Al0.13Si1.93O6\n", "spinel: 1.100424 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.25Cr0.66Ti0.01O4\n", "...Adding the solid phase clinopyroxene to the assemblage.\n", "...Dropping phase clinopyroxene from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1494.070312\n", "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 \n", "Activity of H2O = 0.0132713 Melt fraction = 0.286008\n", "olivine: 46.363698 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 23.776873 g, composition opx Na0.00Ca0.08Fe''0.13Mg1.71Fe'''0.01Ti0.00Al0.14Si1.92O6\n", "spinel: 1.215221 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.30Cr0.61Ti0.01O4\n", "...Adding the solid phase clinopyroxene to the assemblage.\n", "...Dropping phase clinopyroxene from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1484.070312\n", "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 \n", "Activity of H2O = 0.0153017 Melt fraction = 0.26408\n", "olivine: 47.113450 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.01)2SiO4\n", "orthopyroxene: 25.099600 g, composition opx Na0.00Ca0.09Fe''0.13Mg1.70Fe'''0.01Ti0.00Al0.14Si1.92O6\n", "spinel: 1.333875 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.34Cr0.56Ti0.01O4\n", "...Adding the solid phase clinopyroxene to the assemblage.\n", "...Dropping phase clinopyroxene from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1474.070312\n", "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 \n", "Activity of H2O = 0.0175466 Melt fraction = 0.244318\n", "olivine: 47.761057 g, composition (Ca0.00Mg0.91Fe''0.08Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 26.304323 g, composition opx Na0.00Ca0.09Fe''0.13Mg1.69Fe'''0.01Ti0.00Al0.15Si1.92O6\n", "spinel: 1.456306 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.39Cr0.52Ti0.01O4\n", "...Adding the solid phase clinopyroxene to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1464.070312\n", "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 \n", "Activity of H2O = 0.022457 Melt fraction = 0.211448\n", "olivine: 48.588784 g, composition (Ca0.00Mg0.91Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 25.514705 g, composition opx Na0.00Ca0.09Fe''0.14Mg1.69Fe'''0.01Ti0.00Al0.15Si1.92O6\n", "clinopyroxene: 3.049418 g, composition cpx Na0.02Ca0.57Fe''0.11Mg1.18Fe'''0.02Ti0.01Al0.21Si1.89O6\n", "spinel: 1.653782 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.45Cr0.46Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1454.070312\n", "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 \n", "Activity of H2O = 0.0285482 Melt fraction = 0.183011\n", "olivine: 49.290361 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.942667 g, composition opx Na0.00Ca0.09Fe''0.14Mg1.68Fe'''0.01Ti0.00Al0.16Si1.91O6\n", "clinopyroxene: 5.569546 g, composition cpx Na0.02Ca0.58Fe''0.11Mg1.16Fe'''0.02Ti0.01Al0.22Si1.89O6\n", "spinel: 1.846146 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.49Cr0.42Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1444.070312\n", "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 \n", "Activity of H2O = 0.0358609 Melt fraction = 0.158908\n", "olivine: 49.880834 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.566195 g, composition opx Na0.00Ca0.08Fe''0.14Mg1.68Fe'''0.01Ti0.00Al0.17Si1.91O6\n", "clinopyroxene: 7.586043 g, composition cpx Na0.02Ca0.59Fe''0.11Mg1.14Fe'''0.02Ti0.01Al0.22Si1.88O6\n", "spinel: 2.024578 g, composition Fe''0.17Mg0.84Fe'''0.08Al1.53Cr0.38Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1434.070312\n", "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 \n", "Activity of H2O = 0.0444289 Melt fraction = 0.138614\n", "olivine: 50.378410 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.332580 g, composition opx Na0.00Ca0.08Fe''0.15Mg1.68Fe'''0.01Ti0.01Al0.17Si1.91O6\n", "clinopyroxene: 9.189705 g, composition cpx Na0.02Ca0.60Fe''0.12Mg1.12Fe'''0.02Ti0.01Al0.23Si1.88O6\n", "spinel: 2.185259 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.55Cr0.35Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1424.070312\n", "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 \n", "Activity of H2O = 0.0542477 Melt fraction = 0.121566\n", "olivine: 50.798768 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.200332 g, composition opx Na0.00Ca0.08Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.17Si1.90O6\n", "clinopyroxene: 10.463339 g, composition cpx Na0.03Ca0.61Fe''0.11Mg1.11Fe'''0.02Ti0.01Al0.24Si1.88O6\n", "spinel: 2.327302 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.57Cr0.33Ti0.01O4\n", "...Adding the solid phase garnet to the assemblage.\n", "...Dropping phase garnet from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1414.070312\n", "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 \n", "Activity of H2O = 0.0652847 Melt fraction = 0.107232\n", "olivine: 51.155082 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.138395 g, composition opx Na0.01Ca0.08Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6\n", "clinopyroxene: 11.477210 g, composition cpx Na0.03Ca0.62Fe''0.11Mg1.09Fe'''0.02Ti0.01Al0.24Si1.87O6\n", "spinel: 2.451633 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.59Cr0.32Ti0.01O4\n", "...Adding the solid phase garnet to the assemblage.\n", "...Dropping phase garnet from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1404.070312\n", "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 \n", "Activity of H2O = 0.0774864 Melt fraction = 0.0951415\n", "olivine: 51.458212 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.124315 g, composition opx Na0.01Ca0.08Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6\n", "clinopyroxene: 12.288111 g, composition cpx Na0.03Ca0.63Fe''0.11Mg1.08Fe'''0.02Ti0.01Al0.25Si1.87O6\n", "spinel: 2.560042 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.60Cr0.30Ti0.01O4\n", "...Adding the solid phase garnet to the assemblage.\n", "...Dropping phase garnet from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1394.070312\n", "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 \n", "Activity of H2O = 0.0907855 Melt fraction = 0.0848974\n", "olivine: 51.717009 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.142219 g, composition opx Na0.01Ca0.07Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6\n", "clinopyroxene: 12.940628 g, composition cpx Na0.04Ca0.63Fe''0.11Mg1.07Fe'''0.02Ti0.01Al0.25Si1.87O6\n", "spinel: 2.654583 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.62Cr0.29Ti0.01O4\n", "...Adding the solid phase garnet to the assemblage.\n", "...Dropping phase garnet from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1384.070312\n", "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 \n", "Activity of H2O = 0.105106 Melt fraction = 0.076174\n", "olivine: 51.938661 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "orthopyroxene: 24.181036 g, composition opx Na0.01Ca0.07Fe''0.15Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6\n", "clinopyroxene: 13.469233 g, composition cpx Na0.04Ca0.64Fe''0.11Mg1.05Fe'''0.02Ti0.01Al0.26Si1.87O6\n", "spinel: 2.737252 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.62Cr0.28Ti0.01O4\n", "...Adding the solid phase garnet to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1374.070312\n", "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 \n", "Activity of H2O = 0.120429 Melt fraction = 0.0686709\n", "olivine: 52.150643 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 0.089132 g, composition (Ca0.18Fe''0.14Mg0.68)3Al2Si3O12\n", "orthopyroxene: 24.175149 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6\n", "clinopyroxene: 13.870519 g, composition cpx Na0.04Ca0.65Fe''0.11Mg1.04Fe'''0.02Ti0.01Al0.26Si1.87O6\n", "spinel: 2.790367 g, composition Fe''0.17Mg0.84Fe'''0.07Al1.63Cr0.28Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1364.070312\n", "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 \n", "Activity of H2O = 0.137078 Melt fraction = 0.0619795\n", "olivine: 52.510356 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 0.903057 g, composition (Ca0.18Fe''0.14Mg0.68)3Al2Si3O12\n", "orthopyroxene: 23.701151 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6\n", "clinopyroxene: 13.953296 g, composition cpx Na0.04Ca0.65Fe''0.11Mg1.03Fe'''0.02Ti0.01Al0.26Si1.87O6\n", "spinel: 2.675436 g, composition Fe''0.17Mg0.83Fe'''0.07Al1.62Cr0.29Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1354.070312\n", "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 \n", "Activity of H2O = 0.154397 Melt fraction = 0.0562602\n", "olivine: 52.826652 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 1.641180 g, composition (Ca0.18Fe''0.15Mg0.68)3Al2Si3O12\n", "orthopyroxene: 23.273605 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.18Si1.90O6\n", "clinopyroxene: 14.003803 g, composition cpx Na0.04Ca0.66Fe''0.11Mg1.02Fe'''0.02Ti0.01Al0.26Si1.87O6\n", "spinel: 2.568454 g, composition Fe''0.18Mg0.83Fe'''0.07Al1.60Cr0.30Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1344.070312\n", "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 \n", "Activity of H2O = 0.172301 Melt fraction = 0.051337\n", "olivine: 53.106072 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 2.314755 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12\n", "orthopyroxene: 22.884334 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.90O6\n", "clinopyroxene: 14.030766 g, composition cpx Na0.05Ca0.67Fe''0.11Mg1.02Fe'''0.02Ti0.01Al0.26Si1.87O6\n", "spinel: 2.468664 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.59Cr0.32Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1334.070312\n", "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 \n", "Activity of H2O = 0.190714 Melt fraction = 0.0470717\n", "olivine: 53.353971 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 2.932822 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12\n", "orthopyroxene: 22.527017 g, composition opx Na0.01Ca0.07Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.90O6\n", "clinopyroxene: 14.040591 g, composition cpx Na0.05Ca0.67Fe''0.11Mg1.01Fe'''0.02Ti0.01Al0.26Si1.87O6\n", "spinel: 2.375397 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.57Cr0.33Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1324.070312\n", "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 \n", "Activity of H2O = 0.209569 Melt fraction = 0.0433543\n", "olivine: 53.574796 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 3.502746 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12\n", "orthopyroxene: 22.196672 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.91O6\n", "clinopyroxene: 14.038032 g, composition cpx Na0.05Ca0.68Fe''0.11Mg1.00Fe'''0.02Ti0.01Al0.26Si1.87O6\n", "spinel: 2.288062 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.56Cr0.34Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1314.070312\n", "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 \n", "Activity of H2O = 0.22881 Melt fraction = 0.0400961\n", "olivine: 53.772272 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 4.030605 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12\n", "orthopyroxene: 21.889309 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.67Fe'''0.01Ti0.01Al0.17Si1.91O6\n", "clinopyroxene: 14.026655 g, composition cpx Na0.05Ca0.68Fe''0.10Mg0.99Fe'''0.02Ti0.02Al0.26Si1.87O6\n", "spinel: 2.206131 g, composition Fe''0.18Mg0.83Fe'''0.08Al1.54Cr0.36Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1304.070312\n", "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 \n", "Activity of H2O = 0.248388 Melt fraction = 0.037225\n", "olivine: 53.949552 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 4.521471 g, composition (Ca0.17Fe''0.15Mg0.68)3Al2Si3O12\n", "orthopyroxene: 21.601681 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.17Si1.91O6\n", "clinopyroxene: 14.009152 g, composition cpx Na0.06Ca0.69Fe''0.10Mg0.99Fe'''0.02Ti0.02Al0.26Si1.87O6\n", "spinel: 2.129135 g, composition Fe''0.19Mg0.83Fe'''0.08Al1.53Cr0.37Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1294.070312\n", "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 \n", "Activity of H2O = 0.268264 Melt fraction = 0.0346821\n", "olivine: 54.109313 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 4.979616 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12\n", "orthopyroxene: 21.331113 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6\n", "clinopyroxene: 13.987567 g, composition cpx Na0.06Ca0.69Fe''0.10Mg0.98Fe'''0.02Ti0.02Al0.26Si1.88O6\n", "spinel: 2.056649 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.51Cr0.38Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1284.070312\n", "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 \n", "Activity of H2O = 0.288408 Melt fraction = 0.0324187\n", "olivine: 54.253844 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 5.408668 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12\n", "orthopyroxene: 21.075380 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6\n", "clinopyroxene: 13.963453 g, composition cpx Na0.06Ca0.70Fe''0.10Mg0.97Fe'''0.02Ti0.02Al0.26Si1.88O6\n", "spinel: 1.988296 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.50Cr0.40Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1274.070312\n", "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 \n", "Activity of H2O = 0.308796 Melt fraction = 0.0303942\n", "olivine: 54.385101 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 5.811732 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12\n", "orthopyroxene: 20.832617 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6\n", "clinopyroxene: 13.937990 g, composition cpx Na0.06Ca0.70Fe''0.10Mg0.97Fe'''0.02Ti0.02Al0.26Si1.88O6\n", "spinel: 1.923733 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.48Cr0.41Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1264.070312\n", "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 \n", "Activity of H2O = 0.32941 Melt fraction = 0.0285752\n", "olivine: 54.504761 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.191479 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12\n", "orthopyroxene: 20.601247 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.16Si1.91O6\n", "clinopyroxene: 13.912073 g, composition cpx Na0.07Ca0.70Fe''0.10Mg0.96Fe'''0.02Ti0.02Al0.26Si1.88O6\n", "spinel: 1.862651 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.47Cr0.43Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1254.070312\n", "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 \n", "Activity of H2O = 0.350238 Melt fraction = 0.0269334\n", "olivine: 54.614264 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.550223 g, composition (Ca0.17Fe''0.16Mg0.67)3Al2Si3O12\n", "orthopyroxene: 20.379933 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.68Fe'''0.01Ti0.01Al0.15Si1.91O6\n", "clinopyroxene: 13.886375 g, composition cpx Na0.07Ca0.71Fe''0.10Mg0.96Fe'''0.02Ti0.02Al0.26Si1.88O6\n", "spinel: 1.804774 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.45Cr0.44Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1244.070312\n", "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 \n", "Activity of H2O = 0.37127 Melt fraction = 0.0254453\n", "olivine: 54.714847 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.889977 g, composition (Ca0.16Fe''0.16Mg0.67)3Al2Si3O12\n", "orthopyroxene: 20.167534 g, composition opx Na0.01Ca0.06Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.15Si1.91O6\n", "clinopyroxene: 13.861396 g, composition cpx Na0.07Ca0.71Fe''0.10Mg0.95Fe'''0.02Ti0.02Al0.25Si1.88O6\n", "spinel: 1.749849 g, composition Fe''0.19Mg0.82Fe'''0.08Al1.44Cr0.46Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1234.070312\n", "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 \n", "Activity of H2O = 0.392501 Melt fraction = 0.0240909\n", "olivine: 54.807573 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.212501 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12\n", "orthopyroxene: 19.963075 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.15Si1.92O6\n", "clinopyroxene: 13.837503 g, composition cpx Na0.07Ca0.71Fe''0.09Mg0.95Fe'''0.02Ti0.02Al0.25Si1.88O6\n", "spinel: 1.697649 g, composition Fe''0.20Mg0.82Fe'''0.08Al1.42Cr0.47Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1224.070312\n", "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 \n", "Activity of H2O = 0.413926 Melt fraction = 0.0228534\n", "olivine: 54.893359 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.519342 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12\n", "orthopyroxene: 19.765716 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.15Si1.92O6\n", "clinopyroxene: 13.814956 g, composition cpx Na0.08Ca0.72Fe''0.09Mg0.94Fe'''0.02Ti0.02Al0.25Si1.88O6\n", "spinel: 1.647966 g, composition Fe''0.20Mg0.82Fe'''0.08Al1.41Cr0.49Ti0.01O4\n", "alphaMELTS at: P 20000.000000, T 1214.070312\n", "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 \n", "Activity of H2O = 0.435546 Melt fraction = 0.0217186\n", "olivine: 54.972993 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.811862 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12\n", "orthopyroxene: 19.574738 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.14Si1.92O6\n", "clinopyroxene: 13.793936 g, composition cpx Na0.08Ca0.72Fe''0.09Mg0.94Fe'''0.02Ti0.02Al0.25Si1.88O6\n", "spinel: 1.600614 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.39Cr0.50Ti0.01O4\n", "...Adding the solid phase hornblende to the assemblage.\n", "...Dropping phase hornblende from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1204.070312\n", "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 \n", "Activity of H2O = 0.457358 Melt fraction = 0.0206743\n", "olivine: 55.047157 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.091268 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12\n", "orthopyroxene: 19.389518 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.69Fe'''0.01Ti0.01Al0.14Si1.92O6\n", "clinopyroxene: 13.774559 g, composition cpx Na0.08Ca0.72Fe''0.09Mg0.93Fe'''0.02Ti0.02Al0.25Si1.89O6\n", "spinel: 1.555422 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.37Cr0.52Ti0.01O4\n", "...Adding the solid phase hornblende to the assemblage.\n", "...Dropping phase hornblende from the assemblage.\n", "...Adding the solid phase biotite to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1194.070312\n", "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 \n", "Activity of H2O = 0.48032 Melt fraction = 0.0195234\n", "olivine: 55.101579 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.353334 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12\n", "orthopyroxene: 19.217252 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.70Fe'''0.01Ti0.01Al0.14Si1.92O6\n", "clinopyroxene: 13.768639 g, composition cpx Na0.08Ca0.73Fe''0.09Mg0.93Fe'''0.02Ti0.02Al0.25Si1.89O6\n", "biotite: 0.019341 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.512253 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.36Cr0.53Ti0.01O4\n", "...Adding the solid phase hornblende to the assemblage.\n", "...Dropping phase hornblende from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1184.070312\n", "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 \n", "Activity of H2O = 0.504807 Melt fraction = 0.0182712\n", "olivine: 55.138126 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.599245 g, composition (Ca0.16Fe''0.17Mg0.67)3Al2Si3O12\n", "orthopyroxene: 19.055767 g, composition opx Na0.01Ca0.05Fe''0.16Mg1.70Fe'''0.01Ti0.01Al0.14Si1.92O6\n", "clinopyroxene: 13.775658 g, composition cpx Na0.09Ca0.73Fe''0.09Mg0.92Fe'''0.02Ti0.02Al0.25Si1.89O6\n", "biotite: 0.057551 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.470694 g, composition Fe''0.20Mg0.81Fe'''0.08Al1.34Cr0.55Ti0.01O4\n", "...Adding the solid phase hornblende to the assemblage.\n", "...Adding the solid phase apatite to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1174.070312\n", "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 \n", "Activity of H2O = 0.561124 Melt fraction = 0.00688762\n", "olivine: 54.216269 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.034694 g, composition (Ca0.16Fe''0.18Mg0.66)3Al2Si3O12\n", "orthopyroxene: 20.990514 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.14Si1.92O6\n", "clinopyroxene: 11.608657 g, composition cpx Na0.08Ca0.74Fe''0.09Mg0.92Fe'''0.02Ti0.02Al0.24Si1.89O6\n", "hornblende: 3.724794 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.192739 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.433585 g, composition Fe''0.21Mg0.80Fe'''0.09Al1.31Cr0.58Ti0.01O4\n", "apatite: 0.032 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1164.070312\n", "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 \n", "Activity of H2O = 0.618219 Melt fraction = 0.00269047\n", "olivine: 53.845070 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.554240 g, composition (Ca0.16Fe''0.18Mg0.66)3Al2Si3O12\n", "orthopyroxene: 21.717746 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.13Si1.92O6\n", "clinopyroxene: 10.592297 g, composition cpx Na0.08Ca0.75Fe''0.09Mg0.92Fe'''0.02Ti0.02Al0.24Si1.89O6\n", "hornblende: 5.275140 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.234543 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.389359 g, composition Fe''0.22Mg0.80Fe'''0.09Al1.28Cr0.60Ti0.02O4\n", "apatite: 0.043 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1154.070312\n", "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 \n", "Activity of H2O = 0.667748 Melt fraction = 0.00107007\n", "olivine: 53.698077 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.479341 g, composition (Ca0.16Fe''0.18Mg0.66)3Al2Si3O12\n", "orthopyroxene: 21.952849 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.13Si1.92O6\n", "clinopyroxene: 10.087276 g, composition cpx Na0.07Ca0.76Fe''0.09Mg0.92Fe'''0.02Ti0.03Al0.23Si1.89O6\n", "hornblende: 5.952685 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.248069 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.347398 g, composition Fe''0.22Mg0.79Fe'''0.09Al1.25Cr0.63Ti0.02O4\n", "apatite: 0.047 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1144.070312\n", "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 \n", "Activity of H2O = 0.712375 Melt fraction = 0.000273933\n", "olivine: 53.628726 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.530873 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12\n", "orthopyroxene: 22.018594 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.70Fe'''0.01Ti0.01Al0.13Si1.93O6\n", "clinopyroxene: 9.772084 g, composition cpx Na0.07Ca0.76Fe''0.09Mg0.92Fe'''0.01Ti0.03Al0.23Si1.89O6\n", "hornblende: 6.330544 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.253575 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.308516 g, composition Fe''0.23Mg0.79Fe'''0.09Al1.22Cr0.65Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "...Dropping a liquid from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1134.070312\n", "olivine: 53.623437 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.672821 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12\n", "orthopyroxene: 21.948711 g, composition opx Na0.01Ca0.05Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6\n", "clinopyroxene: 9.620627 g, composition cpx Na0.07Ca0.77Fe''0.09Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6\n", "hornblende: 6.475084 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255172 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.273360 g, composition Fe''0.23Mg0.79Fe'''0.10Al1.20Cr0.67Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1124.070312\n", "olivine: 53.658848 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 6.867438 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12\n", "orthopyroxene: 21.794215 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6\n", "clinopyroxene: 9.576887 g, composition cpx Na0.07Ca0.77Fe''0.09Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6\n", "hornblende: 6.474978 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255171 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.241208 g, composition Fe''0.23Mg0.79Fe'''0.10Al1.18Cr0.69Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1114.070312\n", "olivine: 53.692537 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.056088 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12\n", "orthopyroxene: 21.643001 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6\n", "clinopyroxene: 9.536212 g, composition cpx Na0.07Ca0.78Fe''0.08Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6\n", "hornblende: 6.474870 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255170 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.210412 g, composition Fe''0.23Mg0.78Fe'''0.10Al1.16Cr0.71Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1104.070312\n", "olivine: 53.724590 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.239090 g, composition (Ca0.16Fe''0.19Mg0.65)3Al2Si3O12\n", "orthopyroxene: 21.494946 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.12Si1.93O6\n", "clinopyroxene: 9.498400 g, composition cpx Na0.07Ca0.78Fe''0.08Mg0.92Fe'''0.01Ti0.03Al0.22Si1.89O6\n", "hornblende: 6.474759 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255169 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.180899 g, composition Fe''0.24Mg0.78Fe'''0.10Al1.14Cr0.73Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1094.070312\n", "olivine: 53.755083 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.416727 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12\n", "orthopyroxene: 21.349934 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.71Fe'''0.01Ti0.01Al0.11Si1.93O6\n", "clinopyroxene: 9.463270 g, composition cpx Na0.07Ca0.78Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.89O6\n", "hornblende: 6.474645 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255168 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.152600 g, composition Fe''0.24Mg0.78Fe'''0.10Al1.12Cr0.75Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1084.070312\n", "olivine: 53.784087 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.589261 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12\n", "orthopyroxene: 21.207863 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.93O6\n", "clinopyroxene: 9.430657 g, composition cpx Na0.07Ca0.79Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.89O6\n", "hornblende: 6.474529 g, composition NaCa2Mg3.70Fe2+0.30Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255167 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.125454 g, composition Fe''0.24Mg0.77Fe'''0.10Al1.09Cr0.77Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1074.070312\n", "olivine: 53.811662 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.756924 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12\n", "orthopyroxene: 21.068640 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.93O6\n", "clinopyroxene: 9.400412 g, composition cpx Na0.07Ca0.79Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.90O6\n", "hornblende: 6.474411 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255167 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.099406 g, composition Fe''0.25Mg0.77Fe'''0.10Al1.07Cr0.79Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1064.070312\n", "olivine: 53.837863 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 7.919928 g, composition (Ca0.15Fe''0.20Mg0.65)3Al2Si3O12\n", "orthopyroxene: 20.932182 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.94O6\n", "clinopyroxene: 9.372399 g, composition cpx Na0.07Ca0.79Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.90O6\n", "hornblende: 6.474290 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255166 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.074409 g, composition Fe''0.25Mg0.77Fe'''0.10Al1.05Cr0.81Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1054.070312\n", "olivine: 53.862741 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.078465 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12\n", "orthopyroxene: 20.798414 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.11Si1.94O6\n", "clinopyroxene: 9.346495 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.21Si1.90O6\n", "hornblende: 6.474167 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255165 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.050419 g, composition Fe''0.25Mg0.77Fe'''0.10Al1.03Cr0.84Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1044.070312\n", "olivine: 53.886339 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.232706 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12\n", "orthopyroxene: 20.667270 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.72Fe'''0.01Ti0.01Al0.10Si1.94O6\n", "clinopyroxene: 9.322588 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.91Fe'''0.01Ti0.03Al0.20Si1.90O6\n", "hornblende: 6.474041 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255164 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.027398 g, composition Fe''0.26Mg0.76Fe'''0.11Al1.00Cr0.86Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "alphaMELTS at: P 20000.000000, T 1034.070312\n", "olivine: 53.908697 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.382809 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12\n", "orthopyroxene: 20.538689 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.73Fe'''0.01Ti0.01Al0.10Si1.94O6\n", "clinopyroxene: 9.300575 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6\n", "hornblende: 6.473912 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255163 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 1.005312 g, composition Fe''0.26Mg0.76Fe'''0.11Al0.98Cr0.88Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase water to the assemblage.\n", "alphaMELTS at: P 20000.000000, T 1024.070312\n", "olivine: 53.930719 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.530180 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12\n", "orthopyroxene: 20.410844 g, composition opx Na0.01Ca0.04Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.10Si1.94O6\n", "clinopyroxene: 9.282727 g, composition cpx Na0.07Ca0.80Fe''0.08Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6\n", "hornblende: 6.470963 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255163 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.984166 g, composition Fe''0.26Mg0.75Fe'''0.11Al0.95Cr0.90Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.000 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 1014.070313\n", "olivine: 53.992584 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.733804 g, composition (Ca0.15Fe''0.21Mg0.64)3Al2Si3O12\n", "orthopyroxene: 20.201784 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.10Si1.94O6\n", "clinopyroxene: 9.377289 g, composition cpx Na0.08Ca0.81Fe''0.07Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6\n", "hornblende: 6.335583 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255162 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.965428 g, composition Fe''0.27Mg0.75Fe'''0.11Al0.93Cr0.92Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.003 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 1004.070313\n", "olivine: 54.056671 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 8.939458 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12\n", "orthopyroxene: 19.988425 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.10Si1.94O6\n", "clinopyroxene: 9.481035 g, composition cpx Na0.08Ca0.81Fe''0.07Mg0.90Fe'''0.01Ti0.03Al0.20Si1.90O6\n", "hornblende: 6.190129 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255161 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.947424 g, composition Fe''0.27Mg0.75Fe'''0.11Al0.91Cr0.94Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.006 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 994.070313\n", "olivine: 54.123351 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 9.147760 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12\n", "orthopyroxene: 19.770058 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.73Fe'''0.00Ti0.01Al0.09Si1.94O6\n", "clinopyroxene: 9.594762 g, composition cpx Na0.08Ca0.81Fe''0.07Mg0.89Fe'''0.01Ti0.03Al0.20Si1.90O6\n", "hornblende: 6.033525 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255160 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.930133 g, composition Fe''0.27Mg0.74Fe'''0.11Al0.89Cr0.97Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.009 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 984.070313\n", "olivine: 54.193047 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 9.359398 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12\n", "orthopyroxene: 19.545873 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6\n", "clinopyroxene: 9.719410 g, composition cpx Na0.09Ca0.81Fe''0.07Mg0.89Fe'''0.01Ti0.03Al0.20Si1.91O6\n", "hornblende: 5.864521 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255159 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.913535 g, composition Fe''0.28Mg0.74Fe'''0.11Al0.87Cr0.99Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.013 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 974.070313\n", "olivine: 54.266249 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 9.575158 g, composition (Ca0.14Fe''0.22Mg0.64)3Al2Si3O12\n", "orthopyroxene: 19.314933 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6\n", "clinopyroxene: 9.856092 g, composition cpx Na0.09Ca0.80Fe''0.07Mg0.89Fe'''0.01Ti0.03Al0.20Si1.91O6\n", "hornblende: 5.681651 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255159 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.897614 g, composition Fe''0.28Mg0.74Fe'''0.11Al0.85Cr1.01Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.017 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 964.070313\n", "olivine: 54.343534 g, composition (Ca0.00Mg0.89Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 9.795942 g, composition (Ca0.14Fe''0.23Mg0.64)3Al2Si3O12\n", "orthopyroxene: 19.076147 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6\n", "clinopyroxene: 10.006133 g, composition cpx Na0.09Ca0.80Fe''0.07Mg0.88Fe'''0.01Ti0.03Al0.20Si1.91O6\n", "hornblende: 5.483185 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255158 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.882352 g, composition Fe''0.28Mg0.73Fe'''0.11Al0.82Cr1.03Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.021 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 954.070313\n", "olivine: 54.425581 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 10.022804 g, composition (Ca0.14Fe''0.23Mg0.64)3Al2Si3O12\n", "orthopyroxene: 18.828229 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.09Si1.95O6\n", "clinopyroxene: 10.171117 g, composition cpx Na0.10Ca0.80Fe''0.07Mg0.88Fe'''0.01Ti0.03Al0.20Si1.91O6\n", "hornblende: 5.267059 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255157 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.867734 g, composition Fe''0.29Mg0.73Fe'''0.12Al0.80Cr1.05Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.026 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 944.070313\n", "olivine: 54.513207 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 10.256990 g, composition (Ca0.13Fe''0.23Mg0.64)3Al2Si3O12\n", "orthopyroxene: 18.569643 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.08Si1.95O6\n", "clinopyroxene: 10.352961 g, composition cpx Na0.10Ca0.80Fe''0.07Mg0.87Fe'''0.01Ti0.03Al0.20Si1.91O6\n", "hornblende: 5.030797 g, composition NaCa2Mg3.71Fe2+0.29Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255156 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.853744 g, composition Fe''0.29Mg0.72Fe'''0.12Al0.78Cr1.07Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.031 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 934.070313\n", "olivine: 54.607395 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 10.499991 g, composition (Ca0.13Fe''0.23Mg0.64)3Al2Si3O12\n", "orthopyroxene: 18.298542 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.74Fe'''0.00Ti0.01Al0.08Si1.95O6\n", "clinopyroxene: 10.553993 g, composition cpx Na0.10Ca0.80Fe''0.07Mg0.87Fe'''0.01Ti0.03Al0.20Si1.91O6\n", "hornblende: 4.771393 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255156 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.840369 g, composition Fe''0.30Mg0.72Fe'''0.12Al0.76Cr1.09Ti0.02O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.036 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 924.070313\n", "olivine: 54.709348 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 10.753619 g, composition (Ca0.13Fe''0.24Mg0.64)3Al2Si3O12\n", "orthopyroxene: 18.012671 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.95O6\n", "clinopyroxene: 10.777075 g, composition cpx Na0.11Ca0.80Fe''0.07Mg0.87Fe'''0.01Ti0.03Al0.20Si1.91O6\n", "hornblende: 4.485163 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255155 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.827593 g, composition Fe''0.30Mg0.72Fe'''0.12Al0.74Cr1.11Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.042 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 914.070313\n", "olivine: 54.820555 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 11.020101 g, composition (Ca0.13Fe''0.24Mg0.64)3Al2Si3O12\n", "orthopyroxene: 17.709249 g, composition opx Na0.01Ca0.03Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.96O6\n", "clinopyroxene: 11.025755 g, composition cpx Na0.11Ca0.80Fe''0.07Mg0.86Fe'''0.01Ti0.03Al0.21Si1.92O6\n", "hornblende: 4.167540 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255154 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.815404 g, composition Fe''0.30Mg0.71Fe'''0.12Al0.72Cr1.13Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.049 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 904.070313\n", "olivine: 54.942876 g, composition (Ca0.00Mg0.90Fe''0.10Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 11.302216 g, composition (Ca0.12Fe''0.24Mg0.64)3Al2Si3O12\n", "orthopyroxene: 17.384800 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.96O6\n", "clinopyroxene: 11.304489 g, composition cpx Na0.12Ca0.79Fe''0.06Mg0.86Fe'''0.01Ti0.03Al0.21Si1.92O6\n", "hornblende: 3.812798 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255153 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.803787 g, composition Fe''0.31Mg0.71Fe'''0.12Al0.70Cr1.15Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.057 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 894.070313\n", "olivine: 55.078676 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 11.603480 g, composition (Ca0.12Fe''0.24Mg0.64)3Al2Si3O12\n", "orthopyroxene: 17.034920 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.08Si1.96O6\n", "clinopyroxene: 11.618940 g, composition cpx Na0.12Ca0.79Fe''0.06Mg0.85Fe'''0.01Ti0.03Al0.21Si1.92O6\n", "hornblende: 3.413659 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255153 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.792730 g, composition Fe''0.31Mg0.70Fe'''0.12Al0.68Cr1.17Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.065 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 884.070313\n", "olivine: 55.230995 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 11.928412 g, composition (Ca0.12Fe''0.24Mg0.64)3Al2Si3O12\n", "orthopyroxene: 16.653946 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.07Si1.96O6\n", "clinopyroxene: 11.976403 g, composition cpx Na0.13Ca0.79Fe''0.06Mg0.84Fe'''0.01Ti0.03Al0.21Si1.92O6\n", "hornblende: 2.960746 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255152 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.782220 g, composition Fe''0.31Mg0.70Fe'''0.12Al0.66Cr1.19Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.075 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 874.070313\n", "olivine: 55.403811 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 12.282912 g, composition (Ca0.12Fe''0.25Mg0.64)3Al2Si3O12\n", "orthopyroxene: 16.234481 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.75Fe'''0.00Ti0.01Al0.07Si1.96O6\n", "clinopyroxene: 12.386425 g, composition cpx Na0.14Ca0.78Fe''0.06Mg0.84Fe'''0.01Ti0.03Al0.22Si1.92O6\n", "hornblende: 2.441783 g, composition NaCa2Mg3.72Fe2+0.28Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255151 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.772247 g, composition Fe''0.32Mg0.70Fe'''0.12Al0.64Cr1.21Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.086 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 864.070313\n", "olivine: 55.602415 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 12.674821 g, composition (Ca0.11Fe''0.25Mg0.64)3Al2Si3O12\n", "orthopyroxene: 15.766697 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6\n", "clinopyroxene: 12.861706 g, composition cpx Na0.14Ca0.78Fe''0.06Mg0.83Fe'''0.01Ti0.03Al0.22Si1.93O6\n", "hornblende: 1.840431 g, composition NaCa2Mg3.73Fe2+0.27Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255150 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.762800 g, composition Fe''0.32Mg0.69Fe'''0.12Al0.62Cr1.23Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.099 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 854.070313\n", "olivine: 55.833981 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 13.114759 g, composition (Ca0.11Fe''0.25Mg0.64)3Al2Si3O12\n", "orthopyroxene: 15.237290 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6\n", "clinopyroxene: 13.419453 g, composition cpx Na0.15Ca0.77Fe''0.06Mg0.83Fe'''0.01Ti0.03Al0.23Si1.93O6\n", "hornblende: 1.134537 g, composition NaCa2Mg3.73Fe2+0.27Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255150 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.753873 g, composition Fe''0.32Mg0.69Fe'''0.12Al0.61Cr1.25Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.114 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 844.070313\n", "olivine: 56.108425 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 13.617393 g, composition (Ca0.11Fe''0.25Mg0.64)3Al2Si3O12\n", "orthopyroxene: 14.627896 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6\n", "clinopyroxene: 14.083439 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.82Fe'''0.01Ti0.03Al0.23Si1.93O6\n", "hornblende: 0.293464 g, composition NaCa2Mg3.73Fe2+0.27Al1.00Fe3+0.00Al2Si6O22(OH)2\n", "biotite: 0.255149 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.745463 g, composition Fe''0.32Mg0.69Fe'''0.12Al0.59Cr1.27Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.132 g, composition H2O\n", "...Dropping phase hornblende from the assemblage.\n", "alphaMELTS at: P 20000.000000, T 834.070313\n", "olivine: 56.206913 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 13.859489 g, composition (Ca0.11Fe''0.26Mg0.64)3Al2Si3O12\n", "orthopyroxene: 14.359445 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.07Si1.96O6\n", "clinopyroxene: 14.307738 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6\n", "biotite: 0.255148 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.736116 g, composition Fe''0.33Mg0.68Fe'''0.12Al0.57Cr1.29Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.138 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 824.070313\n", "olivine: 56.210501 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 13.960153 g, composition (Ca0.11Fe''0.26Mg0.63)3Al2Si3O12\n", "orthopyroxene: 14.273795 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.76Fe'''0.00Ti0.01Al0.06Si1.97O6\n", "clinopyroxene: 14.298327 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6\n", "biotite: 0.255148 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.726683 g, composition Fe''0.34Mg0.67Fe'''0.12Al0.55Cr1.31Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.138 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 814.070313\n", "olivine: 56.213232 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 14.058245 g, composition (Ca0.11Fe''0.26Mg0.63)3Al2Si3O12\n", "orthopyroxene: 14.189909 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.77Fe'''0.00Ti0.01Al0.06Si1.97O6\n", "clinopyroxene: 14.290043 g, composition cpx Na0.16Ca0.77Fe''0.06Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6\n", "biotite: 0.255147 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.717796 g, composition Fe''0.34Mg0.67Fe'''0.12Al0.53Cr1.33Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.138 g, composition H2O\n", "alphaMELTS at: P 20000.000000, T 804.070313\n", "olivine: 56.215135 g, composition (Ca0.00Mg0.90Fe''0.09Mn0.00Co0.00Ni0.00)2SiO4\n", "garnet: 14.153870 g, composition (Ca0.10Fe''0.27Mg0.63)3Al2Si3O12\n", "orthopyroxene: 14.107733 g, composition opx Na0.01Ca0.02Fe''0.17Mg1.77Fe'''0.00Ti0.01Al0.06Si1.97O6\n", "clinopyroxene: 14.282833 g, composition cpx Na0.16Ca0.77Fe''0.05Mg0.81Fe'''0.01Ti0.03Al0.23Si1.93O6\n", "biotite: 0.255147 g, composition K(Fe''0.00Mg1.00)3AlSi3O10(OH)2\n", "spinel: 0.709425 g, composition Fe''0.35Mg0.66Fe'''0.12Al0.50Cr1.35Ti0.01O4\n", "apatite: 0.049 g, composition Ca5(PO4)3OH\n", "water: 0.138 g, composition H2O\n", "Minimum Temperature reached\n", "Successful return from alphamelts...\n", "\n", "Choose:\n", " 1. Read MELTS file to set composition of system\n", " 2. Twiddle starting or continuation parameters\n", " 3. Single (batch) calculation\n", " 4. Execute (follow path, mineral isograd or melt contour)\n", " 5. Set fO2 buffer\n", " 6. Set H2O (ppm) or aH2O\n", " 7. Impose initial entropy, enthalpy or volume\n", " 8. Adjust solid phase setting(s)\n", " 9. Turn liquid on / off\n", "10. Turn phase diagram mode on / off\n", "11. Update state using restart file\n", "12. Update composition using MELTS file\n", "13. Write out restart file\n", "14. Write out MELTS file\n", "15. Write thermodynamic output for all phases\n", "16. Calculate integrated melt and output file(s)\n", "17. Fit parental melt composition (amoeba)\n", "18. Cumulate Invertor (not yet implemented)\n", "-1. Turn off menu display for options 1-18\n", " 0. QUIT\n", "Your choice: \n", "\n" ] }, { "data": { "text/plain": [ "Dict{Any,Any} with 15 entries:\n", " \"Temperature\" => [1624.07, 1614.07, 1604.07, 1594.07, 1584.07, 1574.07, 1…\n", " \"water\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.…\n", " \"orthopyroxene_0\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 7.96064, 11.4105, 14.2824…\n", " \"elements\" => Any[\"Pressure\", \"Temperature\", \"mass\", \"liquid_0\", \"oliv…\n", " \"clinopyroxene_0\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 11…\n", " \"mass\" => [99.9575, 99.9544, 99.9515, 99.949, 99.9467, 99.9445, 99…\n", " \"olivine_0\" => [0.0, 8.88223, 16.6072, 23.3962, 29.1114, 33.9796, 36.01…\n", " \"hornblende_0\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 3.…\n", " \"garnet_0\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 11…\n", " \"spinel_0\" => [0.0, 0.0, 0.0, 0.134547, 0.235384, 0.321473, 0.464989, …\n", " \"liquid_0\" => [99.9575, 91.0722, 83.3443, 76.4183, 70.5999, 65.6434, 5…\n", " \"apatite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.…\n", " \"Pressure\" => [20000.0, 20000.0, 20000.0, 20000.0, 20000.0, 20000.0, 2…\n", " \"\" => [NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN, NaN … Na…\n", " \"biotite_0\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.…" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## --- # # # # # # # # # # # pMelts equil. batch melting # # # # # # # # # # # #\n", "\n", "# Conditions\n", "P_range = [20000,20000]\n", "T_range = [1700,800]\n", "# Starting composition\n", "elements = [\"SiO2\", \"TiO2\",\"Al2O3\",\"Fe2O3\",\"Cr2O3\", \"FeO\", \"MnO\", \"MgO\", \"NiO\", \"CoO\", \"CaO\", \"Na2O\", \"K2O\", \"P2O5\", \"H2O\",]\n", "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)\n", "# Run simulation\n", "melts_configure(meltspath, scratchdir, composition, elements, T_range, P_range,\n", " batchstring=\"1\\nsc.melts\\n10\\n1\\n3\\n1\\nliquid\\n1\\n1.0\\n0\\n10\\n0\\n4\\n0\\n\",\n", " dT=-10, dP=0, index=1, version=\"pMELTS\",mode=\"isobaric\",fo2path=\"FMQ\")\n", "\n", "# Read results\n", "melt_comp = melts_query_liquid(scratchdir, index=1)\n", "solid_comp = melts_query_solid(scratchdir, index=1)\n", "modes = melts_query_modes(scratchdir, index=1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Plot results" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "0\n", "\n", "\n", "20\n", "\n", "\n", "40\n", "\n", "\n", "60\n", "\n", "\n", "80\n", "\n", "\n", "100\n", "\n", "\n", "0\n", "\n", "\n", "10\n", "\n", "\n", "20\n", "\n", "\n", "30\n", "\n", "\n", "40\n", "\n", "\n", "\n", "\n", "Percent melt\n", "\n", "\n", "Abudance (wt. %) in melt\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "SiO2\n", "\n", "\n", "\n", "Al2O3\n", "\n", "\n", "\n", "CaO\n", "\n", "\n", "\n", "MgO\n", "\n", "\n", "\n", "FeO\n", "\n", "\n", "\n", "Na2O\n", "\n", "\n", "\n", "K2O\n", "\n", "\n" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## --- Plot melt composition\n", "melt_comp = melts_query_liquid(scratchdir, index=1)\n", "\n", "h = plot(xlabel=\"Percent melt\",ylabel=\"Abundance (wt. %) in melt\")\n", "for e in [\"SiO2\",\"Al2O3\",\"CaO\",\"MgO\",\"FeO\",\"Na2O\",\"K2O\"]\n", " plot!(h,melt_comp[\"mass\"],melt_comp[e],label=e)\n", "end\n", "plot!(h,xlims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:right)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "scrolled": false }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "0\n", "\n", "\n", "20\n", "\n", "\n", "40\n", "\n", "\n", "60\n", "\n", "\n", "80\n", "\n", "\n", "100\n", "\n", "\n", "0\n", "\n", "\n", "10\n", "\n", "\n", "20\n", "\n", "\n", "30\n", "\n", "\n", "40\n", "\n", "\n", "50\n", "\n", "\n", "\n", "\n", "Percent melt\n", "\n", "\n", "Abudance (wt. %) in solid\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "SiO2\n", "\n", "\n", "\n", "Al2O3\n", "\n", "\n", "\n", "CaO\n", "\n", "\n", "\n", "MgO\n", "\n", "\n", "\n", "FeO\n", "\n", "\n", "\n", "Na2O\n", "\n", "\n", "\n", "K2O\n", "\n", "\n" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## --- Plot solid composition\n", "solid_comp = melts_query_solid(scratchdir, index=1)\n", "\n", "h = plot(xlabel=\"Percent melt\",ylabel=\"Abundance (wt. %) in solid\")\n", "for e in [\"SiO2\",\"Al2O3\",\"CaO\",\"MgO\",\"FeO\",\"Na2O\",\"K2O\"]\n", " plot!(h,100 .- solid_comp[\"mass\"],solid_comp[e],label=e)\n", "end\n", "plot!(h,xlims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:right)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "scrolled": false }, "outputs": [ { "data": { "image/svg+xml": [ "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", " \n", " \n", " \n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "800\n", "\n", "\n", "1000\n", "\n", "\n", "1200\n", "\n", "\n", "1400\n", "\n", "\n", "1600\n", "\n", "\n", "0\n", "\n", "\n", "20\n", "\n", "\n", "40\n", "\n", "\n", "60\n", "\n", "\n", "80\n", "\n", "\n", "100\n", "\n", "\n", "\n", "\n", "Temperature (C)\n", "\n", "\n", "Abudance (wt. %)\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "liquid_0\n", "\n", "\n", "\n", "olivine_0\n", "\n", "\n", "\n", "garnet_0\n", "\n", "\n", "\n", "orthopyroxene_0\n", "\n", "\n", "\n", "clinopyroxene_0\n", "\n", "\n", "\n", "hornblende_0\n", "\n", "\n", "\n", "biotite_0\n", "\n", "\n", "\n", "spinel_0\n", "\n", "\n", "\n", "apatite\n", "\n", "\n", "\n", "water\n", "\n", "\n" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## --- Plot phase modes\n", "modes = melts_query_modes(scratchdir, index=1)\n", "\n", "h = plot(xlabel=\"Temperature (C)\",ylabel=\"Abundance (wt. %)\")\n", "for m in modes[\"elements\"][4:end]\n", " plot!(h,modes[\"Temperature\"],modes[m],label=m)\n", "end\n", "plot!(h,ylims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:topleft)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Julia 1.7.1", "language": "julia", "name": "julia-1.7" }, "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "1.7.1" } }, "nbformat": 4, "nbformat_minor": 2 }