polishing

src/IO.f90

@@ -783,6 +783,8 @@ subroutine IO_error(error_ID,el,ip,g,instance,ext_msg)
       msg = 'incomplete loadcase'
     case (838)
       msg = 'mixed boundary conditions allow rotation'
+    case (839)
+      msg = 'non-positive restart frequency in grid load case'
     case (841)
       msg = 'missing header length info in grid mesh'
     case (842)

src/grid/DAMASK_grid.f90

@@ -299,7 +299,7 @@ program DAMASK_spectral
      write(6,'(2x,a,i5)')  'increments: ', newLoadCase%incs
      if (newLoadCase%outputfrequency < 1)  errorID = 836                                            ! non-positive result frequency
      write(6,'(2x,a,i5)')  'output  frequency:  ', newLoadCase%outputfrequency
-     if (newLoadCase%restartfrequency < 1)  errorID = 836                                           ! non-positive restart frequency
+     if (newLoadCase%restartfrequency < 1)  errorID = 839                                           ! non-positive restart frequency
      write(6,'(2x,a,i5)')  'restart frequency:  ', newLoadCase%restartfrequency
      if (errorID > 0) call IO_error(error_ID = errorID, ext_msg = loadcase_string)                  ! exit with error message
    endif reportAndCheck

src/plastic_dislotwin.f90

@@ -1239,8 +1239,8 @@ pure subroutine kinetics_trans(Mp,T,dot_gamma_sl,instance,of,&
  enddo
 
  significantStress: where(tau > tol_math_check)
-   StressRatio_s = (dst%tau_hat_tr(:,of)/tau)**prm%s
-   dot_gamma_tr  = dst%V_tr(:,of) * Ndot0*exp(-StressRatio_s)
+   StressRatio_s   = (dst%tau_hat_tr(:,of)/tau)**prm%s
+   dot_gamma_tr    = dst%V_tr(:,of) * Ndot0*exp(-StressRatio_s)
    ddot_gamma_dtau = (dot_gamma_tr*prm%s/tau)*StressRatio_s
  else where significantStress
    dot_gamma_tr  = 0.0_pReal