this depends a bit on how realistic you want your model to be. We’ve made a model with melt transport and a simple stress-limiter rheology for a manuscript that is currently in review, you can find the manuscript here: https://arxiv.org/abs/1810.10105 (the model is described in Section 3.3)
and a description of how to reproduce the models here:
However, this does not include hydrofracturing, so there is no influence of the fluid pressure on the brittle failure mechanism. For a mid-ocean ridge system, this can be a good approximation, but it depends on the process you are interested in and may be more important for other applications. I’ve started to look at implementing hydrofracturing with @jbnaliboff, but this is technically challenging, and I am not sure when we will have this done. If you’re interested in the details, the method we want to implement is described in Keller et al., 2013 (https://academic.oup.com/gji/article/195/3/1406/2874184), and we can discuss this in more detail if you think this is what you need.
On the other hand, if you think that for your application you don’t need hydrofracturing, the model we use in the manuscript could ba a good starting point.