223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

Fix for 3881 because gfortran was complaining

This commit is contained in:
Su Leen Wong 2015-02-24 10:15:07 +00:00
parent fc1a948873
commit 2baac6fc10
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
code/plastic_dislotwin.f90
View File

@ -1871,16 +1871,17 @@ subroutine plastic_dislotwin_dotState(Tstar_v,Temperature,ipc,ip,el)
shear_trans = 0.0_pReal
shearrate_trans = 0.0_pReal
probrate_trans = 0.0_pReal
#ifdef TRANS
j = 0_pInt
do f = 1_pInt,lattice_maxNtransFamily ! loop over all trans families
index_myFamily = sum(lattice_NtransSystem(1:f-1_pInt,ph)) ! at which index starts my family
!* Projection of shear and shear rate on fault band (twin) systems
shear_trans = matmul(plastic_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%state(2_pInt*ns+1_pInt:3_pInt*ns, of)) This is not safe
shearrate_trans = matmul(plastic_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%dotState(2_pInt*ns+1_pInt:3_pInt*ns, of)) This is not safe
if (nr > 0_pInt .and. ns == nr) then
shear_trans = matmul(plastic_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%state(2_pInt*ns+1_pInt:3_pInt*ns, of))
shearrate_trans = matmul(plastic_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%dotState(2_pInt*ns+1_pInt:3_pInt*ns, of))
endif
do i = 1_pInt,plastic_dislotwin_Ntrans(f,instance) ! process each (active) trans system in family
j = j+1_pInt
@ -1930,7 +1931,7 @@ subroutine plastic_dislotwin_dotState(Tstar_v,Temperature,ipc,ip,el)
enddo
enddo
#endif
end subroutine plastic_dislotwin_dotState