bug fixes

Author: wilfonba

src/common/m_helper_basic.fpp

@@ -110,10 +110,10 @@ contains
         res = f_approx_equal(var, real(nint(var), wp))
     end function f_is_integer
 
-    pure subroutine s_configure_coordinate_bounds(weno_polyn, buff_size, idwint, idwbuff, &
+    pure subroutine s_configure_coordinate_bounds(weno_polyn, igr_order, buff_size, idwint, idwbuff, &
                                                   viscous, bubbles_lagrange, m, n, p, num_dims, igr)
 
-        integer, intent(in) :: weno_polyn, m, n, p, num_dims
+        integer, intent(in) :: weno_polyn, m, n, p, num_dims, igr_order
         integer, intent(inout) :: buff_size
         type(int_bounds_info), dimension(3), intent(inout) :: idwint, idwbuff
         logical, intent(in) :: viscous, bubbles_lagrange
@@ -125,9 +125,9 @@ contains
         ! the next one
         if (igr) then
             if (viscous) then
-                buff_size = 6
+                buff_size = (igr_order - 1) / 2 + 2
             else
-                buff_size = 4
+                buff_size = (igr_order - 1) / 2 + 1
             end if
         else
             if (viscous) then

src/common/m_variables_conversion.fpp

@@ -1147,10 +1147,12 @@ contains
                         end do
                     end if
 
-                    $:GPU_LOOP(parallelism='[seq]')
-                    do i = advxb, advxe
-                        qK_prim_vf(i)%sf(j, k, l) = qK_cons_vf(i)%sf(j, k, l)
-                    end do
+                    if (.not. igr .or. num_fluids > 1) then
+                        $:GPU_LOOP(parallelism='[seq]')
+                        do i = advxb, advxe
+                            qK_prim_vf(i)%sf(j, k, l) = qK_cons_vf(i)%sf(j, k, l)
+                        end do
+                    end if
 
                     if (surface_tension) then
                         qK_prim_vf(c_idx)%sf(j, k, l) = qK_cons_vf(c_idx)%sf(j, k, l)

src/post_process/m_global_parameters.fpp

@@ -529,11 +529,15 @@ contains
             adv_idx%end = E_idx + num_fluids
 
             if (igr) then
-                sys_size = adv_idx%end - 1
+                adv_idx%beg = E_idx + 1
+                adv_idx%end = E_idx + num_fluids - 1
             else
-                sys_size = adv_idx%end
+                adv_idx%beg = E_idx + 1
+                adv_idx%end = E_idx + num_fluids
             end if
 
+            sys_size = adv_idx%end
+
             if (bubbles_euler) then
                 alf_idx = adv_idx%end
             else

src/pre_process/m_assign_variables.fpp

@@ -386,7 +386,6 @@ contains
             end do
         end if
 
-
         if (.not. igr .or. num_fluids > 1) then
             ! Volume fraction(s)
             do i = adv_idx%beg, adv_idx%end

src/pre_process/m_global_parameters.fpp

@@ -605,11 +605,15 @@ contains
             adv_idx%end = E_idx + num_fluids
 
             if (igr) then
-                sys_size = adv_idx%end - 1
+                adv_idx%beg = E_idx + 1
+                adv_idx%end = E_idx + num_fluids - 1
             else
-                sys_size = adv_idx%end
+                adv_idx%beg = E_idx + 1
+                adv_idx%end = E_idx + num_fluids
             end if
 
+            sys_size = adv_idx%end
+
             if (bubbles_euler) then
                 alf_idx = adv_idx%end
             else
@@ -870,7 +874,7 @@ contains
         chemxb = species_idx%beg
         chemxe = species_idx%end
 
-        call s_configure_coordinate_bounds(weno_polyn, buff_size, &
+        call s_configure_coordinate_bounds(weno_polyn, igr_order, buff_size, &
                                            idwint, idwbuff, viscous, &
                                            bubbles_lagrange, m, n, p, &
                                            num_dims, igr)

src/pre_process/m_mpi_proxy.fpp

@@ -46,7 +46,7 @@ contains
             & 'weno_order', 'precision', 'perturb_flow_fluid', &
             & 'perturb_sph_fluid', 'num_patches', 'thermal', 'nb', 'dist_type',&
             & 'R0_type', 'relax_model', 'num_ibs', 'n_start', 'elliptic_smoothing_iters', &
-            & 'num_bc_patches', 'mixlayer_perturb_nk' ]
+            & 'num_bc_patches', 'mixlayer_perturb_nk', 'igr_order' ]
             call MPI_BCAST(${VAR}$, 1, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr)
         #:endfor
 

src/simulation/m_global_parameters.fpp

@@ -861,15 +861,17 @@ contains
                 mom_idx%beg = cont_idx%end + 1
                 mom_idx%end = cont_idx%end + num_vels
                 E_idx = mom_idx%end + 1
-                adv_idx%beg = E_idx + 1
-                adv_idx%end = E_idx + num_fluids
 
                 if (igr) then
-                    sys_size = adv_idx%end - 1
+                    adv_idx%beg = E_idx + 1
+                    adv_idx%end = E_idx + num_fluids - 1
                 else
-                    sys_size = adv_idx%end
+                    adv_idx%beg = E_idx + 1
+                    adv_idx%end = E_idx + num_fluids
                 end if
 
+                sys_size = adv_idx%end
+
                 if (bubbles_euler) then
                     alf_idx = adv_idx%end
                 else
@@ -1211,7 +1213,7 @@ contains
             fd_number = max(1, fd_order/2)
         end if
 
-        call s_configure_coordinate_bounds(weno_polyn, buff_size, &
+        call s_configure_coordinate_bounds(weno_polyn, igr_order, buff_size, &
                                            idwint, idwbuff, viscous, &
                                            bubbles_lagrange, m, n, p, &
                                            num_dims, igr)

src/simulation/m_sim_helpers.fpp

@@ -109,7 +109,7 @@ contains
         if (igr) then
             if (num_fluids == 1) then
                 alpha_rho(1) = q_prim_vf(contxb)%sf(j, k, l)
-                alpha(1) = q_prim_vf(advxb)%sf(j, k, l)
+                alpha(1) = 1._wp
             else
                 $:GPU_LOOP(parallelism='[seq]')
                 do i = 1, num_fluids - 1