1. phasemod(nphase,2) (integer*4): An array which contains the Young's modulus, phasemod(i, 1), and Poisson ratio, phasemod(i, 2), of the ith> phase. This array is now initialized in its own subroutine called phasemod_init instead of within MAIN.
2. eigen(nphase,6) (double precision): Thermal strains of each phase: eigen (nphase, 1) = xx, eigen(nphase, 2) = yy, eigen(nphase, 3) = zz, eigen(nphase, 4) = xz, eigen(nphase, 5) = yz, eigen(nphase, 6) = xy.