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

BOUNDARIES: Do not apply curl-free condition to reconnection boundaries.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2016-05-13 14:01:42 -03:00
parent bac3af0617
commit 1576802a48
1 changed files with 0 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
src/boundaries.F90
View File

@ -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) &