From 4a93dbec117cbb0092ddc8f33e1a5cbe00097d2b Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Tue, 19 Apr 2022 13:11:40 -0400 Subject: [PATCH] more descriptive variable names --- src/material.f90 | 30 +++++++++++++++--------------- 1 file changed, 15 insertions(+), 15 deletions(-) diff --git a/src/material.f90 b/src/material.f90 index 7264738e5..2258c40a1 100644 --- a/src/material.f90 +++ b/src/material.f90 @@ -95,9 +95,9 @@ subroutine parse() integer, dimension(:), allocatable :: & counterPhase, & counterHomogenization, & - temp_ho - integer, dimension(:,:), allocatable :: temp_ph - real(pReal), dimension(:,:), allocatable :: temp_v + ho_of + integer, dimension(:,:), allocatable :: ph_of + real(pReal), dimension(:,:), allocatable :: v_of real(pReal) :: v integer :: & @@ -135,9 +135,9 @@ subroutine parse() allocate(material_O_0(materials%length)) allocate(material_F_i_0(materials%length)) - allocate(temp_ho(materials%length)) - allocate(temp_ph(materials%length,homogenization_maxNconstituents),source=-1) - allocate(temp_v(materials%length,homogenization_maxNconstituents),source=0.0_pReal) + allocate(ho_of(materials%length)) + allocate(ph_of(materials%length,homogenization_maxNconstituents),source=-1) + allocate( v_of(materials%length,homogenization_maxNconstituents),source=0.0_pReal) ! parse YAML structure select type(materials) @@ -147,10 +147,10 @@ subroutine parse() item => materials%first do ma = 1, materials%length material => item%node - temp_ho(ma) = homogenizations%getIndex(material%get_asString('homogenization')) + ho_of(ma) = homogenizations%getIndex(material%get_asString('homogenization')) constituents => material%get('constituents') - homogenization => homogenizations%get(temp_ho(ma)) + homogenization => homogenizations%get(ho_of(ma)) if (constituents%length /= homogenization%get_asInt('N_constituents')) call IO_error(148) allocate(material_O_0(ma)%data(constituents%length)) @@ -158,14 +158,14 @@ subroutine parse() do co = 1, constituents%length constituent => constituents%get(co) - temp_v(ma,co) = constituent%get_asFloat('v') - temp_ph(ma,co) = phases%getIndex(constituent%get_asString('phase')) + v_of(ma,co) = constituent%get_asFloat('v') + ph_of(ma,co) = phases%getIndex(constituent%get_asString('phase')) call material_O_0(ma)%data(co)%fromQuaternion(constituent%get_as1dFloat('O',requiredSize=4)) material_F_i_0(ma)%data(1:3,1:3,co) = constituent%get_as2dFloat('F_i',defaultVal=math_I3,requiredShape=[3,3]) end do - if (dNeq(sum(temp_v(ma,:)),1.0_pReal,1.e-9_pReal)) call IO_error(153,ext_msg='constituent') + if (dNeq(sum(v_of(ma,:)),1.0_pReal,1.e-9_pReal)) call IO_error(153,ext_msg='constituent') item => item%next end do @@ -187,7 +187,7 @@ subroutine parse() do el = 1, discretization_Nelems ma = discretization_materialAt(el) - ho = temp_ho(ma) + ho = ho_of(ma) do ip = 1, discretization_nIPs ce = (el-1)*discretization_nIPs + ip @@ -196,10 +196,10 @@ subroutine parse() material_homogenizationEntry(ce) = counterHomogenization(ho) end do - do co = 1, size(temp_ph(ma,:)>0) + do co = 1, size(ph_of(ma,:)>0) - v = temp_v(ma,co) - ph = temp_ph(ma,co) + v = v_of(ma,co) + ph = ph_of(ma,co) do ip = 1, discretization_nIPs ce = (el-1)*discretization_nIPs + ip