Changed cell rotation to read in as double

This commit is contained in:
Bryan Herman 2015-02-19 07:59:15 -05:00
parent ea1df445d2
commit e2f7fce57d

View file

@ -904,7 +904,7 @@ contains
integer :: universe_num
integer :: n_cells_in_univ
integer :: coeffs_reqd
integer :: temp_int_array3(3)
integer :: temp_double_array3(3)
integer, allocatable :: temp_int_array(:)
real(8) :: phi, theta, psi
logical :: file_exists
@ -1052,10 +1052,10 @@ contains
end if
! Copy rotation angles in x,y,z directions
call get_node_array(node_cell, "rotation", temp_int_array3)
phi = -temp_int_array3(1) * PI/180.0_8
theta = -temp_int_array3(2) * PI/180.0_8
psi = -temp_int_array3(3) * PI/180.0_8
call get_node_array(node_cell, "rotation", temp_double_array3)
phi = -temp_double_array3(1) * PI/180.0_8
theta = -temp_double_array3(2) * PI/180.0_8
psi = -temp_double_array3(3) * PI/180.0_8
! Calculate rotation matrix based on angles given
allocate(c % rotation(3,3))