From 1576802a48f4e9d031e396c603dbb17888cc25f4 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Fri, 13 May 2016 14:01:42 -0300 Subject: [PATCH] BOUNDARIES: Do not apply curl-free condition to reconnection boundaries. Signed-off-by: Grzegorz Kowal --- src/boundaries.F90 | 18 ------------------ 1 file changed, 18 deletions(-) diff --git a/src/boundaries.F90 b/src/boundaries.F90 index 3be0c59..cdc3e84 100644 --- a/src/boundaries.F90 +++ b/src/boundaries.F90 @@ -5144,15 +5144,6 @@ module boundaries ! qn(ivx,i,j,k) = min(0.0d+00, qn(ivx,ib,j,k)) -! get the tangential components of magnetic field from the curl-free condition -! - qn(iby,i,j,k) = qn(iby,ip2,j,k) & - + (qn(ibx,ip1,jm1,k) - qn(ibx,ip1,jp1,k)) * dxy -#if NDIMS == 3 - qn(ibz,i,j,k) = qn(ibz,ip2,j,k) & - + (qn(ibx,ip1,j,km1) - qn(ibx,ip1,j,kp1)) * dxz -#endif /* NDIMS == 3 */ - ! update the normal component of magnetic field from divergence-free condition ! qn(ibx,i,j,k) = qn(ibx,ip2,j,k) & @@ -5199,15 +5190,6 @@ module boundaries ! qn(ivx,i,j,k) = max(0.0d+00, qn(ivx,ie,j,k)) -! get the tangential components of magnetic field from the curl-free condition -! - qn(iby,i,j,k) = qn(iby,im2,j,k) & - + (qn(ibx,im1,jp1,k) - qn(ibx,im1,jm1,k)) * dxy -#if NDIMS == 3 - qn(ibz,i,j,k) = qn(ibz,im2,j,k) & - + (qn(ibx,im1,j,kp1) - qn(ibx,im1,j,km1)) * dxz -#endif /* NDIMS == 3 */ - ! update the normal component of magnetic field from divergence-free condition ! qn(ibx,i,j,k) = qn(ibx,im2,j,k) &