Category
Published on
Abstract
A simple Matlab function for calculating the viscosity of magma. It is based on the method of Shaw (1972). Viscosities of magmatic silicate liquids: an empirical method of prediction. American Journal of Science 272, pp. 870-893
It requires the inputs:
SiO2
TiO2
Al2O3
Fe2O3
MnO
MgO
CaO
Na2O
K2O
P2O5
H2O - water content
mTc - magma temperature in celsius
The outputs are:
viscosity_gas = viscosity of the magma with the water content
viscosity_no_gas = viscosity of the magma with no gas
totalweight = total weight of all the inputs
During creation I also used the excel spreadsheet by John D. Winter as a guide (Available at: http://www.whitman.edu/geology/winter/) - which is called "Magma Viscosity spreadsheet". As per my magma density code, the excel spreadsheet is probably more user friendly, however, it is useful to have this file in Matlab. It is also useful to use in combination with modelling work.
If you have any questions please don't hesitate to contact me.
Cite this work
Researchers should cite this work as follows: