From 88376568707bb07eaf13099cddc788b5b0d4ab4f Mon Sep 17 00:00:00 2001 From: Zhuowen Zhao Date: Tue, 9 Jan 2018 15:17:51 -0500 Subject: [PATCH] added "plastic_kinematichardeing.f90" to commercialFEM_fileList Changed outputID type --- src/commercialFEM_fileList.f90 | 1 + src/plastic_kinematichardening.f90 | 4 ++-- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/src/commercialFEM_fileList.f90 b/src/commercialFEM_fileList.f90 index 51848ece5..f57f03467 100644 --- a/src/commercialFEM_fileList.f90 +++ b/src/commercialFEM_fileList.f90 @@ -28,6 +28,7 @@ #include "plastic_none.f90" #include "plastic_isotropic.f90" #include "plastic_phenopowerlaw.f90" +#include "plastic_kinematichardening.f90" #include "plastic_dislotwin.f90" #include "plastic_disloUCLA.f90" #include "plastic_nonlocal.f90" diff --git a/src/plastic_kinematichardening.f90 b/src/plastic_kinematichardening.f90 index 3bd35aaea..04ac60b9e 100644 --- a/src/plastic_kinematichardening.f90 +++ b/src/plastic_kinematichardening.f90 @@ -226,8 +226,8 @@ subroutine plastic_kinehardening_init(fileUnit) Nchunks_SlipFamilies = count(lattice_NslipSystem(:,phase) > 0_pInt) ! maximum number of slip families according to lattice type of current phase Nchunks_SlipSlip = maxval(lattice_interactionSlipSlip(:,:,phase)) Nchunks_nonSchmid = lattice_NnonSchmid(phase) - allocate(param(instance)%outputID(phase_Noutput(phase)), source=0_pInt) ! allocate space for IDs of every requested output - allocate(param(instance)%crss0 (Nchunks_SlipFamilies), source=0.0_pReal) + allocate(param(instance)%outputID(phase_Noutput(phase)), source=0) ! allocate space for IDs of every requested output + allocate(param(instance)%crss0 (Nchunks_SlipFamilies), source=0.0_pReal) allocate(param(instance)%tau1 (Nchunks_SlipFamilies), source=0.0_pReal) allocate(param(instance)%tau1_b (Nchunks_SlipFamilies), source=0.0_pReal) allocate(param(instance)%theta0 (Nchunks_SlipFamilies), source=0.0_pReal)