gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Add two arrays of conserved variables to the block structure.

Browse Source 
In these two arrays we store temporary states during the Runge-Kutta
integration.  The new pointer %u in the block structure points to the
current array of the conservative variables.
This commit is contained in:
Grzegorz Kowal 2012-07-31 14:14:42 -03:00
parent b09f21351a
commit 2e26d47362
1 changed files with 15 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/blocks.F90
View File

@ -134,10 +134,14 @@ module blocks
! !
type(block_meta), pointer :: meta type(block_meta), pointer :: meta
! an allocatable array to store all conserved ! a pointer to the array conserved variables
!
real, dimension(:,:,:,:) , pointer :: u
! an allocatable arrays to store all conserved
! variables ! variables
! !
real, dimension(:,:,:,:) , allocatable :: u real, dimension(:,:,:,:) , allocatable :: u0, u1
! an allocatable array to store all primitive ! an allocatable array to store all primitive
! variables ! variables
@ -1040,12 +1044,17 @@ module blocks
! allocate space for conserved variables ! allocate space for conserved variables
! !
allocate(pdata%u(nvars,nx,ny,nz)) allocate(pdata%u0(nvars,nx,ny,nz))
allocate(pdata%u1(nvars,nx,ny,nz))
! allocate space for primitive variables ! allocate space for primitive variables
! !
allocate(pdata%q(nvars,nx,ny,nz)) allocate(pdata%q(nvars,nx,ny,nz))
! initiate the array pointer
!
pdata%u => pdata%u0
! allocate space for the numerical fluxes ! allocate space for the numerical fluxes
! !
if (nflux .gt. 0) allocate(pdata%f(ndims,nflux,nx,ny,nz)) if (nflux .gt. 0) allocate(pdata%f(ndims,nflux,nx,ny,nz))
@ -1097,7 +1106,8 @@ module blocks
! deallocate conservative variables ! deallocate conservative variables
! !
if (allocated(pdata%u)) deallocate(pdata%u) if (allocated(pdata%u0)) deallocate(pdata%u0)
if (allocated(pdata%u1)) deallocate(pdata%u1)
! deallocate primitive variables ! deallocate primitive variables
! !
@ -1109,6 +1119,7 @@ module blocks
! nullify pointers ! nullify pointers
! !
nullify(pdata%u)
nullify(pdata%next) nullify(pdata%next)
nullify(pdata%prev) nullify(pdata%prev)
nullify(pdata%meta) nullify(pdata%meta)