In case of constitutive_pheno.f90, the computation can be optimized. Mattex file slightly modified
The dislocation based constitutive subroutine can be only use for fcc crystal !!!! The constitutive equations are different for bcc and hcp.
Some syntax corrections