%************************************************
% Define elastic stiffness tensor (GPa)
%************************************************
%
% Reference Elastic constants
% Omphacite Bhagat et al.1992 JGR v.97,6843-6848.
%
% Reference Crystal Structure
% Di38 Jd62 space group C2/c XIIa*,YIIb,ZIIc d=3.327g/cm3
%
% Define density (g/cm3)
rho= 3.3270;
%
% Define Cartesian Tensor crystal symmetry & frame
cs_Tensor = crystalSymmetry('2/m',[  9.5388  8.6976  5.2504],...
 [  90.0000 106.0660  90.0000]*degree,'x||a*','z||c',...
 'mineral','Omphacite');
%
% elastic Cij stiffness tensor (GPa) as matrix M
 M =.... 
 [[  257.30   85.90   76.20    0.00    7.10    0.00];...
 [    85.90  216.20   71.80    0.00   13.30    0.00];...
 [    76.20   71.80  260.20    0.00   33.70    0.00];...
 [     0.00    0.00    0.00   80.20    0.00   10.20];...
 [     7.10   13.30   33.70    0.00   70.60    0.00];...
 [     0.00    0.00    0.00   10.20    0.00   85.80]];
%
% M as stiffness tensor C with MTEX tensor command
C = tensor(M,cs_Tensor);
%
