From 9b3ddcd2c2373493abc46fb1152e37063d5d255a Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 28 Aug 2018 06:27:38 +0200 Subject: [PATCH] more explicit --- src/homogenization_isostrain.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/homogenization_isostrain.f90 b/src/homogenization_isostrain.f90 index af16aaddd..3ac64e6da 100644 --- a/src/homogenization_isostrain.f90 +++ b/src/homogenization_isostrain.f90 @@ -137,8 +137,9 @@ subroutine homogenization_isostrain_partitionDeformation(F,avgF,el) instance = homogenization_typeInstance(mesh_element(3,el)) associate(prm => param(instance)) - F = 0.0_pReal F(1:3,1:3,1:prm%Nconstituents) = spread(avgF,3,prm%Nconstituents) + if (homogenization_maxNgrains > prm%Nconstituents) & + F(1:3,1:3,prm%Nconstituents+1_pInt:homogenization_maxNgrains) = 0.0_pReal end associate end subroutine homogenization_isostrain_partitionDeformation