@@ -178,7 +178,7 @@ contains
178178 !! other procedures that are necessary to setup the module.
179179 subroutine s_initialize_rhs_module () ! ---------------------------------
180180
181- integer :: i, j, k, l !< Generic loop iterators
181+ integer :: i, j, k, l, d !< Generic loop iterators
182182
183183 ! Configuring Coordinate Direction Indexes =========================
184184 ix%beg = - buff_size; iy%beg = 0 ; iz%beg = 0
@@ -308,7 +308,6 @@ contains
308308 iy%beg:iy%end, iz%beg:iz%end, 1 :sys_size))
309309
310310 end if
311-
312311 ! Allocation of dq_prim_ds_qp ======================================
313312
314313 if (any (Re_size > 0 )) then
@@ -448,7 +447,6 @@ contains
448447 end if
449448 end if
450449 end if
451-
452450 ! ==================================================================
453451
454452 ! Allocation of gm_alphaK_n =====================================
@@ -547,7 +545,6 @@ contains
547545
548546 end if
549547 end do
550-
551548 ! END: Allocation/ Association of flux_n, flux_src_n, and flux_gsrc_n ===
552549
553550 if (alt_soundspeed) then
@@ -575,7 +572,6 @@ contains
575572!$acc update device(Res, Re_idx, Re_size)
576573 end if
577574
578-
579575 ! Associating procedural pointer to the subroutine that will be
580576 ! utilized to calculate the solution of a given Riemann problem
581577 if (riemann_solver == 1 ) then
@@ -584,6 +580,7 @@ contains
584580 s_riemann_solver = > s_hllc_riemann_solver
585581 end if
586582
583+
587584 ! Associating the procedural pointer to the appropriate subroutine
588585 ! that will be utilized in the conversion to the mixture variables
589586 if (model_eqns == 1 ) then ! Gamma/ pi_inf model
@@ -597,20 +594,15 @@ contains
597594 s_convert_species_to_mixture_variables
598595 end if
599596
600- !$acc parallel loop collapse(4 ) gang vector default(present)
597+
598+ !$acc parallel loop collapse(5 ) gang vector default(present)
601599 do i = 1 , sys_size
602600 do l = startz, p - startz
603601 do k = starty, n - starty
604602 do j = startx, m - startx
605- flux_gsrc_n(1 )%vf(i)%sf(j, k, l) = 0d0
606-
607- if (n > 0 ) then
608- flux_gsrc_n(2 )%vf(i)%sf(j, k, l) = 0d0
609- end if
610-
611- if (p > 0 ) then
612- flux_gsrc_n(3 )%vf(i)%sf(j, k, l) = 0d0
613- end if
603+ do d = 1 , num_dims
604+ flux_gsrc_n(d)%vf(i)%sf(j, k, l) = 0d0 ;
605+ end do
614606 end do
615607 end do
616608 end do
0 commit comments