how to suppress mineral in py for 2.3.1 and find liquids temperature

[Kai Wen]

Hello,

I tried to use 2.3.1 alphamelts-py to stimulate magmatic processes. How could I set the suppress phase?

Moreover, in the app, we could find liquid temperature for wet or not, but we could not find that function in Python. Does that influence the calculation?

Many thanks.

Best,
Kai

[Kai Wen]

Moreover, sorry, in the app, it seems to adjust my content automatically, if I choose find wet liquids

find wet liquids
Initial alphaMELTS calculation at: P 4000.000000 (bars), T 1113.085938 (C)
liquid:    SiO2 TiO2 Al2O3 Fe2O3 Cr2O3  FeO  MnO  MgO  NiO  CoO  CaO Na2O  K2O P2O5  H2O
98.8228 g 59.87 0.63 17.73  0.98  0.00 4.13 0.13 2.91 0.00 0.00 7.08 2.86 1.67 0.00 2.02
Activity of H2O = 0.21086  Melt fraction = 0.994465
clinopyroxene: 0.550002 g, composition Na0.02Ca0.73Fe''0.23Mg0.90Fe'''0.05Ti0.01Al0.20Si1.88O6

find liquids
Initial alphaMELTS calculation at: P 4000.000000 (bars), T 1113.281250 (C)
liquid:    SiO2 TiO2 Al2O3 Fe2O3 Cr2O3  FeO  MnO  MgO  NiO  CoO  CaO Na2O  K2O P2O5  H2O
99.3730 g 59.82 0.63 17.66  0.99  0.00 4.14 0.13 2.98 0.00 0.00 7.14 2.84 1.66 0.00 2.01
Activity of H2O = 0.209584  Melt fraction = 1

not same

[Paula]

Hi Kai,

To suppress a phase you add a .melts file-like line in the Settings. Something like:

Code Select Expand

ptpath.engine.setSystemProperties(["Suppress: olivine"])
The Find Wet Liquidus result differs slightly because it performs one or more equilibration calculations as part of the algorithm, whereas the original Find Liquidus routine only uses affinities. In this case, the temperatures agree to within 0.2 ^oC because the system is not water-saturated.

The MATLAB and Python versions always use Find Wet Liquidus. The app has both options for backwards compatibility and so it can be more easily compared to other software, such as easyMelts.

Paula