From 7b28dedfe395d2478c263f075068a8e0c758750c Mon Sep 17 00:00:00 2001
From: Martin	Diehl <m.diehl@mpie.de>
Date: Tue, 22 Oct 2013 15:52:58 +0000
Subject: [PATCH] removed temperature calculation and debugging statements

---
 code/DAMASK_spectral_utilities.f90 | 19 ++++++++-----------
 1 file changed, 8 insertions(+), 11 deletions(-)

diff --git a/code/DAMASK_spectral_utilities.f90 b/code/DAMASK_spectral_utilities.f90
index 200ef5b87..186c787ce 100644
--- a/code/DAMASK_spectral_utilities.f90
+++ b/code/DAMASK_spectral_utilities.f90
@@ -1014,16 +1014,16 @@ subroutine utilities_temperatureUpdate()
    materialpoint_heat
   
  implicit none
- integer x,y,z,e
- real(pReal) :: a
- do e=1_pInt, product(grid)
-   if (e==1) print*, materialpoint_heat(1,e)*0.0001_pReal
+ integer :: &
+   x,y,z,e
+ real(pReal) :: & 
+  a
+ forall(e=1_pInt:product(grid)) &
    crystallite_temperature(1,e) = crystallite_temperature(1,e) + materialpoint_heat(1,e)*0.0001_pReal
- enddo
  e = 0_pInt
- z = 0
- y = 0
- x = 0
+ z = 0_pInt
+ y = 0_pInt
+ x = 0_pInt
 
 !< 6 or less neighboring IPs as [element_num, IP_index, neighbor_index that points to me]
  do z = 0_pInt,grid(3)-1_pInt
@@ -1055,9 +1055,6 @@ subroutine utilities_temperatureUpdate()
                                       + y * grid(1) &
                                       + x &
                                       + 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
      enddo
    enddo