Elementelasticconstant(6,6,Nel)

Elementelasticconstant(6,6,Nel) are the elastic constants of each element. The strain tensor is in the form of
${\mathrm{\Sigma }}_{i,j=1}^6{c}_{i,j}$
Nel is the element number. These variables save the elastic constants of each elements.

For example: For wurtzite lattice

elementelasticconstant(1:6,1:6,totalelements) =0.0
do i=1,totalelements
  elementelasticconstant(1,1,i) = c11
  elementelasticconstant(1,2,i) = c12
  elementelasticconstant(1,3,i) = c13
  elementelasticconstant(2,2,i) = c11
  elementelasticconstant(3,3,i) = c33
  elementelasticconstant(2,1,i) = c12
  elementelasticconstant(2,3,i) = c13
  elementelasticconstant(3,1,i) = c13
  elementelasticconstant(3,2,i) = c13
  elementelasticconstant(3,3,i) = c33  
  elementelasticconstant(4,4,i) = c44
  elementelasticconstant(5,5,i) = c44
  elementelasticconstant(6,6,i) = c66
  ...
  ...
  ...
end do