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

BOUNDARIES: Prepare block coordinates in boundaries_specific().

Browse Source 
If specific boundary condition depends on the position, it would be good
to pass the coordinates to them in order to update the data block.

Also, if boundary conditions are time dependent, we should pass time and
the time increment too.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2017-03-07 16:13:28 -03:00
parent e503c4fc0a
commit 88141891c3
1 changed files with 17 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
src/boundaries.F90
View File

@ -1056,6 +1056,7 @@ module boundaries
use blocks , only : list_meta, list_leaf
use blocks , only : ndims, nsides
use coordinates , only : im, jm, km
use coordinates , only : ax, ay, az
use equations , only : nv
#ifdef MPI
use mpitools , only : nproc
@ -1078,6 +1079,12 @@ module boundaries
! local variables
!
integer :: i, j, k, n, m
! local arrays
!
real(kind=8), dimension(im) :: x
real(kind=8), dimension(jm) :: y
real(kind=8), dimension(km) :: z
!
!-------------------------------------------------------------------------------
!
@ -1109,6 +1116,16 @@ module boundaries
if (pmeta%process == nproc) then
#endif /* MPI */
! prepare block coordinates
!
x(1:im) = pmeta%xmin + ax(pmeta%level,1:im)
y(1:jm) = pmeta%ymin + ay(pmeta%level,1:jm)
#if NDIMS == 3
z(1:km) = pmeta%zmin + az(pmeta%level,1:km)
#else /* NDIMS == 3 */
z(1:km) = 0.0d+00
#endif /* NDIMS == 3 */
#if NDIMS == 2
! iterate over all directions
!