223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

reflect change of mesh_element(2,e) in mesh.f90 

-change of  line, and those below, will be ignored--

M    constitutive_pheno.f90
This commit is contained in:
Philip Eisenlohr 2008-03-25 13:59:22 +00:00
parent 3640a48095
commit 7faf4093a5
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
trunk/constitutive_pheno.f90
View File

@ -535,7 +535,7 @@ subroutine constitutive_Assignment()
!********************************************************************* !*********************************************************************
use prec, only: pReal,pInt use prec, only: pReal,pInt
use math, only: math_sampleGaussOri,math_sampleFiberOri,math_sampleRandomOri,math_symmetricEulers,math_EulerToR use math, only: math_sampleGaussOri,math_sampleFiberOri,math_sampleRandomOri,math_symmetricEulers,math_EulerToR
use mesh, only: mesh_NcpElems,FE_Nips,FE_mapElemtype,mesh_maxNips,mesh_element use mesh, only: mesh_NcpElems,FE_Nips,mesh_maxNips,mesh_element
use IO, only: IO_hybridIA use IO, only: IO_hybridIA
use crystal, only: crystal_SlipIntType use crystal, only: crystal_SlipIntType
@ -591,7 +591,7 @@ constitutive_maxNstatevars = material_maxNslip + 0_pInt
allocate(texture_totalNgrains(texture_maxN)) ; texture_totalNgrains=0_pInt allocate(texture_totalNgrains(texture_maxN)) ; texture_totalNgrains=0_pInt
do i=1,mesh_NcpElems do i=1,mesh_NcpElems
texID = mesh_element(4,i) texID = mesh_element(4,i)
texture_totalNgrains(texID) = texture_totalNgrains(texID) + FE_Nips(FE_mapElemtype(mesh_element(2,i)))*texture_Ngrains(texID) texture_totalNgrains(texID) = texture_totalNgrains(texID) + FE_Nips(mesh_element(2,i))*texture_Ngrains(texID)
enddo enddo
! generate hybridIA samplings for ODFfile textures to later draw from these populations ! generate hybridIA samplings for ODFfile textures to later draw from these populations
@ -621,7 +621,7 @@ constitutive_HardeningMatrix=0.0_pReal
do e=1,mesh_NcpElems do e=1,mesh_NcpElems
matID=mesh_element(3,e) matID=mesh_element(3,e)
texID=mesh_element(4,e) texID=mesh_element(4,e)
do i=1,FE_Nips(FE_mapElemtype(mesh_element(2,e))) do i=1,FE_Nips(mesh_element(2,e))
g = 0_pInt ! grain counter g = 0_pInt ! grain counter
do m = 1,multiplicity(texID) do m = 1,multiplicity(texID)
o = 0_pInt ! component counter o = 0_pInt ! component counter