From e698403076f8ed6476b15e314285a05b9bcd7c79 Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Mon, 18 Dec 2017 10:33:30 +0700 Subject: [PATCH] Don't allocate/deallocate p0 from openmc_material_set_densities --- src/material_header.F90 | 7 ++----- 1 file changed, 2 insertions(+), 5 deletions(-) diff --git a/src/material_header.F90 b/src/material_header.F90 index a3c14f8cab..9d93131acf 100644 --- a/src/material_header.F90 +++ b/src/material_header.F90 @@ -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))