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

corrected wrong integer arithmetics in calculation of coordinates in DAMASK_spectral.f90

This commit is contained in:
Krishna Komerla 2012-04-24 16:59:38 +00:00
parent 99f2b8603a
commit e33c34d86b
2 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
code/DAMASK_spectral.f90
View File

@ -535,12 +535,12 @@ program DAMASK_spectral
ielem = ielem + 1_pInt ielem = ielem + 1_pInt
F(i,j,k,1:3,1:3) = math_I3 F(i,j,k,1:3,1:3) = math_I3
F_lastInc(i,j,k,1:3,1:3) = math_I3 F_lastInc(i,j,k,1:3,1:3) = math_I3
coordinates(i,j,k,1:3) = geomdim/real(res * [i,j,k], pReal) - geomdim/real(2_pInt*res,pReal) coordinates(i,j,k,1:3) = geomdim/real(res,pReal)*real([i,j,k],pReal) - geomdim/real(2_pInt*res,pReal)
call CPFEM_general(2_pInt,coordinates(i,j,k,1:3),math_I3,math_I3,temperature(i,j,k),& call CPFEM_general(2_pInt,coordinates(i,j,k,1:3),math_I3,math_I3,temperature(i,j,k),&
0.0_pReal,ielem,1_pInt,sigma,dsde,P_real(i,j,k,1:3,1:3),dPdF) 0.0_pReal,ielem,1_pInt,sigma,dsde,P_real(i,j,k,1:3,1:3),dPdF)
C = C + dPdF C = C + dPdF
enddo; enddo; enddo enddo; enddo; enddo
C_ref = C * wgt ! linear reference material stiffness C_ref = C * wgt
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! calculate the gamma operator ! calculate the gamma operator

1
code/IO.f90
View File

@ -300,6 +300,7 @@ subroutine IO_write_jobBinaryFile(myUnit,newExt,recMultiplier)
open(myUnit,status='replace',form='unformatted',access='direct', & open(myUnit,status='replace',form='unformatted',access='direct', &
recl=pReal,iostat=myStat,file=path) recl=pReal,iostat=myStat,file=path)
endif endif
if (myStat /= 0_pInt) call IO_error(100_pInt,ext_msg=path) if (myStat /= 0_pInt) call IO_error(100_pInt,ext_msg=path)
end subroutine IO_write_jobBinaryFile end subroutine IO_write_jobBinaryFile