{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# MELTS zircon saturation example\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 the StatGeochem package\n", "using StatGeochem\n", "using Plots" ] }, { "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 = abspath(\"/$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/\";" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Run alphamelts" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ALPHAMELTS_DELTAP 0.0\n", "ALPHAMELTS_MODE isobaric\n", "ALPHAMELTS_CELSIUS_OUTPUT true\n", "ALPHAMELTS_DELTAT -7.0\n", "ALPHAMELTS_MAXT 3000.0\n", "ALPHAMELTS_MINT 450.0\n", "ALPHAMELTS_MINP 2.0\n", "ALPHAMELTS_MAXP 90000.0\n", "ALPHAMELTS_VERSION pMELTS\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 61.3266'\n", "Successfully read 'initial composition: tio2 0.3793'\n", "Successfully read 'initial composition: al2o3 14.8372'\n", "Successfully read 'initial composition: fe2o3 3.3769'\n", "Successfully read 'initial composition: cr2o3 0.0'\n", "Successfully read 'initial composition: feo 0.0'\n", "Successfully read 'initial composition: mno 0.0668'\n", "Successfully read 'initial composition: mgo 1.6031'\n", "Successfully read 'initial composition: nio 0.0'\n", "Successfully read 'initial composition: coo 0.0'\n", "Successfully read 'initial composition: cao 3.28'\n", "Successfully read 'initial composition: na2o 3.4288'\n", "Successfully read 'initial composition: k2o 3.4877'\n", "Successfully read 'initial composition: p2o5 0.2275'\n", "Successfully read 'initial composition: h2o 7.9856'\n", "Successfully read 'initial temperature: 1700.0'\n", "Successfully read 'initial pressure: 5800.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 5800.000000, T 1700.000000\n", "liquid: 99.702 g 61.51 0.38 14.88 0.41 2.68 0.07 1.61 3.29 3.44 3.50 0.23 8.01 \n", "Activity of H2O = 0.561749 Melt fraction = 1\n", "...Adding the solid phase biotite to the assemblage.\n", "...Dropping phase biotite from the assemblage.\n", "...Adding the solid phase biotite to the assemblage.\n", "...Dropping phase biotite from the assemblage.\n", "Isograd solution at: P 5800.000000, T 1016.507813\n", "liquid: 99.707 g 61.51 0.38 14.88 0.46 2.64 0.07 1.61 3.29 3.44 3.50 0.23 8.01 \n", "Activity of H2O = 0.776445 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: Initial alphaMELTS calculation at: P 5800.000000, T 1016.507813\n", "liquid: 99.707 g 61.51 0.38 14.88 0.46 2.64 0.07 1.61 3.29 3.44 3.50 0.23 8.01 \n", "Activity of H2O = 0.776445 Melt fraction = 1\n", "...Adding the solid phase biotite to the assemblage.\n", "alphaMELTS at: P 5800.000000, T 1009.507813\n", "liquid: 99.276 g 61.59 0.38 14.89 0.46 2.64 0.07 1.49 3.30 3.45 3.46 0.23 8.03 \n", "Activity of H2O = 0.780322 Melt fraction = 0.995674\n", "biotite: 0.431308 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 1002.507813\n", "liquid: 98.876 g 61.66 0.38 14.90 0.47 2.65 0.07 1.38 3.32 3.47 3.43 0.23 8.04 \n", "Activity of H2O = 0.784258 Melt fraction = 0.991663\n", "biotite: 0.831234 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 995.507813\n", "liquid: 98.503 g 61.73 0.39 14.91 0.47 2.66 0.07 1.28 3.33 3.48 3.40 0.23 8.05 \n", "Activity of H2O = 0.788257 Melt fraction = 0.987921\n", "biotite: 1.204323 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 988.507813\n", "liquid: 98.156 g 61.80 0.39 14.92 0.48 2.66 0.07 1.18 3.34 3.49 3.38 0.23 8.07 \n", "Activity of H2O = 0.792322 Melt fraction = 0.984431\n", "biotite: 1.552335 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 981.507813\n", "liquid: 97.831 g 61.86 0.39 14.93 0.48 2.67 0.07 1.09 3.35 3.50 3.35 0.23 8.08 \n", "Activity of H2O = 0.796458 Melt fraction = 0.981175\n", "biotite: 1.876965 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 974.507813\n", "liquid: 97.529 g 61.92 0.39 14.94 0.48 2.67 0.07 1.00 3.36 3.52 3.32 0.23 8.09 \n", "Activity of H2O = 0.800668 Melt fraction = 0.978138\n", "biotite: 2.179847 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 967.507813\n", "liquid: 97.246 g 61.97 0.39 14.95 0.49 2.67 0.07 0.92 3.37 3.53 3.30 0.23 8.10 \n", "Activity of H2O = 0.804958 Melt fraction = 0.975302\n", "biotite: 2.462567 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 960.507813\n", "liquid: 96.982 g 62.02 0.39 14.96 0.49 2.67 0.07 0.85 3.38 3.54 3.28 0.23 8.11 \n", "Activity of H2O = 0.809331 Melt fraction = 0.972654\n", "biotite: 2.726670 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 953.507813\n", "liquid: 96.735 g 62.07 0.39 14.96 0.49 2.67 0.07 0.78 3.39 3.54 3.26 0.24 8.12 \n", "Activity of H2O = 0.813791 Melt fraction = 0.970177\n", "biotite: 2.973673 g, composition K(Fe''0.01Mg0.99)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 946.507813\n", "liquid: 96.504 g 62.12 0.39 14.97 0.50 2.67 0.07 0.72 3.40 3.55 3.24 0.24 8.13 \n", "Activity of H2O = 0.818342 Melt fraction = 0.967856\n", "biotite: 3.205078 g, composition K(Fe''0.02Mg0.98)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 939.507813\n", "liquid: 96.287 g 62.16 0.39 14.98 0.50 2.67 0.07 0.66 3.41 3.56 3.22 0.24 8.14 \n", "Activity of H2O = 0.822987 Melt fraction = 0.965676\n", "biotite: 3.422385 g, composition K(Fe''0.02Mg0.98)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 932.507813\n", "liquid: 96.083 g 62.20 0.39 14.98 0.51 2.67 0.07 0.60 3.41 3.57 3.21 0.24 8.15 \n", "Activity of H2O = 0.827731 Melt fraction = 0.963623\n", "biotite: 3.627111 g, composition K(Fe''0.02Mg0.98)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 925.507813\n", "liquid: 95.889 g 62.24 0.40 14.99 0.51 2.66 0.07 0.55 3.42 3.58 3.19 0.24 8.16 \n", "Activity of H2O = 0.832576 Melt fraction = 0.961681\n", "biotite: 3.820812 g, composition K(Fe''0.03Mg0.97)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 918.507813\n", "liquid: 95.705 g 62.28 0.40 15.00 0.51 2.65 0.07 0.51 3.43 3.58 3.17 0.24 8.16 \n", "Activity of H2O = 0.837527 Melt fraction = 0.959833\n", "biotite: 4.005108 g, composition K(Fe''0.03Mg0.97)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 911.507813\n", "liquid: 95.529 g 62.32 0.40 15.00 0.52 2.64 0.07 0.46 3.43 3.59 3.16 0.24 8.17 \n", "Activity of H2O = 0.842588 Melt fraction = 0.958061\n", "biotite: 4.181724 g, composition K(Fe''0.03Mg0.97)3AlSi3O10(OH)2\n", "alphaMELTS at: P 5800.000000, T 904.507813\n", "liquid: 95.358 g 62.36 0.40 15.01 0.52 2.63 0.07 0.42 3.44 3.60 3.15 0.24 8.18 \n", "Activity of H2O = 0.847761 Melt fraction = 0.956348\n", "biotite: 4.352531 g, composition K(Fe''0.04Mg0.96)3AlSi3O10(OH)2\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 897.507813\n", "liquid: 95.191 g 62.39 0.40 15.01 0.52 2.62 0.07 0.38 3.45 3.60 3.13 0.24 8.19 \n", "Activity of H2O = 0.853052 Melt fraction = 0.954673\n", "biotite: 4.519606 g, composition K(Fe''0.04Mg0.96)3AlSi3O10(OH)2\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase feldspar to the assemblage.\n", "...Dropping phase feldspar from the assemblage.\n", "...Adding the solid phase feldspar to the assemblage.\n", "...Dropping phase feldspar from the assemblage.\n", "...Adding the solid phase clinopyroxene to the assemblage.\n", "...Adding the solid phase apatite to the assemblage.\n", "alphaMELTS at: P 5800.000000, T 890.507813\n", "liquid: 94.983 g 62.45 0.39 15.02 0.52 2.60 0.07 0.35 3.44 3.61 3.12 0.23 8.20 \n", "Activity of H2O = 0.858716 Melt fraction = 0.952588\n", "clinopyroxene: 0.036804 g, composition cpx Na0.06Ca0.92Fe''0.25Mg0.33Fe'''0.07Ti0.40Al0.74Si1.23O6\n", "biotite: 4.679182 g, composition K(Fe''0.05Mg0.95)3AlSi3O10(OH)2\n", "apatite: 0.012 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase feldspar to the assemblage.\n", "alphaMELTS at: P 5800.000000, T 883.507813\n", "liquid: 93.581 g 62.54 0.36 14.98 0.52 2.59 0.07 0.32 3.41 3.64 3.01 0.24 8.32 \n", "Activity of H2O = 0.869644 Melt fraction = 0.938524\n", "clinopyroxene: 0.343151 g, composition cpx Na0.06Ca0.92Fe''0.26Mg0.32Fe'''0.07Ti0.40Al0.73Si1.24O6\n", "biotite: 4.756278 g, composition K(Fe''0.06Mg0.94)3AlSi3O10(OH)2\n", "feldspar: 1.015215 g, composition K0.83Na0.16Ca0.01Al1.01Si2.99O8\n", "apatite: 0.015 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 876.507813\n", "liquid: 91.556 g 62.61 0.32 14.92 0.53 2.60 0.07 0.31 3.41 3.68 2.81 0.24 8.50 \n", "Activity of H2O = 0.883285 Melt fraction = 0.918208\n", "clinopyroxene: 0.648666 g, composition cpx Na0.07Ca0.92Fe''0.27Mg0.32Fe'''0.07Ti0.39Al0.72Si1.25O6\n", "biotite: 4.807883 g, composition K(Fe''0.06Mg0.94)3AlSi3O10(OH)2\n", "feldspar: 2.684726 g, composition K0.82Na0.17Ca0.01Al1.01Si2.99O8\n", "apatite: 0.014 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 869.507813\n", "liquid: 89.683 g 62.67 0.29 14.85 0.53 2.60 0.07 0.29 3.41 3.72 2.63 0.25 8.67 \n", "Activity of H2O = 0.896833 Melt fraction = 0.89943\n", "clinopyroxene: 0.929735 g, composition cpx Na0.07Ca0.92Fe''0.27Mg0.31Fe'''0.07Ti0.39Al0.71Si1.26O6\n", "biotite: 4.866644 g, composition K(Fe''0.07Mg0.93)3AlSi3O10(OH)2\n", "feldspar: 4.217472 g, composition K0.81Na0.18Ca0.01Al1.01Si2.99O8\n", "apatite: 0.014 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 862.507813\n", "liquid: 87.940 g 62.73 0.26 14.80 0.54 2.60 0.08 0.28 3.41 3.75 2.46 0.25 8.84 \n", "Activity of H2O = 0.910372 Melt fraction = 0.881942\n", "clinopyroxene: 1.190842 g, composition cpx Na0.07Ca0.91Fe''0.28Mg0.31Fe'''0.07Ti0.38Al0.70Si1.27O6\n", "biotite: 4.932307 g, composition K(Fe''0.08Mg0.92)3AlSi3O10(OH)2\n", "feldspar: 5.633933 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "apatite: 0.015 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 855.507813\n", "liquid: 86.302 g 62.79 0.23 14.74 0.54 2.60 0.08 0.26 3.42 3.78 2.30 0.26 9.01 \n", "Activity of H2O = 0.924002 Melt fraction = 0.865522\n", "clinopyroxene: 1.437056 g, composition cpx Na0.08Ca0.91Fe''0.29Mg0.31Fe'''0.07Ti0.37Al0.69Si1.28O6\n", "biotite: 5.004499 g, composition K(Fe''0.09Mg0.91)3AlSi3O10(OH)2\n", "feldspar: 6.951835 g, composition K0.79Na0.20Ca0.01Al1.01Si2.99O8\n", "apatite: 0.016 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase leucite ss to the assemblage.\n", "alphaMELTS at: P 5800.000000, T 848.507813\n", "liquid: 82.903 g 63.07 0.20 14.53 0.54 2.62 0.08 0.25 3.45 3.61 2.17 0.25 9.22 \n", "Activity of H2O = 0.942133 Melt fraction = 0.831431\n", "clinopyroxene: 1.727847 g, composition cpx Na0.08Ca0.91Fe''0.30Mg0.30Fe'''0.08Ti0.36Al0.66Si1.31O6\n", "biotite: 5.123926 g, composition K(Fe''0.10Mg0.90)3AlSi3O10(OH)2\n", "feldspar: 8.216228 g, composition K0.80Na0.20Ca0.01Al1.01Si2.99O8\n", "leucite ss: 1.689022 g, composition K0.01Na0.99AlSi2O5.12(OH)1.76\n", "apatite: 0.051 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 841.507813\n", "liquid: 79.665 g 63.37 0.18 14.32 0.55 2.63 0.08 0.23 3.49 3.43 2.04 0.24 9.43 \n", "Activity of H2O = 0.960138 Melt fraction = 0.798963\n", "clinopyroxene: 2.031540 g, composition cpx Na0.08Ca0.91Fe''0.32Mg0.30Fe'''0.08Ti0.34Al0.63Si1.34O6\n", "biotite: 5.248428 g, composition K(Fe''0.11Mg0.89)3AlSi3O10(OH)2\n", "feldspar: 9.300912 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 3.378858 g, composition K0.01Na0.99AlSi2O5.11(OH)1.78\n", "apatite: 0.086 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...One of the solid phases has undergone phase separation.\n", "alphaMELTS at: P 5800.000000, T 834.507813\n", "liquid: 76.474 g 63.68 0.16 14.17 0.54 2.60 0.09 0.22 3.43 3.26 1.94 0.23 9.68 \n", "Activity of H2O = 0.981413 Melt fraction = 0.766963\n", "clinopyroxene: 2.002338 g, composition cpx Na0.08Ca0.91Fe''0.32Mg0.30Fe'''0.08Ti0.34Al0.63Si1.34O6\n", "clinopyroxene: 0.674243 g, composition cpx Na0.08Ca0.87Fe''0.37Mg0.46Fe'''0.08Ti0.09Al0.26Si1.78O6\n", "biotite: 5.229945 g, composition K(Fe''0.13Mg0.87)3AlSi3O10(OH)2\n", "feldspar: 10.279334 g, composition K0.81Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 4.935999 g, composition K0.01Na0.99AlSi2O5.11(OH)1.79\n", "apatite: 0.114 g, composition Ca5(PO4)3OH\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase water to the assemblage.\n", "alphaMELTS at: P 5800.000000, T 827.507813\n", "liquid: 73.367 g 64.03 0.15 14.05 0.53 2.54 0.09 0.20 3.35 3.09 1.84 0.23 9.90 \n", "Activity of H2O = 1.00016 Melt fraction = 0.735811\n", "clinopyroxene: 1.883447 g, composition cpx Na0.08Ca0.91Fe''0.33Mg0.29Fe'''0.08Ti0.34Al0.63Si1.34O6\n", "clinopyroxene: 1.538316 g, composition cpx Na0.08Ca0.87Fe''0.38Mg0.45Fe'''0.08Ti0.09Al0.25Si1.79O6\n", "biotite: 5.164668 g, composition K(Fe''0.14Mg0.86)3AlSi3O10(OH)2\n", "feldspar: 11.212812 g, composition K0.81Na0.18Ca0.01Al1.01Si2.99O8\n", "leucite ss: 6.368844 g, composition K0.01Na0.99AlSi2O5.10(OH)1.80\n", "apatite: 0.140 g, composition Ca5(PO4)3OH\n", "water: 0.034 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 820.507813\n", "liquid: 69.333 g 64.64 0.14 14.02 0.51 2.45 0.10 0.19 3.13 2.92 1.74 0.23 9.95 \n", "Activity of H2O = 1.00013 Melt fraction = 0.695353\n", "clinopyroxene: 1.582366 g, composition cpx Na0.08Ca0.91Fe''0.33Mg0.28Fe'''0.08Ti0.35Al0.64Si1.34O6\n", "clinopyroxene: 3.144754 g, composition cpx Na0.09Ca0.87Fe''0.39Mg0.45Fe'''0.08Ti0.09Al0.25Si1.80O6\n", "biotite: 4.823611 g, composition K(Fe''0.15Mg0.85)3AlSi3O10(OH)2\n", "feldspar: 12.571806 g, composition K0.81Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 7.780767 g, composition K0.01Na0.99AlSi2O5.10(OH)1.80\n", "apatite: 0.167 g, composition Ca5(PO4)3OH\n", "water: 0.306 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 813.507813\n", "liquid: 65.967 g 65.20 0.14 14.00 0.48 2.35 0.10 0.19 2.93 2.75 1.64 0.22 10.00 \n", "Activity of H2O = 1.00011 Melt fraction = 0.661595\n", "clinopyroxene: 1.373531 g, composition cpx Na0.08Ca0.90Fe''0.34Mg0.28Fe'''0.08Ti0.35Al0.64Si1.33O6\n", "clinopyroxene: 4.434305 g, composition cpx Na0.09Ca0.86Fe''0.40Mg0.45Fe'''0.08Ti0.08Al0.24Si1.80O6\n", "biotite: 4.558973 g, composition K(Fe''0.16Mg0.84)3AlSi3O10(OH)2\n", "feldspar: 13.699524 g, composition K0.81Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 8.955902 g, composition K0.01Na0.99AlSi2O5.10(OH)1.81\n", "apatite: 0.189 g, composition Ca5(PO4)3OH\n", "water: 0.530 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 806.507813\n", "liquid: 63.109 g 65.72 0.13 13.98 0.46 2.26 0.11 0.18 2.75 2.60 1.55 0.22 10.04 \n", "Activity of H2O = 1.00008 Melt fraction = 0.632936\n", "clinopyroxene: 1.230601 g, composition cpx Na0.08Ca0.90Fe''0.34Mg0.27Fe'''0.07Ti0.35Al0.64Si1.33O6\n", "clinopyroxene: 5.487351 g, composition cpx Na0.09Ca0.86Fe''0.41Mg0.45Fe'''0.07Ti0.08Al0.23Si1.81O6\n", "biotite: 4.352237 g, composition K(Fe''0.17Mg0.83)3AlSi3O10(OH)2\n", "feldspar: 14.648592 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 9.953361 g, composition K0.01Na0.99AlSi2O5.09(OH)1.82\n", "apatite: 0.208 g, composition Ca5(PO4)3OH\n", "water: 0.719 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase rhm oxide to the assemblage.\n", "...Dropping phase rhm oxide from the assemblage.\n", "...Adding the solid phase rhm oxide to the assemblage.\n", "...Dropping phase rhm oxide from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 799.507813\n", "liquid: 60.648 g 66.22 0.12 13.96 0.44 2.16 0.11 0.17 2.59 2.45 1.47 0.22 10.09 \n", "Activity of H2O = 1.00006 Melt fraction = 0.608266\n", "clinopyroxene: 1.135641 g, composition cpx Na0.08Ca0.90Fe''0.34Mg0.27Fe'''0.07Ti0.36Al0.65Si1.32O6\n", "clinopyroxene: 6.359241 g, composition cpx Na0.09Ca0.86Fe''0.42Mg0.45Fe'''0.07Ti0.08Al0.23Si1.82O6\n", "biotite: 4.190382 g, composition K(Fe''0.18Mg0.82)3AlSi3O10(OH)2\n", "feldspar: 15.455866 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 10.813941 g, composition K0.01Na0.99AlSi2O5.09(OH)1.82\n", "apatite: 0.224 g, composition Ca5(PO4)3OH\n", "water: 0.879 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase rhm oxide to the assemblage.\n", "...Dropping phase rhm oxide from the assemblage.\n", "...Adding the solid phase rhm oxide to the assemblage.\n", "...Dropping phase rhm oxide from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 792.507813\n", "liquid: 58.506 g 66.69 0.11 13.94 0.42 2.06 0.11 0.16 2.44 2.32 1.39 0.22 10.13 \n", "Activity of H2O = 1.00005 Melt fraction = 0.58679\n", "clinopyroxene: 1.076013 g, composition cpx Na0.08Ca0.90Fe''0.35Mg0.26Fe'''0.07Ti0.36Al0.65Si1.32O6\n", "clinopyroxene: 7.089433 g, composition cpx Na0.09Ca0.86Fe''0.42Mg0.44Fe'''0.07Ti0.07Al0.22Si1.82O6\n", "biotite: 4.063905 g, composition K(Fe''0.19Mg0.81)3AlSi3O10(OH)2\n", "feldspar: 16.148100 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 11.566786 g, composition K0.01Na0.99AlSi2O5.09(OH)1.83\n", "apatite: 0.238 g, composition Ca5(PO4)3OH\n", "water: 1.017 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase rhm oxide to the assemblage.\n", "alphaMELTS at: P 5800.000000, T 785.507813\n", "liquid: 56.512 g 67.07 0.10 14.02 0.40 1.93 0.12 0.15 2.34 2.18 1.32 0.21 10.16 \n", "Activity of H2O = 1.00003 Melt fraction = 0.566796\n", "clinopyroxene: 0.416188 g, composition cpx Na0.08Ca0.90Fe''0.35Mg0.26Fe'''0.07Ti0.37Al0.66Si1.31O6\n", "clinopyroxene: 8.254124 g, composition cpx Na0.09Ca0.85Fe''0.43Mg0.45Fe'''0.07Ti0.07Al0.22Si1.83O6\n", "biotite: 3.807323 g, composition K(Fe''0.19Mg0.81)3AlSi3O10(OH)2\n", "feldspar: 16.885056 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 12.286933 g, composition K0.01Na0.99AlSi2O5.08(OH)1.83\n", "rhm oxide: 0.137389 g, composition Mn0.01Fe''0.88Mg0.04Fe'''0.15Ti0.93O3\n", "apatite: 0.252 g, composition Ca5(PO4)3OH\n", "water: 1.154 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Dropping phase clinopyroxene from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 778.507813\n", "liquid: 54.779 g 67.44 0.10 14.08 0.37 1.80 0.12 0.15 2.24 2.05 1.25 0.21 10.19 \n", "Activity of H2O = 1.00002 Melt fraction = 0.549411\n", "clinopyroxene: 9.129142 g, composition cpx Na0.09Ca0.85Fe''0.43Mg0.45Fe'''0.07Ti0.07Al0.21Si1.84O6\n", "biotite: 3.611659 g, composition K(Fe''0.19Mg0.81)3AlSi3O10(OH)2\n", "feldspar: 17.503284 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 12.914241 g, composition K0.01Na0.99AlSi2O5.08(OH)1.84\n", "rhm oxide: 0.233720 g, composition Mn0.01Fe''0.88Mg0.04Fe'''0.14Ti0.93O3\n", "apatite: 0.263 g, composition Ca5(PO4)3OH\n", "water: 1.271 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 771.507813\n", "liquid: 53.287 g 67.83 0.09 14.09 0.36 1.70 0.12 0.14 2.13 1.93 1.18 0.21 10.22 \n", "Activity of H2O = 1.00001 Melt fraction = 0.534457\n", "clinopyroxene: 9.569315 g, composition cpx Na0.09Ca0.85Fe''0.44Mg0.45Fe'''0.06Ti0.06Al0.21Si1.85O6\n", "biotite: 3.519006 g, composition K(Fe''0.20Mg0.80)3AlSi3O10(OH)2\n", "feldspar: 17.979324 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 13.455435 g, composition K0.01Na0.99AlSi2O5.08(OH)1.84\n", "rhm oxide: 0.255811 g, composition Mn0.01Fe''0.89Mg0.04Fe'''0.14Ti0.93O3\n", "apatite: 0.272 g, composition Ca5(PO4)3OH\n", "water: 1.365 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 764.507813\n", "liquid: 51.950 g 68.19 0.08 14.09 0.34 1.60 0.13 0.13 2.03 1.82 1.12 0.21 10.26 \n", "Activity of H2O = 1 Melt fraction = 0.521049\n", "clinopyroxene: 9.963522 g, composition cpx Na0.09Ca0.85Fe''0.44Mg0.45Fe'''0.06Ti0.06Al0.20Si1.85O6\n", "biotite: 3.440781 g, composition K(Fe''0.21Mg0.79)3AlSi3O10(OH)2\n", "feldspar: 18.396548 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 13.947181 g, composition K0.00Na1.00AlSi2O5.07(OH)1.85\n", "rhm oxide: 0.275492 g, composition Mn0.01Fe''0.89Mg0.03Fe'''0.14Ti0.93O3\n", "apatite: 0.280 g, composition Ca5(PO4)3OH\n", "water: 1.449 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 757.507813\n", "liquid: 50.745 g 68.54 0.07 14.10 0.32 1.50 0.13 0.12 1.94 1.72 1.06 0.21 10.28 \n", "Activity of H2O = 0.999994 Melt fraction = 0.508975\n", "clinopyroxene: 10.317670 g, composition cpx Na0.09Ca0.85Fe''0.45Mg0.45Fe'''0.06Ti0.06Al0.20Si1.86O6\n", "biotite: 3.374627 g, composition K(Fe''0.22Mg0.78)3AlSi3O10(OH)2\n", "feldspar: 18.762584 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 14.397376 g, composition K0.00Na1.00AlSi2O5.07(OH)1.86\n", "rhm oxide: 0.293191 g, composition Mn0.01Fe''0.89Mg0.03Fe'''0.14Ti0.93O3\n", "apatite: 0.286 g, composition Ca5(PO4)3OH\n", "water: 1.524 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 750.507813\n", "liquid: 49.656 g 68.87 0.07 14.10 0.30 1.41 0.13 0.12 1.85 1.62 1.01 0.21 10.31 \n", "Activity of H2O = 0.999988 Melt fraction = 0.498062\n", "clinopyroxene: 10.636691 g, composition cpx Na0.09Ca0.84Fe''0.46Mg0.45Fe'''0.06Ti0.05Al0.19Si1.87O6\n", "biotite: 3.318685 g, composition K(Fe''0.22Mg0.78)3AlSi3O10(OH)2\n", "feldspar: 19.083673 g, composition K0.80Na0.19Ca0.01Al1.01Si2.99O8\n", "leucite ss: 14.812332 g, composition K0.00Na1.00AlSi2O5.07(OH)1.86\n", "rhm oxide: 0.309198 g, composition Mn0.01Fe''0.89Mg0.03Fe'''0.13Ti0.93O3\n", "apatite: 0.292 g, composition Ca5(PO4)3OH\n", "water: 1.590 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 743.507813\n", "liquid: 48.669 g 69.19 0.06 14.11 0.28 1.32 0.13 0.11 1.77 1.52 0.96 0.21 10.33 \n", "Activity of H2O = 0.999983 Melt fraction = 0.488166\n", "clinopyroxene: 10.924728 g, composition cpx Na0.09Ca0.84Fe''0.46Mg0.45Fe'''0.06Ti0.05Al0.19Si1.87O6\n", "biotite: 3.271454 g, composition K(Fe''0.23Mg0.77)3AlSi3O10(OH)2\n", "feldspar: 19.365018 g, composition K0.80Na0.19Ca0.00Al1.00Si3.00O8\n", "leucite ss: 15.197149 g, composition K0.00Na1.00AlSi2O5.07(OH)1.87\n", "rhm oxide: 0.323729 g, composition Mn0.01Fe''0.89Mg0.03Fe'''0.13Ti0.93O3\n", "apatite: 0.297 g, composition Ca5(PO4)3OH\n", "water: 1.650 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 736.507813\n", "liquid: 47.772 g 69.49 0.06 14.12 0.27 1.23 0.14 0.10 1.70 1.43 0.91 0.21 10.35 \n", "Activity of H2O = 0.999979 Melt fraction = 0.479169\n", "clinopyroxene: 11.185283 g, composition cpx Na0.09Ca0.84Fe''0.47Mg0.45Fe'''0.05Ti0.05Al0.18Si1.88O6\n", "biotite: 3.231711 g, composition K(Fe''0.24Mg0.76)3AlSi3O10(OH)2\n", "feldspar: 19.611030 g, composition K0.80Na0.19Ca0.00Al1.00Si3.00O8\n", "leucite ss: 15.555971 g, composition K0.00Na1.00AlSi2O5.06(OH)1.87\n", "rhm oxide: 0.336954 g, composition Mn0.01Fe''0.89Mg0.03Fe'''0.13Ti0.94O3\n", "apatite: 0.300 g, composition Ca5(PO4)3OH\n", "water: 1.704 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 729.507813\n", "liquid: 46.954 g 69.77 0.05 14.13 0.25 1.14 0.14 0.10 1.63 1.35 0.86 0.21 10.37 \n", "Activity of H2O = 0.999976 Melt fraction = 0.470969\n", "clinopyroxene: 11.421340 g, composition cpx Na0.08Ca0.84Fe''0.47Mg0.45Fe'''0.05Ti0.05Al0.18Si1.88O6\n", "biotite: 3.198444 g, composition K(Fe''0.24Mg0.76)3AlSi3O10(OH)2\n", "feldspar: 19.825511 g, composition K0.81Na0.19Ca0.00Al1.00Si3.00O8\n", "leucite ss: 15.892180 g, composition K0.00Na1.00AlSi2O5.06(OH)1.88\n", "rhm oxide: 0.349012 g, composition Mn0.01Fe''0.90Mg0.03Fe'''0.13Ti0.94O3\n", "apatite: 0.303 g, composition Ca5(PO4)3OH\n", "water: 1.752 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 722.507813\n", "liquid: 46.206 g 70.04 0.05 14.13 0.24 1.06 0.14 0.09 1.56 1.27 0.82 0.21 10.39 \n", "Activity of H2O = 0.999973 Melt fraction = 0.463481\n", "clinopyroxene: 11.635463 g, composition cpx Na0.08Ca0.84Fe''0.47Mg0.45Fe'''0.05Ti0.04Al0.17Si1.89O6\n", "biotite: 3.170812 g, composition K(Fe''0.25Mg0.75)3AlSi3O10(OH)2\n", "feldspar: 20.011786 g, composition K0.81Na0.19Ca0.00Al1.00Si3.00O8\n", "leucite ss: 16.208538 g, composition K0.00Na1.00AlSi2O5.06(OH)1.88\n", "rhm oxide: 0.360026 g, composition Mn0.01Fe''0.90Mg0.03Fe'''0.12Ti0.94O3\n", "apatite: 0.306 g, composition Ca5(PO4)3OH\n", "water: 1.796 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 715.507813\n", "liquid: 45.523 g 70.29 0.05 14.14 0.22 0.98 0.14 0.09 1.50 1.19 0.78 0.21 10.40 \n", "Activity of H2O = 0.999972 Melt fraction = 0.45663\n", "clinopyroxene: 11.829868 g, composition cpx Na0.08Ca0.84Fe''0.48Mg0.45Fe'''0.05Ti0.04Al0.17Si1.89O6\n", "biotite: 3.148103 g, composition K(Fe''0.25Mg0.75)3AlSi3O10(OH)2\n", "feldspar: 20.172809 g, composition K0.81Na0.19Ca0.00Al1.00Si3.00O8\n", "leucite ss: 16.507296 g, composition K0.00Na1.00AlSi2O5.06(OH)1.89\n", "rhm oxide: 0.370101 g, composition Mn0.01Fe''0.90Mg0.03Fe'''0.12Ti0.94O3\n", "apatite: 0.307 g, composition Ca5(PO4)3OH\n", "water: 1.835 g, composition H2O\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "...Adding the solid phase spinel to the assemblage.\n", "...Dropping phase spinel from the assemblage.\n", "alphaMELTS at: P 5800.000000, T 708.507813\n", "liquid: 44.896 g 70.54 0.04 14.15 0.21 0.91 0.14 0.08 1.45 1.12 0.74 0.22 10.42 \n", "\n", "\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "sh: line 1: 7599 Killed: 9 /usr/local/bin/doalarm 300 /usr/local/bin/alphamelts < batch.txt\n", "RUN_ALPHAMELTS.COMMAND WARNING: alphamelts may have crashed!\n" ] }, { "data": { "text/plain": [ "Dict{String, Union{Vector{Float64}, Vector{String}}} with 23 entries:\n", " \"Temperature\" => [1016.51, 1009.51, 1002.51, 995.51, 988.51, 981.51, 974.51…\n", " \"water\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 2.31…\n", " \"elements\" => [\"Pressure\", \"Temperature\", \"mass\", \"solids\", \"liquid\", \"a…\n", " \"albite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 3.38…\n", " \"clinopyroxene\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 6.89…\n", " \"mass\" => [99.707, 99.7072, 99.7074, 99.7077, 99.7079, 99.7081, 99.7…\n", " \"rhm\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0,…\n", " \"ilmenite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0,…\n", " \"muscovite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0,…\n", " \"anorthite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.12…\n", " \"biotite\" => [0.0, 0.431308, 0.831234, 1.20432, 1.55234, 1.87696, 2.179…\n", " \"feldspar\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 23.2…\n", " \"quartz\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0,…\n", " \"leucite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 19.3…\n", " \"magnetite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0,…\n", " \"hematite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.0,…\n", " \"garnet\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.26…\n", " \"apatite\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 0.31…\n", " \"Pressure\" => [5800.0, 5800.0, 5800.0, 5800.0, 5800.0, 5800.0, 5800.0, 5…\n", " \"liquid\" => [99.707, 99.2759, 98.8762, 98.5033, 98.1556, 97.8312, 97.5…\n", " \"amphibole\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 9.22…\n", " \"solids\" => [0.0, 0.431308, 0.831234, 1.20432, 1.55234, 1.87696, 2.179…\n", " \"orthoclase\" => [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 … 19.7…" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## --- # # # # # # # # # # # equilibrium batch crystallization # # # # # # # # # # # #\n", "\n", "# Conditions\n", "P_range = (5800,5800) # bar\n", "T_range = (1700,450) # C\n", "\n", "# Elements to include in simulation \n", "elements = [\"SiO2\",\"TiO2\",\"Al2O3\",\"Fe2O3\",\"Cr2O3\", \"FeO\", \"MnO\", \"MgO\", \"CaO\",\"Na2O\", \"K2O\",\"P2O5\", \"H2O\",]\n", "# Starting composition (must correspond directly with elements above)\n", "composition = [61.437, 0.380, 14.864, 3.383, 0, 0, 0.067, 1.606, 3.286, 3.435, 3.494, 0.228, 8.0,]\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=-7, dP=0, index=1, version=\"pMELTS\",mode=\"isobaric\",fo2path=\"FMQ\")\n", "\n", "# Read results\n", "melt = melts_query_liquid(scratchdir, index=1)\n", "solid = melts_query_solid(scratchdir, index=1)\n", "modes = melts_clean_modes(scratchdir, index=1)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Plot results" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdd3hUVfoH8Pfc6ZNJnzSSAAGSEHoRAlJWAbFRpSmooAiIoGBbAVHX37KrIMjiirIWggqCirQVRaUIivSlt5AEAunJJJNkernn98fFIYQBApiZkHw/j4/PzJmbmXfGMd+ce09hnHMCAABoqAR/FwAAAOBPCEIAAGjQEIQAANCgIQgBAKBBQxACAECDhiAEAIAGDUEIAAANGoIQAAAaNAQhAAA0aAhCAABo0OT+LsC7I0eOPP7440lJSf4u5DbGOWeM+bsKaNDwJbxF+ABvUWVlZXh4+PLly699WB0NwoqKCovFMmLECH8XchuzWq0ajcbfVdzG3G63KIoKhcLfhdzG8CW8RXa7XalUIgtv2t69ew8ePHjdw+poEMrl8rCwMAThraisrAwMDPR3Fbcxp9MpiqJKpfJ3IbcxfAlvkdls1mq1CMKbJpfL09PTr3sYrhECAECDhiAEAIAGDUEIAAANGoIQAAAaNAQhAAA0aAhCAABo0BCEAADQoNXReYQA0EAUWWllprjklBggp57RrFcU6xktRGEWPvgQghAAfM3N6bcC/t158cdcfsHEBzYWPuklI6JfC3hauvjUr+4oDesZzXpFs55RrHkQppND7UIQAoDvOEX6xyH3ByfEeB0b2Jh90kvWWc9kfyRdjyhG7QWR07EyvqOA/3CBz9onElGv6Iu52DaUCYhF+LMhCAHAR04a+djt7gg17RksTwi8aqAJjNqFsXZhbGorIqKsSv5rAf+1gC8+IeZbeLdI1lTHIjQUoWZ6NUWqWaSG9GrSq5kSYx7gpiAIAaBGXCL9mMuPl/EQJWnkpJZRiJKpZaSVU5CCVDIKVDCdghRXpJHRQevOiZ+dEU8a+WsdZc+0urFOXbNA1iyQjU0kIiqy0t5ift7ES+yUUcF3F1GBRSyxU7GVl9hJJ6dIDYtQk17N9GqK1ly8oVezaA3p1RShZirZn/V5QP2BIASA69hbzL/MFFdlii2CWPdIlllBVhfZ3FTmEO1usrio3EF2N5lc3OQkN6dgJWnlTCVQkFwpk7kyKnjfRsKzrYWBjYUrY/KGRGpoQGNG5D1JS+1UZOUlNiqx8yIrFVnpnInvK6Zim1hkpWIbldi4SkaRGhb5R1hGaihSzaSMjNBQpJoi1EyD34sNDP6DA4B3J438qyzxy0zOiEY3F3YOlNdk3IqbU4WDLC5uc1Oe0a7SKFNCWKBPNrMKU1GYSqrwqnWWO6hQCksbL7ZRoZUumPlBA5XYxMI/wlLGKELDojSkV5FezSLUFKVlEWqK0bLkYGqsw2YQ9Q2CEAAuk1nBvz7Lv84Si200MoGtuEvWJeIGfvPLGIWqKFTFiCiSiYFXvxboF8FKClaypGC6RlianFRo5cU2KrbxEhsVWSnfwo+WUr5FPGUkg50nB7OWISwlhCUHU3IISw5mapxxvZ0hCAGAiOh0OV9zjq8+K+aZ+bAEYVF3Wc+oBjpEU6cgnYI1DyKvYWly0ulyfsrITxr5N2fplFHMrOQxGtYyhFJCmCcj9Wqf1w03C0EI0KAdNPB158Q153iZg4Y2Ye+mynpFN9D8qyGdgjrrWWf9pc/IzelsJT9p5KeMtKeYLzsjnjJyRnSx1xjCWgazlBBKCMQHW0chCAEaHKdIOwr4hmxxfTZXymhIE/ZRL1lqBH5N3yQZoxZBrEUQG9j4UmOxjU4a+SkjP13Ot+aJp4xUbONtQ1n7cNYxnHUIZ21DMSqnrsB/B4CGwmCnTRfEDef5z7liUjAb3ET4/j6hVQjSr1ZEqCkimvWOvvTxljvocCk/ZOC7i/h/ToknjTxBx9qHsw5/RGMEzqb6CYIQoD7jRIcM/IcLfOMF8Vgp79NIGNCYvdddgcU8fS9YSb2rRKNTpBNGfsjADxv4WxfEQ6VcI2MdwqnDH9GIteV8BkEIUA8Z7PRzjvhTLv8xh+sU9EA8+1snWe9oTCevQxQCtQ9j7cMYJV5sOVfJD5Xywwb6MoP/da943sSDFIowlTtYSUFKClayYAV5bocoKURJQUoWrKRgJQUpKFjpo2kq9Q+CEKCecIj0eyHfkif+lMNPl/O7YoR749irHQR0LG4XTQNZ00A2pMnFu25O+UaLXaYud7AKJ5U7eLmDKpxU7qASG8+sIKODyh1iuYMqHFTupHIHt7goWEkhl9LxspgMVV3eqKRgJQtCdiIIAW5rIqejZXxzLt+cK/5eyFuGsH6xbF5X2Z1R7BbXcAG/kzEKVXKtljF2nVUCPERO5Q4qkyLTQeUOXu4k6XapnWdVktFBFQ6x3EF/ZCo3OSlEVb27KWVniJKFqChYUb3fGays7bfuawhCgNvP8TK+LZ//ks+354vhKtY3lk1oKXx5txCq8ndl4FdCldUMiKjm2Wl08D/S0XODyhw8q1K6/Ud2OqjcwSucFKqq3t2UQjTUE5lXnMityxCEALcBkdMJI9+ez38t4L/kizoFuyuGDW3C3usub6TFmU+4eTeRnZzIaL+su+k5bVvm4GcrqdxZvd9Z7vByRVO6feWJXOm2L/+qQxAC1FFOkQ4a+K8FfEcB/61AjFCzXtHswcZsXld5Yx3CD/yGUbXspJrEZ5n9su5m+R/nb8scPNskBapYNVmNjovr4QUpqGovM0RFIZcaWZCSkoPZLY6CRhAC1CFldtpVxH8vFH8r5P8r4c0CWe8Y9mgL9p+eimhMeIDb2U1kp9FBFY7q2SmdyD1vunTO9v54YWb7W7okjiAE8CeR0+lyvqeI7yriOwv5eRPvGsF6RLEZ7WV3RmFEHzRoIUoKUd5Ydt4cBCGArxXbaE8R31ss7i7ie4t5hJqlRrLUCDappdAujMkx2hPAtxCEALWLE52t5EdL+dFSOlLK/2fgBhvvGsFSI9m0NrLUCGxTAOBnCEKAP1mZnY6W8aOl/EgpP1rKj5XxUBVrG0ptw9hDTdn/dRaSQ7CzK0AdgiAEuCXnTfx0OZ008pNGftrITxi5zU2tQ1nbUNY+jD3aQmgbVtcnUQE0cAhCgBtwrpIfLuUnjHSyjJ808tPlPEjJWgZTyxDWNpQNTxBSQggT+wBuLwhCgKuyuuhIBT9k4EdK+ZFSfriUBypYuzBqG8b6NGJTWwvJwfVwuSmAhgZBCHARJ8qu5MeNdKRU2hyHzpuFlBB3uzDWLowNbSq0D2fhWMMMoN5BEELDdcHMT5TRsTJ+oowfK+MnjTxEyVqFUvswNqgxm9mWWuhEnQbRB1DPIQihociz8ONldLyMnzDyY6X8hJEHyFmrUGoTyrpFsqdaCq1CLjvP6XS6RdF/5QKAryAIoR6yuSm9nJ8u5+nldMrIT5fz00aullObUNYqhHXWs8dbCK1DfbqqLwDUWQhCuO2V2SmrkktdvawKOl7Gsyp5Iy1rFkStQlivaDY2UWgTxrBWJwB4hSCE20mRlTIqeEYFz6zgGRWUXs7Ty7lKRikhLDmEJQWznlGsZQg11WGhMgCoKQQh1EWcKMfMMyvIk3mZFTyjgisFahHMmgeyFkF0Xzx7trWQHIwznABwSxCE4GcOkc6beFYlZZTzzEqeUUEZ5TyrkoepWIsgah7EWgSx4QnUPEhoEYQlWgDgz4cgBB+xu+m8iZ8z0blKnv3HjXMmKrHxWC1LCKQWQax5EOsZdTHztPhuAoBP4JcN/MkcIuWYeVYF5Vl4voWyKnlW5cW7jbQsRkvSMJbe0WxcotAsiBoH4HoeAPgTghBuhtlFeWZeaKVCK8+zUJ6FZ1/s4fFSOzXWsSY6aqpjTXSsVzR7rIXQNJAaaZmANTgBoO5BEIIXLpGKbLzAQvlWKrLyXDMV2Xi+hQqsvNBKeWZORDFaFq2lKA1rpKUoDXswnpoGCk111CgAewwBwO0EQdjgOEQqsfESGxXbqMjKS2xUYuPFNiqyUbF0187L7BShZlEaaqSlSA1rpKXEINYrmmI0QqSGYgNYAL44AFBf4PdZvVLuoCIrN9ipxEa55TIzicU2XmSlEvvF8Cu0crOL9CqK0DC9iqK0TK8ivZq1DqW7NRSpFvRq0qtZpIbQqwOABgJBeHswOqjMLvXeyGDnBhuV2LnBRkU2KrFxg/RvO6llFKFmERoKV1GwTIjWcb2aNQ+iCDXp1YJeTRFqFoZZdwAAVSAI/cbsojI7L7NTqZ3KHH/csPMyx8UbpXYq++NGkJJCVSxMJUUaC1eRXs3ahFKEhvRqQa+icDXp1UxZZfhlZaU1MFDtv/cHAHB7QBBeVYWT3CI5RTK5OBEZ7cSJLC6yu8nFqdJZvdEhktnFOSejg4io0kkukaxubnOTw01mF7lEqnQSJzI6uMVFAqMwFQtVUpiKQlUsVEVhKgpVslYhUotQ9SEZzlQCANSOBhGElU76Xwk/Vc6Ndip38AonVTqp0kkVDm50kEMks5OIqMLJ3ZycIpmcRESBCpILpBBIJ2dEFKIiRqSVk0pGckaBCkZEwUoS2MVGhUChSkZEzYOIiHQKUgikkQlqGSllFCAnGaMgJRFRqJJp5KSW+enjAACAKuphELpEyqjgx8qkf+iQgRdYeLsw1jqUhaspSMEaBVCgggIVFKQQQlWkEEinICIKUjAZu3QXAAAagts+CJ0inangJ8v4SSMdN/ITZfxMBY/VsrZhrHUojUxgczoLicEMpxYBAMCr2ywITU46Xc5PGvlJIz9lpJNGfs7E4wNY61DWMpgGxLO/thNSQhjOOgIAQA3V6SAss9P+En66nP+xyTgZ7Dw5mLUMYa1C2CPNKSVESAy+bKgkAADADanTQZiWLm68ICYFs5bBbFATITmYGuuwfBcAAPyZ6nQQvtBWeKEtunsAAFCLajEI3W73jh078vLymjZtmpqaKpfLiYhz/ssvv+Tn5/fo0aNJkya19+oAAAA1UVtBWFZWdv/991ut1nbt2p05c2bRokWpqalENHLkyNOnT3fu3Pm5555bvnz5fffdV0sFAAAA1ERtBeFLL73UpEmTlStXCsKlc5s7d+789ddf09PTg4KCli1bNnPmTAQhAAD4V61cgRNFceXKlX/961+PHj26b98+h8MhtX/33Xf33ntvUFAQEQ0bNuzw4cMXLlyojQIAAABqqFZ6hHl5eVardcaMGUqlsrS0tLKycsuWLVFRUbm5ufHx8dIxgYGBwcHBVVuqstvtRUVFS5Ys8bTcddddiYmJtVFtfeV2u91ut7+ruI253W5RFPEZ3gp8CW+R9AEyDJa/WaIo1uSwWglCq9VKRHfddderr77KOR86dOjcuXPfffddl8tV9UypXC53Op1en8Fms5nN5v3793takpKSmjZtWhvV1ldOp/NqHy/UhNPpFEWx6jcWbhS+hLdI+gARhDfN5XJxzq97WK0EYUxMDBHdfffdRMQY69Onz4YNG6T24uJi6Rin01lWVtaoUSOvzxAcHJyQkPDJJ5/URnkNhNPpVKuxDdPNk8lkoiiqVNi/8ebhS3iL3G63Wq1GEN40pVJZk0+vVv7a1el0Xbp0ycrKku5mZGTExsYSUe/evbdu3Sr1Vbdu3RobG4tOHgAA+Jf3HuEXX3zRv3//qKioqo2FhYXr16+fOHFiTZ739ddfnzhxYmlpqcFgWL58+bZt24jowQcffOONNx5++OHevXvPnz9/xowZMhlWBQUAAH/y3iN8+eWXMzMzqzVmZWVNmjSphs87YMCANWvW5ObmqtXqffv2tW/fnojkcvn27dtTU1PPnTv34YcfTp48+VZKBwAAuHU3cI3QZDLpdLqaH9+tW7du3bpVawwODn7xxRdr/iQAAAC16rIgPHr06K5du4jIarWuX7/+2LFjnocsFsvKlStTUlJ8XSAAAEBtuiwIN2/e/MILL0i3582bV/UhxlhycvJ//vMf35UGAABQ+y4LwsmTJ48bN46IkpOTv/jii65du3oe0mq1GEcOAAD1z2VBqFarpUk/+/bti4qKwgQgAACo97wPlsEGSQAA0EBcCsKvv/7ac4HwGnJycmqzHgAAAJ+6FIRNmjQZOnSoH0sBAADwvUtBmJqaKu2dCwAA0HBgZX0AAGjQrhqE//3vf//yl7/o9fq4uDipZe7cuf/61798VRgAAIAveA/CL774YvDgwTKZbMCAAZ7G6Ojot99+u4b7HAIAANwWvAQh53zmzJlTp07dunXrE0884Wnv2bNnYWEhRo0CAEB94iUICwsLc3Nzx48fX609OjqaiIqKinxRFwAAgE94CUKlUklENputWnt2djYRBQcH+6AsAAAA3/AShGFhYSkpKR9++CEReTa555zPnTs3Li6uRYsWPi0QAACgNnlfYu3tt98eOnRobm5uq1atbDbb4sWLv/nmm+3bt3/++eeeaAQAAKgHvI8aHTRo0Lp1686fP//ee+8ZDIapU6eeOXNm2bJljz32mI/rAwAAqFVX3aF+4MCBAwcOPHv2bElJSWBgYFJSkiBg9j0AANQ3Vw1CSUJCQkJCgm9KAQAA8L2rBmFxcfHmzZvz8vJcLlfV9ldeeaX2qwIAAPAR70H47bffjh071mw2X/kQghAAAOoT7yvLTJgwoXv37idPnnQ6nfxyvi8RAACg9njpERYXF5eVlc2ZM6dly5a+LwgAAMCXvPQIw8PDw8PDjUaj76sBAADwMS9BKJPJ/v73v7/22muFhYW+LwgAAMCXvA+Wefrpp/fv39+sWbN27drpdLqqD/38888+KQwAAMAXvM+Rnzlz5tKlSwMCAjCJHgAA6jcvPUKr1fruu+8+99xzCxYskMuvM+MeAADgtualw1dRUeF0OseOHYsUBACAes9LEEZGRjZv3vzUqVO+rwYAAMDHvPT5GGNpaWkTJ04MCwvr06ePtE8vAABAveT95OdTTz11/vz5+++/Xy6XBwYGVn2otLTUJ4UBAAD4gvcgfOyxx0wmk49LAQAA8D3vQTh79mwf1wEAAOAXmCYIAAANGoIQAAAaNAQhAAA0aAhCAABo0BCEAADQoCEIAQCgQbvqaqKVlZU7duzIz88XRbFq+8SJE2u/KgAAAB/xHoQbNmx48sknDQbDlQ8hCAEAoD7xHoSTJ09u2rTpf//736SkJGxJCAAA9ZiXICwpKcnLy1u5cmX37t19XxAAAIAveentBQQEqFQqdAQBAKAh8JJ2Go3mqaee+s9//sM5931BAAAAvuT9GmHz5s3/8Y9/3Hnnnf3799dqtVUfeuWVV3xSGAAAgC94D8K5c+caDAaDwbB79+5qDyEIAQCgPvEehKdPn642fRAAAKBe8h6EwcHBPq4DAADALzA0FAAAGrRLQbh8+XKdTrdgwQIiatasme4q/FcqAADAn+/SqdHWrVs///zzXbt2JaJJkyaZTCb/VQUAAOAjl4KwY8eOHTt2lG5jaCgAADQQuEYIAFB3Mcb8XUL9hyAEAIAGDUEIAAANGoIQAAAaNAQhAAA0aAhCAABo0LwvsUZExcXFx44dKyoqUiqVERERHTt2DAgI8GVlAAAAPlA9CA0Gw2efffb5558fOXKk6n6ECoWiR48e48ePHz58uFqt9m2RAAAAteVSEFqt1oULF7799tuCINx///1PPvlkcnJyWFiYy+UqLS09dOjQzp07x48fP2vWrLfeemv06NGY3QIA0ABVVlZ++eWX586d0+l0nTp1uueee+Ry+YIFC/R6/dixY6VjDh8+vGHDBrPZnJKSMnLkSI1GQ0Sc8/3792/btq2ioqJ169YjRoyQy696VtKXLl0jXL169ZdffrlkyZL8/PyVK1c+99xz9957b5cuXbp37/7ggw+++uqr33//fV5e3rRp01588cUzZ874sWgAAPCL0tLSdu3abd26NT4+3u12//3vf8/NzSWixo0bx8TESMfMmzevd+/eFRUVUVFRn3/+efv27S9cuEBE6enpw4cPLywsDAoKmjdv3sCBA6ued/Qj5qmjtLQ0JCREEK4/fMZsNguCICV8Ldm9e/f06dOv3BYYaq6ysjIwMNDfVdzGnE6nKIoqlcrfhdzG8CW8RWazua4Nzli+fPn8+fMPHTpUrf3MmTNKpbJJkyZ79+7t0aPHrl277rjjDiISRXHUqFFms/n77793OByCIEi9wOLi4ujo6JMnTyYlJdVetWvXrl22bNn69euvfdil2AsLC6tJChJRQEBAraYgAADUTWq1Oj8/Pysrq1r7v/71r+XLlxPRmjVr7r33XikFiUgQhJkzZ27atKm4uFipVHrOhbrdbs65Vqv1ZfFXc/3zs7t27dq1axcR9ejRIzU1tfZLAgCAq/oyU9ye76MzioEKequLTFGlizRkyJDVq1cnJiYmJyf/5S9/mThxome3BsmZM2dSUlKqtkh3MzIyIiIiPI3Tp09/5JFH4uLiavcN1Mx1gvCll15atGhRmzZtiOiVV1554YUX5s6d65PCAADAi4RAZnL66LVUMpJdPixSLpevWrXq3Xff/eWXXzZs2JCamrpt27YePXp4DlAoFE7nZfW5XC7OuUwm87TMnDnzzJkzW7dureXya+paQZiXl/fxxx8fPXq0ZcuWRHTq1KmuXbu+8MILUVFRvioPAAAu0z2SdY/086D9Ro0ajR49evTo0YMHD165cmXVIGzZsuWWLVuqHnzgwAGFQuHpJr755psbN27cunVrcHCwT4u+uks93oKCghEjRuTk5HhaDAZDZGRkcnKydDc5OTkyMrKkpMTXNQIAQN2Qk5NjtVql25xzo9EYGhpa9YBHH310//79q1evlu6aTKbZs2c/9thj0rCp+fPnr1y58qefftLr9T6u/BouGywTGhraqlWrefPmSR1bqSM4YsSI5cuXL1++fPjw4YIgeHIRAAAaGmnixPDhw6dMmZKampqfnz9lypSqB7Ro0WLFihUTJkx44IEHxo8f37ZtW51Ot3DhQiI6cuTIyy+/rFarH3vssXvuueeee+7Zv3+/n97HZVi1aRx79uyZMmWKxWJ5//33+/Tpk5GRMWPGjN9//52IevTo8dZbb7Vo0cIHZWH6xK3DyPVbhOkTtw5fwltUB6dPEFFOTs7//ve/srKyxo0b9+rVSxoImpWVpVQqPYNfysrKtm/fXl5e3qpVqy5dukiNZrP51KlTVZ8qMTExKCio9kqt4fSJ6tcIU1NT9+7du3z58lGjRnXr1m3x4sWeHi4AAEBcXNyVoz2bNWtW9W5oaOiQIUOqHRMQENC5c+faLe6meJk4KAjC448/fvr06WbNmrVp0+Zvf/ubw+HwfWUAAAA+UD0IzWbzzp07f/75Z7vdvmjRok2bNq1fv75z5847duzwS30AAAC16rJTo1u2bHnkkUeKi4vlcjlj7PXXX589e/b+/fsXL148ePDgAQMGvPPOO9HR0f6qFQAA6ojc3FzpIlrTpk2lll9++eWLL7749NNPXS7X9u3bpWuEnTt3fuSRRxQKhXQM53zjxo3btm1zuVypqakjR46ULjHa7fZVq1YdO3bM4XC0a9duzJgxvtzm6LIe4dixY2fOnGm1Wu12+08//bRw4cJjx47JZLLnnnvu5MmTjLGUlJTTp0/7rDgAAKibPvroo7/97W+LFy/2tJSXl2dkZBDRwYMHp0+fLpPJEhMTFy5c+NBDD0kHuN3u4cOHT58+PSoqqnnz5gsXLuzVq5fJZCKiioqKn3/+OTY2NjEx8aOPPho4cKBP3wz/Q2VlpVqtltZ/kwwcOHD16tW8iu3bt585c4bXvl27dqWmpvrgheqxiooKf5dwe3M4HDabzd9V3N7wJbxFJpPJ3yV453a7mzRp8uGHH0ZFRTkcDqlx3bp1vXv35pw7HA5RFKVGad+JvLw8zvkHH3wQHR1tMBikh2w2W8eOHV966aVqT37+/Hki8hx2K9asWTNo0KDrHnapR6jT6Zo1azZ9+vRjx45lZmZ++umn27dvr7aIXO/evX0zfQIAAOqszZs3K5XKSZMmRUdH//DDD9UeVSgUng1rTSaTTCbT6XREtHbt2okTJ4aFhUkPqVSq6dOnr1ixotqPb926NSEhwZfrzlx2jXDVqlWPP/74v//9byKKiIhYvHhxtRGxN8pms82ZM6ddu3YjR46UWrKzs995552CgoI+ffo8/fTTNdzvAgAAJJZ9m+1nT/jmtZggCx46icmqT7RLS0sbO3YsY+zxxx9PS0sbNGiQ1x93u91Tp06dOnWqNJ00IyNj3LhxVQ9ITk7Oz8/3TJfs2LHjmTNnNBrNpk2bqq5NWtsue3tt27Y9ePBgcXGxzWaLjY299ZR6880309LS7rvvPikIbTZb7969hw8fPnbs2NmzZ5eUlLz++uu3+BIAAA2KPLoJd/poShtTaZhQPZDKy8s3bNjwz3/+k4geffTR2bNnFxQUXDmOknM+YcIEInr77belFrVabbPZqh5jt9sZY569mX755ReLxbJs2bLBgwefOHGiVufaV+Vl0e2qO2XcioMHD+7YsePRRx/1LE+6evXq0NDQBQsWEFFYWNjQoUNnzJihVCr/lJcDAGgIlPGJyvhEPxawfPlyi8VS9XzhihUrXnzxxarHcM6nTJmSkZHxww8/eMZ/pqSk7Nu378knn/Qctnfv3pSUFM/6TcHBwcHBwTNnzly8ePHvv/9+33331f67IfI6of5P4XK5Jk2a9P7773uinoj27dvXs2dP6Xb37t3LysrOnj1bSwUAAEBtSEtL++STTzwjTVatWpWWllbtmFdeeeXgwYMbN26sukTcpEmTPvvss71790p3MzMz33333WeffZaIpGFBnvbi4mLPrAwfuP7GvDdnzpw599xzT8eOHb/66itPY2FhoWfNbkEQwsLCCgoKvK7iXVJScvLkyT59+nhannnmGZ/9dVA/mM1mz/VquAnSWqPVdlaDG4Iv4S2yWCwqlapqd8Lvjh07duzYsWHDhnlaBg8ePHHiRE+8EdGOHTveeeed5OTku+++W2pJS0tr27Zt//7933777Xvvvbd790/0UiMAACAASURBVO4qlWrHjh1PPvnkpEmTiGjDhg2vv/56+/btHQ7Hb7/99sorr0i7PlQliqLFYrmhaq1WK+fX38S4Vj7fo0ePrl69+splxbVard1u99y1Wq1XW082NDQ0NjZ21qxZnpb27dtL446ghjjn+MRuBRbdvnX4Et6iqtfP6oioqKj//e9/ISEhnha1Wr13796QkJCWLVu2bt2aiDp06FDt97/nPOpzzz03evTo3bt3u1yu9957Lz4+XmofPXp0x44d09PTFQrFhx9+6HXnekEQbvTrpNFoavKnWK18xFu3bs3MzGzUqBERWa1WURSPHj164MCBuLi49PR06RiDwWAymby+WyKSyWRBQUH9+vWrjfIAAODmREREXDmOxHNiTxreEhQUdI3FtfV6/YABA65sT0lJ8Wze62O1co3w6aefzsvLy8zMzMzMnDRp0kMPPbR161YiGj58+I8//pifn09Ey5Yt69mzJxZsAwAA/6qVHqFKpfKcUFKr1SqVSpoa2a5du/Hjx3fq1Ck5OTk9Pf26e0QBAADUtpoGocvlkuY5iKJ4Qy/wxhtvVP2R+fPnP/vss3l5ee3bt9dqtTf0VAAAUBfk5OR89tlnmZmZer2+Z8+eAwcOvMalOFEU16xZs2XLFofD0bVr18cff1yj0fiy2uuq6alRxtigQYOutnzANWg0mmojYpo0adK9e3ekIADA7Wjr1q2tWrXKyMjo1atXo0aN3n333b///e9XO9jtdg8dOvS1115r3br1nXfeuXLlytTU1LKyMl8WfF017RHKZLJ169bVaikAAFDH2Wy2MWPGvPHGG54Z9NOmTcvOziaioqKiLVu25OTkNGnSZMiQIdJJxI8//njPnj2nTp2SBpqOHTv2rrvumjVr1ocffujHd1GN9x7hkSNHpK0xqjKbzQcOHKj9kgAAoI7atm1bRUWFNAtewhiTJr8vXLjwwIEDMpls1apVffr0kS6Kfffdd+PHj/dMt5DL5c8999zXX39dk+l9PuO9R9i/f/81a9bceeedVRuPHDly55131qnqAQAamg1nNh0oOOyb11LKFH9NfVYhU3hazp4927RpU69LY7711lvSjeeff75Dhw579+7t1q1bRkbGiBEjqh6WlJRUWlpqNBpDQ0Nrtfiau4FRow6HA5OLAQD8q01ESpAq0DevpZQp5JdvPaFUKq+2vMuOHTtmz55dVFQUEBCQnZ2dnZ3drVu3gICAasdLd+vUKtOXvcPy8nKDwUBEbrc7Ly8vKyvL81BpaenSpUsbN27s6wIBAKCKZiFNmoU08derd+rU6cKFC7m5ubGxsVXbXS7X4MGD161b95e//IWIevXq5XK5iCglJWX37t2TJ0/2HLlnz57k5OSrLSvmF5cF4dKlS1944QXpdrXOLBEpFIrFixf7qC4AAKh7OnXq1Lt37ylTpqxYsUIKs/Pnzx85cqRHjx4mk0laUCYrK8uzxNrkyZP79u07duxYae3o9PT0d955p65twHdZED744IPSmmcTJkx46aWXqi6Hrdfrk5OTpVXTAACgwfrqq6/GjBnTrFmzzp07V1RUZGdnL1q0KDQ0dODAgb169Wrfvv2xY8c8S2b36NFjyZIlw4cPb9OmjVar3bNnz+TJk6WFtuuOy4IwKSkpKSmJiGw2W//+/aOiovxUFQAA1FERERE//fTTmTNnMjIywsPDW7duLXUNv/322927d5tMph49elRUVHgWyB43btzw4cP37t3rcDg6deoUGRnp1/K98D5YRqlUWq1WH5cCAAC3i8TExMTEy/YHZox1795dul1tyRSdTld1W726xvs8wmnTpiUkJDRq1GjkyJEfffRRbm6uj8sCAADwDe89wqNHj27btm3Lli1btmz55ptvGGNt27bt27dv3759H3zwQR+XCAAAUHu89wgjIiJGjhz5n//8JyMj49y5c5988klMTMzChQu9biIFAAANxG+//TZjxoyq+9EfOHBgxowZ27dvv/YPGgyGf/7zn4888siTTz757bff1qm1Wa41od5isfz6669Sv/DQoUPBwcHSBBEAAGiY9u3b9/77758+fXrt2rVSy9y5c7///vvQ0NBrBMTZs2d79ep15513jhw50mAwvPzyy+vXr//ss89qsn28D3gPwvnz52/cuHHXrl3Sxc+HHnpo8eLFd9xxh1xeK/sXAgDA7aJ379579+4tLi6OiIgoLS397bff7rrrLs+jubm533zzDef80UcfXbt27dixY1Uq1ZQpU3r27Llq1SrpmMGDByclJa1Zs2bYsGH+eQ+X835qdP78+b///vuECRPOnj27devWV199tVu3bkhBAABQKpUjRoz48ssviejLL78cNmyYWq2WHsrPz+/SpUt2djZjbNSoUU8//bTFYnE6nZs3b37uuec8zxAdHT1y5MhvvvnGP2/gCt6zbfbs2Rs3bkxLS1uyZMkdd9zRt2/ffv36de/eHWuNAgD4V94OQ+mJSt+8lkwlJD8WL8irn8AcO3bs+PHjp02btmzZsiVLlrz99ttS++LFi++9996FCxcSUYcOHbZt20ZE2dnZTqezRYsWVZ+hRYsWdWdrP+9BOHXq1KlTpzocjj179kjXCOfNm6dQKHr27Pnjjz/6uEQAAPAIaxWojfZRn4TJ2JUpSEQdO3bknK9YscJms91xxx2e9hMnTtxzzz3S7dTUVOlGYGAgEZnN5qrPYDKZPP1Iv7vW2U6lUtmrVy+9Xh8eHq5Wq7ds2fLTTz/5rDIAALiSWq9U6/2/dcO4ceOefvrpv/3tb1UbAwMDKysv9lYrKiqkG1FRUXq9fufOnQkJCZ4jf//996oJ6l/egzA3N3fz5s1SXzAvL08mk3Xu3HnGjBl9+/b1cX0AAFAHPfroozab7fHHH6/a2L9//7lz506ePDkwMPCDDz7wtE+bNu2NN97o2bOntIXvsmXL9uzZ8/HHH/u45qvxHoSdO3cuLCxs1qzZgAED+vXr169fv7qzgyIAAPhdeHj4K6+8Uq3xkUce+e2331q0aBEeHj5o0CDGmHT+c+bMmQaDoUOHDp06dTIYDCUlJWvXrm3evLk/CveCeZ3VuHbt2m7dusXExPi+IMnu3bunT5++e/dufxVQD1RWVkqn5uHmOJ1OURQxQOxW4Et4i8xmc53at4+IysvLrVZrdHR01cbCwkKVShUSEiLdtVqtnPOTJ08OGjSo6gqdJSUlx48fDwgI6NChg2+mIaxdu3bZsmXr16+/9mHeSxk6dGgtlAQAALe34ODg4ODgao1Vtyp69tlnU1NTLRbLggULXnzxxaqH6fX6urkqy6UgPHHihFarlU7gXtvu3bsTExPDw8NrsS4AALgN9ezZ88iRIzKZ7IMPPrhdhpVcCsLMzMwRI0aMHj36iSee6Nmz55Ur39jt9o0bN3788ce//PJLeno6ghAAAKoZNWrUqFGj/F3FjbkUhAMHDty+ffsrr7zSu3fvyMjI1NTUpKSk8PBwl8tlMBgOHz584MABq9U6atSoU6dOxcfH+7FoAADwi19++WXNmjWeu/3797/GZgx79+5NS0s7f/58XFzc448/3qNHD5/UeMMuu0aYmpr6yy+/HD58+PPPP9+2bdsPP/zgcrmISKPRdOnS5a9//esTTzwRGxvrp1IBAMDPDh48uHnz5qlTp0p3IyIirnbk8uXLn3766VdffXXo0KGHDx++//77FyxYMGHCBF9VegO8DJZp3779ggULpNsGg0GpVGLcFwAASOLi4p555pmqLRUVFStXrjx//nznzp0feughIiotLZ08efInn3zy8MMPE1H//v2Tk5NHjRo1ePDgyMhI/9R9dd4X3fYIDw9HCgIAgIcoipY/cM4NBkPnzp3T09NbtGjx/vvvv/DCC0S0a9cutVo9cuRIz08NHDhQr9d///33/iv8qrChBADA7eTMgSUF57b55rUEmaL7wKWC7LIV3bZv3+6ZZX7gwIFPP/20f//+0nnE4cOHx8XFvfbaaxkZGQkJCYJwqa/FGGvWrFl2drZvKr8hCEIAgNtJbOKDwRGtfPNacqWuWgoS0d1331113emDBw+eP3/es9a2y+XKzMwMCQnxrDXqUVFRodVqa7Xgm4MgBAC4nWiD4rVBdWjcvlarfeKJJ5566ilPS2BgIGMsKysrJycnLi5OaiwqKjp9+nTXrl39VOa1XOcaIQAAwDXcf//9q1atUiqVoaGhoaGhDodDLpd37ty5R48ekydPtlgsRGSz2aZMmdK+ffvevXv7u14vEIQAAHDzxo8f37Nnz8TExEGDBnXv3t0Tdd98840oivHx8b169YqPjzcajWvXrr1yqZa64FqnRo8cOXL8+HGj0Th58mQiysnJ0Wg0WFAGAKDBeuaZZ6qeBSUiQRAWLVr05ptvZmVlBQUFefaU0Ov1GzduzM3NzcvLi46OrsvLsHgPwvLy8pEjR0qXQ2NjY6UgnDNnTkZGxubNm31aIAAA1BkqlcrrliwhISGdOnW6sj02NrbuL8Pi/dTo5MmTjxw5snbt2nXr1nkax4wZs337dpPJ5KvaAAAAap2XILRard9+++2iRYuGDBni2V+KiJKTk10u14ULF3xYHgAAQO3yEoSlpaUOh6Ndu3bV2qV9FKUhQAAA0ABt3rx50qRJR48ele46nc5JkyaVl5df7fhTp07NnDlzwIABDz30UFpamiiKnodOnz797LPP3nPPPUOHDv3888+rPuRjXoIwPDxcpVJ53qfHb7/9JghCTTYsBACAeuno0aOffvrpyy+/LN11uVwfffTRNTpI0lYVzz777KhRo9544w3PQta7du3q0qWLVqudPXv2yJEj/+///m/cuHG1X753XgbLqNXqIUOGzJw5s02bNp6hrocOHXr++ef79++PUaMAAA1Z3759jx8/vnXr1j59+lRtz8rKWrNmzblz52JiYp588klpGbZZs2Z5DjAYDCtWrJBCdMKECVOmTHnrrbekh3r16pWYmDhmzJh7773Xh2/lIu+DZd577z2VStWmTZsxY8aUlJSkpKR06tTJ5XItWbLEx/UBAECdotFoXn311ZkzZ3LOq7Zv2bKFMXb33XfbbLYuXbpcucRaenp648aNiaiwsPD48eOTJk3yPBQXFzdgwICqwzN9yfv0icjIyH379i1dunTz5s05OTlhYWGPPfbY5MmTQ0NDfVwfAABU9Y+Tp9fm5vnmtQLl8k29e6iE6l2mp5566l//+tf69eurduCkvQbtdnufPn3279//448/jhgxwvPozp0709LS9u7dS0SZmZkymazazMKEhIQTJ07U4pu5uqtOqNdqtVOnTvXsvggAAHXB+IQm90VH+ea1VIJwZQoSkVwuf+ONN2bMmNG3b19P408//fT8888LgqBSqbKzs6tm5KFDhx566KEvv/wyOTmZiMLCwtxut8lkCg4O9hxjNBoDAgJq891clfcgPHHihM1mqzY78vjx406ns0OHDj4pDAAAvIhWq6PVan9XQQ8//PD8+fNXrFjhaXniiSdWrVrVq1cvIhoyZIhnFOixY8ceeOCBDz/88MEHH5RamjVrFhAQsG3btiFDhkgtbrd7x44dTz/9tG/fxEXerxGOHDnyhx9+qNa4f//+/v37u1yu2q8KAADqNMbYnDlz/v73v0t3Oefl5eV6vZ6Iqq5Bdvr06fvuu2/+/PnSzvUSpVI5ffr0v/71r5mZmdLPvvrqqxUVFWPHjvX5+yDy2iM0m83Hjx/v169ftfZ+/foVFxefO3euRYsWPqkNAADqrvvvvz85OTkvL4+IGGOzZs266667WrduXVlZ2aNHD+mY+fPnFxQUeC60xcfHHz58mIjefPNNl8vVsWPHhISE/Pz8xo0b//zzz/4ahsKqDfshory8vNjY2KNHj7Zp06Zqe1lZWVhY2N69e7t06VLbZe3evXv69Om7d++u7ReqxyorKwMDA/1dxW3M6XSKouh1WUWoIXwJb5HZbPbXZbObU1RUVFpampSUJHi7sngli8WSlZUVGhpaS+uRrl27dtmyZevXr7/2YV56hHq9XqvV/vrrr9WCcMeOHURUl1cQBwAAP4qMjIyMjKz58VqttlrQ+IWX0FYqlcOGDZs5c+Z///tfT+OOHTumTp169913R0dH+7A8AACA2uV91OjChQsPHTo0aNCgiIiIuLi4/Pz8goKChISEpUuX+rg+AACAWuX9NG54ePiePXuWLFnSq1evkJCQ1NTUhQsXHj58GAuNAgA0ZBs3bpwzZ450m3M+Z86cN954QxTFlStXPvzww3feeeeoUaOqDu/gnH/xxRcDBw7s2rXrmDFjfv/9dz8Vfi1XnVCv0WgmTZpUdQkcAABo4NLT03/99VcicrvdkydPPnDgwKZNmwRB+PbbbwcNGpSUlPTbb7/17dt3//79KSkpRDRt2rT169e/8847LVq02L59e//+/T/77LNhw4b5+31c5qpBCAAA4JXL5Ro/fvzZs2e3bt0qrQ6zevVq6aFu3bqtX79+y5YtKSkp+/bt++CDD44ePSqFYqdOnTQazTPPPDNgwIA6NR7b+6lRk8n05ptvdujQITIyMuxyPq4PAADqFKvVOmTIkOLi4k2bNlVdI83zaEZGhjTdfPv27XfeeaeUgpJx48YZDIa6NjXOe49wzJgx33333QMPPNCvXz9pP14AAKgLbOv2OQ6e9c1rMZkQOGsoyWVVG3ft2iWXy0+dOqXVaq/8kWnTprVt21ZaaDQrK6tJkyZVH1Wr1VFRUdIc/LrDS8jZ7fbvv/9+wYIF06dP931BAABwDap+bRUdmlz/uD/nxRTVUpCI7r777sTExP79+2/btq1Ro0ZVH5o9e/a+ffu2bt0q7WWr1+uzs7OrHsA5Lysr0+l0tV34DfEShBUVFS6X6y9/+YvvqwEAgGtjOrVM589Ftxlj77//vjSzvGoWzpkzZ926ddu2bfOslNa+ffsPPvjAZDJ5km/nzp1Op7Nr167+Kf0qvFwj1Ov1KSkphw4d8n01AABQ90lZ2K9fvz59+kjnORcsWPDFF1/8/PPPERERnsMGDhzYqFGjSZMm2Ww2Irpw4cKUKVMmTpwYFeWjbaRqyEsQMsbS0tLmzp27YcMGu93u+5oAAKCOk7Kwb9++ffr0yc/Pf+2119LT0xs1asQYY4zNnj2biJRK5aZNm8rKyiIiIlq2bNmqVas+ffr861//8nft1XlZdJuIkpOTz58/b7PZBEGoNiiotLTUB2Vh0e1bh/WObxEW3b51+BLeottu0e2rKS8vNxgMsbGxPv4f6uYX3Saixx57zGQy1UJVAADQ4AQHB1850aLu8B6EUq8WAACg3qvRllEAAAD11VUny5eXl//0009nz56tdlHw7bffrv2qAAAAfMR7EO7bt++BBx4oKSlRKBQymUwa+apQKHQ6HYIQAADqE++nRp955pmEhIQLFy6MGTPmxRdftFgsX3/9daNGjbAfIQAA1DNeeoROp/PQoUPff/99XFwcEbndbo1GM2LECI1GM27cuAceeECpVPq8TgAAgFrhpUdoMBhcLpe0UmpQUJDRaJTa7777boPBcOrUKZ8WCAAAUJu8BGFERIRCoSgoKCCixo0b79q1S5p0f+bMGSLC/GIAAJ/xuuYJ/Lm8BKFMJuvevfvmzZuJaNSoUSdOnBg4cODrr78+dOjQFi1aNG/e3OdFAgAA1Bbvg2Xee++9u+66i4ji4uK++OKLs2fPzp8/Py4u7ttvv8X2hAAAUJ94T7X27dt7bo8aNWrUqFG+qgcAAMCnsLIMAAA0aJd6hPv371+9evV1fwAT6gEAoD65FITp6ekfffSR567FYpE2IwwICDCbzYSVZQAAoD66dGp09OjRpX9YtWqVXq9PS0urqKgwmUwWi2X16tVxcXH//ve//VgrAADAn87LYBlRFCdOnPjOO+888sgjUotGoxk2bJherx8wYMCQIUPqx0aRAAAA5HWwTHFxcXZ2dqdOnaq1d+rUyWQynT592ieFAQAA+IKXIAwICJDL5T/99FO19p9//pmIQkNDfVEXAACAT3gJQp1ON3r06JdeemnmzJn79u27cOHCwYMH33rrrXHjxt1zzz0JCQnXfVKLxTJ69OjIyEi1Wt26devvvvvO89DKlStjYmJ0Ol3fvn3z8/P/zLcCAABw47zPI1yyZMkjjzwyb968rl27Nm7cuFOnTrNmzerTp89XX31Vkyd1Op0dO3Y8fPiw1WqdOXPmyJEjpZVL8/LyJkyYsHr1aqPR2Lx582nTpv2ZbwUAAODGeV9ZRqPRLFu27B//+Mfhw4dzc3Ojo6PbtGlTk76gJDg4+OWXX5ZuP/roo88++2x6enp0dPSKFSt69erVo0cPIpo5c2ZycrLRaAwJCflT3gkAAMBNuNbCobGxsbGxsbf4Ajt37hRFsW3btkR05syZNm3aSO0JCQlyuTw7O/tqQeh2u8vKyjx3g4KCZDLZLRYDAABQjfcgPHLkiNPp9PpQ586da/7s+fn5jz766MKFC6UhNkajsVGjRp5Hg4KCSktLvf5gUVHRoUOHmjVr5mmZO3fu6NGja/7SYDabGWP+ruI25nQ6RVG82v8IUBP4Et4ii8UiiiI+w5tmtVprso+V9yDs379/YWGh14dqvjlWcXFxv379nnrqqSeffFJq0ev1FRUVngOMRmNERITXn42MjOzcufPu3btr+FpwJc65TqfzdxW3MSkIsQHnrcCX8BYxxrRaLYLwpmk0mpp8et6D8Ouvv3Y4HJ675eXl27dvX7ly5fz582v48kaj8f777x82bNirr77qaWzZsuW6deuk26dOnRIEoWnTpjV8QgAAgNrgPQh79+5drWXYsGHNmjX74IMPxo4de90nNZvNffv2jYmJGTx48IEDB4ioWbNmoaGho0ePfu2111avXn3XXXfNnj374Ycfxl+LAADgXzewDdPw4cP37t175syZ6x5pNBrDwsLsdvuMPxw+fJiI9Hr9mjVr5s2b165du4CAgIULF9584QAAAH+GG9huvri4mIhqMnYgNjZWWobmSn379t27d2/NXxQAAKBW1WjUqNvtzszMnDNnjl6vT0pK8lVtAAAAte4GRo0mJiZ+/fXXcvkNdCIBAADquBqNGpXJZLGxsc2bN8eUdgAAqGdqOmoUAACgXrrWec5z585Ja41GRUW1adMmOTnZZ2UBAAD4hvcgtFgsEyZMWLlyZdV1ZB544IHPP/88PDzcV7UBAADUOu9BKG2WNGvWrOHDh0dFRRkMho0bN/7zn/8cMWLE1q1bfVwiAEAD4TaWOPPPOvPOOvPOilYzD4ngjRIUUfHyyDhZMDohtcVLEJpMpq+++mrRokVTpkyRWmJiYqRTo0OHDs3Kyqq6FjYAANwc7rA5C7KduVnO/HPOvLPOvCwmVyhiEhSxzdQtOwvaQHPuWWdelvXQdldhDnfa5ZFx8sg4RVRjeWSsPDJeHhHL5Ap/v4n6wEsQms1mt9vdr1+/au1Si9Fo9EVdAAD1jrvc4Mw54yw47yzIdl444zIUyMOjFfGJiugm6uSOivgkWVDYZcc3beNZdFu0mV3FeW5DgbMg23Jwh9tQ4Cw8L6gD5NFN5PoYRVRjRUwTWXi0PCyasEj3DfIShJGRkU2bNt2/f3+10TH79u0LDAxs2bKlr2oDALhdcafDVZTjKs51Fue4Ci+4inKchRdkumBFTFNFTIKmVWpQv4flEbEk1HSdS0EdoIxPpPhEDfW62CS6XaWFrsILzsILzvxz1sO/OYsukNstj4yTR8UrIuPkkXHyyHh5RCMmw/zva/Hy6TDGPv3007Fjx1osloceeig8PLy8vPyHH36YNWvWxx9/rNVqfV8lAEBd5jYWu4pynUU5rqILrqIcV1GOu7JMHh4tj4yXR8apEjvoeg6SRzcW1H/q709BJtc3kusbqVunetpES6WrKMdZeN5VlGvZt8VZeN5dViQLiZBHxknXGuWRcYrIeEEX/GdWcpu7FIRffvnl1KlTPXctFsvEiRMnTpyoUCik5dYUCsXkyZNHjRrlhzIBAOoG7nK6SvJcBeddhnyXocCZn+3My2IymTw8Rh7dWBHdRNWsjTy6sSIyjgQ/rEAiaAOVTVOUTVMuNYluV1mR21DgzM92nE+3/G+7qyCbOx1yfYwsPFoR3UQR3US6wRRK3xdcF1wKwqSkpIkTJ/qxFAAA/xNFt7lcNFdc/MdU7jaXi6YK0VLhLje4inJFc7k8Mk4eEauIiFM1axPQ7T55ZKygDvB33VcnyOThMfLwGFVSR0+baK5wFl6Q+q/mfVtcRRfcxmJZSIQiurGn1yiPjBMCgvxYuM9cCsI77rjjjjvu8GMpAAC1SrRUStkmWircJinnjKK5wm2qEC0VoqlcNFeIVpMQECwEBMkCggTpn8AQeWiEENdcCAyVR8bJQyPrwWgUISBI1ay1qllrTwt3u1wlea7CHFfRBUfmMfPvP7iKc4gJiqh46QSvIipOHhkvD4/2S0+3VtX0Cqrb7d60aVNaWtrq1atrtSAAgBvFnQ7RUilaTaLFJFaUussNotVUvcVcQYIgaHSCVicEhcuCwwSNTtAGKoL1Kq1OFhwuBIUJGp0sMKT+/aKvCSaTK6IaK6IaV20ULSaXId9VkO0sOG/et8VVkO0yFMiCwmXh0YroxoroJnJ9jDymqSww1F9l/ymuH4Tp6elpaWmff/55Xl5eYGCgD2oCAPBCdDuyT9szjrgrSkVzhdtULprLRVO5aKkkxi524wJDBG2QoAsStEHyqHiZLlhqFwKChYBADJ68UYJWp9QmKuMTPS3c5XQV50oDguxnj5t3b3IV5ZJcoYiOl6Y2KqIayyPj5GFRNR8Q63dX/VpUVlZ+8803S5cu3blzpyAIDz744Lvvvjto0CBfFsedDneFgYjI5RIdNiIil5M77UTEnXbudBKR6LCRy0lEOIwfeAAAIABJREFUos1CXCTORav54gFSu6WSiLjdyt1uQakiuYKIBJWGZPKLNwQZE2RMrSEiptYyJpBcISjVRMQ0AYwxkiuZQskYY5oAImJKNZMpmEzGVBpffhoADZOz4Lw9/aDt9EFH5lGZPlqd2EEeESs0TbkYewFBQkBwgx3l4XtMrlDENFXENK3a6K4ovThFpCjHnn7QVZTjriiV6xtdutwYFV+Xr6RWD0LO+a+//pqWlvbNN9+YzeZOnToNGjRo69atGzZs8H1xpl/WmHf/SEQkl0vJRHIFU6iIiCmU0ldfUKovZptaS0wgxgStjoiYIlxac0HQBhIRU6mZTC7abeR2EZFot166Ibq52y1aTETES4s4F8nllHKXW82cc+50kMvBucitFiLiDht3O7nbze3Wi68rCEymYBez82KUSnX+UYzqYrUaHRExpYopVIJay9RaQaVhKg1TaQSNTlBrmFLDlCqffsoAdY+73GBPP2hLP2hPP8TkClViB+0dd4c98jwG/ddNsqAwWVCYKrG9p4U7HVLH0VmUYzu137VjvbPogqDSyKPi5ZFxCmmCo9RxrAMXXC8LwoULF37wwQcZGRkxMTGTJ08eO3ZsmzZtli5d6q/1RQPveTjwnof98tI1J9osJIrc7eRVspNcDu50ENHFfHXaL961moiIO2yiyegqyeNWs2i3coeV222ixSTardxu4S6noAlgaq2g1DC1hik1glZ3KS/VWkGjYyqNoNYwVYCg1jBNgKAOYCoNFluCuoJz0VzuNlW4ivKt3EGcmFLF5AoipohvIf05eCXRZrFnHLGnH7SnH3RXGlUt2quTOwbd+6hcH+Pj8uHWMYVS0ShB0Sih6nkzt7HEVZQjrTBgO7nfWXRBrDReHIIbGSePaixlpO9Ptl0WhIsWLSosLFy8ePGkSZOwB28N/ckzZIlIFEWbRbSZud3K7VbRbuVWs2i3cLtNtFtEi8llKOA2i2i3cptFtFu41SLazNxmISKm1grqAKYJENQat0zpDAhkaq2g1pIgEzQBTDrHq9YymUzQ6KRuq6DSkEwmaAJxshdqiNut7soysdLoNleIJqO7oky6VueuKBNNRtFc7jaVCxqdoAshjU4MDiPGuMPOXU4SXc68s5oOvZWNk0RTubui1GUocBWeFy2VgkbnNpUrm7RUJ3cMHfOyMq5FXegowJ9LFqKXhehVSR08LdxhdxXnOItyXIUXbMd3O7fluopyBK2uaq9RERkvC42o1e/DZUE4cODApUuXTp069auvvho7duzw4cODghrEJJK6RRAErU46p3pDuMvJ7VbRZhYtJm63WowGFSNut4o2C7ldosXEXQ7udHCbmbvdotUsXXCVzhKLVhN3ubjDxpRqJpcztZYxGVNrmCBjKg0Jsot5qVRL8cnkCunstHSbiEgmF66Xo9IJ6mu+dzlTqa/1DLLrHnDxHPVVD1Cq/tyus3T5WbSZiYvcauZc5FYzSb1/fvEaNhcvnkiXzhDU9JltFs7FGh4sXQu/6lM5L56iqO6Pi+7Vn81qIn5ZC3faucshPRUxQRYYIgSGCAEhMl2wEBQqD4sSGicJuhBZYKigC5bpgqWBl5WVldVG2Lkryyy7f3ScOyUEBMnDY1SJHRTRjYWAYNFqkgWF4VJfQ8OUKkVsc0Vs80tNnLuNxc6iHFdRjqvwgvXYbldRjmipkEfEKSLjqpxZjbv2/+k35LLfSv/+97//8Y9/fPXVV2lpaePHj586derQoUN1uhv+jQx+weQKJlcIAUEUTkTkqKzU3vgoX+6wcZfr4q91m4WLIrdbSXSLdiu5pX+7uMPG3S7usIsWE3capEFJ5HaJdut1ntxu427XtY4QXdxuu9YzuK97wMVz1Fc9wHFxFJVX1w1aQRNIRNxp4y6n9HaYSsNksot/OmgCGGNMoyMiQRNAjDGlmskVF/+e+OOacQ2xsEjGarwQpfZa/60919Sr++NidvVn0+jo8r+/mULF5EqSvma3cBlbFhjq9XrHTfzlB/UTY7LQSFloJCV38rRxu9VVnCt1HK1Hf68synUV5wgBwVIiatrdqUrscI2nvK7qf54HBQVNmDBhwoQJJ0+eTEtL++KLLwoKCmQy2fPPPz9mzBjMuK/3mFLNlA33t1LVoHU6naIoqlSX/dIXrZV0cbCx/PodXAD4MzCVRhHXQhHX4lIT566yImkWh3jNP45r4qp/b6akpMybN+/ChQvr168fMGDA4sWLu3Tpgq0noH5jMrl0Xvpq/0hLVQm6EEGrQwoC+A1j8rAodcvOut6DNW263eKTXefEi1wuHzRo0Lp16y5cuPD/7d17lBTVnQfw371169mveQ8DAwzgKxgSFEUxREkEH3uMbtBwTHysyRFH92zWnMPJYvboye7qbiLrru5xMfGRNSbqmhiSHCXBAD6zEFFJjGgUkKc8h2Fm+lmPW/fe/aN6mmaYGYaH0zP273M4c6qrb1fdLrr7W/fWrar77ruPMfzmI4QQ+kQZ6hGI5ubmRYsWvfvuux9rbRBCCKFhNmougYMQQgh9HDAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVwyBECCFU1TAIEUIIVTUMQoQQQlUNgxAhhFBVY5WuwCiW5lyq4nRBhL6U0bRUkOa8VMyXshCGhz0UIsGYQWlK1y2N2pqWZLpBaVLH/w6EEBpuI/uXV0hQCpjWd7ZSGR4CQDcPACDLw1CpghC+EIGUeSFCqbJhqED1BBwAejhXxWISAEKlsryYTHkhgt4Ay3AuVDHZejiPJhWonrJU6wl4b/ZBStcpKU47GjNpsXlNCaR0vfQSk1KHscMealo2DH0pM5y7QnhCpjkPpMyGoa1plkZrdN2gNM5YjDGD0lpdNyiNMS3OmKNpccZSum5rmsO02t6JlK7HNOYwLcFG9v8pQgiNMCP6R/PFx1+YtmabIJAxaI+upXXo1rUuRnp06pokbVDX0KLpjKE8k+YMJaliEColDRKGUpg05CI0ICSgdCIogK7pGiEGEQBACUkxXYNitE2kipLitGUJIAAABMAislQli4oTfVcKIAQ39EIZAgFgAAzABADwhe8KJRQJFJOCckG5p/ky3COIUJogGgcqlBYCC5TmSxCgcTAEUAEsBD0EJhQVRNcU10hIVciAUwg1EEz5FIRBlQFCA0FBGMAdjWlEGMA1IkxCErrGQDAidRBxIwYAlNCY7gCARjWH2QCga7qlmQBgMsPQDACwmcUo04jm6Laj26ZmWsyM6zFTM0xmxo0YATLIxkAIoYob0UG4aWZhScP2JsLinNVJmuJGPKT1whjrk1io2z61MtQJiOEpkxPdlbqnJCPC1oTNhMOI40iHaXFbObq0dWVoytZdU3JLk7YmHV1QleeFirw1i5k61fvMNDVT1/r+j5iaaWh9S2qEOrrTz2I1U9f0DA8LIjyQyXDdKAjhCdHNuStET+D2BNwVwhUizcMM9z0hs0KmQ+EKkQvDQig9KfNCOho1KXE0YlEwKbEomESalOhEmlQYBKgKTerrBJTKmhASEFJ4QhRC6Svh5oMcyIIv/GyQMzTD0sy4ETOZaWpGTHdsZpma4eiOrduWZtq6FdOdKEETRtxipsUsm1lxI2Zrln7Ee0cIoZNrRAfhbZ+55rbPXHNML1FuoPK+Kvgq58m8p/J+8eEBTwW+yqejOdFTRNdIzCJxiyZtErdJwqIJi8RtmrRJwiZxiyZskrCAjLI2TVJnSZ3FOE8kEse3hGwYRtGY4aErRD4M0zx0hSgI0cOLUdodBK6QnhC7ow5eKaJ+5kzIBVHdglNCUnHGCIlpmk7BpMQkoCllhYpwaQZSSmHSUEpO1QGAkEguZV7IQImCFL4X5qTIKxUmGEvqRowZFrMc3Y7pTpSUMd1xdNtmdhSlhx4yK27EYrpDCY4FQwgd3YgOQt/tKmR2AUAYZJQq9k8qpbifObKwCF0puW4ki48ZaPWm1mwxI64bjYZda9p1cHg3nfIClfNVzpVZT2VdlfVkzpU7DvCsq3Keyroy66mcR+IGiVkkbkDcogkLEhaNWxA3SNwiCYskLIiZxOh7IHMgupkghCoZ8iA/SLHytzwU3M8opcrnFAqFsFBsNWrM1JjV5yXMSBByWLUJpbqRAIAEYwnGmkxz6BU4klQqzUOuZC4MPSGjQI0OhYZK9QRcKJXmnCuVC0NXCE+INC8WKIShL2VP8dApz4eCc5XUKfOJRcCkoIGyNElVSCFvkB5QoZI+lZ6SvpS+FDkuCgaBmEbjjCV0q0ZnCUNPMMvRnYQRjzpyHd1xmB2lZumho9sn8q4RQqPOiA7Cj95ftvODXwAAM5Kkd++eEKKbySMLa8ymVOfBoYwUoSdCPwyy3M8EXrfvdh3Dum0AG6AJCBBL1ZkqqQtbL5gs7RjcZqGth5bObT20dW7poaUAuO5yvcB1l2teqLsB87heCHWPMzfQvVB3Q80HAO5nlVKEaroRG2T95W95KHQzSQ5vuQohNE3r3RS+CL0+LwmDrFKHHfJUUvIgW3pINZ0d6oCl5Zu9fHVMjxNa/CAZVqpYmupa72uZEaeUmQAmKQYtAGjMKmWzbiSAUOjdSwAA3UxF5/b0LpDoZqqbhxwgoFYuDH0h05y7QnpSdPf295Zaqz0BLwiRD3lPEOwVYUGE2bxwM0oqZTOwKeggGXENkmOKE+AgA6o8JV0pXCm8OCNxpsU1La6xWtuq1c0600rqTsyIxXXH0R2bmY7u2MwyNN3RHUszDc2ID/ofihAasUZ0EJ5y9sJTzl5Y6VoMifK4yroyU1BZT+U8mXVVxpWllmWnqzKuCkKasEncIgmbxEwiCbENIIQ4JhAgjgmEEKd3DgCJmSSaA4TETCBA7N6SQ5DNZo+7azQiBQ97j6EqJcKyjCxvgIY8p2RxFG7gpXtfG4jQLRYIcjIqoCT3iwW8/P5SNvMgC0pCcS8hmkgDqNIClZJhkIGyqCaURXsSMeYkNINQTTfi0LsDUWwBR7lLgcVilOpUM5Rme8TOa7ZPmE+sgmJ5arhK85SekyQPtCBJQZIuznvCcK/v58IwlxfZUBaE0ghYRBmkYJIsI6EOXFOBBgGRHkiPSleInE2VRVVK1+MaSTAtZUR5aVrMiusxi1lRgsaNWMKIR925UXs0GnmE0KgglRz66Io8L8ihdW75IghEMJSSQkq39+cFACYmWxuc+iHWp18jOghHEWLpxNJpYz9N1UNCIXOeynoq68qCD1IpNwClVN4HAFXwQSnZke6dBlXwlVKqEIBSxTlucT4AkJgJUAzFYojaBqEELJ1oFExdSlGI24RpwDRi6gBAbAMoiZ4tf0gsHQDAMohGgRJiGQBATEaZbmipUvVNu+5j2nrHQUoeBgUAEGFBikDKMAzy0NulHHUGgCqmZpTEUgQySDPREeeuI/yojAg9KbzStBC+DN2Qu1IEupEAohlWimoG0+2QJTwtxvWaPLF9FndZwqWxglZTACsPZg5YRrKsgoxQm7nIhiJdEGFBxTQSo+BQsDRpgTCIx0iGSR+UC6KgRE6EWR5mGAQ1OksZTkyPxY1Y3IjFixPxRPS3/KEewzFEI5lUKs/zAJANcgBQ4AWhpBf6oQyj4eLlf73Q4yL0hMdF6As/ELz0l4tQSklp356hiudQaTz5UAz9aL2pGUPcI9QotdmhQxgXjJt5zRlfGmJ9+oVBOIyYRmtiUHNyOtAOxWdviCo3UFKBx5WQ4PMwX6BEg1BCKFTeAwDZmQGpomcBQLlB9FB5HADAC5SQIJXyAgBQHgchSw1QojPQNSiF7qH41IlGAQ61U6MGbqkkABBLh6iMWZpg0emhxGBE1wAAdEYMdlgZxxhkmBKlem+vaWqgMieokDsoBKfgSxGE3A15XkkeeGklecgLgdcTBp08yIZBlvtZHuTCIBP4mdDP8iALoHQjCWZNYNb7ZoPLUi5NuVqiQGN56mRJKgdNWcUySssLyAjVzUUPD01Ok6DFJcS4cqg0SWhAoEMnlTuJLCiZEzwT8rTPuwxK4kY8UfwX651IJIx4wowlowkjnjQTCSNBR9tor5NIKFHgLgCEUkRtiCiBACAKHigLAC/0ueQAvWc3Abi8OFEIXSGFApUL8gCQ43mllBu6QkpP+FzwQAS+CKKFU0JiegwAotOHHN3WiGYxk1EWDRfv89cyrQat3tD0aIh46a+uMdd1Lcvqc9TjY8ohQzPMau2ZGNFB2PWXrLvfT0xy4q02ZdX7Ze5XsUUYG7Cb1M9mrRPrGgUAiJqkAIqHwAUUA1ipUCo/ik+uhITeSO4tANAbz9AbwACgfA5RYZ9DKAFABaHiAgCAhyoIDytT8EFBsbsYADRabNcaDHQGUStco4f6inWNRPOjJC61g6MYLrV0o3Yw04jJin3OAGAwwjQopTsAAOhmUpPSPK4RQ1IE3M/wKCP9Hh7keG9G8mBf6Gd5WWqGQZYHOcELodXoWU2u0eAadQVWU2DJnBbP0liW1GRoUw+wNKM9oLqIMCitNbRanSUoxJS0A2FxbuR9pg4QsV3JrORpznsyQSbr5xzdTpqJpJlIGomkEU+ayaQZTxrJpBkvzux99jje6UlR4O7Wnh3dXo9QImoMlRo9eV6QUkZZFSrhchfKWjk5no9OA47aXqWgKr380LmwRIuGQTHKbGYBgK4xMzoptrchYjLToDqUnd1k65atWwDQHGvUqEYA4kYcetPFZrZGaXTOkkF1k5mlhZ8s+XzecRxSxbsyw2NEB6ERZz0bcwd+udft8GPjrOQkJzkplmhz9NhQh2iiE0WKUUGgNw8ah7cCvUkMohi9KgiBhxAdlxWy1FcMvJipxSQ+vB18qKUbtYN5qAIB0NsvHaVvb7oXo9TUgRLf1EFnhFFi6qVULuuRBmIZxZBmlJgMNI0YDHRNMxjTY7ZRQ+w2SDHCtFLLeIA3Krif4V468DPcT/Pev4G3N/B7uJcOvO7AS3M/HXhdaZ97sRbXbCmYDXmjLqfV5FgiTWP7Id6jaruFdjCU6VA2mEZDwqzTNdPUPE1RKggJPeJ1+65W6FDyw7yfyQbZtJ9N+5kCLySNZNJMpMxEyowmkjVWKmUkkmYymlNjJpNmotRyOlKpLRXJF/JZkhdSFEJXKlngbjbIRX+7vZ4ur2d//sCe7N5MkJuUmtDg1GlEixpDpQCL6Q7VaYNmGJpRCrOojRU9G4VEwogDAAHS51oQCA3FiA7C+AQ7PsEGABHI7A43sy2/d83BjU99ZKb0RJuTmhxLtDl2Y5W25atFbxIDAIFhOrEhSkeeK8hQ6IpAKIot4Oiw7mE90lGySnkwC0IqPwQhlB9CKIqNXZ+rUERBG/U2E8sAjRDHBEqIbRY7mXWNGAxMnWjUckxL04ndQlgrmIyYOiQoaTaLrVhTJ7oGFgtknodROnYHbrfvdfmFHYHb5XtdgdvluwcLbldPSAtOa8FqyVrNOaMuq6W20USPsruU1SHszqA+xvQxltmcssZa1umWUa+TGBFJFprgWVDwgkyX2729Z2faz0RhmfYzeZ43NOPIy0FEoisNlR5Gh7iiWCKElAYHxY14jZmclJrQ6DSMTYxpchqruf8WVdxwB6FSau3atfv27Tv//PPHjRs3xFdpBq05NVZzagwAlFSFfX5ma757Y3bHiv1KqOQkJzk5lmxzYq0Wofh1Qieq2Ag2KJFSO7GTKY+kvACEUsXRUn7xGC0XKgjB50pIlfdBSJX3JBcQhMrjEArlBioUyg/BC1QoletDECouNMd0dM0xGLFrid5ILB0sneiMWDqpNRSjUgtCLeAqx0U+EFlf7PBFt8+78v6eAt9/MCjknNa83Zqxxub0+s2stofEOpTVKbQ9gaDEbHUmjLWscY7d2mB/2rbG206DaRiUwuHX1LU1zaLFfppa41BGnvjQZYSGwXAH4XXXXff2229Pnz69vb39mWeemTt37rEugVASG2vFxlots+sBwO/mma35zLbC/nXdXlcQH28nJ8eSE+1Em8Ns7EFFI07xUOXAB3ePiSr4KsrLgq+4UD4Hj6sgVD5XbkCCkPiKppUexMG3lJtQXoPyQ+VzVfCVx4EQMDSwqDJAsFDQvKAHA5ILIBeI7qzKZZ1E1kl1O3WdTmqjkVhtxHsI62aaIiAVZISKzvt0pfSlAgIApIcXT041KbE1SoA0W3aLbY237baY02yZLZbVatunxOPlkYlQBQ1rEL7++uurV6/etGlTTU3NY489dscdd7z11lsnuEyzVm+cUdM4owYAQldktxcy2wu7X+nMfuSaKT0xwYmPt+PjLWeMpZl4wS30SUMc84T6QEKhPK7coPjP48oLlBtEMakKfpjLiGxW7CmovA9ehvhpUErjGqjDL9JElGShAgVlVzeSRISUAwDXlKuTgmmkHStjOD2G8b5hPERkj80SDammhtSEMfWfSibqjOJhjjhjpRu76ITGy26oktSZhp2o6GQb1iB87rnnLrvsspqaGgBYsGDBLbfcsmvXrtbW1pO1fGZrtZ9K1H4qAb09qNmdbm5nYf8b3YX9fmncqWbR8h5UzaKDDMqiOqH6gAlK2cDPEhioSUooaOaxtVY165gvnMklz8f6H9HQL2Zrx3qjCM2gRPvYf5WoTql+MtdCKBnKXlEYhlJKZfSz6uOu0mAfmIpgGolrJH7CAx17j56WU6HId6djsZjwPL9zj3dgl9u5O+jeHaQPys48C00GzRpP6L5uBqrHYgElAJDT6cHyhQAIBQAQUPAZCRX0uZRgH5SQPkmZY0QMcMREI4QO+UOvU5I3tH7XrRPaZw2CEk+nAGDS/r8hGiU6IQDgG6y8ehoQo2xElRSCah9vz5ZBKevvB9C3RuIIEkaJQftukPrJLWedfeoJLfZEXnysdu3aNX78+Gg6mUymUqmBgjAIgo6Ojh/+8IfRQ0rpxRdf3NbWdkyrs5p1q1lvPLd4knvoFntshCeVPPR5Fr6Cgb9akivJBzwdVYYDP6sOrbHvMxKEf2y3cxK+OJYrjwIAePkg7BrsJ6OP0BVwDMUBAEQglTjG1xw7yaXkJ3MtSioZHH2B0Q9uvztJx12lwT4wRyLArCH9CB51J4ZoZPDk7rN3WI7qhLJjDu9ogZxzXY/ebyNAI7XOphNolLqhTLvhtkDu8PwdvrtT5A+A5ADAhF5qbirpKih+UzSpxWU/v1dKhVId2uELiRaSw4rVCkOD/keQhpQJoABAiA5ksE1NgHLq1In+NyJXWnTbNgosWg4NSa1PASAg/WQnIVQCC4ECgCMILSuhADgc2trl250A/ThuAOAClf19epxj2IsePiEQ/4jaFrp3fuazk/stL+WQvm7DGoSc8/KrJDDGgqD/CxkUCoV8Pv/mm2+W5rS1tbW0tJzQ6ktrdg7bkCwadz0gcvincdTI5Wg8Hq90LUYxzrk83vMITw4Fwh/S1zgsHG2/SoIIBluU9FX53uFhT3Elw2PcCyst0A/Nss4PGSrVuygGCU18xvQ/nWAAMYCGY13D0KlQ5I5SAlylBvvhVyoU0h2kQG8xX0Hfi/oeScpAKv+oxeDwKwYrCJQc0qs+6fp+UGvGj/H9/rcM53zwLoTIsAZhS0tLZ2dnNM057+7uHjt2bL8la2pqJk2a9KMf/WgYa/dJI4RwHDyV6vhVPggBYIiXIRpBF787DI4aPUF4Qv0JMk1zKFtvWNs6s2fPfvnll6O26iuvvNLS0jJp0qThrABCCCHUx7AG4ZVXXskYu+GGGx5++OH29vZ/+Id/0D7m48AIIYTQ4IY1CBljv//976dNm/buu+/ef//93/zmNwcq6bpuJtPP3XfR0L300ksDHYJFQ7Ft27a//OUvla7FKKaU+t3vflfpWoxur7/++sGDB49eDg0gm812d3cftdhwDwOpra294447HnzwwauuumqQYlu2bPnoo4+GrVafSIsXL96yZUulazGK/epXv3ryyScrXYtRrKurq729vdK1GN0eeOCBNWvWVLoWo9j777+/efPmoxYbleMhERoGQxlshtDHDT+HJ2KIWw+DECGEUFXDIEQIIVTVRuJFdAAgl8t5njdv3rxKV2QU27Vr180334ynEh63HTt2BEHwzjvvVLoioxXnPJvN4rf4RPz5z3/+6KOPli5dWumKjFZbtmxx3aNfCWGEBuHXv/71QqEwc+bMSldkFNu+ffvEiRPxVNzjlk6nwzCsr6+vdEVGsW3btuG5widiz5499fX1Fb6qw2iWyWSGEoQEj8QihBCqZniMECGEUFXDIEQIIVTVMAgRQghVNQxChBBCVU37p3/6pwquftu2bevXr6+trbVtuzQzk8k899xz77777vjx4y3r0L2z//znP69YscLzvNLdfVG5rVu3rlixYvPmzc3NzeXbc82aNS+++CJjrLm5uYLVGxV27Njx/PPPf/DBBzU1Nclk8ZbOu3bt+vWvf7179+62tja8TPxQvPbaa7lcrqmpKXoopVy9evXvf//7RCJRVzdSbxk1AoRhuGbNmldffbWjo2PixInld2/94x//+MILLwRB0O+dzFG/giBYsWLFunXrGhsbj3I7MFU5TU1NNTU1jLGXX365NHP//v1tbW1XXnnl1Vdf3draumvXrmj+gw8+2NLS0t7ePnny5MWLF1emxiPYE088UV9ff9ttt1177bX19fXvvPNONP+2CufWAAAOHklEQVRb3/rWKaec0t7ePmbMmIcffriylRzhHnnkkfr6+q985SvXXnvtjTfeGM38v//7v7q6um984xuzZs2aM2dOdJ9PNIhnn32WMVbagEqpL3/5y9OnT1+4cGF9ff3y5csrWLcR7uyzzz733HNvuumm6dOnf/azn02n09H8++67b9y4ce3t7W1tbXfddVdlKzla+L5/wQUXfO5zn/v6179eV1e3bt26QQpXMgg//PBDKWVzc3N5EN51113z58+Ppq+//vpvf/vbSqlCoVBbW/uHP/xBKbVz507btnfv3l2JKo9cZ5111qOPPhpN33jjjX//93+vlNq+fbtt29HOxGuvvdbY2Oh5XiVrOYJt3LgxFott2LChz/yLL774vvvuU0r5vn/GGWf86le/qkTtRo3Ozs6pU6d+4xvfKN+TGDNmTCaTUUr95Cc/+exnP1vRCo5oH374YTQRhuG0adMefPBBpVQmk0kmk+vXr1dKbdmyxbbtjo6OStZylHjmmWc+/elPB0GglPre9753+eWXD1K4kscIp0yZcuTp3suXL7/mmmui6auvvvr5558HgDVr1jiOc/755wPA+PHjzzrrrBdeeGGYazvC1dXVFQqFaDqfz0enga9YsWLmzJnjxo0DgNmzZ1NK161bV8lajmDLli275JJLmpqaVq9evXPnzmhmPp9/6aWXrr76agAwDONLX/rS8uXLK1rNke7222//zne+U+oUBYDly5dfeumlUcfU/PnzN2zYUNq8qI8pU6ZEE5qmNTY2RrdRe+211xoaGs4++2wAmDx58tSpU1euXFnJWo4Sy5cvv/LKK3VdB4Brrrlm5cqVnPOBCo+4wTK7d++OfrgBYNy4cbt27YpmlveMjxs3bvfu3ZWp30j18MMPP/3001ddddWFF17IGFu0aBEcvt0IIWPHjsXtNpAtW7bs3bv3i1/84mOPPXbOOefcf//9ALB3716lVPkHEjfgIH7zm9/s3bv3uuuuK59Z/iGMxWKpVAq34VG9+uqrf/rTnxYsWAD463e8+kSJEGLv3r0DFR5xl1gTQpQOEWuaJoSIZpa3HRlj0XxU8rOf/YxzvmDBgnQ6/e///u9r166dN2/ekdstDMMKVnIk8zxv69atmzdvTiaTb7/99nnnnXf99ddHG7C0DTVNww04kHQ6vWjRot/+9rd9unnwQ3isPvjgg69+9auPPPJIlH+4AY9PnygBgEG224gLwpaWlo6Ojmh6//79Y8eO7TMTAPbt23fRRRdVpn4jUhAE//zP/7xu3brp06cDgJTy7rvvnjdvXktLy/r160vFStsTHamlpeWss86KRopOnz49Fott3Lhx2rRpSqnOzs4xY8YAwP79+1taWipd0xHqqaeeklLee++9APDGG2/k8/k777zznnvuKf/yBkHQ1dWFH8JBbNmy5ZJLLvne975XOkLU59cPv8VD1CdKojkDFR5xXaNz5sz53e9+F02vXLnyC1/4AgCcd955e/fu/fDDDwEgk8msW7duzpw5FazkSEMpJYRERxQAwPf9aA9ozpw5a9euzeVyAPD+++93dXWde+65lazoCPbFL35x69atUkoA2L9/fzqdbm1tTaVSZ599dnRIRim1atWq6AOJjjRv3rx//dd/nTt37ty5c9va2lpaWj73uc8BwJw5c1588cWoC+fFF18cP348XoZ7IDt37pw3b94//uM//s3f/E1p5gUXXLBt27YdO3YAQFdX1/r167EZMBRz5swpHUxduXLl+eefX35SWV/DMn6nf//yL/9yyy232LZ9xRVX3HLLLXv27FFKbdq0qaamZvHixXfeeWcqlXrvvfeiwosWLZo2bdoDDzwwe/bsBQsWVLDaI9Pf/u3fnnbaaf/1X/91991319TU/O///m80/6qrrrrooovuv//+M8888zvf+U5lKzmSCSFmzZr1la98ZenSpeedd95NN90UzV+2bFlDQ8OSJUtuuOGGU089NZ/PV7aeo8Idd9xRGjUahuGMGTP++q//+j/+4z9aW1vxHJ5BTJs2bcKECbf0evrpp6P5f/d3fzd9+vQHHnhg1qxZN9xwQ2UrOVpkMplJkybddNNN9957b319/fPPPz9I4UrefWLZsmUHDx4sPVywYEFNTQ0AbNmyJfoEXHvttaeddlr0rFLqF7/4xfr1608//fQbbriBsRHXqVtZSqnf/OY3b775pmmal1566YwZM6L5QRD89Kc//fDDD88555z58+fjXZkG4bruE088sWvXrrPOOqt8W61Zs2bFihW1tbXRCUmVreSoEHWNllrP2Wz2xz/+8b59+77whS/MnTu3snUbyZ588snS2G8AOPPMM6NWtZTy5z//+dtvvz116tSvfe1r+Os3RAcOHHjiiSfS6fQVV1xx3nnnDVISb8OEEEKoqo24Y4QIIYTQcMIgRAghVNUwCBFCCFU1DEKEEEJVDYMQIYRQVcMgRAghVNUwCBFCx+P111//5S9/WelaIHQSYBAidHQzZsyYMmXKlClTTjnllM9//vO33Xbbpk2bKl2pw7iu+8gjj2zcuHHY1vjjH/84uslJRdaO0EmEQYjQ0W3fvt1xnFtuuWXhwoXTpk376U9/OnPmzM2bN1e6XodkMpn29vY1a9ZU4doROkEYhAgNyWmnnbZ48eLFixc/9NBDy5YtS6fTjz76aHmBAwcOdHV1DfRyIcS+fft83y/NKRQK+/bt6/fWMPl8ft++fYPca+zgwYPlizoOHR0d5UuQUh44cGCgwt3d3eU3QEDoEwaDEKFjNnfuXE3TSj2Bjz766MSJE5uamurr688888yXXnqpVPLSSy+97rrrfvCDH4wZM6alpeUHP/gBALzzzjsXX3xxMplsaWmxLGvu3Lmle2evX7/+wgsvjJ5qaGj4t3/7t9JFEJcuXVpXV/fHP/7xnHPOaWhoiMfjl19+eXd3NwBs2rTp9NNPB4BvfvObdXV1dXV1TzzxxJHVvuuuu0499dRXX331lFNOaW5uTqVS3/3udwHgoYceampqampqam5uXr58eflLfv3rX3/qU5+qq6trbm6eOHHiz3/+8yMXu3HjxqGsHaGRaxiuAo7QaFdXVzd//vzSw+3btwPAwoULlVL/+Z//SSldvHjxW2+9tW7dui9/+cuWZW3YsCEqOWvWrKampjPOOONnP/vZ2rVrN2zYsHHjxmQyedpppz377LPvvffe6tWrb731Vs/zlFLvvfdePB6fM2fO6tWrN2zYcO+99+q6/v3vfz9a1H333QcAp59++tKlS996662lS5cahvGtb31LKZXL5aKIWrRo0apVq1atWrVr164j38WiRYts2z7jjDOeeuqpN95449ZbbwWAm2++eebMmS+88MKaNWsuuuiiZDLZ3d0dlV+2bBkh5KabbvrDH/6wfv36W2+9lRCycuXK6Nn29va2trahrx2hEQuDEKGjq6uru+KKK7q6ujo6OtauXTt79mwAWLZsWSaTSSQSt912W6mk7/uTJk2KMlIpNWvWLF3Xt27dWirw1a9+NRaL7d69+8i1XHPNNRMnTiy/09Ptt99eX18vpVS9Qfj444+Xnr3++usnTpwYTe/btw8AfvSjHw3yLqKxLatWrSpVtba2NpFIdHR0RHP+9Kc/Re9LKSWlnDRp0mWXXVa+hNmzZ19yySXRdCkIh7h2hEYsvJ0HQkOyfPny0j2YUqnUkiVL5s+fv3Llymw2O378+NWrV5dKtrW1vfvuu6WH06dPL78V7erVqy+77LJ+bzK+cuXKCy64YO3ataU5NTU1Bw8e3Lt3b6n85ZdfXnp26tSpzzzzjJSS0qEe4zAMo3R3JMMwJk2aFI/HGxsbozlRD+dHH30EAJs3b962bdsVV1xR/tbGjx//yiuvDHFdCI0WGIQIDcns2bO/+93vMsYmTJgwYcKE6J5w+/fvB4Dvf//7mqaVFx4/fnxpurm5uTStlOrs7Cx/tsR13Uwm8+qrr65bt658fm1tbUdHRykIa2trS08ZhhGG4TEFYTKZLK+qYRh9FggA0SCa6K09/vjjTz75ZJ+FKKXwxpbokwSDEKEhaWpqOvKmsqlUCgCWLVs2yP1my1OKEFJTUxN1JPZhWZZpmjfeeOMPf/jDk1TlExK9tfvvv//mm2+udF0Q+njhqFGEjt+sWbMMw3j22WeH/pILL7xw1apVPT09feYTQj7/+c//9re/zefzx1GTeDwOAJ7nHcdr+zV16tTGxsYhvrWTvnaEhhMGIULHr7GxcdGiRY8++uidd965bds213U3bdq0dOnS//mf/xnoJXfddVc2m73yyivfeuutQqGwdevWJUuWBEEAAHffffe+ffuuuuqqN99803Xd3bt3P/fcc7fffvtQahKLxdra2p566qmXX355/fr1nZ2dJ/jWGGN33333ypUrFy5cuHHjRtd1t27d+pOf/GTJkiUDrf3pp58+WWtHaDhhECJ0Qu6555577rnnv//7vydPnuw4zumnn75kyRLHcQYqP2PGjOXLl+/evfvcc8+NxWJTpkwppeb555+/YsWKPXv2zJw503Gc1tbWr33ta0NvZj322GOFQuGv/uqvzjnnnJNyFdD29vZHHnnk+eefP+OMMxzHmTJlyre//W3TNAdaez6fP4lrR2jYENV7ui5C6Lhxzt9//33P88aNGzdu3LijlldKbdq0KZ1OjxkzZsKECX2e3bp1a2dnZ01NTVtbWzSApYKklB988EEul2tubm5tbS2NtYnGnQ99nA5CIxYGIUIIoaqGe3MIIYSqGgYhQgihqoZBiBBCqKphECKEEKpqGIQIIYSqGgYhQgihqvb/IeFoEp6AZZIAAAAASUVORK5CYII=", "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", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n" ], "text/html": [ "\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", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "## --- Plot melt composition\n", "h = plot(xlabel=\"Percent melt\",ylabel=\"Abudance (wt. %) in melt\")\n", "for e in [\"SiO2\",\"Al2O3\",\"CaO\",\"MgO\",\"FeO\",\"Na2O\",\"K2O\"]\n", " plot!(h,melt[\"mass\"],melt[e],label=e)\n", "end\n", "plot!(h,xlims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:right,xflip=true)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Zircon saturation calculations\n", "#### Trace element partitioning" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Dict{String, Union{Float64, String, Vector{String}}} with 2 entries:\n", " \"Zr\" => 170.68\n", " \"elements\" => [\"Zr\"]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Enter trace element data\n", "# Can include any elements you want, the only required one is zirconium (Zr)\n", "trace = elementify([170.68,],[\"Zr\",],skipstart=0,importas=:Dict)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "scrolled": false }, "outputs": [], "source": [ "# Since the mineral/melt partition coefficients we're using are averaged from GERM \n", "# as a function of silica (from 40-80), first figure out which coefficients we want\n", "si_index = round.(Int, melt[\"SiO2\"] .|> x -> isnan(x) ? 80 : x) .- 39\n", "si_index[si_index.<1] .= 1\n", "si_index[si_index.>40] .= 40\n", "\n", "# Melt fraction\n", "F = modes[\"liquid\"] ./ (modes[\"liquid\"] + modes[\"solids\"])\n", "\n", "# Use GERM partition coeffs to calculate bulk solid/melt d for each element\n", "germ_kd[\"Whitlockite\"] = germ_kd[\"Apatite\"]\n", "germ_kd[\"minerals\"] = unique([germ_kd[\"minerals\"]; \"Whitlockite\"])\n", "d = Dict()\n", "for e in trace[\"elements\"]\n", " # Calculate bulk partition coeff.\n", " d[e] = zeros(size(modes[\"solids\"]))\n", " for m in germ_kd[\"minerals\"]\n", " for k in filter(x -> containsi(x, m), keys(modes))\n", " # Note that minerals that we don't have data for end up being\n", " # treated like all elements are incompatible in them.\n", " # Note, geometric mean = log average\n", " d[e] = nansum([d[e] modes[k]./modes[\"solids\"] .* (10.0.^germ_kd[m][e][si_index]) ], dim=2)\n", " end\n", " end\n", "end\n", "\n", "# Calculate trace elements in melt as a function of melt fraction (equilibrium)\n", "calculated = Dict()\n", "calculated[\"elements\"] = trace[\"elements\"]\n", "for e in calculated[\"elements\"]\n", " calculated[e] = trace[e] ./ (d[e].*(1.0.-F) + F)\n", "end" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Calculate mass of zircon saturated, using Boehnke et al zircon saturation equations" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "scrolled": false }, "outputs": [], "source": [ "# Some quantities we'll need to calculate M\n", "melt[\"FeOT\"] = melt[\"FeO\"] + melt[\"Fe2O3\"]/1.11\n", "melt[\"MnO\"] = zeros(size(melt[\"FeOT\"]))\n", "\n", "# The `tzircZr` function calculates how much zirconium is required to saturate zircon \n", "# for a given temperature and bulk composition\n", "ZrSat = tzircZr(([\"SiO2\", \"TiO2\", \"Al2O3\", \"FeOT\", \"MnO\", \"MgO\", \"CaO\", \"Na2O\", \"K2O\", \"P2O5\", \"Temperature\"] .|> x -> melt[x])...)\n", "modes[\"zircon\"] = modes[\"liquid\"] .* max.(calculated[\"Zr\"] - ZrSat, 0)*2.009/1_000_000\n", "modes[\"elements\"] = unique([modes[\"elements\"]; \"zircon\"])\n", "modes[\"solids\"] .+= modes[\"zircon\"];\n", "modes[\"liquid\"] .-= modes[\"zircon\"];" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Plot results" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAIAAAD9V4nPAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nOzdZ1gU19sG8DMzu0vvvVsARelW7BV7MBY0iYVoYow9GmNssZeYpvkbNTEaNdFYo8YQBUUFu6g0FRuCShMQkLbA7s68HzYvIq4IyNa5f5cf2OHs7MO6Z585Zc6hOI4jAAAAfEWrOwAAAAB1QiIEAABeQyIEAABeQyIEAABeQyIEAABeQyIEAABeQyIEAABeQyIEAABeQyIEAABeQyIEAABeE6g7gEY2d+7cGzduWFlZqTsQLcNxHEVR6o5C++B9awC8aQ2D960BiouLrays/vjjj9qL6VoizM3NdXV1HThwoLoD0TJisdjAwEDdUWif8vJyfX19dUehZSorKwUCAU2jO6p+UEkb4OrVq3FxcW8spmuJ0MTExN3dfeTIkeoORMsUFxebmJioOwrtU1JSYmxsrO4otIxYLBaJRAzDqDsQLYNK2gACgeDevXtvLIaLMgAA4DUkQgAA4DUkQgAA4DUkQgAA4DUkQgAA4DUkQgAA4DUkQgAA4DVdu48QAHRAJUv+eMD+7xYrYUkHGyrIjupgQ7WyoBisrAJKgEQIAJrlbBY35YLM2Yh824Gx0iOXcrjz2dx3SWxmKdfWhgqypTrYUh1taRss6QONBIkQAN7WnynsukQ2wIpqY00FWlP+lpRBg75aHpVwX1xlr+ZyP3Skh7r9N3Djb0V96kUIIQUV5HIOdyWX/ek2O+6szEqf6mhLdbChOtpS/laUEOM80FBIhADwViLSuc8uy37tKsgo427kcb/dY28Xcu6mVBvr//75WVKGr/+mkbLkeh73rIKcz2a33mWnt2Z2dGMU5lELPTLAhRrgwhBCOELuFHJXcrjLOdyvd9lsMTfFi5nWmrbSU9rfCboLiRAAGi42lxsXLT3cR9DJjiLkvxG8SpYk5XPX87jredyOe+ztQq65yYu82MKQiESEEFIkIYdS2VXxrJGAOBkRd1Mq/l2Bk1GdhgEpQrzMKS9zKsyTEELuPee+TWI990vGuNOfedNNTDCWCPWARAgADZRWzA09KdvWVZ4FXxDRRJ7z5A8rWXLz//Pirvvs7UJBU2PWSp+Lf8Z1tae2d2O62b9t3vI0o37pwixvw6y/KWt3VNrPmZ7rS/tZIh1CnSARAkBDFFaSQRGyBf70YNc35BsRTQKtqUBr6mNCCCFFpeJ7paK8SrqrPWXUqN9A9gZkbTtmgT/zczI7KELmY0Hm+TE9HJAO4Q2QCAGg3iQsGRklDXampraq9xwVIU0CrAijtDshTIVkri8905v+4wH76QWZqZB84Uu/24SmkRDhNTDRCgDq7bPLMn2GfNtBc/cUFNFkgid9a7hgvj/9XRLb8qD0lztshUzdYYFGQiIEgPrZeZ89lcH90UOg+be30xQZ6kZffEewvStz7DHrtley9IYsv0LdYYGGQSIEgHq4lMPNvSL7qy9jJlJ3KPXRxZ46Fiw4N1hQUEFaHJDMvCR7UsqpOyjQFEiEAFBX2WIyMkq2s7uglbnGNwYV8TCjNgQxN94VMBTx/0s6MUaWXIh0CEiEAFA3Mo6MPSv9uAU9wEUrs2AVFyPq+45M6mihryXV97hsSKT0wlOkQ15DIgSAOll0TcZxZFGAjnxpmArJTG86JVQwsik9MUbW5Zj0QCrLIiHyko58pgFAqcKfcHtSuD97acEEmXrRY8g4D/r2CME8P/q7JLbFAemGm2w5JpfyDBIhALzB3efcxBjpnz0ZXd3wgabIEFf68juCzZ2Z8Ces537pDzfZEom6wwJVQSIEgNrkV5AhkbLlbZga66jppD5OVOQAwbF+TFwe57ZXMvOSLLMMvaW6D4kQAF5LvoLM8CbUpJY8+q7ws6R29WCuDxUQQrwPScedld19jnSoy3j04QaA+pp8XmYmola11dwVZJSniQm1IYi5O1LYxIR0PSYddVp2PQ/pUDchEQKAYrvus7F53O89GD6v0mmjT5a3YR6OEgbZUsNOyfr8Kz2ZgXSoa5AIAUCBzDJu7lXZb92Yxt0gQksZC8ksb/pBqOATL3pBrMz/L+mu+6yUVXdY0EiQCAGgJo6Qj87JprViqvYUBEKIkCYjm9JXhwpWtWO23WVbHpQee4xkqAuQCAGgpm132Rwxme+H7wcFKEIGuVDRgwW/dGHmXWVDTsrSitFZqt3wQQeAl6SXcouuyX7rxgjw9VCrXo5UwjBBLweqw9/SpTdk2ONJe+GTDgAvmRgj+8yb8bFEp+ibCWky05uODREkPCN+f0lPYR6NdsI4OAC88McDNltMZvvgErkeXI2pw32ZqExu2kWZuyn5qRPjaozLCG2CjzsA/Ce/gsy9Ivu5CyPEF0P99XakbgwVtLGmAg9Ll96QVWIajfbA5x0A/vP5Fdl7zemOtmjNNJCBgCwNZK6ECK7mcu2OSM9no6dUO6BrFAAIISQ6izudyd0cju+Et9XclPq3n+DYY3ZstKyrHfVtB8bWQN0xQa3QIgQAUiEjn16Q/RhEGwvVHYquGOJK3xouaGZKfP+SbLjJytA41GBIhABAVsTJvC2od9zwhdCYDAVkaSATNVDwVxob9LcUS5VqLHzuAfguMZ/75Q77Q0d8GyhFawsqerBgri89JFI67qzsWYW6A4JXKHE8oLi4ODY21tzcPDAwUGEBlmVjY2PFYnHHjh319f/b8VMikSQkJJSUlHh5ednZ2ckPlpWVZWdnVz3R3t7e0NBQeZED8IeUJRNiZN92YJyMMEdGiUY2pfs60Uuuy1odlKxow3zcks8rmWscZSXCmzdv9u3b18/PLzU11cfHZ//+/TT90vVmRUVF//798/PzraysUlNTY2JiXFxc7t+/HxQU5OLiYmNjc+XKlWXLls2aNYsQEhUVFRoa6ujoKH/utm3bevTooaTIAXjl2yTWWp+M80BzUOnMRWRDEDPOg55yQbb/IftnL4GNvrpjAkKI8rpGFy9ePHHixBMnTly/fv3GjRuRkZE1CuzZs6eoqOjatWunT5/u06fPmjVrCCEWFhYXLlyIi4uLjIz8559/Pv/888LCQnn5Nm3apPw/ZEGARnH3Ofd9kuznLnzcblBd2lhTl94RBNlR7Y9K455h1FAjKCURVlZW/vPPP2PHjiWEGBsbDx069K+//qpR5tChQ6NGjRIKhYSQMWPGHDp0iBBibW3dokULeQFvb2+ZTFZcXCx/KJVKExMT09PTlREwAA+xHPn4nGxpIOOGZVBUi6bIijbM/4KY/ieku+7jxnv1U0rXaFZWllQqdXNzkz90c3OLiIioUSY9Pb16gdzc3PLy8qqRQkLIhg0bOnfu7OLiIn94//79CRMmpKWltWzZ8sCBAw4ODgpfuqioKC4ubv/+/fKHQqFwyJAhNXpl4VUsy7IsKmS9afX79lMyRwiZ1IKo+E+Qv2kUxffsO9CZRPajh0WxtwrYlW1o5k3vh1Z/2NSlju+YUhJheXk5IUQg+O/kenp6YrH41TLy5iAhRCQScRxXUVFRlQgPHjz4888/x8TEyB/27NkzJyeHYRixWDxq1KjPP/989+7dCl86Nzc3IyOjqh0pEAg6d+5sbGzcqH+fDqqoqKj674C6Ky8vr/qca5dnFdTyG8ypYFllhap758rLy1mWZRj0xxIPQxLdjxp7jh4cIdvRmTUT1vZ/gUraABKJhOPe/AlXSh2WN9eePXsmn/aZm5v7agPO3t7+2bNn8p/z8vIMDQ3NzMzkD48dOzZt2rQTJ054eHjIj1RlMgMDg0mTJk2fPv11L928efN+/frNnDmzUf8g3SeTyTARtwFYltXS921+guw9d+JnJ1L9S1MUJRKJkAjlDA3JqcFk4TVZtwjuSF+mtcVrG4aopA2gp6dXl74HpfQZmpqatmrVKjo6Wv4wOjo6KCioRpmgoCCFBSIjIz/66KO///7b399f4ckfPHhga2urjLABeCK1mNv9gF3oj1SkERiKrG3HLG9D9/pXejgNnZ9qoKxenblz586ZM0cikdy8efPWrVsHDx4khNy9e9fX1zcnJ8fMzGzy5Mn+/v5r1qxxcHBYtmzZ77//Li8QEhIycODAM2fOnDlzhhDywQcfODs7L1u2TE9Pz9XVNTk5ef369du3b1dS2AB88GUsO9uHscMCmJrkveZ0S3Nq2CnZlVxudVsGtxmqkrJmkYSFhW3YsOHkyZMSieTixYvm5uaEECsrqwULFsgHAt3c3M6fP5+VlXXu3Lk///xzwIABhBADA4OlS5e2b9++xtl69OiRk5Nz/PjxioqKM2fOjBw5UklhA+i8q7ncxafcLG/MINM4AVbU1RDBlRwu5KS0SKLuaPiEqstAohaZPn26u7s7xgjrq7i42MTERN1RaJ+SkhKtm4rVK1w61oP+0FNtiVAsFmOMsBaVLJl6QXYumzsazLQwe9EwRCVtgMOHD+/YsePo0aO1F8NVIQCPHH3E5lVgHRmNJqLJ1q7MbB+6xz/Skxk61VDRWKgPAHwh48jCa+y69swbb1kDtZvUkj7YR/BhjOzrBEyfUTokQgC+2H6XtTcg/Z2RBrVDZzvq4hBmfyo7PlpWLlN3NDoNiRCAFypkZGU8u6otRua0iasxdW6woFxGJsYgEyoREiEAL/x6l/W1JB1s0RzUMoYCsr0bE5XJJj/H/52yIBEC6L5yGVmbwH4VgOagVjISkBmtmR+StXIxP62ARAig+7Yks+1sqHY2aFJoq2mt6ZPZ9P3nmESqFEiEADpOLCXfJrFLAlHZtZipkEzykH2diBmkSoG6AaDjNt5mO9tRfpZoDmq3qZ6yvx+xacVoFDY+JEIAXVYqJd8nyRb5o6ZrPVMh91ELeh0ahUqA6gGgyzbcZHs70T5oDuqEOb7M/odsRikahY0MiRBAZxVJyPqbssUBqOY6wkqPhHnS3yShUdjIUEMAdNb6m+xAF7r6ws2g7T73ZX6/zz4VqzsO3YJECKCbSqVk023ZAowO6hZ7A/KBO/1dEhaaaUyoJAC66Zc7bDd72hPNQZ3zhS+97S6bW67uOHTIa5cqePLkSUJCQl5enkAgsLOza9u2rYWFhSojA4AGk7Bk/U32rz5YSkYHORtRI5vSP96SrWiD/9/GUTMRZmVlbdu2bdeuXffv369+nKbptm3bTpw48f3339e6nUgB+GbXfdbLnLSxRnNQNy3wp9sckc72Ziz01B2KTnjRNVpSUrJo0aLmzZv/+OOPXbt2/fXXX2NiYm7evBkfH3/q1KnvvvvOxcVl1qxZ7u7uW7dulcnQQw2goViOfJfEzvdDc0FnuRpTQ1zpjbcxfbRxvGgRHj58+MSJE3v27Bk0aJBQKKxRrnfv3rNmzSoqKtq9e/eKFSu6d+/u6emp2lABoE4OpbGmItLdAc1BXbbQn+50TDrLmzap+W0N9faiRRgSEhIbGzt06NBXs2AVU1PTTz/99P79+66urioJDwDqbV0Cu8AP8+B0XHNTqq8TvQmNwsbworaYmppSVJ0uIfX09PT19ZUWEgA0XGQGVyolg12RCHXfQn/6+5uyEom649B+b97gKioq6tKlS4SQTp069erVS/khAUDDrYmXLfCnaXSL8oCXOdXNnt56l/3MG9c9b+UNb9+UKVMGDRp0/Pjx48ePDxw4cMaMGaoJCwAa4Gou97CYjGqGr0W+WBpIf5MoE0vVHYeWe1FhOK7mQq7p6el79uxJTk6+cOHChQsXkpOTd+7cmZWVpdoIAaCuVsWz83xpIfIgb7S2oNrZ0NvvYaTwrbyoMU+fPh00aFBKSkrVkcLCQmtrazc3N/lDNzc3KyurgoICVccIAHVwu5C7msN+6Ik0yC+LA+h1iWwlUuFbeFFnLC0tPTw8/P39ly5dKhaLCSEtW7YUiUTvvPPOtm3btm3bNmTIEAMDA9w1AaCZvklkZ3gzBm8e9wed0taaamVOdqJR+BZeJEKRSLR+/frz58+fOnXK29s7PDxcIBAcP37c0tJyzZo1a9assbKy+vfffwUC1DMAjZNRyh17xH7qheYgHy0JZNYmsFKkwoaqmdX8/PzOnTv3+++/T5w4sV27dj/++OOuXbvUEhkA1N0PN9kwT9pcpO44QB062lJNTMjuFHa8B66EGkLBu0ZR1Lhx4+7cudOsWTNfX9+lS5eWl2OdcwDNVSQhO++zM1rjS5C/FgcwK+NYGfaub5CaNaegoODkyZN///13YWHhhg0boqOjIyIifH19IyIi1BIfALzRT7fZQS60qzFuHuSvHg6UgyE58BDdow3xUiI8fvx48+bNBw8e/P7777u7uy9atCgwMPDChQtffPHFmDFjhg8f/vjxY3UFCgAKVcjIxlvs575oDvLdogBmRRzLolFYfy9VngkTJqxZs6asrKykpOTixYs///xzQkICTdMfffTR3bt3bWxsvL29k5OT1RUrALxq53020Jp4W6A5yHfBTpS5HjnyCI3CentpG6bCwsKPP/6YYRhCSPv27YOCgh48eCD/raWl5ZYtW6KiovT0sP8VgKZgObL+JjvXFzsuASGEfOlHr4hDm7DeXiRCY2PjFi1aTJo0KTY29tatW//73//Onj3bpk2b6qXbtWvXrFkzlQcJAIodecSaCEk3ezQHgRBChrjSDEXCHyMV1s9Lt0/s3bt34sSJ7du3J4S4uLhs3769SZMm6okLAOrg20T2S+y4BNXM96OX3pANchXg4qjuXkqELVu2vHDhQmlpaVlZmY2NjbpiAoC6iM7inlWQEDckQnhhWFN6WRx7MoMLdkIqrCsFVcjIyAhZEEDzfZMo+9wHOy7BSyh5o/C6TN2BaBNcSwJopeRC7sYzbixWEoFXhDaj8yvI2SyMFNYVahGAVtpwk53ixehjuii8gqHIfH96RRwahXWFRAigfUok5EAqO6EFekVBsTHu9JNScj4bjcI6QSIE0D57UtgeDrSjIRIhKMZQZK4PvSoejcI6QSIE0D6/3GEntUTlhdqEedJ3npPYXDQK36wedam8vBzbUACoXfwzLrec9MHkeKiVkCaf+9Cr4rHi2pvVNRFKJBIDAwMDAwOlRgMAb7Q5mf2kJc0gD8KbfNSCvpbH3chDo/AN6rrdPE3T7733nlJDAYA3KpGQg6ls0vC61lzgMz2GzPGh1yaw+3tjenFt6lqdGIbZs2ePUkMBgDfajWkyUB+ftKS/TpDcLKCxP0ktFHeNJiYmlpSU1DhYWlp6/fp15YcEAK+1FdNkoD4MBWSWN7M2ASOFtVFco4KDgxMTE2scTExMbNu2bR3PW1pa+tVXX73zzjtz587Nz89XWGb37t3Dhg0LCwu7ceOG/EhWVta33347evTo995779dff5XJXsz9/fvvv0eOHPnBBx+cO3eujjEA6JhreVx+BemLaTJQH9Na0acy2HvPMVL4WvW4tKyoqKj7ZoQTJ068du3atGnTMjMz33333VcL7Nmz54svvhg/fnxAQECvXr3S09MJIUeOHElOTh46dOg777zz9ddfz58/X144MjLyww8/DA0N7dmz5+DBg2/dulX3sAF0xi932I9bYnFRqB9jIZnSCo3C2rw0Rvj06VN5QpJIJHfv3q2e9srKyrZs2VLHXZkeP358+PDh9PR0GxubHj162NraXrt2rUZr8ocffli1alVISAgh5NKlS1u3bl22bNmnn35aVcDAwGDWrFnr1q2TF54/f/7IkSMJIQkJCT/99NOmTZsa/DcDaKMSCTmEaTLQIDNb0x77JanFdFMTXEYp8FKl2rNnz+zZs+U/T5gwoUZRPT29X375pS4njYuL8/DwkG9hIRKJ2rVrFxsbWz0RSqXSuLi4Ll26yB926dIlPDy8xknS0tIcHR3lP8fGxn711VdVhb/99tu6hAGgS3ansD0dMU0GGsJMRCa1pL9OYLd0wfRRBV5KhMOHD/fx8SGEhIaGLl26tFWrVlW/MjIy8vT0tLKyqstJs7OzLS0tqx5aW1tnZ2dXL5CbmyuTyarKWFlZ1Sjw8OHD5cuXHzx4kBAikUjy8/OrF87KynrdS9+9ezciIuLo0aP//XkCwY4dO0xNTesSNp+VlpZSFL5h6620tFRlr7Xltmipr7SkROs7uMRisUgkYhh8I9fPW1bSSc0o/3+EMzzErkaNGJSmE4vFHPfmwdGXEqGrq6urqyshZOfOnZ06dapj2nuVkZFRRUVF1cOysjIjI6MaBQghVWXEYrGxsXHVbzMyMoKDg5ctW9arVy9CiEAg0NPTe13hGlxcXDw8PKpGJQ0NDR0cHPAV/0Ycx9XyrkItVPO+xeZyzyWyIc0NdWCAkGEYJMIGeMtKakzIRy1lWx4a/BjEo3fewMCgLt//iifL5OXlPX36tMGv7eLi8ujRo6o8nJaWJs+vVUxNTc3MzNLS0qoKuLi4yH/Ozs7u3bv3J598Mn36dPkRiqKcnJwUFn6VoaGhp6dnn//XqVMnZEHQAZgmA29vjg+z+wGbWYbpozUpToQLFy5s3bq1o6PjmDFjfvvtt8ePH9frpJ07d2YY5p9//iGEXLt2LTU1deDAgYSQ5OTkv/76S14mNDR027ZthJDS0tJ9+/bJJ8Lk5ub26dNnzJgxc+fOrX7CUaNG/fbbbxzHSSSSP/74Q14YgCdKJORQGjveA2kQ3oqdARnrQf+QpPW9641OcSK8e/fusWPHRo0alZiYOHHiRDc3Nw8Pj8mTJx84cKAuJxUIBD/99NOHH37Yo0eP/v37//DDD/JRurNnzy5fvlxeZvHixTExMe3bt/f29vb393/nnXcIIZs2bbp9+/b3339vaWlpaWlpa2srLzx79uy0tLTAwEBvb29zc/Nx48Y1wp8OoCX2pLC9ME0GGsNcH/q3e2wudk94GfXGgcScnJzTp0/v2rXr+PHjhJC6DDzKFRYW3rt3r2nTpvLpo4SQioqK8vJyMzMz+UOpVJqYmGhqauru7i4/IhaLa2xwYWFhIf+BZdmkpCSRSOTl5VXLi06fPt3d3X3mzJl1DBLkiouLTUxM1B2F9ikpKVHBGGHHv6VLApgBLjqSCDFZpmEaq5JOuSCz0COr2vLi/T98+PCOHTuqpk++Tm33JBUWFp45cyYqKioqKurOnTu2trby2St1ZG5u3r59++pH9PT0qt+bKBAIAgMDqxeoZYMLmqb9/Pzq/uoAuuHucy69lAQ760gWBLWb708HHpZ+7sNY1HV9FN2nuGt0xYoVHTp0sLa2Hj9+fFpa2qRJkxISErKzs//8808VxwfAc7/eYce5U9h0CRqLixEV4kb/eAsjhS8obhH+9NNPBQUFn3766fz586vuagcAFZOyZHcKe3ogVpOBxrTAn25/RDrTmzYXqTsUzaC4Rbhq1arBgwfv2bPH2dnZz89v9uzZ4eHhr+5HAQBK9e8TtpkJ1dIc7UFoTM1MqAEu9OZkNAr/ozgRTpw48dChQ7m5ubGxsR988MGtW7dCQ0MtLS27du2q4vgA+Oy3e9yHnth0CRrfQn96/U1ZiUTdcWiG2rpcaJoOCAiQ/1xRUREdHX3+/HmVRAUA5FkFic5md/UQqjsQ0EEtzakeDvTPd9g5PrjSek0ivH///qlTp6Kios6cOZOfn6+npxcUFLRixYrevXurOD4A3tp5j33HlTZBHgTl+CqA7ntcOsWLNuD9GLTiN6Br1665ubkBAQFjxozp0qVLv379sG41gIrtvM/yallIULHWFlQHG/rXu+z01nxvFCpOhDt37mzfvn3VzewAoGLX8rhiCenmgGkyoERfBdIhkbJJLWk9fl9xKb4Q6NevH7IggBr9do+d4IlFtkG5AqwoH0uy4z7fp4++SISJiYl37typy3NOnz6dm5urtJAA+K5cRvalsOOwyjYo31cBzOp4tpLfqfBFIszIyPD19R09evSJEydkMtmrRUtLS//444/u3buHhIRUVlaqMEgAfjmSxraxplyNkQhB6TrYUp5mZPcDXmfCF2OEAwYMiI2NnT9//sCBA01MTNq1a9eiRQsrKyupVJqTk3Pz5s34+HhCyLhx4/bu3evg4KC+mAF03G/3WNw+CCqzLJAZFy0b604L+Pqhe2myjJ+f37///nvv3r3ff//9zJkzu3btkq8mY21t3a5du9WrV48dO7ZqKwkAUIb0Uu5aHnekL1+/k0DlOtlRToZkfyr7fnOefuoUzBr19PRcsWKF/GexWCwQCIRC3MoEoCI77nGjm+PWLlCpxQHMjEuy0c1ofk7QekP+NzAwQBYEUBmOkJ330S8KqtbHibLUI3+l8XSkEPUNQIPEZHH6DGlrzcvLclCr+f7Myji2rhuv6xYkQgAN8ts9dkIL1EpQg0EulJAmxx7xsVGIKgegKUok5O/H7Ad8nbAAarcwgF4Rh0QIAOrz92O2sx1la6DuOICvQtzoCpYcf8K7/lEkQgBNsTeFG90MVRLUhiJkoT+9PE7Bgiq67bW1juO4y5cvb9u2bdOmTfIjGRkZz549U1VgAPxSWElistkhbkiEoE4jm9LFEnI6k1+NQsU3K+Xn54eEhMi34XVycpoyZQohZM2aNbdu3Tpz5oxKAwTgh7/S2L5OtCluVgK1oikyz49eESfr5cijW1kVX35+8sknaWlpJ06cOHHiRNXB995779y5c8XFxaqKDYBH9qWwo5rhrglQv/eb0xll5Fw2jxqFChJhWVnZ0aNHN2zY0K9fP319/arjLVq0kMlkT548UWF4ALyQV06u5HIDXNAvCurHUGSeL70qnkcjhQoqXkFBgUQi8fLyqnGcoihCiFgsVkVcAHxyMJUd5EIb8agvCjTaOA/67nNy8SlfGoUKEqG1tbW+vn5CQkKN49HR0TRNN2vWTCWBAfDIvofoFwUNIqTJF7702gS+3FOoIBHq6emNGDFi3rx5N27ckLcCCSGXL1/+7LPPBg0ahJ3rARpXtpgk5nP9nNEvChpkgicd94y7nseLRqHiurd+/Xpra+s2bdqMGDEiNzfXzc0tKChIJBJt3gp2wZ4AACAASURBVLxZxfEB6Lx9KWyIG63HqDsOgGr0GPK5D72GH41CxYnQysrq4sWL27Zt69mzZ+fOnf39/b///vu4uDgnJycVxweg8/Y9ZEfhPnrQPJ940Zeeckn5ut8ofO3ovJ6e3oQJEyZMmKDKaAD45kkp96CI6+WIAULQOPoMmelNr0lg9/TU8f4KxdehsbGx0dHRNQ5evXpVfos9ADSWvSncsCa0EA1C0EhTW9GnM9m7z3W8Uai4/o0fP/7y5cs1Dt6/f3/YsGFSqVT5UQHwxb6H7ChsNwGaykhAprVi1sTr+EihghpYUlKSnJzcq1evGsd79uyZm5ubmpqqksAAdF9KEZdZxnWzR78oaK4Z3nT4E/a+TjcKFSTCoqIiQkj1NWXk5EcKCwtVEBYAH/yZwo1sSjPIg6DBTIXkUy/6myRdbhQqSIQ2NjZGRkavLq595swZiqJcXV1VEhiA7sN8UdAKs7yZQ6nsoxKdbRQqqIRCoXDUqFELFy7cv38/x/33l0dEREybNq1v3752dnaqjRBAN90p5J5XkiA7tAdB01nqkY9b0t8k6myjUPHV6Hfffefl5TVq1CgTE5OWLVuam5v379/fzMxs27ZtKo4PQFftSWFHNaOQBkErfO7D7E1hM8t0s1Go+D5Cc3Pz8+fP79279/Tp01lZWe3atevatevYsWMNDAxUHB+ArjqYyu3qoeO3Z4HOsNYn4zzo75LY7zro4If2tTfUi0SicePGjRs3TpXRAPBE3DOuQkbaWKNBCFrjCz/G+6Bkni9jq3MNIgzUA6jBvofse83RLwraxN6AjG5O/3BTB/cpVJwInz9/Pm/evBYtWhgYGFAvU3F8ADrpYCqH+aKgdeb50b/cYfPK1R1HY1PcNTp69OioqKihQ4e+//77r95QCABv43IOJ6SJjyUuK0HLuBhRw5vS/7slW9ZGp0YKFSTC8vLyU6dOrV+/furUqaoPCEDn7XvIvodl1UA7felHtz8i/cyHMRepO5TGo6A2FhcXS6XSTp06qT4aAJ3HcuRAKjeyKZqDoJWamVCDXOmfbuvUPYWKV5bx8fG5du2a6qMB0HmXczgrPeJljkQI2mqhP/3jLVmJRN1xNB7F/TM7d+78/vvv9+/fj5VFARrXkUfsu02QBUGLeZpRvRzpzcm60yhUnAgHDBhw586dUaNGWVhYvM2s0fT0dPkS3q9TUFCQnZ1dr3MCaLWjj7gQNwwQgnZbHED/cFNWpiub8imeNbpo0aLS0tK3OW9GRsaQIUPy8/OfP38+ZcqUVatW1SjAsuzkyZMPHTpkYGDg7u5+5MgRc3NzQsi8efOOHTt29+7dH3/8sWq2zrFjx955552q54aHhw8cOPBtwgNQi1sFXIWM+FuhRQjarZU5FWRLb73DzvTWhas6xYlw2rRpb3neRYsWBQQEbNu2LSsry8/Pb8iQIR07dqxe4OjRo6dPn3748KGJiUlISMi6detWr15NCGnbtm1ISMi8efNqnLBz587nz59/y6gA1OvII25oE9yNC7pgcQA9MEI6qSVt8NoFyrSGUpK5VCrdt2+fPJs6ODgMHz58z549Ncrs3r177NixZmZmNE1PnTp19+7d8uMjR47s1KmTnp7eq6cVi8XKiBZAZY4+YtEvCrrB34pqY03tuK8LI4WvTeW5ubknTpxISUkpL39pFYG1a9e+8aRPnz4Vi8Wenp7yhx4eHmfPnq1RJjU1ddiwYVUF0tPTJRKJUCh83TmvXr1qb2/PsmxoaOiGDRuMjY0VFqusrHzy5Mn169flDw0MDFq1avXGgAFUIKOUSy3mumLfJdAViwOYEadkE1vQIi2/ulOcCC9cuDBo0KDi4mKBQEDTtEwmk2cpY2PjuiRC+QSZqq0qjIyMnj9/XqNMSUlJ9QIsy5aUlFhYWCg8Yfv27Z88eWJnZ/fkyZPhw4fPnz//f//7n8KS9+/fv3fvXlRUlPwhwzBHjhwxMzN7Y8w8V1JSou4QtFJpaWnVnp1vtPc+E2xPi0t1bn2qehKLxSKRiGF0amkSFdDASuqlTzyMhb/erBzbVEMXIBWLxXWpoYoT4fTp01u3bn306NG5c+c6OTktXbo0MjJy2rRpa9asqctryzfvLSwstLS0JIQUFBS8up2vra1t1b0Z+fn5enp68skytZyQEOLi4rJgwYK5c+e+rmTr1q1DQkJmzpxZlzihOhMTE3WHoH0oinpd58SrIp5KP/WiTUz4vmahQCBAImwYDaykKztwY8/KJnkLBBrZKJQvl/3GYgpir6ysTExMXLFihbW1NSFEKpUKBIKBAwdu3rx58uTJFRUVbzyphYWFi4vLlStX5A+vXLni5+dXo4yvr+/Vq1flP1+9etXX17eO92YUFBQYGRnVpSSA5nheSa7mcsHOGvltAdBQQbaUixHZ+1C7RwoVtAgLCgpkMpmzszMhxNTUtKrd1rVr18LCwjt37rya1WqgKGrq1Knz58+3sbFJSko6c+bMTz/9RAhJTU0NCQm5cOGCiYnJ5MmTu3Tp0rdvX3t7+xUrVixZskT+3CtXrjx+/Pjp06dxcXEHDhwICgpydnbeunWrtbW1m5vb7du3FyxY8OWXXzbmewCgfOFP2B4OtJH2z68DqGFRADPtouz95jSttcPfCuqltbW1SCTKysry9PRs0qTJ9u3bWZalafr27duk2shf7T7//HOWZWfMmGFpaRkeHu7o6EgI0dPT8/b2lneJ+Pj4/Pnnnxs2bBCLxbNnz67aAfjKlSvnz5/38vIqKio6cOCAi4uLs7OzkZHR9u3bs7OznZycvvvuu/fee6/R3gAAlTj6iAtx09rvCYDX6+1IWeuTg6lsqNbuLEYpHEjs1atXx44dV69enZWV1bRp086dO/v7++/bt8/ExOTmzZua3Lk/ffp0d3d3jBHWV3FxsQYOP2i+kpKSuowRVsiIwx7J3ZFCG76PDxKCyTINpcmV9PgTbu5VWeIwgaY1Cg8fPrxjx46jR4/WXkxxAl+/fn3v3r0JIQ4ODvv27SsuLv79999btWp1+PBhfHwB6isqk/O1pJAFQVcNcKEMBeTvx9o6Uqh4yMLX17fq55CQkJCQEFXFA6CDcB896LyF/vSyG2yIm6a1CesElRNAuViOHHvMYoAQdNs7bjRFkeNP6npbrUZ50SK8ePHizp073/iEn3/+WZnxAOiayzmcjT7VzASJEHQZRciXfvTyONlAF+2bG/0i4szMzFOnTlU9zM3NLS4upmna0tKysLBQKpUaGBg4ODioI0gALXb0ETsUzUHggeFN6OU32FMZXB8nLfvAv+gaHTFiRMr/W79+vZmZ2b59+8rLy3Nzc8vLy//9918nJ6cFCxaoMVYAbYQNCIEnaIp86UcvuaGhy63VQkH9ZFl2ypQp33//fWhoqHwVbIZhBgwYsHPnzhkzZmjgencAGut2ISeWkQBrLbtABmiY0c3pvHISk61lI4UKEmFOTk56enr1iaNyPj4+ZWVl9+7dU0lgALrgSBo31A0bEAJfMBT50o9eGadljUIFidDExEQoFP7zzz81joeHhxNCrKysVBEXgE7AjRPAN2Pc6QdF5MJTbWoUKpjeY2RkNH78+C+//PLJkyfvvvuug4NDbm7u8ePH169fP2jQIDc3N9VHCaCNMsu4lCKumz0ahMAjQprM86PXxMv+6ac100cVB7px40ahULhp06YNGzbIj9A0PWrUqC1btqgwNgDtdiSNG+RKa+b2NADKM8GTXhPPXsvj2mrJ6LjiRKinp7dp06bly5ffvHkzIyPD3t6+VatWuHcCoF6OPmI/8UIaBN4R0mS2D706nv2rj3YsyVlb09Xa2rpHjx6qigRApzyvJFdyub/6IhECH01qSa9LlCbmc76WWtAofG0ilEqlSUlJT548qbET78iRI5UfFYDWC3/CdscGhMBX+gyZ5U2vjmf39tKCRqHianr+/PmwsLCUlJRXf6Vw2yYAqAEbEALPfepFN0+U3CqgW1toekVQnAjHjx/PMMzBgwdbtWqlp6en4pgAtF2FjJzMYDd2Eqo7EAC1MRKQGa2ZrxPYXT00vVGoIBEWFBQ8fPgwIiIiODhY9QEB6ICoTM7HAhsQAt9Na0033ye5/5z2MNPoRqGCkXyRSCQQCExNTVUfDYBuwH30AIQQUyGZ2or+OlHTN+xVUFeNjIxGjRq1Z88e1UcDoAPkGxAOcdXoS2AA1fjMm/n7EZtWrNGTSxSPEQ4aNOizzz5LT0/v37+/hYVF9V9h1ihA7a7kcjb6lIb3BQGohpmIfNSCXpfIbuqsuSOFihPhZ5999vTp08OHDx8+fLjGrzBrFKB2/6A5CFDNbB+m5QHJQn/ayUhD64XiRHjlyhWZTMuWDwfQEBHp3PqOmnvxC6Bi1vokzJP+Non9QVPrheJEiJW1ARomr5ykFHEdbDX0yhdALeb6Mq0OSr7wZRwM1R2KIpjYBtCYIjPYno60EBULoBo7A/KBO73+poZ2NCpuEbZo0SI3N1fhr/Lz85UZD4B2i0jn+jmjOQhQ0xe+tN9f0s99GQ28v1ZxIhw7dmxJSUnVw7y8vHPnzuXm5k6cOFFVgQFoH46QkxnskkAsMApQk7MRNaIp/eMt2Yo2GjdSqLjGLlq0qMYRqVQaFhb2/Plz5YcEoK0S8zkTIdXMBC1CAAUW+tNtjkhnezMWGrZwZ12HMgQCwaJFi7Zt2/b06VOlBgSgvdAvClALV2NqsCu98bbGLTRTjzF9IyMjlmUzMzOVFw2AVotIZ/s5Y54MwGst8qc33pYVS9Qdx8vqWmkLCwsXLFggEAg8PDyUGhCAliqTkthcrrs9WoQAr9XclOrjSG/SsEZhXWeNFhQU0DS9fPlyY2NjlQQGoGXOZnHtbChj7LwEUKtFAXSPcOnUVrTmVJY6zRolhDg4OPTp06d169YqiQpA+0Sks8FO6BcFeAMvc6qbPb31LvuZt6bUl7rOGgWA2kWkc3t7aUrFBtBkiwLoASekk1vSBppxq9Eb6m1FRUVaWppYLFZNNABa6lEJV1jJ+VlhgBDgzfwsqXY29PZ7mjJS+NpEuHv37tatWxsYGDRt2tTQ0NDd3f2nn37C1hMACkWkc8FONNIgQB0tCaDXJbKVmpEKFSfC9evXjxkzRk9Pb8WKFVu3bl27dq2Tk9O0adPmzZun4vgAtALuIASol0BrqpU52XVfIzKhgg7aysrKJUuWTJw48ddff606OG/evCVLlqxevfqLL76wtrZWYYQAmk7GkegsdnNnjZkDB6ANlgQyH5yRhXnQAnWPrSt4/by8vKKiohkzZtQ4PmPGDKlUmpaWpoq4ALTH5RzOzYSyNVB3HABapaMt1cSE7ElRf6NQQSK0trY2MTEpLCyscVx+K6Grq6tKAgPQGpHpbD8n9IsC1NuiAGZFHCtT9+QTBYlQJBLNnz9/5syZjx49qjqYk5MzefLkyZMn29raqjA8AC0QkcFhZTWABujpQDkYkgMP1dwofDFGePbs2U2bNlU9fPz4sbu7e5s2bRwcHHJzc69fv07TtIWFhTqCBNBcBRUkuYDriC3pARpkoT8z+7IstBmtxlnXLxJhcXHxw4cPqx42bdqUECKVSp88eUIIka8pk5qaqvIIATTayQy2uwOtp3E7rAFoh37OlLkeOfKIHdZEbd0qLxLhkCFDhgwZoq44ALQUbpwAeEtf+tFfXWffbaK2NiEGNgDeyskMJEKAtzLElWYoEv5YbXNmalvoraysLDU1tby8vPrBNm3aKDkkAK2R/JwS0MTdFIkQ4K3M96NXxMkGu6pn7VHFr5qfnz9hwoS///771TXVsMoaQJVTWXR/NAcB3tqwpvSyODYygwtWx51IirtGJ0+efP78+S1btgQHB7///vsHDhz48MMPjY2Nt2zZouL4ADRZVDaNflGAt0cRMt+PXnpdppZXV9Ai5DguPDx806ZN48ePv3TpkqOj44gRI0aMGNGiRYstW7Z88skndTkvy7J//vlnYmKiu7t7WFiYUKhg9amrV68eOXLE2Ng4LCzM0dFRfjA9Pf3atWs5OTkjR46sfrdGUlLSvn37RCLR2LFj5TNaAdSrXEau5lE9HDDQDtAIQpvRy26w0VlcdwdVX1wqqMM5OTllZWUdOnQghOjp6RUXF8uPT5w4MT4+PiUlpS7nnT179jfffOPs7Lxr166xY8e+WuDUqVPBwcFmZmbp6ent2rXLz88nhBQUFLRu3Xr9+vWTJ0/OzMysKnz9+vXOnTvLg2nbtq38jg4A9YrJ5nwtODORuuMA0AkMReb70yvi1NAoVNAiNDMzoyhKnv8cHR0vXrwoPy6VSgkhVXmxFnl5eT///POtW7eaNWs2duxYR0fHe/fueXp6Vi+zdu3apUuXzpo1ixDy8OHD3377bc6cOebm5vKF3ASClwL75ptvZsyYsXjxYkJIdnb2pk2b1qxZ08C/GKCRRKSzvR0wZA7QaMa406vi2fPZXBd7lTYKFbQI9fX1W7RocePGDUJI3759o6Kifvrpp8uXL0+ZMsXY2NjDw+ONJ71y5Yqrq2uzZs0IIebm5u3bt4+Ojq5egGXZmJiYfv36yR8GBwefPXuWEEJRFE0rCOns2bPBwcE1CgOoV0Q619tePUMaADqJochcH3p1gqqrleJZo7Nnz87JySGEBAUFjRs3btq0aYQQPT29zZs3GxkZvfGkWVlZNjY2VQ/t7Oyq93MSQvLy8iQSSdWypXZ2dllZWa87m1Qqzc3NrV64xtmqe/To0bVr15KSkuQPhULhypUr6xIzz5WXlyscx4XXySwj2WV0SwNxeTkWlamf8vJylmUZBu9b/fCkko52JSvjmAsZFW2sGqG7pbKysi53OihOhB9//HHVz9u2bVu8eHFaWpqXl5ednV1dXlsgEMhkL1K6RCKp8f8n7/mU97XKf6jlP5imaYZh6ljY2NhYX1+/bdu2VS9kZGTEh0/PWxIKhXiX6uV0DtfXieiJ8L7Vm7wKIxHWF08qqZCQOT7culvMX70boXdUIBBQ1JvPU6e7F5s0adKkSZO6v7ajo2P1RltmZqaTk1P1AhYWFgYGBpmZmfLMmpmZWTVr9FU0Tdvb22dmZnp7e7+xsJWVlbu7++TJk+seLRBCGIbBF1O9nMyUDXSh8L41APP/1B2IluHPmzbJi6xLkiYU0IHWb5sLFY61KSj2li+jUJcuXYqKiq5evUoISUtLS0hIkA8HZmVlyYceKYoaPHjwwYMHCSEymezw4cODBw+u5YSDBw8+dOiQ/OeDBw9iTVRQLxlHojLYvtiDEEAJ9Bky25tem6C6vZmUsp6NoaHh8uXLhw4dOmTIkMjIyDlz5tjb2xNCjhw58vPPP8fHxxNCFi1a1KtXrydPnsh3PRw9erT8uaNHj3727BnLsh9//LGRkdHBgwfNzMzmzp3buXPnwsLCoqKiR48effTRR8oIG6COruVyTkaUoyFVUqLuUAB00WQvel2i5GYB7W2histNZS3sNn369O7duyckJEycOLF9+/byg8OGDevUqZP8Z19f39u3b0dFRZmamvbp00dPT09+fMGCBRKJpOo88nkuTZs2vX37dmRkpEgkCg4ONjQ0VFLYAHURgYW2AZTJUEBmeTNrE9g/eqiiN1iJK5z6+vr6+vpWP2JnZ1d9uo2tre1777336rMUns3c3Dw0NLTRgwRogMh0dlkbXozWAKjLtFa0+37Jvee0p5nSLzqxOhRA/RRJyM0CrosdWoQASmQsJFNaMaoZKaytRVhRUZGVlcWyL8Uhv00egLdOZbCd7ChsSQ+gbDNb0x77JanFdFMT5V53Kk6EKSkpM2bM+Pfff1/9FbZhAp47mcEFO6ErBUDpzERkUkv66wR2SxflXngqToShoaFpaWnLly/39PSs430YADwRncVN9kKlAFCFz3wYz/2SBf60q7ESG4UKEmFxcXFcXNzevXsxOQWghtxyklXGqWZKNwBY6ZGJLejvktgNQUpsFCq4sJXJZBzHNW/eXHmvCqClorPYrvY0gzwIoCpzfJg/HrCZZUoclVOQCM3NzXv37h0REaG8VwXQUjHZXDeV7xoKwGd2BmSsB/1DkhKnjyoeI1y+fPmYMWPKysr69+9fY2FPzBoFPovJ4rZ2xQAhgErN82V8Dkm+8GNs9JVyfsWJcNiwYU+fPl21atWqVatq/AqzRoG3CipIajEXYIUWIYBKORiSkc3o9Tdlq9oqZaRQcSLcunVreXm5Ml4PQHudy2aD7CgBGoQAKjffj25zRPq5D2Oh1/gnV5wIsb0DwKtisrlu9kiDAGrgaky940r/eItdEtj4dbC2M4rF4sTExH///TcuLq4Ey+wD78Vkc90xUwZATRb40/+7JSusbPwzK06EHMctW7bM1tbWz89v0KBBgYGB1tbWs2bNqqxUQggA2qBEQu4Ucm3feqdQAGiY5qbUABd6c3LjTx9V3DW6evXqpUuXDhkyJDQ01MHBITc399ixYxs3bpRKpRs3bmz0IAA03/mnXFtrLDEKoE4L/enu4dLprWhjYWOeVkEilMlk33///axZs3744Yeqg6NHjw4ICJg/f/6qVavMzMwaMwQAbRCTxeIOQgD1amlO9XCgf77DzvFpzJFCBefKycnJz88PCwurcTwsLEwqlT548KARXx5AW8Rkc90dMFMGQM2+CqC/S5KJpY15TgUV29DQkKKozMzMGsflR4yNjRvz9QG0gVhKEvK5DjZoEQKoWWsLqoMNve1eY44UKkiEZmZmQUFBM2fOTE5Orjr4+PHjTz75pFmzZp6eno348gBa4VIO52dJGda2fScAqMhXgfS6BLZC1mgnVFyzN23a1KtXL29vb29vb0dHx9zc3MTERKFQ+M8//1AULoqBd2KyWdw4AaAhAqwoH0uy8z47qWXjjFYoPoufn19SUtK8efNMTU0fPnwoEommTZuWlJTUs2fPRnlVAO0Sk4Vb6QE0yFcBzKp4trKR+kdf29fj6Oi4evXqxnkRAG1WyZLreVwnO7QIATRFB1vK04zsecCGeTbCFariUyQnJ1+7dq3GwZs3byYkJLz9SwJol6s5XAtzyqRR71sCgLe0yJ9ZGc9KG6NRqDgRjh49Ojw8vMbBy5cv9+/fXyZrvAFKAG0Qk811t0dzEECzdHegnAzJ/tRGyIQKEqF8idHg4OAax/v375+dnZ2Wlvb2rwqgRWKycSs9gCZaFMCsjGPZt94bUEEiLCgoIISYmprWOG5iYkIIyc/Pf9vXBNAeUpZczuE622GmDIDG6etEWeqRv9LetlGooHpbW1vr6+tfunSpxvELFy4QQpycnN7yJQG0SNwzrokxZamELdAA4O3N92dWxikhEYpEopCQkC+//PLMmTNVB69fvz5t2rQuXbo4Ojq+5UsCaJHobA79ogAaa6AL9W4T+i2nzCju8Fm/fr2VlVWvXr2aNm3arVs3Dw+Pdu3aVVRUbNu27a1eDUDbxGRhpgyA5qIIWRJIC95u7ELxs+3t7a9fv/7tt9+2atWqrKysadOmK1euTEpKwvpqwCssRy4+ZbviVnoAnabghvqioqKPP/545syZc+bMmTNnjupjAtAQifmcrQFla6DuOABAmRRc6paXl+/fv59hsAMp8F1MNoclRgF0noJEaGNj07Rp08TERNVHA6BRYrK5bhggBNB1ChIhRVGbN29euXLlP//8g3VkgLc4Qs5ls12QCAF0neJFt7/66qvCwsIhQ4YIhUIHBweB4EWxlJQUVcUGoE7JhZyJkHIxQiIE0HGKE2G3bt38/f1VHAqARsGNEwA8oTgRfvPNNyqOA0DTxGRz/ZyRCAF0H26QAlAMM2UAeOJFi/DOnTunT5/u2LFjYGDgjh07ysrKFD5hypQpqooNQG0eFHEUIU1NkAgBdN+LRHjp0qWpU6euXr06MDBw3rx5OTk5Cp+ARAh8EJPN9cAdhAD88CIRjh07dsSIEXp6eoSQhw8fsmxj7PsLoJ1istAvCsAXLxKhQCCQ7zhICDEyMlJTPAAaITqbm++PEXQAXnipqm/evPnw4cOvK3r+/PnQ0FDlhwSgZumlXLmU8zRDixCAF15KhF9//fWwYcM++eQTiUTyatHHjx8fOHBAVYEBqM3ZLK6bA400CMATNTt/PDw8tm7d2q9fv7y8PLUEBKB2uHECgFdqJsIPP/xw7969V69eDQgIuHbtmlpiAlCvaMyUAeATBdMBQkNDL168KBQKu3btumvXLtXHBKBG2WKSV861tkAiBOALxfPifH19Y2Nju3TpMn78+JkzZ2IPCuCPmCy2mz2NEUIA/njtBHErK6sTJ07Mmzfvxx9/HDx4cEFBQb3OW1xcPHr0aDMzMzc3tx07digs89133zk4OFhaWlafnpOent63b19TU1MvL6+TJ0/KD54+fbp5NWfPnq1XMAB1F5PNdcOt9AB8UtudUgzDrF27dseOHWfPnu3YsWNycnLdz7tkyZKioqLMzMxDhw7NnDnz1edGR0evW7cuOjr64cOHCQkJ69evlx//5JNPWrRokZeXt2bNmtDQ0KKiIkJIaWmphYXFyf/Xvn37ev6ZAHWFAUIAvnnzLcPjx4+/ePFiRUXFypUr63hSlmV37NixcOFCIyOjtm3bhoSE7Ny5s0aZbdu2jR8/3tPT09zcfO7cudu3byeEZGZmRkZGLl68WCQSDR06tEWLFocOHZKX19fXb/b/DA0N6/M3AtRVfgVJL+X8rZAIAXjkpUTo6upqYWHxaqGAgIDLly937ty5jifNyckpKCjw8fGRP/Tx8bl3716NMnfv3q1e4MGDBzKZ7MGDBzY2NnZ2dq8+MSkpqXnz5m3atPn666+lUmkdIwGol5hstpMdxSAPAvDJS/sRxsTEvK6cvb39+fPn63jS/Px8iqKMjY3lD83MzF69K7GgoKBqRTdTU1OpVFpUVFRQUFD1LPlx+RNbMV5oAQAAGRlJREFUt2597Nix5s2b37lzZ9KkSRUVFV999ZXCl46Li9u4ceOsWbPkD01MTJKTk83MzOoYOW+VlpZSFL7+SdRjQUdLUlJSXsfypaWlSo1HJ4nFYpFIxDCMugPRMqikDSAWizmOe2MxxRvzviVra2uO44qLi+UZqLCw0MbGpkYZKysr+fgfIeT58+dCodDMzMzKyqq4uLiqTGFhobx1KO8RJYQ4OTmtWLFi1apVr0uEAQEBI0eOnDlzpjL+Lh3GcVz1SxDeuvhMurETY2ysX/en4H2rL4ZhkAgbAJW0AQwMDOpy9aCUZYWtra2trKwSExPlD5OSkjw9PWuU8fT0TEpKqirg4eFB03Tz5s2fPXuWnZ1dyxOFQiF2xgBlKJKQ+8+5Nta46AbgF6UkQpqmw8LCVq1aVVRUdOnSpb///jssLIwQkpGRMXLkSPmWvx999NHOnTtv376dl5e3bt26CRMmEEIcHBz69eu3bNkysVh88ODBBw8ejBgxghASHh5++/btkpKS2NjYJUuWvPvuu8oIG3jufDbX3oYSYc8JAJ5RStcoIWTZsmVTpkxp3ry5mZnZ5s2bW7RoQQiRSCQpKSny9lzXrl0XL148YMCA8vLyUaNGVXVm/vLLLx9//LGTk5Ozs/Nff/0l7wq4c+fOzJkz8/Ly7O3tQ0NDFy5cqKSwgc9OZbC9nZAGAXiHqstAohaZPn26u7s7xgjrq7i4uGruEm95H5Lu7M7Uq2u0pKQEwzb1hckyDYNK2gCHDx/esWPH0aNHay+G618AQgjJKOWeirkA3EEIwD9IhACEEBKZwfV1whKjAHyERAhACCEnM7i+TkiDAHyERAhAWI5EZbJ9kAgBeAmJEIDEPeNs9CkXIyRCAD5CIgQgkRlcMJqDAHyFRAhATqazfXEHIQBfofID35VKybU8bMYLwF9IhMB3Z7O49jaUkbIWWQIATYdECHx3MgP9ogC8hvoPfBeZzgU7o18UgL+QCIHX0ku5vHLOzxKJEIC/kAiB106kc8HOWFkNgNeQCIHXsLIaACARAn+xHDmTySIRAvAcEiHw17U8zt6QcjREIgTgNSRC4K/IdKysBgBIhMBjuIMQAAgSIfBWsYTEPeO62qNFCMB3SITAU2cy2Y62lCFWVgPgPSRC4KmTGRz6RQGAIBECb0VmYGU1ACAEiRD46VEJ97yS88XKagCARAj8FJHOBTthYTUAIASJEPjpZAbXF/2iAEAIQSIEHpJx5Ewm28cRH34AIASJEHgoNpdzNqIcDNUdBwBoBiRC4B3MFwWA6pAIgXewshoAVIevA+CXYglJeMZ1sUOLEAD+g0QI/BKVyXayowywshoA/D8kQuCXkxlcMPpFAaAafCMAv0SmY6YMALwEiRB4JK2YK5FwrS2QCAHgBSRC4JET6Vx/F6ysBgAvQSIEHjmZwfV1Qh4EgJcgEQJfyDgSnYWV1QCgJnwpAF9cyeHcTChbA3XHAQAaBokQ+CIygw1GvygAvAKJEPjiZAaHldUA4FX4XgBeKKwkN/O5zlhZDQBegUQIvHA6k+1iT+kx6o4DADQPEiHwAvpFAeB18NUAvICV1QDgdZAIQfc9KOIqWNLKHIkQABRAIgTdF5nO9cONEwDwGkiEoPsisbIaALyeEvcnvX79+o0bN7y8vLp06aKwQGZm5smTJ42NjQcOHGhg8GLBj9OnTz98+LB9+/a+vr5VB/Py8k6cOCEUCgcOHGhiYqK8sEHHSFkSncVu7SpUdyAAoKGU1SL8/vvvQ0JCEhISwsLC5s2b92qB+Ph4b2/v6OjoLVu2dO7cWSwWy49Pnjx56tSp8fHx/fr1+/XXX+UHU1JSWrdufeLEiT179gQGBubn5yspbNA9l3I4d1PKRl/dcQCAplJKIiwtLV2+fPmRI0c2btx4+vTpH3/8MTMzs0aZlStXTp06dfv27RERESKRaM+ePYSQBw8e/P7779HR0Rs3btyzZ8/ixYsrKysJId98882wYcP++OOPo0ePenh4/PLLL8oIG3TSyQwW80UBoBZKSYTnz583Nzdv27YtIcTV1dXf3z8yMrJ6AY7jwsPDhw8fTgihaXro0KHh4eGEkH///bdz5862traEkB49ekgkkuvXrxNCwsPDhw0bJn/u8OHD5YUB6iISdxACQK2UMkaYmZnp5ORU9dDJySkjI6N6gfz8/PLy8qoyjo6O8iZj9SdSFOXg4JCRkcGybHZ29quFFcrNzU1NTT3NelZQIkKIUChs3749w2BBkTeoqKD09CTqjqLxsYTcKSTtLKUS5fxxEolEoqRT6y6JREJRFMuy6g5Ey+DD1gBSqbQuxZSSCGUyGUW96IyiaVomk9UoQAipKsMwjDxchU9kWZZl2VcLK1RRUVFaJq6oIDKaI4TQUumzcplAgJ6xN6is5EREN7+Y1rUhDMe9/AFsNDKZTKakU+suvGkNg/etAViW5TjujcWUkggdHBxycnKqHj59+jQ4OLh6AWtra6FQmJOTY21tLS/g6Ogof2JycnL1Jzo6OgoEAltb25ycHC8vL0JIdna2vLBCzs7OPXq4z5w5oNH/KN1WXFxpYqKn7ii0j1Qq1dfHPJz64ThOJBKhn6a+JBIJPmz1JRKJqjeuXkcpYydBQUHp6ekPHjwghBQUFMTGxvbo0YMQIhaL5RM+aZru0aPHiRMn5OVPnDjRq1cvQkjPnj3PnTtXWlpKCImPjy8vLw8MDJQfj4iIkBeOiIiQFwYAAHh7SmkRWlpaTpky5d133w0LCzt48OCwYcPc3d0JITt27Pj555/j4+MJIfPnz3/33XdLS0sfPXp09+7dffv2EUICAgK6d+8+aNCgIUOG/PLLL3PmzDEyMiKEzJ07t0ePHgzDFBUVnTx58saNG8oIGwAAeEhZs+m++eabpUuX5ufnT5s2befOnfKDvXr1Wr58ufznnj17nj17lmXZVq1axcbGWlhYyI8fPHgwLCwsPz//22+/Xbx4sfxgQEDAlStX9PT0HB0db9y44ezs/LrXzcvLw12GDRAREVGXnnSo7vnz5xcvXlR3FNonPj6+lvlu8DpVXWhQd8XFxQUFBW8sRunY15+fn5+Hh8fBgwfVHYiWsba2vnPnjnzIFuro5MmTa9eujYqKUncgWiYsLKxr164TJ05UdyBaxtDQ8NmzZ9UX4YI3mj9//o4dO7KysmovpoP3V+lYagfQPaikoBp1/KTpYCIEAACoOyRCAADgNV0bI2zSpEl5ebmPj4+6A9EyMTExQUFBQiG2aKiH/Pz8tLQ0+R0+UHe3b982Nzev5W5gUOjMmTPdu3enabRe6iElJSU/P7+wsLD2YrqWCMPDw/Pz8x0cHNQdiJZJTU1t2rSpuqPQMhKJ5OnTp7XMYQaFnj59amxsLL8zCuoOlbQBioqKxGLxBx98UHsxXUuEAAAA9YJWNgAA8BoSIQAA8BoSIQAA8BoSIQAA8BqzdOlSdcdQb5WVlfHx8SkpKW5ubtWP371799ixY4WFhU2aNKnaekMikZw4ceLy5cvW1tampqbqiFcT3blz5/jx42lpac7OziKRSH6Q47gzZ86cPXvWyMjIyspKvRFqpvz8/GPHjt24cUMkEtnY2MgPFhYWHjlyJDk52dXVVU8Pu1kp9ujRo2vXrjk7O1dtwPTo0aOjR49mZ2c3bdoUdwXUUFRUFBUVdenSJZlMVn0avEwmO3ny5IULF8zNzc3NzdUYoearrKw8fvz4lStXbGxsTExMaivK/V979x7UxNU2APw0gggkkCzBkKIz3EeQaRVE5R5HkCK3UmJxoFwUG2+jpSoy9jJOB61VLjodLx07U6U6HS2ZwkCqsOBEFKkgWKZAi8WiSBMSQlhBCddkvz/O17z5oPDytWiCPL+/dp895Dw5E/aB3ZM9c41EIlm4cCGXy3V0dDSMf//991wuVyQSeXt7v/feezg4NjYWEhLi7++/detWgiBqa2uNkbLJyc/P5/F4e/bsiY+Pd3BwePz4MY4nJyd7e3uLRCIul1tUVGTcJE2QVCq1s7OLjo5OTU0NDAzEwa6uLkdHx4SEhNjYWBcXF6VSadwkTZP+270KhQJHbty4QRBERkaGn59fRESEVqs1boYmRSaTsVisDRs2pKen8/l8kUiE4zqdLjIy0tfXNyMjw87OrrKy0rh5mrKRkZGAgIDAwMAtW7YQBFFXVzdN47lXCNVqtVqtrqqqMiyEOp3Ozc1NLBbTNE1RFEEQP//8M03TYrF42bJlIyMjNE3n5uaGh4cbK22T4ujoWFJSgrcjIyNzcnJomr5//z6Hw6EoiqZpsVjs7u6OF3cGmEajcXBwuHr16oT4/v37U1JS8HZ8fPzhw4dfdmZzwaFDh7Kzsw0LYUBAwJkzZ2ia1mg0Tk5O5eXlRk3QtAwODsrlcrzd2dnJYDAePHhA03RlZeXSpUsHBwdpmj5//vyaNWuMmaVpu3Llire39+joKE3Tx44di4yMnKbx3LscQRAEQRATgr/99ltXV1d0dDRCiM1mr1+/XiKRIIQkEklsbCy+9CcUCquqqoaHh19+zqbGzs5Oo9Hg7cHBQbzohEQiWb9+Pb7YEhMT09XV1dbWZswsTYxUKrWysoqMjJRKpb/++qs+XlZWlpCQgLeFQiH+4AFDTU1NlZWVmZmZ+ohara6trcXjZmlpGR0dDeNmyMrKSn851N7e3tzcfGRkBCEkkUg2btxoZWWFEBIKhXV1dT09PcZM1IThkz9+WpZQKCRJcmxsbKrGc68Q/i25XG5vb6+/PePo6CiTyRBCMpnM0dFRH8QtjZWk6bh06VJubm58fPyaNWs8PT23bt2KEJLJZPqHpOCLz3gMAfbHH3+Ym5sHBgaeP38+KioKDxpCSC6XG37GYNAmGB8fF4lEX331leED/ORyubm5+eLFi/EujNs0jh8/vnLlyuXLl6P/e0LjcDiWlpYwblOZcPLXarXTLMb0Qlaof/m0Wq1+dgxCaMGCBePj4ziuvwmPG+D4PHfx4kVra+vExMTu7u5Tp041NTWtXr16whiamZnBWBkaHh5ub29vbm728vKiKMrV1TUtLS00NNRw3PQfPKB39OhRgUDg6+urVqv1QTxoMG7/1ZUrV77++uvq6mp8HjM8oSH4JZ2W4Vjh+VnTjNUrUgj5fL5ardbpdPidK5VK/FA+Pp+vv3SgUqlomoZH/SqVylOnTuEpHgghiqJOnDghFov5fH5HRwduo9PpVCoVjJUhPp9vb2/v5eWFEOJwOG+++WZzc3NoaKiDg4NKpcJtlEolDNoEeXl5b7/99vbt2/FdiaysrMzMzNdff310dLS/v9/W1hYhpFQq4fnAkxUXF3/44YckSbq5ueGI4QlNo9E8e/YMPm9TMRwrpVKJI1M1fkUujS5btozFYt2+fRshNDY2JpVK161bhxASCAQkSdI0jRAiSXLVqlVMJtPIuRqbmZkZTdOjo6N4d2RkBP+5JBAIpFIpvoxeU1PDYrE8PT2NmaiJCQ0NHRgY6OvrQwhptdpHjx4tXboUIbRu3TqSJHEbkiQFAoERkzRBFy9ejI2NDQsLCw0NRQiFhIQsXryYx+N5enricaNpurKyEv/CAr3y8vKdO3dKJBLDtXQEAkFVVZVOp0MIVVRUeHh4QCGcCj75422SJNeuXWtpaTll65cxfWdWyWQykUgUFRVlbW0tEomOHDmC4/n5+c7OzgUFBVFRUUFBQXjG4/Pnz11dXVNSUk6cOMHlcouLi42au6nYtGnTihUrTp8+/fHHHzOZTDwJW6fTBQYGbty4saCgwNnZ+eTJk8ZO0+Rs3749MDDwzJkz8fHxfn5+eEJaS0uLra3tJ598kp2dzeFwHj58aOw0TVRvby8ymDV6+fJlHo+Xl5eXmJjo5eWFp3YDrLOz08LCIiAgQPSXpqYmmqZHR0e9vb2FQmFeXh6fzy8sLDR2pqZrYGDA2dk5PT39+PHjdnZ2ZWVl0zSee6tP9PX1icVi/S6Xy33nnXfw9rVr12pqapYsWZKeno4nViGE1Gr1xYsXKYrauHFjQECAETI2PVqttri4+JdffrG2to6JicGX+xBCGo3mwoULMpksODg4MjLSuEmaIJ1OV1RU1NTU5O7unpSUtGjRIhx/8ODB1atXGQxGUlKSi4uLcZM0WcPDw99++21KSor+D/Pq6mqSJO3t7dPT0+G74YYmnOUQQhEREfj5If39/RcuXFCpVOHh4XD5YXoqlaqwsLC/vz86OnrNmjXTtJx7hRAAAACYRa/IPUIAAADgn4FCCAAAYF6DQggAAGBeg0IIAABgXoNCCAAAYF6DQggAAGBeg0IIAJhl3d3dhYWF+Ck8AJg+KIQA/CsfffSR6198fHyEQmFpaamxk5qotLS0vLz8pXXX0tKSnp7e2dlplN4B+P96RR66DYCx9Pb2dnR0fPHFFwihgYGBoqKiuLi4s2fP7ty509ip/Udubq6tre1bb71lrN7ZbLaxegfgv4L/CAGYBdnZ2dnZ2UePHm1sbOTz+fn5+YZHBwcHFQqFVqud6scpiqIoSr87Pj6uUCiGhoYmtxwdHZ3qEDY0NPT06dN/9Cb+1/Pnzye8Ql9f31SLmg4NDSkUCv0z3AGYi6AQAjCbWCxWQEDAo0eP8JLijY2NwcHBNjY2eBWnzz//XP9Qw9OnTxMEUVdXt3LlSoIg4uLiEELPnj3btWsXQRB8Pt/KysrT0/Onn37C7fv7+7dt28bhcPh8vo2NTUJCgn75pydPnhAEcenSpfT0dBsbGw6H4+7uXlNTg4+uXbv27t27JEkSBEEQhFAonJx2fX09QRAlJSVRUVH4FcLDwymKamxsfOONN+zs7JhM5r59+/C6B9jDhw9xYz6fTxBEZmbm35ZD3HtFRQXufdOmTbM21gDMErg0CsAs6+zsZLPZFhYWra2tAoHAz8+PJEkej/fjjz9++umnZmZmBw8eRAgNDw9TFLV58+Zdu3adPXt2aGhIq9VGRUU1NDR89tlnYWFho6Ojt27d0mg0CKHx8fHIyMjHjx+fO3du1apVv//+e2ZmZkxMTG1tLYPB0Ol0FEUdPHgwPj6+urp6YGDggw8+SEtLa2trMzc3P3nypEgksra2PnLkCELIzs5ucs7j4+MURe3evXvbtm2HDx9ubW3dvXt3RkZGa2trVlbWihUrfvjhh2PHjvn7++NKplAoQkJCeDxeSUmJs7Pz7du3Dxw4MDIycu7cuQmvjHtnMpk5OTlT9Q6Akb2cFTEAeFW9//77CKG+vr6+vr62trYDBw4ghJKSkmiaTkhIcHJyGhwc1Dfeu3cvl8vF27m5uQihL7/8Un8ULzjwzTffTO7lu+++QwjdunVLH6murkYISaVSmqYfPXqEEIqOjtYfLS4uRgg1NDTg3aCgoKioqGnexZ07dxBCO3bs0EdSU1MRQpcvX8a7Op3Ozc0tOTkZ7+7fv9/Gxka/phJN0wUFBWZmZn19fTRN43Xg7t+/r+/dMDcATA38RwjALCAIAm8wGIyEhIQzZ84ghEiSDAoKqq2t1Tdjs9m9vb3d3d36xbLxFVGssrLS3Nw8JSVl8uuTJMlkMkdGRqqqqnBEq9UuWLCgpaVFvxaP4cpZeFHlrq4uX1/fmb+LiIgI/baHhwdCaMOGDXj3tdde8/Dw6OrqwrsVFRVubm7Nzc3Nzc04YmFhMT4+3tbW5u/vP/MeATAFUAgBmAWVlZUIIR6P5+LiYm1tjRDSaDTPnj27efPm3bt3DVtyOJyenh59IeTxePpDvb29Dg4OZmZ/81upVCqHhobeffddw6CNjY3hFBsOh6PftrCwQAjh+5QzZ/gKCxcunBzR3wXs6enp7++fkA+Hw8Gr7wIwt0AhBGAWhIWFTYgsWrRo4cKFaWlpk2+bGWIw/jNhjc1mq1QqnU5nGMRsbW0Jgujp6ZmthP8lGxsbHx+f69evGzsRAGYBzBoF4IVgMBjBwcHXrl0bHByc4Y+EhIQMDw+XlZVNPhQaGqpSqW7evPnPkmEymdN84+IfCA0NvXPnjlwuN0rvAMwuKIQAvCg5OTnd3d1xcXH37t0bGhqSyWSlpaWZmZlTtU9MTPT29haJREVFRU+fPu3p6RGLxfjKampqqqenZ3JyclFRkVqtpiiqvr5+37597e3tM8lk+fLl9+7dKyoqamhomOGPTO/QoUMMBiMmJqa6uvr58+cKhaKqqmrLli1T9V5fXy8Wi2erdwBmFxRCAF4Uf3//69evy+Xy1atXW1lZLVmyJDk5eXh4eKr2FhYWJEn6+PgkJiZyOBwejycSifA/lFZWVjdu3PDz89u8eTOXyyUIwt/fv66uztLSciaZZGVlrVu3bseOHX5+fnv27Pn3b83V1VUqlZqZmQkEAhaLxefzo6Ojp7pBmJWVJRAIRCKRn5/f3r17/33vAMyu1+i/vt4LAHhBOjo6ent72Wy2k5MTnoQyPYVC8eTJExaL5erqOqE9RVHt7e24rLLZ7BeW8kz9+eefcrmcxWI5OTkZVmU8qdWIiQEwc1AIAQAAzGtwaRQAAMC8BoUQAADAvAaFEAAAwLwGhRAAAMC8BoUQAADAvAaFEAAAwLz2PyueX78u6MQqAAAAAElFTkSuQmCC", "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" ], "text/html": [ "\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" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h = plot(xlabel=\"Percent melt\", ylabel=\"Zircon abundance (wt. %)\")\n", "plot!(h, modes[\"liquid\"], modes[\"zircon\"], label=\"\")\n", "plot!(h,xlims=(0,100),framestyle=:box,fg_color_legend=:white,legend=:right,xflip=true)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Julia 1.9.0-beta4", "language": "julia", "name": "julia-1.9" }, "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "1.9.0" } }, "nbformat": 4, "nbformat_minor": 2 }