From 7a84209a14aba85770c6e4a9f711c125e4d84763 Mon Sep 17 00:00:00 2001
From: Pratheek Shanthraj
Date: Thu, 24 Jan 2013 13:12:51 +0000
Subject: [PATCH] bug fix in microstructure -This line, and those below, will
be ignored--
M constitutive_dislotwin.f90
---
code/constitutive_dislotwin.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/code/constitutive_dislotwin.f90 b/code/constitutive_dislotwin.f90
index 48e089f46..1bb400a90 100644
--- a/code/constitutive_dislotwin.f90
+++ b/code/constitutive_dislotwin.f90
@@ -856,7 +856,7 @@ forall (s = 1_pInt:ns) &
!* 1/mean free distance between 2 twin stacks from different systems seen by a moving dislocation
!$OMP CRITICAL (evilmatmul)
state(g,ip,el)%p((3_pInt*ns+nt+1_pInt):(4_pInt*ns+nt)) = 0.0_pReal
-if (nt > 0_pInt) &
+if (nt > 0_pInt .and. ns > 0_pInt) &
state(g,ip,el)%p((3_pInt*ns+nt+1):(4_pInt*ns+nt)) = &
matmul(constitutive_dislotwin_interactionMatrix_SlipTwin(1:ns,1:nt,myInstance),fOverStacksize(1:nt))/(1.0_pReal-sumf)
!$OMP END CRITICAL (evilmatmul)