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

INTERPOLATIONS: Make variable initialization OpenMP conformant.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2022-01-08 10:33:43 -03:00
parent 9b5198a5ff
commit 9913468203
1 changed files with 20 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
sources/interpolations.F90
View File

@ -852,7 +852,7 @@ module interpolations
real(kind=8), dimension(:,:,:) , intent(in) :: q
real(kind=8), dimension(:,:,:,:,:), intent(out) :: qi
integer :: i , j , k = 1
integer :: i , j , k
integer :: im1, jm1, ip1, jp1
#if NDIMS == 3
integer :: km1, kp1
@ -866,6 +866,10 @@ module interpolations
!-------------------------------------------------------------------------------
!
#if NDIMS == 2
k = 1
#endif /* NDIMS == 2 */
! copy ghost zones
!
do j = 1, NDIMS
@ -982,10 +986,14 @@ module interpolations
real(kind=8), dimension(:,:,:) , intent(in) :: q
real(kind=8), dimension(:,:,:,:,:), intent(out) :: qi
integer :: i, j, k = 1
integer :: i, j, k
!-------------------------------------------------------------------------------
!
#if NDIMS == 2
k = 1
#endif /* NDIMS == 2 */
! copy ghost zones
!
do j = 1, NDIMS
@ -1085,7 +1093,7 @@ module interpolations
logical :: flag
integer :: i, il, iu, im1, ip1
integer :: j, jl, ju, jm1, jp1
integer :: k = 1
integer :: k
#if NDIMS == 3
integer :: kl, ku, km1, kp1
#endif /* NDIMS == 3 */
@ -1099,6 +1107,10 @@ module interpolations
!-------------------------------------------------------------------------------
!
#if NDIMS == 2
k = 1
#endif /* NDIMS == 2 */
! copy ghost zones
!
do j = 1, NDIMS
@ -1268,7 +1280,7 @@ module interpolations
integer :: i, im1, ip1
integer :: j, jm1, jp1
integer :: k = 1
integer :: k
#if NDIMS == 3
integer :: km1, kp1
#endif /* NDIMS == 3 */
@ -1287,6 +1299,10 @@ module interpolations
!-------------------------------------------------------------------------------
!
#if NDIMS == 2
k = 1
#endif /* NDIMS == 2 */
#if NDIMS == 3
do k = 1, nn
km1 = max( 1, k - 1)