Algorithmic modifications to the MELTS software package are presented in order that calculations of heterogeneous phase equilibria can be performed in the subsolidus. Methods are presented for: (1) selecting an ‘‘initial guess assemblage’’ that satisfies the bulk composition constraints; (2) detecting saturation of new phases (including liquid) in an assemblage; (3) adding and removing phases f...
