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

an allocation statement was deleted from rev 518 to rev 520!!!

This commit is contained in:
Christoph Kords 2010-03-30 13:15:18 +00:00
parent 67c523e4b1
commit 1bf740eb9d
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
code/math.f90
View File

@ -168,7 +168,6 @@ subroutine math_misorientation(dQ, Q1, Q2, symmetryType)
integer(pInt), dimension(2), parameter :: NsymOperations = (/24,12/) ! number of possible symmetry operations
real(pReal), dimension(:,:), allocatable :: mySymOperations ! symmetry Operations for my crystal symmetry
dQ = 0.0_pReal
if (symmetryType < 1_pInt .or. symmetryType > 2_pInt) then
@ -176,6 +175,8 @@ subroutine math_misorientation(dQ, Q1, Q2, symmetryType)
return
endif
allocate(mySymOperations(4,NsymOperations(symmetryType)))
mySymOperations = symOperations(:,sum(NsymOperations(1:symmetryType-1))+1:sum(NsymOperations(1:symmetryType))) ! choose symmetry operations according to crystal symmetry
dQ(1) = -1.0_pReal ! start with maximum misorientation angle
do s = 1,NsymOperations(symmetryType) ! loop ver symmetry operations