From 9f8c490ce02538caa33a56f4ee0c804dcdd53a67 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Thu, 24 Jul 2014 12:46:23 -0300 Subject: [PATCH] REFINEMENT: Extend the region for refinement criterion. Now we calculate the refinement criterion over a region which is larger than the block domain by one cell. It should allow to detect incoming gradients earlier and result in gradents entering the already refined block. Signed-off-by: Grzegorz Kowal --- src/refinement.F90 | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/src/refinement.F90 b/src/refinement.F90 index 2ba9510..02ac767 100644 --- a/src/refinement.F90 +++ b/src/refinement.F90 @@ -322,7 +322,8 @@ module refinement ! import external procedures and variables ! use blocks , only : block_data - use coordinates , only : ib, jb, kb, ie, je, ke + use coordinates , only : ibl, jbl, kbl + use coordinates , only : ieu, jeu, keu ! local variables are not implicit by default ! @@ -360,15 +361,15 @@ module refinement ! find the maximum smoothness indicator over all cells ! - do k = kb, ke + do k = kbl, keu #if NDIMS == 3 km1 = k - 1 kp1 = k + 1 #endif /* NDIMS == 3 */ - do j = jb, je + do j = jbl, jeu jm1 = j - 1 jp1 = j + 1 - do i = ib, ie + do i = ibl, ieu im1 = i - 1 ip1 = i + 1