JGonzalez/vlaplex
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Reorder parallel computation inside Z loop

Browse source
This commit is contained in:
JHendrikx 2025-02-04 13:23:40 +01:00
commit a961066dd4
1 changed files with 4 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

6
vlaplex.f90
View file

@ -268,8 +268,9 @@ program VlaPlEx
f_i_old(:,:,nv) = 0.0_dp f_i_old(:,:,nv) = 0.0_dp
sum_ni = 0.0_dp sum_ni = 0.0_dp
! Advect in the r direction ! Advect in the r direction
!$omp parallel do
do iz = 1, nz do iz = 1, nz
!$omp parallel do
do i = 1, nr do i = 1, nr
! Advect negative velocity ! Advect negative velocity
if (i < nr) then if (i < nr) then
@ -296,9 +297,10 @@ program VlaPlEx
end if end if
end do end do
!$omp end parallel do
sum_ni = sum_ni + Zave * n_i(iz,:) sum_ni = sum_ni + Zave * n_i(iz,:)
end do end do
!$omp end parallel do
! Assume quasi-neutrality to start iterating ! Assume quasi-neutrality to start iterating
n_e = 1.0_dp/nz * sum_ni n_e = 1.0_dp/nz * sum_ni