COMAGMAT 5
Calculations of equilibrium and fractional crystallization of Ssaturated and Sundersaturated magmas, including changes in the Fe/Ni ratio in silicate melts, femic minerals, and coexisting sulfides, as well as sulfidesilicate (±FeTi oxides) proportions for multiplesaturated mineral assemblages.
COMAGMAT 5
Calculations of equilibrium and fractional crystallization of Ssaturated and Sundersaturated magmas, including changes in the Fe/Ni ratio in silicate melts, femic minerals, and coexisting sulfides, as well as sulfidesilicate (±FeTi oxides) proportions for multiplesaturated mineral assemblages.
COMAGMAT 3
Calculations of equilibrium and fractional crystallization for dry and hydrous natural magmas in the range of pressures from 1 atm to 1012 kbar and including both open (12 oxygen buffers) and closed system fractionation with respect to oxygen. Simulating formation of layered intrusions.
Features
COMAGMAT 5.x 
COMAGMAT 3.x 


Calculations of equilibrium and fractional crystallization of Ssaturated and Sundersaturated magmas  
Improved precision of calculations at low contents of melt components  
Variable models for melt oxidation state (Fe^{2+}/Fe^{3+} ratio)  
Calculations of trace elements partition coefficients  
Modeling crystallization at atmospheric pressure  
Modeling crystallization at elevated pressures (up to 1012 kbar)  
Modeling crystallization in hydrous systems  
Simulation of open (with respect to oxygen) systems, using fO2buffers  
Simulation of closed to oxygen systems at given Fe^{2+}/Fe^{3+} ratio  
Simulating formation of layered intrusions  
Correction for mineralmelt liquidus temperatures 
COMAGMAT 5 New
The basic core of the COMAGMAT5 program includes the algorithm that has been used in previous versions of the COMAGMAT model (ver. 3.03.5) designed to simulate mafic magma crystallization processes, but limited to atmospheric pressure and include numerical models for simulating crystallization of Ssaturated and Sundersaturated magmas.
Sulfide version of the COMAGMAT program was developed as a part of a cooperative AMIRA project “NiPGE potential of mafic and ultramafic magmas – a combined melt inclusion and numerical modelling approach” (P962) and Russian Science Foundation project (No 161710129).
Download
Recent changes
 ver. 5.2.2.1  fix import of input files with "comma" decimal separator
 ver. 5.2.2  "resizable" user interface, fix minor issues on Windows 10
 ver. 5.2.1  FeNiS sulfide solubility/stability model
 ver. 5.2.0  FeS sulfide solubility/stability model, silicate minerals recalibration
Compatibility
 Tested on Windows XP, Windows 7 and Windows 10 operation systems.
References
 Ariskin A.A., Danyushevsky L.V., Bychkov K.A., McNeill A.W., Barmina G.S. & Nikolaev G.S. (2013). Modeling Solubility of FeNi Sulfides in Basaltic Magmas: The Effect of Nickel. Economic Geology 108, 1983–2003.
 Ariskin A.A., Bychkov K.A., Nikolaev G.S. & Barmina G.S. (2018). The COMAGMAT5: Modeling the Effect of FeNi Sulfide Immiscibility in Crystallizing Magmas and Cumulates. Journal of Petrology. V. 59 (2). 283298.
Applications
 Ariskin A.A., Danyushevsky L.V., Nikolaev G.S., Kislov E.V., Fiorentini M.L., McNeill A.W., Kostitsyn Y., Goemann K., Feig S. & Malyshev A. (2018) The Dovyren Intrusive Complex (Southern Siberia, Russia): Insights into dynamics of an open magma chamber with implications for parental magma origin, composition, and CuNiPGE fertility. Lithos. 302. 10.1016/j.lithos.2018.01.001.
 Ariskin A.A., Bychkov K.A. & Nikolaev G.S. (2017). Modeling of traceelement composition of sulfide liquid in a crystallizing basalt magma: Development of the Rfactor concept. Geochemistry International 55, 465–473.
 Gongalsky B.I., Krivolutskaya N.A., Ariskin A.A. & Nikolaev G.S. (2016). The Chineysky gabbronoriteanorthosite layered massif (NorthernTransbaikalia, Russia): its structure, FeTiV and CuPGE deposits, and parental magma composition. Mineralium Deposita. Mineralium Deposita 51, 1013–1034.
 Ariskin A.A., Kislov E.V., Danyushevsky L.V., Nikolaev G.S., Fiorentini M.L., Gilbert S., Goemann K. & Malyshev A. (2016). Cu–Ni–PGE fertility of the YokoDovyren layered massif (northern Transbaikalia, Russia): thermodynamic modeling of sulfide compositions in low mineralized dunite based on quantitative sulfide mineralogy. Mineralium Deposita. Mineralium Deposita 51, 993–1011.
