Hello ASPECT community!
I’m working on melting within lithospheric delaminations in the East African Rift and wanted to ask if there is currently a way to implement two different melting functions for different materials/ compositions, for a 2D visco plastic model? I am looking for lithospheric drips in which melt fraction increases as a heavy body ‘drip’ descends using modifications of Katz, 2003 melting parameterization. Getting an increasing melt signature on a downgoing body has proven difficult and I think that including heterogeneity, i.e. a pyroxenite body overlying a peridotite body, may increase melt fraction as the drip falls and pulls the easier-to-melt pyroxenite material down.