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

Speed up calculations

Browse source
This commit is contained in:
JHendrikx 2025-04-03 10:32:11 +02:00
commit 57b111c24d
1 changed files with 9 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

9
vlaplex.f90
View file

@ -279,6 +279,9 @@ program VlaPlEx
sum_ni = 0.0_dp sum_ni = 0.0_dp
! Advect in the r direction ! Advect in the r direction
do iz = 1, nz do iz = 1, nz
if (all(n_i(iz,:) < 1.0e-16_dp) .and. iz .NE. z_inj) then
cycle
end if
!$omp parallel do !$omp parallel do
do i = 1, nr do i = 1, nr
! Advect negative velocity ! Advect negative velocity
@ -383,6 +386,9 @@ program VlaPlEx
f_i_old = f_i f_i_old = f_i
do iz = 1, nz do iz = 1, nz
if (all(n_i(iz,:) < 1.0e-16_dp) .and. iz .NE. z_inj) then
cycle
end if
! Advect in the v direction ! Advect in the v direction
! i = 1, v<0 ! i = 1, v<0
i = 1 i = 1
@ -417,6 +423,9 @@ program VlaPlEx
! Reset values for next iteration ! Reset values for next iteration
f_i_old = f_i f_i_old = f_i
do iz = 1, nz do iz = 1, nz
if (all(n_i(iz,:) < 1.0e-16_dp) .and. iz .NE. z_inj) then
cycle
end if
fCum_i(iz,:) = fCum_i(iz,:) + f_i_old(iz,rCum_index,:) fCum_i(iz,:) = fCum_i(iz,:) + f_i_old(iz,rCum_index,:)
end do end do
! Write output ! Write output