crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 21:55:41 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Don't allocate/deallocate p0 from openmc_material_set_densities

Browse source
This commit is contained in:
Paul Romano 2017-12-18 10:33:30 +07:00
parent 9aaa927222
commit e698403076
1 changed files with 2 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

7
src/material_header.F90
View file

@ -455,8 +455,8 @@ contains
associate (m => materials(index))
! If nuclide/density arrays are not correct size, reallocate
if (n /= m % n_nuclides) then
deallocate(m % nuclide, m % atom_density, m % p0, STAT=stat)
allocate(m % nuclide(n), m % atom_density(n), m % p0(n))
deallocate(m % nuclide, m % atom_density, STAT=stat)
allocate(m % nuclide(n), m % atom_density(n))
end if
do i = 1, n
@ -474,9 +474,6 @@ contains
end do
m % n_nuclides = n
! Set isotropic flags to flags
m % p0(:) = .false.
! Set total density to the sum of the vector
err = m % set_density(sum(density))