We have uploaded Brillouin spectra, XRD data, code for modelling calculation and figure data in this paper.

Brillouin spectra was named as 'original data for velocity-intensity'. You can use the calibrition file (named 'velocity-calibrition') as X axis, and every P-T data can be used as Y axis data to fitting the value of velocity.

XRD data is the original data used to detemine the orientation of each platenet of aragonite. You can download the commercial software in Rigaku to open it.

Code is matlab script can be divided into two parts. At first, we use the elastic data to calculate the bulk modulus, shear modulus and density of single mineral along a geotherm. 
You can input 
K0 (bulk modulus at ambient conditions)
G0(shear modulus at ambient conditions)
rho0(density at ambient conditions)
K'(pressure derivative of bulk modulus)
G'(pressure derivative of shear modulus)
dK/dT(temperature derivative of bulk modulus)
dG/dT(temperature derivative of shear modulus)
alpha (coefficient of thermal expansion)
gama(Grüneisen parameter)
to calculate the K, G and rho at any P-T condition.

Second part, using the obtained elasticity of single mineral, you can calculate the subducted oceanic crust as a multiphase assemblage. 

The calculation is based on  3rd order BM finite strain theory.