diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index 5e15bfc..3970de0 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -695,7 +695,7 @@ module user_problem if (.not. periodic(2)) then - if (pmeta%ymin < ymin) then + if (pmeta%ymin - dh(2) < ymin) then #if NDIMS == 3 vv(1,:) = 5.0d-01 * sum(pdata%q(ivx,i,nbm:nb,nb:ne), 1) @@ -735,7 +735,7 @@ module user_problem end if ! resistivity end if - if (ymax < pmeta%ymax) then + if (pmeta%ymax + dh(2) > ymax) then #if NDIMS == 3 vv(1,:) = 5.0d-01 * sum(pdata%q(ivx,i,ne:nep,nb:ne), 1) @@ -779,7 +779,7 @@ module user_problem #if NDIMS == 3 if (.not. periodic(3)) then - if (pmeta%zmin < zmin) then + if (pmeta%zmin - dh(3) < zmin) then vv(1,:) = 5.0d-01 * sum(pdata%q(ivx,i,nb:ne,nbm:nb), 1) vv(3,:) = 5.0d-01 * sum(pdata%q(ivz,i,nb:ne,nbm:nb), 1) @@ -807,7 +807,7 @@ module user_problem end if ! resistivity end if - if (pmeta%zmax > zmax) then + if (pmeta%zmax + dh(3) > zmax) then vv(1,:) = 5.0d-01 * sum(pdata%q(ivx,i,nb:ne,ne:nep), 1) vv(3,:) = 5.0d-01 * sum(pdata%q(ivz,i,nb:ne,ne:nep), 1)