Fixed bugs in the boundary condition of the moment array and correcte…

…d the OpenMP clause
WumingCode · Aug 8, 2024 · f9fe56e · f9fe56e
1 parent aa5eb35
commit f9fe56e
Showing 1 changed file with 8 additions and 8 deletions.
diff --git a/proj/reconnection/boundary_reconnection.f90 b/proj/reconnection/boundary_reconnection.f90
@@ -136,12 +136,12 @@ subroutine boundary_reconnection__particle_y(up,np2)
 
 !$OMP PARALLEL
 
-!$OMP WORKSHARE
+!$OMP PARALLEL WORKSHARE
        cnt(nys-1:nye+1) = 0
-!$OMP END WORKSHARE
-!$OMP WORKSHARE
+!$OMP END PARALLEL WORKSHARE
+!$OMP PARALLEL WORKSHARE
        cnt2(nys:nye) = 0
-!$OMP END WORKSHARE
+!$OMP END PARALLEL WORKSHARE
 
 !$OMP DO PRIVATE(ii,j,idim,jpos)
        do j=nys,nye
@@ -596,7 +596,7 @@ subroutine boundary_reconnection__mom(mom)
 
     do isp = 1,nsp
       !send to rank-1
-!$OMP PARALLEL DO PRIVATE(i,ii)
+!$OMP PARALLEL DO PRIVATE(i,ii,k)
       do i=nxgs-1,nxge+1
          ii = nk*(i-(nxgs-1))
          do k = 1, nk
@@ -610,7 +610,7 @@ subroutine boundary_reconnection__mom(mom)
                         ncomw,nstat,nerr)
 
 !$OMP PARALLEL
-!$OMP DO PRIVATE(i,ii)
+!$OMP DO PRIVATE(i,ii,k)
       do i=nxgs-1,nxge+1
          ii = nk*(i-(nxgs-1))
          do k = 1, nk
@@ -620,7 +620,7 @@ subroutine boundary_reconnection__mom(mom)
 !$OMP END DO NOWAIT
 
       !send to rank+1
-!$OMP DO PRIVATE(i,ii)
+!$OMP DO PRIVATE(i,ii,k)
       do i=nxgs-1,nxge+1
          ii = nk*(i-(nxgs-1))
          do k = 1, nk
@@ -634,7 +634,7 @@ subroutine boundary_reconnection__mom(mom)
                         bff_rcv(1),nk*(nxge-nxgs+3),mnpr,ndown,201, &
                         ncomw,nstat,nerr)
 
-!$OMP PARALLEL DO PRIVATE(i,ii)
+!$OMP PARALLEL DO PRIVATE(i,ii,k)
       do i=nxgs-1,nxge+1
          ii = nk*(i-(nxgs-1))
          do k = 1, nk