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

removed temperature calculation and debugging statements

This commit is contained in:
Martin Diehl 2013-10-22 15:52:58 +00:00
parent ea664688f8
commit 7b28dedfe3
1 changed files with 8 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
code/DAMASK_spectral_utilities.f90
View File

@ -1014,16 +1014,16 @@ subroutine utilities_temperatureUpdate()
materialpoint_heat materialpoint_heat
implicit none implicit none
integer x,y,z,e integer :: &
real(pReal) :: a x,y,z,e
do e=1_pInt, product(grid) real(pReal) :: &
if (e==1) print*, materialpoint_heat(1,e)*0.0001_pReal a
forall(e=1_pInt:product(grid)) &
crystallite_temperature(1,e) = crystallite_temperature(1,e) + materialpoint_heat(1,e)*0.0001_pReal crystallite_temperature(1,e) = crystallite_temperature(1,e) + materialpoint_heat(1,e)*0.0001_pReal
enddo
e = 0_pInt e = 0_pInt
z = 0 z = 0_pInt
y = 0 y = 0_pInt
x = 0 x = 0_pInt
!< 6 or less neighboring IPs as [element_num, IP_index, neighbor_index that points to me] !< 6 or less neighboring IPs as [element_num, IP_index, neighbor_index that points to me]
do z = 0_pInt,grid(3)-1_pInt do z = 0_pInt,grid(3)-1_pInt
@ -1055,9 +1055,6 @@ subroutine utilities_temperatureUpdate()
+ y * grid(1) & + y * grid(1) &
+ x & + x &
+ 1_pInt) + 1_pInt)
if (e==1) print*, a
if (e==1) print*, crystallite_temperature(1,e)
if (e==1) print*,(crystallite_temperature(1,e)+a)/7.0_pReal
crystallite_temperature(1,e) = (crystallite_temperature(1,e)+a)/7.0_pReal crystallite_temperature(1,e) = (crystallite_temperature(1,e)+a)/7.0_pReal
enddo enddo
enddo enddo