%************************************************
% Define elastic stiffness tensor (GPa)
%************************************************
%
% Reference Elastic constants
% Antigorite Bezacier et al.20010 EPSL doi:10.1016/j.eps X=a* Y=b Z=c
% Earth and Planetary Science Letters 289, 198–208.
%
% Reference Crystal Structure
%  M.Mellini 2004 Am.Min.vol.89,pp147-158 (SG P1m1)
%
% Define density (g/cm3)
 rho= 2.6200;
%
% Define Cartesian Tensor crystal symmetry & frame
cs_Tensor = crystalSymmetry('2/m',[ 43.5050  9.2510  7.2630],...
 [  90.0000  91.3200  90.0000]*degree,'x||a*','z||c',...
 'mineral','Antigorite 2010 a* c');
%
% elastic Cij stiffness tensor (GPa) as matrix M
 M =.... 
 [[  208.40   66.20   15.90    0.00    2.40    0.00];...
 [    66.20  201.60    5.00    0.00   -4.40    0.00];...
 [    15.90    5.00   96.70    0.00    2.50    0.00];...
 [     0.00    0.00    0.00   17.40    0.00  -13.10];...
 [     2.40   -4.40    2.50    0.00   18.30    0.00];...
 [     0.00    0.00    0.00  -13.10    0.00   65.00]];
%
% M as stiffness tensor C with MTEX tensor command
C = tensor(M,cs_Tensor);
%
