From 2969a53338d11efec1f026835afc1afdac911172 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 2 Dec 2019 12:54:45 +0100 Subject: [PATCH] not used/needed anymore --- src/plastic_nonlocal.f90 | 137 --------------------------------------- 1 file changed, 137 deletions(-) diff --git a/src/plastic_nonlocal.f90 b/src/plastic_nonlocal.f90 index c2400494c..23bfb50aa 100644 --- a/src/plastic_nonlocal.f90 +++ b/src/plastic_nonlocal.f90 @@ -215,7 +215,6 @@ module plastic_nonlocal plastic_nonlocal_dotState, & plastic_nonlocal_deltaState, & plastic_nonlocal_updateCompatibility, & - plastic_nonlocal_postResults, & plastic_nonlocal_results private :: & @@ -1956,142 +1955,6 @@ subroutine plastic_nonlocal_updateCompatibility(orientation,i,e) end subroutine plastic_nonlocal_updateCompatibility -!-------------------------------------------------------------------------------------------------- -!> @brief return array of constitutive results -!-------------------------------------------------------------------------------------------------- -function plastic_nonlocal_postResults(ph,instance,of) result(postResults) - - integer, intent(in) :: & - ph, & - instance, & - of - - real(pReal), dimension(sum(plastic_nonlocal_sizePostResult(:,instance))) :: & - postResults - - integer :: & - ns, & !< short notation for the total number of active slip systems - c, & !< character of dislocation - cs, & !< constitutive result index - o, & !< index of current output - t, & !< type of dislocation - s !< index of my current slip system - - real(pReal), dimension(param(instance)%totalNslip,8) :: & - rhoSgl, & - rhoDotSgl !< evolution rate of single dislocation densities (positive/negative screw and edge without dipoles) - real(pReal), dimension(param(instance)%totalNslip,2) :: & - rhoDotDip !< evolution rate of dipole dislocation densities (screw and edge dipoles) - - -ns = param(instance)%totalNslip - -cs = 0 - -associate(prm => param(instance),dst => microstructure(instance),stt=>state(instance),dot => dotState(instance)) - -forall (s = 1:ns, t = 1:4) - rhoSgl(s,t+4) = plasticState(ph)%State(iRhoB(s,t,instance),of) - rhoDotSgl(s,t+4) = plasticState(ph)%dotState(iRhoB(s,t,instance),of) -endforall -forall (s = 1:ns, c = 1:2) - rhoDotDip(s,c) = plasticState(ph)%dotState(iRhoD(s,c,instance),of) -endforall - -outputsLoop: do o = 1,size(param(instance)%outputID) - select case(param(instance)%outputID(o)) - - case (rho_sgl_mob_edg_pos_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_mob_edg_pos(:,of) - cs = cs + ns - - case (rho_sgl_imm_edg_pos_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_imm_edg_pos(:,of) - cs = cs + ns - - case (rho_sgl_mob_edg_neg_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_mob_edg_neg(:,of) - cs = cs + ns - - case (rho_sgl_imm_edg_neg_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_imm_edg_neg(:,of) - cs = cs + ns - - case (rho_dip_edg_ID) - postResults(cs+1:cs+ns) = stt%rho_dip_edg(:,of) - cs = cs + ns - - case (rho_sgl_mob_scr_pos_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_mob_scr_pos(:,of) - cs = cs + ns - - case (rho_sgl_imm_scr_pos_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_imm_scr_pos(:,of) - cs = cs + ns - - case (rho_sgl_mob_scr_neg_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_mob_scr_neg(:,of) - cs = cs + ns - - case (rho_sgl_imm_scr_neg_ID) - postResults(cs+1:cs+ns) = stt%rho_sgl_imm_scr_neg(:,of) - cs = cs + ns - - case (rho_dip_scr_ID) - postResults(cs+1:cs+ns) = stt%rho_dip_scr(:,of) - cs = cs + ns - - case (rho_forest_ID) - postResults(cs+1:cs+ns) = stt%rho_forest(:,of) - cs = cs + ns - - case (resolvedstress_back_ID) - postResults(cs+1:cs+ns) = dst%tau_back(:,of) - cs = cs + ns - - case (tau_pass_ID) - postResults(cs+1:cs+ns) = dst%tau_pass(:,of) - cs = cs + ns - - case (rho_dot_sgl_ID) - postResults(cs+1:cs+ns) = sum(rhoDotSgl(1:ns,1:4),2) & - + sum(rhoDotSgl(1:ns,5:8)*sign(1.0_pReal,rhoSgl(1:ns,5:8)),2) - cs = cs + ns - - case (rho_dot_sgl_mobile_ID) - postResults(cs+1:cs+ns) = sum(rhoDotSgl(1:ns,1:4),2) - cs = cs + ns - - case (rho_dot_dip_ID) - postResults(cs+1:cs+ns) = sum(rhoDotDip,2) - cs = cs + ns - - case (v_edg_pos_ID) - postResults(cs+1:cs+ns) = stt%v_edg_pos(:,of) - cs = cs + ns - - case (v_edg_neg_ID) - postResults(cs+1:cs+ns) = stt%v_edg_neg(:,of) - cs = cs + ns - - case (v_scr_pos_ID) - postResults(cs+1:cs+ns) = stt%v_scr_pos(:,of) - cs = cs + ns - - case (v_scr_neg_ID) - postResults(cs+1:cs+ns) = stt%v_scr_neg(:,of) - cs = cs + ns - - case(gamma_ID) - postResults(cs+1:cs+ns) = stt%gamma(:,of) - cs = cs + ns - - end select -enddo outputsLoop -end associate -end function plastic_nonlocal_postResults - - !-------------------------------------------------------------------------------------------------- !> @brief returns copy of current dislocation densities from state !> @details raw values is rectified