https://github.com/charlesll/water-speciation-magma
Model for prediction of water speciation in magmatic liquids
Science Score: 13.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
○CITATION.cff file
-
○codemeta.json file
-
○.zenodo.json file
-
✓DOI references
Found 2 DOI reference(s) in README -
○Academic publication links
-
○Committers with academic emails
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (10.3%) to scientific vocabulary
Keywords
Repository
Model for prediction of water speciation in magmatic liquids
Basic Info
Statistics
- Stars: 0
- Watchers: 2
- Forks: 0
- Open Issues: 0
- Releases: 0
Topics
Metadata Files
README.md
water-speciation-magma
Model for prediction of water speciation in magmatic liquids
Citation
This software was part of the publication:
Moretti, R., Le Losq, C., and Neuville, D.R. (2014) The amphoteric behavior of water in silicate melts from the point of view of their ionic-polymeric constitution. Chemical Geology, 367, 23–33. https://doi.org/10.1016/j.chemgeo.2013.12.012
This should be cited in case of use.
Licence
MIT licence, see Licence file.
Contributors
Roberto Moretti (IPGP), moretti@ipgp.fr
Charles Le Losq (IPGP), lelosq@ipgp.fr
Dependencies
A working fortran compiler. We suggest using gfortran, tested on Mac and Linux. It works well with this software!
How to use
Download the repository, and use the provided example input file. It first requires compilation of the FORTRAN source, then running the compilated software.
Compilation
To create the program, with gfortran, in the terminal on Linux or MacOS:
$ gfortran waspecCLL.f90 -o waspec.o
Running the software
The software takes an input file, INPUT.txt, which contains the compositions of interest.
It returns an output file, OUTPUT.txt
Run in the terminal, after compilation, run the command:
$ ./waspec.o
How to use the input file INPUT.txt
This is a comma-separated text file.
1) first line is an integer (number of compositions) followed by another integer, the KYMPA variable (1 or 0); this is explained later on
4) write the composition line in the form of:
T(°C) P(bar) SiO2 TiO2 Al2O3 Fe2O3 Cr2O3 FeO Mno MgO Cao Na2O K2O P2O5 H2Otot H2Omol
NOTES:
- Oxides are given in wt%.
- H2Omol is the molecular water in wt% (otherwise put a number, such as 0.1, in any case lower than H2Otot)
- T(°C) is the temperature of the liquid (greater than or equal to Tg)
About OUTPUT.txt:
1) If you have set to 0 the second integer, right after the number of composition, waspec considers input H2Omol as an experimental constraint. If you have set it to 1, waspec iterates until convergence. Therefore, results can be different between the "0" and "1" KYMPA options.
2) The output lists many parameters:
- H2Ototwt%: total water in wt%, given in input.
- nH2Otot: total moles of water per 100g of melt (inout value).
- nhtot: total moles of hydrogen per 100g of melt (input value).
- nH2Omolwt%EXP: molecular water in wt% (input value).
- nH2OmolEXP: moles of molecular water per 100g of melt (input value).
- 1/T: the reciprocal of temperature (in Kelvin).
- nfreeOH: calculated moles of OH- per 100g of melt.
- nfreeH: calculated moles of H+ per 100g of melt.
- nOH-IR: the sum of calculated moles of H+ and OH- per 100g of melt.
- nH2OmolEXP: moles of molecular water per 100g of melt (input value) (same as above).
- nH2Omolcalc: calculated moles of molecular water per 100g of melt.
- nHtot_after: calculated moles of total hydrogen per 100g of melt (CAN BE DIFFERENT THAN nhtot, IF KYMPA = 0)
- DISPERSION: the difference between nhtot_after and nhtot
- XfreeOHcalc: mol fraction of calculated free OH
- XfreeHcalc: mol fraction of calculated H
- XH2Omolcalc: mol fraction of calculated H2Omol
- YfreeOHcalc: relative proportion of calculated free OH
- YfreeHcalc: relative proportion of calculated free H
- YH2Omolcalc: relative proportion of calculated H2Omol
Owner
- Name: Charles Le Losq
- Login: charlesll
- Kind: user
- Location: France
- Company: Institut de physique du globe de Paris
- Website: https://charlesll.github.io/
- Twitter: Charles_LeLosq
- Repositories: 49
- Profile: https://github.com/charlesll
Geoscientist, Assistant Professor, willing to improve data reduction protocols through programming and data science
GitHub Events
Total
Last Year
Committers
Last synced: about 2 years ago
Top Committers
| Name | Commits | |
|---|---|---|
| Charles Le Losq | l****q@i****r | 6 |
| Charles Le Losq | c****q@g****m | 1 |
Committer Domains (Top 20 + Academic)
Issues and Pull Requests
Last synced: 11 months ago
All Time
- Total issues: 0
- Total pull requests: 0
- Average time to close issues: N/A
- Average time to close pull requests: N/A
- Total issue authors: 0
- Total pull request authors: 0
- Average comments per issue: 0
- Average comments per pull request: 0
- Merged pull requests: 0
- Bot issues: 0
- Bot pull requests: 0
Past Year
- Issues: 0
- Pull requests: 0
- Average time to close issues: N/A
- Average time to close pull requests: N/A
- Issue authors: 0
- Pull request authors: 0
- Average comments per issue: 0
- Average comments per pull request: 0
- Merged pull requests: 0
- Bot issues: 0
- Bot pull requests: 0