In case of constitutive_pheno.f90, the computation can be optimized. Mattex file slightly modified
Some syntax corrections