From 9249b7f4af0b93d2e54f4feacb91a1e80afddbc1 Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Tue, 14 May 2019 19:52:48 +0200 Subject: [PATCH] improved column alignment of indices such as (g,i,e) --- src/homogenization.f90 | 32 ++++++++++++++++---------------- 1 file changed, 16 insertions(+), 16 deletions(-) diff --git a/src/homogenization.f90 b/src/homogenization.f90 index 831e62718..70b3ff0ab 100644 --- a/src/homogenization.f90 +++ b/src/homogenization.f90 @@ -456,45 +456,45 @@ subroutine materialpoint_stressAndItsTangent(updateJaco,dt) steppingNeeded: if (materialpoint_subStep(i,e) > subStepMinHomog) then ! wind forward grain starting point of... - crystallite_partionedF0(1:3,1:3,1:myNgrains,i,e) = & + crystallite_partionedF0 (1:3,1:3,1:myNgrains,i,e) = & crystallite_partionedF(1:3,1:3,1:myNgrains,i,e) ! ...def grads - crystallite_partionedFp0(1:3,1:3,1:myNgrains,i,e) = & - crystallite_Fp(1:3,1:3,1:myNgrains,i,e) ! ...plastic def grads + crystallite_partionedFp0 (1:3,1:3,1:myNgrains,i,e) = & + crystallite_Fp (1:3,1:3,1:myNgrains,i,e) ! ...plastic def grads - crystallite_partionedLp0(1:3,1:3,1:myNgrains,i,e) = & - crystallite_Lp(1:3,1:3,1:myNgrains,i,e) ! ...plastic velocity grads + crystallite_partionedLp0 (1:3,1:3,1:myNgrains,i,e) = & + crystallite_Lp (1:3,1:3,1:myNgrains,i,e) ! ...plastic velocity grads - crystallite_partionedFi0(1:3,1:3,1:myNgrains,i,e) = & - crystallite_Fi(1:3,1:3,1:myNgrains,i,e) ! ...intermediate def grads + crystallite_partionedFi0 (1:3,1:3,1:myNgrains,i,e) = & + crystallite_Fi (1:3,1:3,1:myNgrains,i,e) ! ...intermediate def grads - crystallite_partionedLi0(1:3,1:3,1:myNgrains,i,e) = & - crystallite_Li(1:3,1:3,1:myNgrains,i,e) ! ...intermediate velocity grads + crystallite_partionedLi0 (1:3,1:3,1:myNgrains,i,e) = & + crystallite_Li (1:3,1:3,1:myNgrains,i,e) ! ...intermediate velocity grads - crystallite_partionedS0(1:3,1:3,1:myNgrains,i,e) = & - crystallite_S(1:3,1:3,1:myNgrains,i,e) ! ...2nd PK stress + crystallite_partionedS0 (1:3,1:3,1:myNgrains,i,e) = & + crystallite_S (1:3,1:3,1:myNgrains,i,e) ! ...2nd PK stress do g = 1,myNgrains plasticState (phaseAt(g,i,e))%partionedState0(:,phasememberAt(g,i,e)) = & - plasticState (phaseAt(g,i,e))%state( :,phasememberAt(g,i,e)) + plasticState (phaseAt(g,i,e))%state (:,phasememberAt(g,i,e)) do mySource = 1, phase_Nsources(phaseAt(g,i,e)) sourceState(phaseAt(g,i,e))%p(mySource)%partionedState0(:,phasememberAt(g,i,e)) = & - sourceState(phaseAt(g,i,e))%p(mySource)%state( :,phasememberAt(g,i,e)) + sourceState(phaseAt(g,i,e))%p(mySource)%state (:,phasememberAt(g,i,e)) enddo enddo forall(i = FEsolving_execIP(1,e):FEsolving_execIP(2,e), & homogState(material_homogenizationAt(e))%sizeState > 0) & homogState(material_homogenizationAt(e))%subState0(:,mappingHomogenization(1,i,e)) = & - homogState(material_homogenizationAt(e))%State( :,mappingHomogenization(1,i,e)) ! ...internal homogenization state + homogState(material_homogenizationAt(e))%State (:,mappingHomogenization(1,i,e)) ! ...internal homogenization state forall(i = FEsolving_execIP(1,e):FEsolving_execIP(2,e), & thermalState(material_homogenizationAt(e))%sizeState > 0) & thermalState(material_homogenizationAt(e))%subState0(:,mappingHomogenization(1,i,e)) = & - thermalState(material_homogenizationAt(e))%State( :,mappingHomogenization(1,i,e)) ! ...internal thermal state + thermalState(material_homogenizationAt(e))%State (:,mappingHomogenization(1,i,e)) ! ...internal thermal state forall(i = FEsolving_execIP(1,e):FEsolving_execIP(2,e), & damageState(material_homogenizationAt(e))%sizeState > 0) & damageState(material_homogenizationAt(e))%subState0(:,mappingHomogenization(1,i,e)) = & - damageState(material_homogenizationAt(e))%State( :,mappingHomogenization(1,i,e)) ! ...internal damage state + damageState(material_homogenizationAt(e))%State (:,mappingHomogenization(1,i,e)) ! ...internal damage state materialpoint_subF0(1:3,1:3,i,e) = materialpoint_subF(1:3,1:3,i,e) ! ...def grad endif steppingNeeded