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

Reading of euler angles from the loadcase corrected.

This commit is contained in:
Mahesh Balasubramaniam 2012-08-02 17:52:07 +00:00
parent 740357172a
commit 3ffba037d4
1 changed files with 0 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
code/DAMASK_spectral.f90
View File

@ -361,8 +361,6 @@ program DAMASK_spectral
forall(k = 1_pInt:3_pInt) temp33_Real(k,1) = IO_floatValue(line,positions,j+l+k) forall(k = 1_pInt:3_pInt) temp33_Real(k,1) = IO_floatValue(line,positions,j+l+k)
if (p==1_pInt) temp33_Real = temp33_Real * inRad if (p==1_pInt) temp33_Real = temp33_Real * inRad
bc(loadcase)%rotation = math_EulerToR(temp33_Real(:,1)) bc(loadcase)%rotation = math_EulerToR(temp33_Real(:,1))
print *,"euler angles", temp33_Real(:,1)
print *,"Rotation Matrix", bc(loadcase)%rotation
case('rotation','rot') ! assign values for the rotation of loadcase matrix case('rotation','rot') ! assign values for the rotation of loadcase matrix
temp_valueVector = 0.0_pReal temp_valueVector = 0.0_pReal
forall (k = 1_pInt:9_pInt) temp_valueVector(k) = IO_floatValue(line,positions,j+k) forall (k = 1_pInt:9_pInt) temp_valueVector(k) = IO_floatValue(line,positions,j+k)