Skip to content
/ VESIcal Public

A generalized python library for calculating and plotting various things related to mixed volatile (H2O-CO2) solubility in silicate melts.


Notifications You must be signed in to change notification settings


Repository files navigation


A generalized python library for calculating and plotting various things related to mixed volatile (H2O-CO2) solubility in silicate melts.

Documentation Status DOI


Check here first for all your VESIcal questions! And be sure to read the manuscripts.

Interactive versions of manuscripts

  • Direct link to interactive VESIcal Part I manuscript: Manuscript on Binder
  • Jupyter Notebook hub with VESIcal: Manuscript on Binder

PDF versions of manuscripts

Installation and online use

In general, VESIcal can be accessed and used in a variety of ways. From most flexible (advanced) to least flexible (novice), these are:

VESIcal on the ENKI server

VESIcal requires installation of not only the VESIcal library but also some other python libraries, one of which is a bit tricky to install (ENKI/thermoengine aka the engine behind MELTS). But, we have a solution! All dependencies and the latest version of VESIcal are all installed on the ENKI server, within a Jupyter Notebook Hub. Steps to use VESIcal on the ENKI server are:

  1. Create a (free) GitLab account, which you'll use to sign into ENKI here:
  2. Email ENKI PI Mark Ghiorso at with your GitLab username and requet access to the ENKI server.
  3. Access the ENKI Production Server by going to and clicking "SERVERS" > "PRODUCTION SERVER"
  4. Sign in with your GitLab credentials: You are now in your own jupyter notebook workspace! You can upload and create files here. They won't be accessible to anyone else.
  5. Click the green "CLOSE THIS SCREEN" button
  6. Create a new notebook by clicking the blue plus button and then selecting Python3 under Notebook. Or select from the menu File > New > Notebook
  7. Be sure to import VESIcal as v at the top of your file, and now you are ready to get to work!

See video tutorials on our ReadTheDocs page for more.

Installing locally

Important! Thermoengine must be installed! Please see below for details on how to install thermoengine, the python implementation of MELTS/MagmaSat

First, obtain Python3.x if you do not already have it installed. If you are new to python, we recommend installing it via anaconda3. VESIcal can be installed with one line. Open a terminal and type the following:

pip install VESIcal

Check that the installation worked by entering the following lines into a terminal:

import VESIcal as v

If no output is returned, VESIcal has installed properly! You will very likely, however, see a warning telling you that no module named 'thermoengine' could be found. The installation you performed via pip attempts to install all dependencies (other libraries that VESIcal requires), but thermoengine is not available via pip and so must be manually installed.

Dependencies that should automatically be installed for you are:

  • pandas
  • numpy
  • matplotlib
  • cycler
  • abc
  • scipy
  • sys
  • sympy
  • copy

If any warnings related to these libraries appear, try installing them as you did VESIcal: with 'pip install [package]'.

Installing thermoengine

Thermoengine is the ENKI implementation of MELTS (MagmaSat), which is the default solubility model implemented in VESIcal. VESIcal cannot be run without thermoengine at this time, however a VESIcal-lite that does not include MagmaSat is planned. To install thermoengine, please refer to the ENKI documentation here:

In almost all cases you will need to install thermoengine using docker. The thermoengine devs have kindly put together a docker image for you. We suggest you follow those instructions here:


To upgrade to the most recent version of VESIcal, type the following into terminal:

pip install VESIcal --upgrade


Issues are tracked on GitHub.

Patches may be submitted via a Github pull request. All changes should include tests (VESIcal uses python's unittest library) and pass flake8.


A generalized python library for calculating and plotting various things related to mixed volatile (H2O-CO2) solubility in silicate melts.