updating physics in pheno+

This commit is contained in:
zhangc43 2016-04-18 11:16:59 -04:00
parent 1e406169ec
commit 20bae5f91b
1 changed files with 1 additions and 1 deletions

View File

@ -981,7 +981,7 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
ENDDO LOOPMYNEIGHBORS ENDDO LOOPMYNEIGHBORS
!*******check if one of the neighbor already can provide a kick for this slip system !*******check if one of the neighbor already can provide a kick for this slip system
IF ( maxval(ne_mprimes) > mprime_cut ) THEN IF ( maxval(ne_mprimes) > mprime_cut ) THEN
plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal + 0.2_pReal*maxval(ne_mprimes) plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal + 0.5_pReal*maxval(ne_mprimes)
ELSE ELSE
plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal - 0.1_pReal*maxval(ne_mprimes) plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal - 0.1_pReal*maxval(ne_mprimes)
ENDIF ENDIF