montesi/MeltMigrator 1.1.0

MATLAB codes that simulate melt migration out of a planetary mantle according to a simplified 3-step process described in Montési et al. (2011). The code itself is described at length in Bai et al. (2016)

MeltMigrator is based on MATLAB and functions for examples with Matlab2015b. It uses the fsolve function, which is available with MATLAB's optimization toolbox. In version 1.0.0, MeltMigrator assumes the availability of a mantle flow and temperature model build with COMSOL Multiphysics and was tested with version 5.0 and 5.1.

You should modify SetParameters.m to enter your preferred parameters and model descriptions

Cartoon diagram describing the three steps of melt migration used by MeltMigrator. See also Gregg et al. (2012) for a slightly different version of the diagram.

New in version 1.1.0:

• The input model no longer need to be COMSOL. You can choose instead to use a ASCII table containing the x,y,z positions and associated temperature and upward velocity.
• all calls to fsolve have been replaced by calls to fzero, which does not require a custom toolbox.

References:

• Montési, L.G.J., M.D. Behn, L.B. Hebert, J. Lin, and J.L, Barry, (2011), The importance of plate-driven flow in generating crustal variations along the Southwest Indian Ridge 10°-16°E, Journal of Geophysical Research, 116, B10102, doi:10.1029/2011JB008259
• Bai, H., Montési, L. G. J. and Behn, M. D. (2016), MeltMigrator: A MATLAB-based software for modeling three-dimensional melt migration and crustal thickness variations at mid-ocean ridges following a rules-based approach. Geochem. Geophys. Geosyst.. Accepted Author Manuscript. doi:10.1002/2016GC006686
• Gregg, P.M., L.B. Hebert, L.G.J. Montési, and R.F. Katz, (2012), Geodynamic models of melt generation and extraction at mid-ocean ridges, Oceanography, 25(1): 78-88, doi:10.5670/oceanog.2012.05