crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-29 06:35:48 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Use move_alloc to expand array

Browse source
This commit is contained in:
Sterling Harper 2013-10-04 10:17:30 -04:00
parent be6a8d22ec
commit beb39a44bb
1 changed files with 4 additions and 9 deletions
Download patch file Download diff file Expand all files Collapse all files

13
src/track_output.F90
View file

@ -33,18 +33,13 @@ contains
subroutine write_particle_track(p)
type(Particle), intent(in) :: p
real(8) :: old_coords(3, n_tracks)
! Save the coordinates gathered in previous calls.
old_coords = coords
real(8), allocatable :: new_coords(:, :)
! Add another column to coords.
n_tracks = n_tracks + 1
deallocate(coords)
allocate(coords(3, n_tracks))
! Put the old coordinates back into the array.
coords(:, 1:n_tracks-1) = old_coords
allocate(new_coords(3, n_tracks))
new_coords(:, 1:n_tracks-1) = coords
call move_alloc(FROM=new_coords, TO=coords)
! Write current coordinates into the newest column.
coords(:, n_tracks) = p % coord0 % xyz