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

real must be pReal for MSC.Marc

This commit is contained in:
Martin Diehl 2018-09-13 05:49:57 +02:00
parent 4266872965
commit 8ecb019566
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/plastic_phenopowerlaw.f90
View File

@ -675,16 +675,16 @@ subroutine kinetics_slip(prm,stt,of,S,gdot_slip_pos,gdot_slip_neg, &
stt
integer(pInt), intent(in) :: &
of
real, dimension(prm%totalNslip), intent(out) :: &
real(pReal), dimension(prm%totalNslip), intent(out) :: &
gdot_slip_pos, &
gdot_slip_neg
real, dimension(prm%totalNslip), optional, intent(out) :: &
real(pReal), dimension(prm%totalNslip), optional, intent(out) :: &
dgdot_dtau_slip_pos, &
dgdot_dtau_slip_neg
real(pReal), dimension(3,3), intent(in) :: &
S
real, dimension(prm%totalNslip) :: &
real(pReal), dimension(prm%totalNslip) :: &
tau_slip_pos, &
tau_slip_neg
@ -741,14 +741,15 @@ subroutine kinetics_twin(prm,stt,of,S,gdot_twin,dgdot_dtau_twin)
of
real(pReal), dimension(3,3), intent(in) :: &
S
real, dimension(prm%totalNtwin), intent(out) :: &
real(pReal), dimension(prm%totalNtwin), intent(out) :: &
gdot_twin
real, dimension(prm%totalNtwin), optional, intent(out) :: &
real(pReal), dimension(prm%totalNtwin), optional, intent(out) :: &
dgdot_dtau_twin
real, dimension(prm%totalNtwin) :: &
real(pReal), dimension(prm%totalNtwin) :: &
tau_twin
integer(pInt) :: i
do i = 1_pInt, prm%totalNtwin
tau_twin(i) = math_mul33xx33(S,prm%Schmid_twin(1:3,1:3,i))
enddo