should not be handled by the individual model

avg(P) and avg(F) exist independently of RGC
This commit is contained in:
Martin Diehl 2018-08-28 13:02:44 +02:00
parent 7f00082d60
commit f9214c8e1b
1 changed files with 2 additions and 7 deletions

View File

@ -162,7 +162,7 @@ subroutine homogenization_RGC_init(fileUnit)
homogenization_RGC_orientation = spread(spread(math_I3,3,mesh_maxNips),4,mesh_NcpElems) ! initialize to identity homogenization_RGC_orientation = spread(spread(math_I3,3,mesh_maxNips),4,mesh_NcpElems) ! initialize to identity
rewind(fileUnit) rewind(fileUnit)
do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>'))/=material_partHomogenization) ! wind forward to <homogenization> do while (trim(line) /= IO_EOF .and. IO_lc(IO_getTag(line,'<','>'))/=material_partHomogenization) ! wind forward to <homogenization>
line = IO_read(fileUnit) line = IO_read(fileUnit)
enddo enddo
@ -200,12 +200,7 @@ subroutine homogenization_RGC_init(fileUnit)
homogenization_RGC_outputID(homogenization_RGC_Noutput(i),i) = maximumrelaxrate_ID homogenization_RGC_outputID(homogenization_RGC_Noutput(i),i) = maximumrelaxrate_ID
case('magnitudemismatch') case('magnitudemismatch')
homogenization_RGC_outputID(homogenization_RGC_Noutput(i),i) = magnitudemismatch_ID homogenization_RGC_outputID(homogenization_RGC_Noutput(i),i) = magnitudemismatch_ID
case('ipcoords')
homogenization_RGC_outputID(homogenization_RGC_Noutput(i),i) = ipcoords_ID
case('avgdefgrad','avgf')
homogenization_RGC_outputID(homogenization_RGC_Noutput(i),i) = avgdefgrad_ID
case('avgp','avgfirstpiola','avg1stpiola')
homogenization_RGC_outputID(homogenization_RGC_Noutput(i),i) = avgfirstpiola_ID
case default case default
homogenization_RGC_Noutput(i) = homogenization_RGC_Noutput(i) -1_pInt ! correct for invalid homogenization_RGC_Noutput(i) = homogenization_RGC_Noutput(i) -1_pInt ! correct for invalid