Note

Go to the end to download the full example code.

MagmaSat: Saturation Pressures#

Converts measured CO2-H2O concentrations in a glass (in a melt inclusion, for example) to a fluid-saturation pressure, given the composition of the glass.

This uses the VESIcal python library and is a preview of a forthcoming version of VESIcal compatible with the new version of thermoengine.

Open this code in an executable MyBinder instance (MyBinder links may be slow to load– please be patient!):

https://mybinder.org/badge_logo.svg

To install the preview version of VESIcal that works with the new pip installable version of thermoengine, run this command (uncomment first):

%pip install kaylai/VESIcal@thermoenginelite

import VESIcal as v

A single glass composition#

Specify the composition of the liquid, including the measured H2O and CO2 concentrations. Here a basaltic composition is given as an example:

oxide_comp = {
        'SiO2':49.91,
        'TiO2': 1.47,
        'Al2O3': 17.91,
        'Fe2O3': 2.45,
        'FeO':  7.02,
        'MnO':  0.16,
        'MgO':  6.62,
        'CaO':  10.02,
        'Na2O':  3.02,
        'K2O':  0.64,
        'P2O5':  0.2,
        'H2O': 1.0,
        'CO2': 0.3
        }

Now convert this to a VESIcal Sample object (this allows the code to deal with unit conversions, normalisation, etc)

sample = v.Sample(oxide_comp)

Now calculate the fluid-saturation pressure:

v.calculate_saturation_pressure(sample=sample, temperature=1200.0).result

/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

5050.0

Run a batch calculation#

A more common use case is converting a suite of melt inclusion compositions into a distribution of saturation pressures. This is exactly what VESIcal is designed to do!

First, import a file containing the compositions. This has to be formatted in the way that VESIcal expects, you can use the example (from the VESIcal manuscript) here as a template:

myfile = v.BatchFile('input_files/magmasat_example_data.xlsx')

myfile.data.head()

/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/batchfile.py:230: FutureWarning: ChainedAssignmentError: behaviour will change in pandas 3.0!
You are setting values through chained assignment. Currently this works in certain cases, but when using Copy-on-Write (which will become the default behaviour in pandas 3.0) this will never work to update the original DataFrame or Series, because the intermediate object on which we are setting values will behave as a copy.
A typical example is when you are setting values in a column of a DataFrame, like:

df["col"][row_indexer] = value

Use `df.loc[row_indexer, "col"] = values` instead, to perform the assignment in a single step and ensure this keeps updating the original `df`.

See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy

  data[column][data[column] < 0] = 0
SiO2 TiO2 Al2O3 Fe2O3 Cr2O3 FeO MnO MgO NiO CoO CaO Na2O K2O P2O5 H2O CO2 Press Temp
BT-ex 77.50 0.0800 12.50 0.207 0 0.4730 0.0000 0.0300 0 0 0.4300 3.98 4.88 0.0 5.500000 0.050000 500 900
TVZMa-ex 78.37 0.1300 11.94 0.000 0 0.9900 0.0400 0.0500 0 0 0.5300 3.80 4.14 0.0 4.060000 0.005000 600 800
TVZOh-ex 77.90 0.0800 12.15 0.000 0 0.9500 0.0500 0.0600 0 0 0.5500 4.05 4.12 0.0 4.630000 0.005000 50 900
Oh48-FTIR1-MI1-a 78.27 0.0298 12.02 0.000 0 0.9828 0.0336 0.0515 0 0 0.4772 4.05 4.09 0.0 4.214912 0.004566 250 950
Oh48-FTIR1-MI1-b 78.27 0.0298 12.02 0.000 0 0.9828 0.0336 0.0515 0 0 0.4772 4.05 4.09 0.0 4.005816 0.004448 500 1025


Now run the saturation pressure calculation. Notice that the temperature is specified here (but could be specified in the excel file- see VESIcal docs for more information).

satPs = myfile.calculate_saturation_pressure(temperature=1200.0)

[=                   ] 7%  Working on sample BT-ex                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[==                  ] 14%  Working on sample TVZMa-ex                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[====                ] 21%  Working on sample TVZOh-ex                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[=====               ] 29%  Working on sample Oh48-FTIR1-MI1-a                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[=======             ] 36%  Working on sample Oh48-FTIR1-MI1-b                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[========            ] 43%  Working on sample Oh48-FTIR1-MI1-IRc                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[==========          ] 50%  Working on sample Oh50-4.1                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[===========         ] 57%  Working on sample Oh50-4.2                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[============        ] 64%  Working on sample Oh49-4.1                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[==============      ] 71%  Working on sample Oh49-4.2                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[===============     ] 79%  Working on sample Ma55-5a.1                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[=================   ] 86%  Working on sample Ma57-3b.2                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[==================  ] 93%  Working on sample Ma57-3c.1                            /home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

[====================] 100%  Working on sample Ma57-3c.2
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()
/home/vscode/.local/lib/python3.10/site-packages/VESIcal/models/magmasat.py:208: RuntimeWarning: invalid value encountered in divide
  fluid_comp = fluid_comp/fluid_comp.sum()

We can look at the calculated Saturation pressures (in bar):

satPs.SaturationP_bars_VESIcal

BT-ex                 2460.0
TVZMa-ex              1120.0
TVZOh-ex              1370.0
Oh48-FTIR1-MI1-a      1180.0
Oh48-FTIR1-MI1-b      1090.0
Oh48-FTIR1-MI1-IRc    1040.0
Oh50-4.1              1370.0
Oh50-4.2              1260.0
Oh49-4.1              1210.0
Oh49-4.2              1190.0
Ma55-5a.1             1300.0
Ma57-3b.2              970.0
Ma57-3c.1             1290.0
Ma57-3c.2             1220.0
Name: SaturationP_bars_VESIcal, dtype: float64

And save them to a csv or excel file:

satPs.to_csv('my_saturation_pressures.csv')

We can make a plot of the saturation pressure distribution using matplotlib, which we should first import:

import matplotlib.pyplot as plt

Here’s some example code for making a plot:

fig, ax = plt.subplots(1,2, figsize=(6,3), sharey='row')

ax[0].scatter(satPs.H2O, satPs.SaturationP_bars_VESIcal,
              lw=1, edgecolor='k')
ax[1].scatter(satPs.CO2, satPs.SaturationP_bars_VESIcal,
              lw=1, edgecolor='k')

fig.tight_layout()

ax[0].set_xlabel('H$_2$O (wt%)',fontsize=12)
ax[1].set_xlabel('CO$_2$ (wt%)',fontsize=12)
ax[0].set_ylabel('Saturation Pressure (bar)', fontsize=12)

for a in ax:
    a.tick_params(axis='both', which='both', top=True, right=True, direction='out')

plt.show
plot magmasat saturation pressures
<function show at 0xffff934f1090>

Total running time of the script: (0 minutes 15.097 seconds)

Gallery generated by Sphinx-Gallery