diff --git a/src/grid/grid_thermal_spectral.f90 b/src/grid/grid_thermal_spectral.f90
index c0a7e3940..bfc9fa52a 100644
--- a/src/grid/grid_thermal_spectral.f90
+++ b/src/grid/grid_thermal_spectral.f90
@@ -38,8 +38,8 @@ module grid_thermal_spectral
   type(tSolutionParams) :: params
 !--------------------------------------------------------------------------------------------------
 ! PETSc data
-  SNES     :: thermal_snes
-  Vec      :: solution_vec
+  SNES :: thermal_snes
+  Vec :: solution_vec
   real(pReal), dimension(:,:,:), allocatable :: &
     T_current, &                                                                                    !< field of current temperature
     T_lastInc, &                                                                                    !< field of previous temperature
@@ -99,6 +99,18 @@ subroutine grid_thermal_spectral_init(T_0)
  call PetscOptionsInsertString(PETSC_NULL_OPTIONS,num_grid%get_asString('petsc_options',defaultVal=''),err_PETSc)
  CHKERRQ(err_PETSc)
 
+!--------------------------------------------------------------------------------------------------
+! init fields
+  allocate(T_current(grid(1),grid(2),grid3), source=T_0)
+  allocate(T_lastInc(grid(1),grid(2),grid3), source=T_0)
+  allocate(T_stagInc(grid(1),grid(2),grid3), source=T_0)
+
+  ce = 0
+  do k = 1, grid3; do j = 1, grid(2); do i = 1,grid(1)
+    ce = ce + 1
+    call homogenization_thermal_setField(T_0,0.0_pReal,ce)
+  end do; end do; end do
+
 !--------------------------------------------------------------------------------------------------
 ! initialize solver specific parts of PETSc
   call SNESCreate(PETSC_COMM_WORLD,thermal_snes,err_PETSc); CHKERRQ(err_PETSc)
@@ -126,20 +138,6 @@ subroutine grid_thermal_spectral_init(T_0)
   call SNESSetFromOptions(thermal_snes,err_PETSc); CHKERRQ(err_PETSc)                               ! pull it all together with additional CLI arguments
 
 !--------------------------------------------------------------------------------------------------
-! init fields
-  allocate(T_current(grid(1),grid(2),grid3), source=0.0_pReal)
-  allocate(T_lastInc(grid(1),grid(2),grid3), source=0.0_pReal)
-  allocate(T_stagInc(grid(1),grid(2),grid3), source=0.0_pReal)
-
-  ce = 0
-  do k = 1, grid3; do j = 1, grid(2); do i = 1,grid(1)
-    ce = ce + 1
-    T_current(i,j,k) = T_0
-    T_lastInc(i,j,k) = T_current(i,j,k)
-    T_stagInc(i,j,k) = T_current(i,j,k)
-    call homogenization_thermal_setField(T_0,0.0_pReal,ce)
-  end do; end do; end do
-
   call DMDAVecGetArrayF90(thermal_grid,solution_vec,T_PETSc,err_PETSc)
   CHKERRQ(err_PETSc)
   T_PETSc = T_current
diff --git a/src/parallelization.f90 b/src/parallelization.f90
index d1e1576bb..29deaf724 100644
--- a/src/parallelization.f90
+++ b/src/parallelization.f90
@@ -114,7 +114,7 @@ subroutine parallelization_init
   if (worldrank /= 0) then
     close(OUTPUT_UNIT)                                                                              ! disable output
     write(rank_str,'(i4.4)') worldrank                                                              ! use for MPI debug filenames
-    open(OUTPUT_UNIT,file='./log.'//rank_str,status='replace')                                             ! close() alone will leave some temp files in cwd
+    open(OUTPUT_UNIT,file='/dev/null',status='replace')                                             ! close() alone will leave some temp files in cwd
   endif
 
 !$ call get_environment_variable(name='OMP_NUM_THREADS',value=NumThreadsString,STATUS=got_env)