From 7ad343c9b61e40a6b29d698ea11c82f4a9a3bd79 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sat, 29 Jun 2019 15:06:47 -0700 Subject: [PATCH] polishing --- src/IO.f90 | 2 ++ src/grid/DAMASK_grid.f90 | 2 +- src/plastic_dislotwin.f90 | 4 ++-- 3 files changed, 5 insertions(+), 3 deletions(-) diff --git a/src/IO.f90 b/src/IO.f90 index d245cccb3..940fb2da9 100644 --- a/src/IO.f90 +++ b/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) diff --git a/src/grid/DAMASK_grid.f90 b/src/grid/DAMASK_grid.f90 index 5beab767d..ef863e017 100644 --- a/src/grid/DAMASK_grid.f90 +++ b/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 diff --git a/src/plastic_dislotwin.f90 b/src/plastic_dislotwin.f90 index 30af6ed5d..82498b0bf 100644 --- a/src/plastic_dislotwin.f90 +++ b/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