COMAGMAT 3
The COMAGMAT model is a programs developed to calculate phase equilibria for dry and hydrous natural magmas crystallizing in the range of pressures from 1 atm to 1012 kbar and including both open (12 oxygen buffers) and closed system fractionation with respect to oxygen. The modeling process may be calculated for systems ranging from basalts to dacites, with modeled major elements including Si  Ti  Al  Fe  Mg  Ca  Na K and P. Moreover, COMAGMAT allows the user to simulate behaviour of 20 trace elements, including Mn, Ni, Co, Cr, V, Sc, Sr, Ba, Rb, Cu, and REE. The modeled minerals include olivine (FoFa solution), plagioclase (AnAb), 3 pyroxenes (augite, pigeonite, and orthopyroxene: EnFsWo solutions plus Al and Ti), ilmenite (IlmHem), and magnetite (MtUlv).
Download
Recent changes
 ver. 3.75  Fix result output for "Simulating formation of layered intrusions", fix Excel report generation.
 ver. 3.74  Fix intermediate output for Pigeonite/Opx
 ver. 3.73  Fix performance issues
 ver. 3.72  Update water effect on liquidus temperature
 ver. 3.57  First release similar to 3.52 with Ilm subroutine changed to escape some computative problems
 ver. 3.65  Specially calibrated to simulate crystallization of FeTi enriched basaltic liquids
 ver. 3.52  Ilmenite code was changed to correctly calculate Mg contents in Ilm; K2O in Pl was also corrected.
 ver. 3.50  More accurate Mtmelt and Ilmmodels were integrated into the program which now allow one to calculate crystallization at watersaturated conditions starting with 0.5 wt% H2O in the melt. Other improvements include corrections of mineralmelt distribution coefficients for TiO2 in clino and orthopyroxenes, as well K2O in Pl. In addition, OXIDES file is now printed out to see Mt and Ilm compositions calculated on single cation basis.
 ver. 3.30  This is INTERNET analogue of the COMAGMAT3.0 release.
 ver. 3.00  This is a basic DOSversion distributed since 1992. See publications for more details.
Compatibility
 Tested on Windows XP, Windows 7 and Windows 10 operation systems.
References
 Ariskin A.A. & Barmina G.S. (2004). COMAGMAT: Development of a magma crystallization model and its petrological applications. Geochemistry International 42, s1–s157.
 Almeev, R. R., Holtz, F., Koepke, J. & Parat, F. (2012). Experimental calibration of the effect of H2O on plagioclase crystallization in basaltic melt at 200 MPa. American Mineralogist 97 (7), 1234–1240.
 Almeev R.R., Holtz F., Koepke J., Parat F. & Botcharnikov R. E. (2007). The effect of H2O on olivine crystallization in MORB: Experimental calibration at 200 MPa. American Mineralogist 92, 670–674.
 Ariskin A.A. (1999). Phase equilibria modeling in igneous petrology: use of COMAGMAT model for simulating fractionation of ferrobasaltic magmas and the genesis of highalumina basalt. Journal of Volcanology and Geothermal Research 90, 115–162.
 Ariskin A.A. & Barmina G.S. (1999). An empirical model for the calculation of spinelmelt equilibria in mafic igneous systems at atmospheric pressure: 2. FeTi oxides. Contributions to Mineralogy and Petrology 134, 251–263.
 Ariskin, A. A., Frenkel, M. Y., Barmina, G. S. & Nielsen, R. L. (1993). Comagmat: a Fortran program to model magma differentiation processes. Computers and Geosciences 19, 1155–1170.
Applications
 Ariskin A.A. (2003). The compositional evolution of differentiated liquids from the Skaergaard Layered Series as determined by geochemical thermometry. Russian Journal of Earth Sciences 5, 1–29.
 Ariskin A.A. (2002). Geochemical thermometry of the layered series rocks of the Skaergaard intrusion. Petrology 10, 495–518.