1) All the codes are implemented into the FEM package Comsol (version 3.5a). There are totally four, corresponding to four sets of parameters. 
2) All the data are in the two-column '.txt' format. The names of the data files explain the x and y variables. For instance,  file "dis vs Phi" means x = dispalcement and y = porosity.  
3) In each case (Result I - IV), we also include the matlab code (.m) that can produce the contour maps of different quantities (i.e., temperature, porosity, grain size and dominant deformation mechanism).