From 53cff7458c572e46966dc24090cc796ee26af0f4 Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Fri, 24 Jul 2009 15:23:45 +0000 Subject: [PATCH] polishing of constitutive output size determination in _init() --- code/constitutive_phenopowerlaw.f90 | 24 +++++++++++------------- 1 file changed, 11 insertions(+), 13 deletions(-) diff --git a/code/constitutive_phenopowerlaw.f90 b/code/constitutive_phenopowerlaw.f90 index 1b0c38552..d55fee913 100644 --- a/code/constitutive_phenopowerlaw.f90 +++ b/code/constitutive_phenopowerlaw.f90 @@ -339,21 +339,19 @@ subroutine constitutive_phenopowerlaw_init(file) do i = 1,maxNinstance do j = 1,maxval(phase_Noutput) select case(constitutive_phenopowerlaw_output(j,i)) - case('resistance_slip') + case('resistance_slip', & + 'shearrate_slip', & + 'resolvedstress_slip' & + ) mySize = constitutive_phenopowerlaw_totalNslip(i) - case('shearrate_slip') - mySize = constitutive_phenopowerlaw_totalNslip(i) - case('resolvedstress_slip') - mySize = constitutive_phenopowerlaw_totalNslip(i) - case('totalshear') - mySize = 1_pInt - case('resistance_twin') + case('resistance_twin', & + 'shearrate_twin', & + 'resolvedstress_twin' & + ) mySize = constitutive_phenopowerlaw_totalNtwin(i) - case('shearrate_twin') - mySize = constitutive_phenopowerlaw_totalNtwin(i) - case('resolvedstress_twin') - mySize = constitutive_phenopowerlaw_totalNtwin(i) - case('totalvolfrac') + case('totalshear', & + 'totalvolfrac' & + ) mySize = 1_pInt case default mySize = 0_pInt