From d953e6bedf687a0f4d56079b48bb2a33519250b4 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 24 Jul 2023 18:29:36 +0200 Subject: [PATCH] unified style --- src/homogenization_damage.f90 | 2 +- src/homogenization_thermal.f90 | 10 ++++++---- 2 files changed, 7 insertions(+), 5 deletions(-) diff --git a/src/homogenization_damage.f90 b/src/homogenization_damage.f90 index 233425ebe..7e83c34d8 100644 --- a/src/homogenization_damage.f90 +++ b/src/homogenization_damage.f90 @@ -158,7 +158,7 @@ module subroutine homogenization_set_phi(phi) integer :: ho, en, ce - do ce=1, ubound(phi,1) + do ce=lbound(phi,1), ubound(phi,1) ho = material_ID_homogenization(ce) en = material_entry_homogenization(ce) damagestate_h(ho)%state(1,en) = phi(ce) diff --git a/src/homogenization_thermal.f90 b/src/homogenization_thermal.f90 index 791286912..069c402eb 100644 --- a/src/homogenization_thermal.f90 +++ b/src/homogenization_thermal.f90 @@ -177,12 +177,14 @@ module subroutine homogenization_thermal_setField(T,dot_T) real(pREAL), dimension(:), intent(in) :: T, dot_T - integer :: ce + integer :: ho, en, ce - do ce=1, min(ubound(T,1),ubound(dot_T,1)) - current(material_ID_homogenization(ce))%T(material_entry_homogenization(ce)) = T(ce) - current(material_ID_homogenization(ce))%dot_T(material_entry_homogenization(ce)) = dot_T(ce) + do ce=max(lbound(T,1),lbound(dot_T,1)), min(ubound(T,1),ubound(dot_T,1)) + ho = material_ID_homogenization(ce) + en = material_entry_homogenization(ce) + current(ho)%T(en) = T(ce) + current(ho)%dot_T(en) = dot_T(ce) call thermal_partition(ce) end do