Elementlatticeconstant(Nel,3)

Elementlatticeconstant(Nel,1:3) are the the variables saving the lattice constants.

   do j=1,totalelements!mshel%tetraelenumstart,mshel%n
       elementlatticeconstant(j,1)=a_x   ! lattice constant of x-axis
       elementlatticeconstant(j,2)=a_y   ! lattice constant of y-axis
       elementlatticeconstant(j,3)=c_z   ! lattice constant of z-axis
      
   end do