From 5d1e648d17e9a7151b27e6238d6fb13f767da08d Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Tue, 2 Oct 2018 15:25:40 +0200 Subject: [PATCH] drag coefficient B introduced and read in --- src/plastic_dislotwin.f90 | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/src/plastic_dislotwin.f90 b/src/plastic_dislotwin.f90 index 7caff4bdd..d1cdf15db 100644 --- a/src/plastic_dislotwin.f90 +++ b/src/plastic_dislotwin.f90 @@ -109,7 +109,8 @@ module plastic_dislotwin q, & !< q-exponent in glide velocity r, & !< r-exponent in twin nucleation rate s, & !< s-exponent in trans nucleation rate - shear_twin !< characteristic shear for twins + shear_twin, & !< characteristic shear for twins + B !< drag coefficient real(pReal), dimension(:,:), allocatable, private :: & interaction_SlipSlip, & !< coefficients for slip-slip interaction for each interaction type and instance interaction_SlipTwin, & !< coefficients for slip-twin interaction for each interaction type and instance @@ -326,6 +327,9 @@ subroutine plastic_dislotwin_init(fileUnit) if (size(prm%burgers_slip) /= size(prm%Nslip)) call IO_error(150_pInt,ext_msg='slipburgers') prm%burgers_slip = math_expand(prm%burgers_slip,prm%Nslip) + prm%B = config_phase(p)%getFloats('B',defaultVal=[(0.0_pReal, i=1,size(prm%Nslip))]) + prm%B = math_expand(prm%B,prm%Nslip) + prm%Qedge = config_phase(p)%getFloats('qedge') prm%Qedge = math_expand(prm%Qedge,prm%Nslip)