montesi/MeltMigrator: Initial release of MeltMigrator

Initial release of 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.

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., L.G.J. Montési, and M.D. Behn, MeltMigrator: a MATLAB-based Software for Modeling Three-dimensional Melt Migration and Crustal Thickness Variations at Mid-Ocean Ridges Following a Rules-Based Approach, submitted to Geochemistry, Geophysics, Geosystems, 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