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

BOUNDARIES: Copy primitive variables instead of conservative.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2014-04-09 09:45:59 -03:00
parent 722d623e0b
commit 10072f5c04
1 changed files with 21 additions and 21 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

42
src/boundaries.F90
View File

@ -1390,27 +1390,27 @@ module boundaries
case(1) case(1)
if (iside == 1) then if (iside == 1) then
call boundary_copy(pdata & call boundary_copy(pdata &
, pneigh%data%u(:,iel:ie,:,:), idir, iside) , pneigh%data%q(:,iel:ie,:,:), idir, iside)
else else
call boundary_copy(pdata & call boundary_copy(pdata &
, pneigh%data%u(:,ib:ibu,:,:), idir, iside) , pneigh%data%q(:,ib:ibu,:,:), idir, iside)
end if end if
case(2) case(2)
if (iside == 1) then if (iside == 1) then
call boundary_copy(pdata & call boundary_copy(pdata &
, pneigh%data%u(:,:,jel:je,:), idir, iside) , pneigh%data%q(:,:,jel:je,:), idir, iside)
else else
call boundary_copy(pdata & call boundary_copy(pdata &
, pneigh%data%u(:,:,jb:jbu,:), idir, iside) , pneigh%data%q(:,:,jb:jbu,:), idir, iside)
end if end if
#if NDIMS == 3 #if NDIMS == 3
case(3) case(3)
if (iside == 1) then if (iside == 1) then
call boundary_copy(pdata & call boundary_copy(pdata &
, pneigh%data%u(:,:,:,kel:ke), idir, iside) , pneigh%data%q(:,:,:,kel:ke), idir, iside)
else else
call boundary_copy(pdata & call boundary_copy(pdata &
, pneigh%data%u(:,:,:,kb:kbu), idir, iside) , pneigh%data%q(:,:,:,kb:kbu), idir, iside)
end if end if
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
end select end select
@ -1539,9 +1539,9 @@ module boundaries
! fill the buffer with data from the current block (depending on the side) ! fill the buffer with data from the current block (depending on the side)
! !
if (pinfo%side == 1) then if (pinfo%side == 1) then
rbuf(l,:,:,:,:) = pinfo%neigh%data%u(:,iel:ie,:,:) rbuf(l,:,:,:,:) = pinfo%neigh%data%q(:,iel:ie,:,:)
else else
rbuf(l,:,:,:,:) = pinfo%neigh%data%u(:,ib:ibu,:,:) rbuf(l,:,:,:,:) = pinfo%neigh%data%q(:,ib:ibu,:,:)
end if end if
! associate the pointer with the next block ! associate the pointer with the next block
@ -1567,9 +1567,9 @@ module boundaries
! fill the buffer with data from the current block (depending on the side) ! fill the buffer with data from the current block (depending on the side)
! !
if (pinfo%side == 1) then if (pinfo%side == 1) then
rbuf(l,:,:,:,:) = pinfo%neigh%data%u(:,:,jel:je,:) rbuf(l,:,:,:,:) = pinfo%neigh%data%q(:,:,jel:je,:)
else else
rbuf(l,:,:,:,:) = pinfo%neigh%data%u(:,:,jb:jbu,:) rbuf(l,:,:,:,:) = pinfo%neigh%data%q(:,:,jb:jbu,:)
end if end if
! associate the pointer with the next block ! associate the pointer with the next block
@ -1596,9 +1596,9 @@ module boundaries
! fill the buffer with data from the current block (depending on the side) ! fill the buffer with data from the current block (depending on the side)
! !
if (pinfo%side == 1) then if (pinfo%side == 1) then
rbuf(l,:,:,:,:) = pinfo%neigh%data%u(:,:,:,kel:ke) rbuf(l,:,:,:,:) = pinfo%neigh%data%q(:,:,:,kel:ke)
else else
rbuf(l,:,:,:,:) = pinfo%neigh%data%u(:,:,:,kb:kbu) rbuf(l,:,:,:,:) = pinfo%neigh%data%q(:,:,:,kb:kbu)
end if end if
! associate the pointer with the next block ! associate the pointer with the next block
@ -2984,13 +2984,13 @@ module boundaries
! Arguments: ! Arguments:
! !
! pdata - the pointer to modified data block; ! pdata - the pointer to modified data block;
! u - the variable array from which boundaries are updated; ! q - the variable array from which boundaries are updated;
! idir - the direction to be processed; ! idir - the direction to be processed;
! iside - the side to be processed; ! iside - the side to be processed;
! !
!=============================================================================== !===============================================================================
! !
subroutine boundary_copy(pdata, u, idir, iside) subroutine boundary_copy(pdata, q, idir, iside)
! import external procedures and variables ! import external procedures and variables
! !
@ -3005,7 +3005,7 @@ module boundaries
! subroutine arguments ! subroutine arguments
! !
type(block_data), pointer , intent(inout) :: pdata type(block_data), pointer , intent(inout) :: pdata
real , dimension(:,:,:,:), intent(in) :: u real , dimension(:,:,:,:), intent(in) :: q
integer , intent(in) :: idir, iside integer , intent(in) :: idir, iside
! !
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
@ -3017,26 +3017,26 @@ module boundaries
case(1) case(1)
if (iside == 1) then if (iside == 1) then
pdata%u(1:nv, 1:ibl,1:jm,1:km) = u(1:nv,1:ng,1:jm,1:km) pdata%q(1:nv, 1:ibl,1:jm,1:km) = q(1:nv,1:ng,1:jm,1:km)
else else
pdata%u(1:nv,ieu:im ,1:jm,1:km) = u(1:nv,1:ng,1:jm,1:km) pdata%q(1:nv,ieu:im ,1:jm,1:km) = q(1:nv,1:ng,1:jm,1:km)
end if end if
case(2) case(2)
if (iside == 1) then if (iside == 1) then
pdata%u(1:nv,1:im, 1:jbl,1:km) = u(1:nv,1:im,1:ng,1:km) pdata%q(1:nv,1:im, 1:jbl,1:km) = q(1:nv,1:im,1:ng,1:km)
else else
pdata%u(1:nv,1:im,jeu:jm ,1:km) = u(1:nv,1:im,1:ng,1:km) pdata%q(1:nv,1:im,jeu:jm ,1:km) = q(1:nv,1:im,1:ng,1:km)
end if end if
#if NDIMS == 3 #if NDIMS == 3
case(3) case(3)
if (iside == 1) then if (iside == 1) then
pdata%u(1:nv,1:im,1:jm, 1:kbl) = u(1:nv,1:im,1:jm,1:ng) pdata%q(1:nv,1:im,1:jm, 1:kbl) = q(1:nv,1:im,1:jm,1:ng)
else else
pdata%u(1:nv,1:im,1:jm,keu:km ) = u(1:nv,1:im,1:jm,1:ng) pdata%q(1:nv,1:im,1:jm,keu:km ) = q(1:nv,1:im,1:jm,1:ng)
end if end if
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */