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

MESH: Rewrite prolong_block() and restrict_block().

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2013-12-27 15:14:31 -02:00
parent c1db19d55f
commit 24b644aa6a
1 changed files with 22 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
src/mesh.F90
View File

@ -1527,18 +1527,32 @@ module mesh
!
!===============================================================================
!
! prolong_block: subroutine expands the block data and copy them to children
! subroutine PROLONG_BLOCK:
! ------------------------
!
! Subroutine prolongs variables from a data blocks linked to the input
! meta block and copy the resulting array of variables to
! the newly created children data block. The process of data restriction
! conserves stored variables.
!
! Arguments:
!
! pblock - the input meta block;
!
!===============================================================================
!
subroutine prolong_block(pblock)
! import external procedures and variables
!
use blocks , only : block_meta, block_data, nchild
use coordinates , only : ng, nh, in, jn, kn, im, jm, km
use coordinates , only : ib, ie, jb, je, kb, ke
use equations , only : nv
use interpolations, only : limiter
! local variables are not implicit by default
!
implicit none
! input arguments
@ -1552,6 +1566,11 @@ module mesh
integer :: ic, jc, kc, ip, jp, kp
real :: dul, dur, dux, duy, duz
! local pointers
!
type(block_meta), pointer :: pchild
type(block_data), pointer :: pdata
! local arrays
!
integer, dimension(3) :: dm
@ -1560,11 +1579,6 @@ module mesh
!
real, dimension(:,:,:,:), allocatable :: u
! local pointers
!
type(block_meta), pointer :: pchild
type(block_data), pointer :: pdata
!-------------------------------------------------------------------------------
!
#ifdef PROFILE
@ -1694,7 +1708,7 @@ module mesh
pchild%data%u(1:nv,1:im,1:jm,1:km) = u(1:nv,il:iu,jl:ju,kl:ku)
#endif /* NDIMS == 3 */
end do
end do ! nchild
! deallocate local arrays
!
@ -1728,7 +1742,7 @@ module mesh
!
subroutine restrict_block(pblock)
! variables and subroutines imported from other modules
! import external procedures and variables
!
use blocks , only : ndims
use blocks , only : block_meta, block_data, nchild