1826 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf
1827 type(scalar_field),
dimension(sys_size),
intent(inout) ::
q_cons_vf
1828 integer,
intent(in) :: idir
1830 real(wp) :: rho_l, gamma_l, pi_inf_l, e_l, mu_l, f_l, pres_l
1831 real(wp) :: rho_r, gamma_r, pi_inf_r, e_r, mu_r, f_r, pres_r
1832 real(wp),
dimension(3) :: vflux_l_arr, vflux_r_arr
1833 real(wp),
dimension(-1:1) :: rho_sf_small
1834# 416 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1835 real(wp),
dimension(num_fluids_max) :: alpha_rho_l, alpha_l, alpha_r, alpha_rho_r
1836 real(wp),
dimension(3) :: vel_l, vel_r
1837 real(wp),
dimension(3, 3) :: dvel
1838 real(wp),
dimension(3) :: dvel_small
1839# 426 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1843# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1845# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1847# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1848#if defined(MFC_OpenACC)
1849# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1851# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1852#elif defined(MFC_OpenMP)
1853# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1855# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1857# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1859# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1861# 430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1863# 433 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1870# 446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1872# 446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1873#if defined(MFC_OpenACC)
1874# 446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1876# 446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1877#elif defined(MFC_OpenMP)
1878# 446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1880# 446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1885# 449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1886#if defined(MFC_OpenACC)
1887# 449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1889# 449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1890#elif defined(MFC_OpenMP)
1891# 449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1893# 449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1898# 452 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1899#if defined(MFC_OpenACC)
1900# 452 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1902# 452 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1903#elif defined(MFC_OpenMP)
1904# 452 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1906# 452 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1908 do r = 1, num_fluids
1911 rho_sf_small(
i) = rho_l
1916 &
l)/rho_sf_small(-1))
1922# 466 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1923#if defined(MFC_OpenACC)
1924# 466 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1926# 466 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1927#elif defined(MFC_OpenMP)
1928# 466 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1930# 466 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1933 dvel(
i, 1) = dvel_small(
i)
1938 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
1939 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(1))/3._wp
1942 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
1943 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(1))/3._wp
1948# 482 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1949#if defined(MFC_OpenACC)
1950# 482 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1952# 482 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1953#elif defined(MFC_OpenMP)
1954# 482 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1956# 482 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1961# 485 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1962#if defined(MFC_OpenACC)
1963# 485 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1965# 485 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1966#elif defined(MFC_OpenMP)
1967# 485 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1969# 485 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1971 do r = 1, num_fluids
1974 rho_sf_small(
i) = rho_l
1984# 498 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1985#if defined(MFC_OpenACC)
1986# 498 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1988# 498 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1989#elif defined(MFC_OpenMP)
1990# 498 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1992# 498 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1995 dvel(
i, 2) = dvel_small(
i)
2000 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
2001 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(2))/3._wp
2004 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
2005 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(2))/3._wp
2010 & 1)**2._wp + dvel(2, 2)**2._wp + (dvel(1, 1) + dvel(2, 2))**2._wp), kind=stp)
2015# 519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2016#if defined(MFC_OpenACC)
2017# 519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2019# 519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2020#elif defined(MFC_OpenMP)
2021# 519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2023# 519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2025 do i = 1, num_fluids
2026 alpha_rho_l(
i) = 0._wp
2027 alpha_rho_r(
i) = 0._wp
2032# 526 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2033#if defined(MFC_OpenACC)
2034# 526 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2036# 526 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2037#elif defined(MFC_OpenMP)
2038# 526 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2040# 526 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2048# 532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2049#if defined(MFC_OpenACC)
2050# 532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2052# 532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2053#elif defined(MFC_OpenMP)
2054# 532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2056# 532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2060# 534 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2061#if defined(MFC_OpenACC)
2062# 534 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2064# 534 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2065#elif defined(MFC_OpenMP)
2066# 534 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2068# 534 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2070 do i = 1, num_fluids
2074 if (num_fluids > 1)
then
2076# 540 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2077#if defined(MFC_OpenACC)
2078# 540 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2080# 540 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2081#elif defined(MFC_OpenMP)
2082# 540 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2084# 540 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2086 do i = 1, num_fluids - 1
2094# 548 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2095#if defined(MFC_OpenACC)
2096# 548 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2098# 548 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2099#elif defined(MFC_OpenMP)
2100# 548 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2102# 548 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2110# 554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2111#if defined(MFC_OpenACC)
2112# 554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2114# 554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2115#elif defined(MFC_OpenMP)
2116# 554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2118# 554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2122# 556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2123#if defined(MFC_OpenACC)
2124# 556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2126# 556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2127#elif defined(MFC_OpenMP)
2128# 556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2130# 556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2132 do i = 1, num_fluids
2136 if (num_fluids > 1)
then
2138# 562 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2139#if defined(MFC_OpenACC)
2140# 562 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2142# 562 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2143#elif defined(MFC_OpenMP)
2144# 562 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2146# 562 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2148 do i = 1, num_fluids - 1
2156# 570 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2157#if defined(MFC_OpenACC)
2158# 570 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2160# 570 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2161#elif defined(MFC_OpenMP)
2162# 570 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2164# 570 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2171 if (num_fluids > 1)
then
2172 alpha_l(num_fluids) = 1._wp
2173 alpha_r(num_fluids) = 1._wp
2176# 580 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2177#if defined(MFC_OpenACC)
2178# 580 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2180# 580 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2181#elif defined(MFC_OpenMP)
2182# 580 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2184# 580 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2186 do i = 1, num_fluids - 1
2187 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
2188 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
2192 rho_l = 0._wp; rho_r = 0._wp
2193 gamma_l = 0._wp; gamma_r = 0._wp
2194 pi_inf_l = 0._wp; pi_inf_r = 0._wp
2197# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2198#if defined(MFC_OpenACC)
2199# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2201# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2202#elif defined(MFC_OpenMP)
2203# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2205# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2207 do i = 1, num_fluids
2208 rho_l = rho_l + alpha_rho_l(
i)
2209 gamma_l = gamma_l + alpha_l(
i)*gammas(
i)
2210 pi_inf_l = pi_inf_l + alpha_l(
i)*pi_infs(
i)
2212 rho_r = rho_r + alpha_rho_r(
i)
2213 gamma_r = gamma_r + alpha_r(
i)*gammas(
i)
2214 pi_inf_r = pi_inf_r + alpha_r(
i)*pi_infs(
i)
2218# 602 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2219#if defined(MFC_OpenACC)
2220# 602 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2222# 602 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2223#elif defined(MFC_OpenMP)
2224# 602 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2226# 602 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2229 vel_l(
i) = vel_l(
i)/rho_l
2230 vel_r(
i) = vel_r(
i)/rho_r
2234 mu_l = 0._wp; mu_r = 0._wp
2236# 610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2237#if defined(MFC_OpenACC)
2238# 610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2240# 610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2241#elif defined(MFC_OpenMP)
2242# 610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2244# 610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2246 do i = 1, num_fluids
2252# 616 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2253#if defined(MFC_OpenACC)
2254# 616 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2256# 616 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2257#elif defined(MFC_OpenMP)
2258# 616 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2260# 616 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2263 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
2265# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2266#if defined(MFC_OpenACC)
2267# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2269# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2270#elif defined(MFC_OpenMP)
2271# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2273# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2276 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j + 1)), kind=stp)
2279# 623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2280#if defined(MFC_OpenACC)
2281# 623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2283# 623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2284#elif defined(MFC_OpenMP)
2285# 623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2287# 623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2290 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j)), kind=stp)
2292# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2293#if defined(MFC_OpenACC)
2294# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2296# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2297#elif defined(MFC_OpenMP)
2298# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2300# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2303 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j)), kind=stp)
2306# 630 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2307#if defined(MFC_OpenACC)
2308# 630 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2310# 630 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2311#elif defined(MFC_OpenMP)
2312# 630 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2314# 630 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2317 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
2319# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2320#if defined(MFC_OpenACC)
2321# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2323# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2324#elif defined(MFC_OpenMP)
2325# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2327# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2330 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j + 1)), kind=stp)
2333# 637 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2334#if defined(MFC_OpenACC)
2335# 637 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2337# 637 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2338#elif defined(MFC_OpenMP)
2339# 637 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2341# 637 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2344 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j)), kind=stp)
2346# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2347#if defined(MFC_OpenACC)
2348# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2350# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2351#elif defined(MFC_OpenMP)
2352# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2354# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2357 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j)), kind=stp)
2360# 644 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2361#if defined(MFC_OpenACC)
2362# 644 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2364# 644 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2365#elif defined(MFC_OpenMP)
2366# 644 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2368# 644 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2371 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
2373# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2374#if defined(MFC_OpenACC)
2375# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2377# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2378#elif defined(MFC_OpenMP)
2379# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2381# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2384 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j + 1)), kind=stp)
2387# 651 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2388#if defined(MFC_OpenACC)
2389# 651 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2391# 651 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2392#elif defined(MFC_OpenMP)
2393# 651 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2395# 651 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2398 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j)), kind=stp)
2400# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2401#if defined(MFC_OpenACC)
2402# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2404# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2405#elif defined(MFC_OpenMP)
2406# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2408# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2411 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j)), kind=stp)
2414# 658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2415#if defined(MFC_OpenACC)
2416# 658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2418# 658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2419#elif defined(MFC_OpenMP)
2420# 658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2422# 658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2425 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
2427# 661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2428#if defined(MFC_OpenACC)
2429# 661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2431# 661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2432#elif defined(MFC_OpenMP)
2433# 661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2435# 661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2438 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j + 1)), kind=stp)
2441# 665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2442#if defined(MFC_OpenACC)
2443# 665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2445# 665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2446#elif defined(MFC_OpenMP)
2447# 665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2449# 665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2452 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j)), kind=stp)
2454# 668 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2455#if defined(MFC_OpenACC)
2456# 668 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2458# 668 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2459#elif defined(MFC_OpenMP)
2460# 668 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2462# 668 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2465 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j)), kind=stp)
2468 e_l = 0._wp; e_r = 0._wp
2471# 675 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2472#if defined(MFC_OpenACC)
2473# 675 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2475# 675 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2476#elif defined(MFC_OpenMP)
2477# 675 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2479# 675 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2486# 680 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2487#if defined(MFC_OpenACC)
2488# 680 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2490# 680 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2491#elif defined(MFC_OpenMP)
2492# 680 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2494# 680 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2500 call s_get_derived_states(e_l, gamma_l, pi_inf_l, rho_l, vel_l, e_r, gamma_r, pi_inf_r, rho_r, &
2501 & vel_r, pres_l, pres_r, cfl)
2503 do i = 1, num_fluids
2505# 689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2506#if defined(MFC_OpenACC)
2507# 689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2509# 689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2510#elif defined(MFC_OpenMP)
2511# 689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2513# 689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2515 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
2516 &
l) + real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
2517 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
2520# 694 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2521#if defined(MFC_OpenACC)
2522# 694 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2524# 694 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2525#elif defined(MFC_OpenMP)
2526# 694 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2528# 694 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2530 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
2531 &
l) - real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
2532 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j))), kind=stp)
2535 if (num_fluids > 1)
then
2537# 701 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2538#if defined(MFC_OpenACC)
2539# 701 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2541# 701 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2542#elif defined(MFC_OpenMP)
2543# 701 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2545# 701 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2547 do i = 1, num_fluids - 1
2549# 703 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2550#if defined(MFC_OpenACC)
2551# 703 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2553# 703 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2554#elif defined(MFC_OpenMP)
2555# 703 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2557# 703 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2560 &
l) + real((0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
2561 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
2564# 708 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2565#if defined(MFC_OpenACC)
2566# 708 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2568# 708 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2569#elif defined(MFC_OpenMP)
2570# 708 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2572# 708 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2576 &
l)*vel_l(1)*(1._wp/dx(
j + 1))), kind=stp)
2579# 713 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2580#if defined(MFC_OpenACC)
2581# 713 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2583# 713 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2584#elif defined(MFC_OpenMP)
2585# 713 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2587# 713 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2590 &
l) - real((0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
2591 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j))), kind=stp)
2594# 718 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2595#if defined(MFC_OpenACC)
2596# 718 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2598# 718 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2599#elif defined(MFC_OpenMP)
2600# 718 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2602# 718 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2606 &
l)*vel_l(1)*(1._wp/dx(
j))), kind=stp)
2611# 725 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2612#if defined(MFC_OpenACC)
2613# 725 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2615# 725 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2616#elif defined(MFC_OpenMP)
2617# 725 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2619# 725 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2622 &
l) + real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j + 1)) &
2623 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j + 1))), kind=stp)
2626# 730 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2627#if defined(MFC_OpenACC)
2628# 730 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2630# 730 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2631#elif defined(MFC_OpenMP)
2632# 730 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2634# 730 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2637 &
l) + real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j + 1)) &
2638 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j + 1))), kind=stp)
2641# 735 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2642#if defined(MFC_OpenACC)
2643# 735 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2645# 735 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2646#elif defined(MFC_OpenMP)
2647# 735 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2649# 735 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2652 &
l) + real((0.5_wp*dt*(vel_l(1)*(e_l + pres_l))*(1._wp/dx(
j + 1)) - 0.5_wp*dt*cfl*(e_l) &
2653 & *(1._wp/dx(
j + 1))), kind=stp)
2656# 740 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2657#if defined(MFC_OpenACC)
2658# 740 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2660# 740 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2661#elif defined(MFC_OpenMP)
2662# 740 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2664# 740 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2667 &
l) - real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j)) &
2668 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j))), kind=stp)
2671# 745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2672#if defined(MFC_OpenACC)
2673# 745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2675# 745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2676#elif defined(MFC_OpenMP)
2677# 745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2679# 745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2682 &
l) - real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j)) &
2683 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j))), kind=stp)
2686# 750 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2687#if defined(MFC_OpenACC)
2688# 750 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2690# 750 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2691#elif defined(MFC_OpenMP)
2692# 750 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2694# 750 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2697 &
l) - real((0.5_wp*dt*(vel_l(1)*(e_l + pres_l))*(1._wp/dx(
j)) - 0.5_wp*dt*cfl*(e_l) &
2698 & *(1._wp/dx(
j))), kind=stp)
2701# 755 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2702#if defined(MFC_OpenACC)
2703# 755 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2705# 755 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2706#elif defined(MFC_OpenMP)
2707# 755 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2709# 755 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2711 do i = 1, num_fluids
2713# 757 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2714#if defined(MFC_OpenACC)
2715# 757 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2717# 757 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2718#elif defined(MFC_OpenMP)
2719# 757 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2721# 757 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2723 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
2724 &
l) + real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
2725 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
2728# 762 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2729#if defined(MFC_OpenACC)
2730# 762 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2732# 762 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2733#elif defined(MFC_OpenMP)
2734# 762 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2736# 762 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2738 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
2739 &
l) - real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
2740 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j))), kind=stp)
2743 if (num_fluids > 1)
then
2745# 769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2746#if defined(MFC_OpenACC)
2747# 769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2749# 769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2750#elif defined(MFC_OpenMP)
2751# 769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2753# 769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2755 do i = 1, num_fluids - 1
2757# 771 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2758#if defined(MFC_OpenACC)
2759# 771 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2761# 771 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2762#elif defined(MFC_OpenMP)
2763# 771 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2765# 771 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2768 &
l) + real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
2769 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
2772# 776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2773#if defined(MFC_OpenACC)
2774# 776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2776# 776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2777#elif defined(MFC_OpenMP)
2778# 776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2780# 776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2784 &
l)*vel_r(1)*(1._wp/dx(
j + 1))), kind=stp)
2787# 781 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2788#if defined(MFC_OpenACC)
2789# 781 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2791# 781 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2792#elif defined(MFC_OpenMP)
2793# 781 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2795# 781 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2798 &
l) - real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
2799 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j))), kind=stp)
2802# 786 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2803#if defined(MFC_OpenACC)
2804# 786 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2806# 786 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2807#elif defined(MFC_OpenMP)
2808# 786 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2810# 786 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2814 &
l)*vel_r(1)*(1._wp/dx(
j))), kind=stp)
2819# 793 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2820#if defined(MFC_OpenACC)
2821# 793 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2823# 793 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2824#elif defined(MFC_OpenMP)
2825# 793 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2827# 793 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2830 &
l) + real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j + 1)) &
2831 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j + 1))), kind=stp)
2834# 798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2835#if defined(MFC_OpenACC)
2836# 798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2838# 798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2839#elif defined(MFC_OpenMP)
2840# 798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2842# 798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2845 &
l) + real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j + 1)) &
2846 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j + 1))), kind=stp)
2849# 803 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2850#if defined(MFC_OpenACC)
2851# 803 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2853# 803 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2854#elif defined(MFC_OpenMP)
2855# 803 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2857# 803 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2860 &
l) + real((0.5_wp*dt*(vel_r(1)*(e_r + pres_r))*(1._wp/dx(
j + 1)) + 0.5_wp*dt*cfl*(e_r) &
2861 & *(1._wp/dx(
j + 1))), kind=stp)
2864# 808 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2865#if defined(MFC_OpenACC)
2866# 808 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2868# 808 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2869#elif defined(MFC_OpenMP)
2870# 808 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2872# 808 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2875 &
l) - real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j)) &
2876 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j))), kind=stp)
2879# 813 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2880#if defined(MFC_OpenACC)
2881# 813 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2883# 813 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2884#elif defined(MFC_OpenMP)
2885# 813 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2887# 813 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2890 &
l) - real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j)) &
2891 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j))), kind=stp)
2894# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2895#if defined(MFC_OpenACC)
2896# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2898# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2899#elif defined(MFC_OpenMP)
2900# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2902# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2905 &
l) - real((0.5_wp*dt*(vel_r(1)*(e_r + pres_r))*(1._wp/dx(
j)) + 0.5_wp*dt*cfl*(e_r) &
2906 & *(1._wp/dx(
j))), kind=stp)
2911# 825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2912#if defined(MFC_OpenACC)
2913# 825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2915# 825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2916#elif defined(MFC_OpenMP)
2917# 825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2919# 825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2921# 825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2923# 827 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2925# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2927# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2929# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2930#if defined(MFC_OpenACC)
2931# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2933# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2934#elif defined(MFC_OpenMP)
2935# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2937# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2939# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2941# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2943# 829 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2945# 832 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2952# 845 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2954# 845 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2955#if defined(MFC_OpenACC)
2956# 845 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2958# 845 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2959#elif defined(MFC_OpenMP)
2960# 845 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2962# 845 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2967# 848 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2968#if defined(MFC_OpenACC)
2969# 848 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2971# 848 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2972#elif defined(MFC_OpenMP)
2973# 848 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2975# 848 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2980# 851 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2981#if defined(MFC_OpenACC)
2982# 851 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2984# 851 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2985#elif defined(MFC_OpenMP)
2986# 851 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2988# 851 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2990 do r = 1, num_fluids
2993 rho_sf_small(
i) = rho_l
3005# 866 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3006#if defined(MFC_OpenACC)
3007# 866 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3009# 866 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3010#elif defined(MFC_OpenMP)
3011# 866 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3013# 866 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3016 dvel(
i, 1) = dvel_small(
i)
3021 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
3022 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(3))
3023 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(1))/3._wp
3026 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
3027 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(3))
3028 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(1))/3._wp
3033# 884 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3034#if defined(MFC_OpenACC)
3035# 884 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3037# 884 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3038#elif defined(MFC_OpenMP)
3039# 884 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3041# 884 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3046# 887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3047#if defined(MFC_OpenACC)
3048# 887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3050# 887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3051#elif defined(MFC_OpenMP)
3052# 887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3054# 887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3056 do r = 1, num_fluids
3059 rho_sf_small(
i) = rho_l
3070# 901 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3071#if defined(MFC_OpenACC)
3072# 901 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3074# 901 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3075#elif defined(MFC_OpenMP)
3076# 901 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3078# 901 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3081 dvel(
i, 2) = dvel_small(
i)
3086 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
3087 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(2))/3._wp
3090 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
3091 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(2))/3._wp
3096# 917 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3097#if defined(MFC_OpenACC)
3098# 917 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3100# 917 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3101#elif defined(MFC_OpenMP)
3102# 917 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3104# 917 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3109# 920 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3110#if defined(MFC_OpenACC)
3111# 920 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3113# 920 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3114#elif defined(MFC_OpenMP)
3115# 920 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3117# 920 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3119 do r = 1, num_fluids
3122 rho_sf_small(
i) = rho_l
3131 &
l - 1)/rho_sf_small(-1))
3134# 935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3135#if defined(MFC_OpenACC)
3136# 935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3138# 935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3139#elif defined(MFC_OpenMP)
3140# 935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3142# 935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3145 dvel(
i, 3) = dvel_small(
i)
3150 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(1))
3151 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(3))/3._wp
3154 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(1))
3155 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(3))/3._wp
3159 jac_rhs(
j,
k,
l) = real(
alf_igr*(2._wp*(dvel(1, 2)*dvel(2, 1) + dvel(1, 3)*dvel(3, &
3160 & 1) + dvel(2, 3)*dvel(3, 2)) + dvel(1, 1)**2._wp + dvel(2, 2)**2._wp + dvel(3, &
3161 & 3)**2._wp + (dvel(1, 1) + dvel(2, 2) + dvel(3, 3))**2._wp), kind=stp)
3166# 957 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3167#if defined(MFC_OpenACC)
3168# 957 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3170# 957 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3171#elif defined(MFC_OpenMP)
3172# 957 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3174# 957 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3176 do i = 1, num_fluids
3177 alpha_rho_l(
i) = 0._wp
3178 alpha_rho_r(
i) = 0._wp
3183# 964 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3184#if defined(MFC_OpenACC)
3185# 964 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3187# 964 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3188#elif defined(MFC_OpenMP)
3189# 964 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3191# 964 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3199# 970 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3200#if defined(MFC_OpenACC)
3201# 970 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3203# 970 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3204#elif defined(MFC_OpenMP)
3205# 970 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3207# 970 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3211# 972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3212#if defined(MFC_OpenACC)
3213# 972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3215# 972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3216#elif defined(MFC_OpenMP)
3217# 972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3219# 972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3221 do i = 1, num_fluids
3225 if (num_fluids > 1)
then
3227# 978 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3228#if defined(MFC_OpenACC)
3229# 978 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3231# 978 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3232#elif defined(MFC_OpenMP)
3233# 978 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3235# 978 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3237 do i = 1, num_fluids - 1
3245# 986 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3246#if defined(MFC_OpenACC)
3247# 986 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3249# 986 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3250#elif defined(MFC_OpenMP)
3251# 986 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3253# 986 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3261# 992 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3262#if defined(MFC_OpenACC)
3263# 992 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3265# 992 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3266#elif defined(MFC_OpenMP)
3267# 992 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3269# 992 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3273# 994 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3274#if defined(MFC_OpenACC)
3275# 994 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3277# 994 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3278#elif defined(MFC_OpenMP)
3279# 994 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3281# 994 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3283 do i = 1, num_fluids
3287 if (num_fluids > 1)
then
3289# 1000 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3290#if defined(MFC_OpenACC)
3291# 1000 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3293# 1000 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3294#elif defined(MFC_OpenMP)
3295# 1000 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3297# 1000 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3299 do i = 1, num_fluids - 1
3307# 1008 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3308#if defined(MFC_OpenACC)
3309# 1008 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3311# 1008 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3312#elif defined(MFC_OpenMP)
3313# 1008 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3315# 1008 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3322 if (num_fluids > 1)
then
3323 alpha_l(num_fluids) = 1._wp
3324 alpha_r(num_fluids) = 1._wp
3327# 1018 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3328#if defined(MFC_OpenACC)
3329# 1018 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3331# 1018 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3332#elif defined(MFC_OpenMP)
3333# 1018 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3335# 1018 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3337 do i = 1, num_fluids - 1
3338 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
3339 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
3343 rho_l = 0._wp; rho_r = 0._wp
3344 gamma_l = 0._wp; gamma_r = 0._wp
3345 pi_inf_l = 0._wp; pi_inf_r = 0._wp
3348# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3349#if defined(MFC_OpenACC)
3350# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3352# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3353#elif defined(MFC_OpenMP)
3354# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3356# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3358 do i = 1, num_fluids
3359 rho_l = rho_l + alpha_rho_l(
i)
3360 gamma_l = gamma_l + alpha_l(
i)*gammas(
i)
3361 pi_inf_l = pi_inf_l + alpha_l(
i)*pi_infs(
i)
3363 rho_r = rho_r + alpha_rho_r(
i)
3364 gamma_r = gamma_r + alpha_r(
i)*gammas(
i)
3365 pi_inf_r = pi_inf_r + alpha_r(
i)*pi_infs(
i)
3369# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3370#if defined(MFC_OpenACC)
3371# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3373# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3374#elif defined(MFC_OpenMP)
3375# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3377# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3380 vel_l(
i) = vel_l(
i)/rho_l
3381 vel_r(
i) = vel_r(
i)/rho_r
3388# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3389#if defined(MFC_OpenACC)
3390# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3392# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3393#elif defined(MFC_OpenMP)
3394# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3396# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3398 do i = 1, num_fluids
3404# 1055 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3405#if defined(MFC_OpenACC)
3406# 1055 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3408# 1055 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3409#elif defined(MFC_OpenMP)
3410# 1055 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3412# 1055 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3415 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
3417# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3418#if defined(MFC_OpenACC)
3419# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3421# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3422#elif defined(MFC_OpenMP)
3423# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3425# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3428 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j + 1)), kind=stp)
3431# 1062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3432#if defined(MFC_OpenACC)
3433# 1062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3435# 1062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3436#elif defined(MFC_OpenMP)
3437# 1062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3439# 1062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3442 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j)), kind=stp)
3444# 1065 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3445#if defined(MFC_OpenACC)
3446# 1065 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3448# 1065 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3449#elif defined(MFC_OpenMP)
3450# 1065 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3452# 1065 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3455 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j)), kind=stp)
3458# 1069 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3459#if defined(MFC_OpenACC)
3460# 1069 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3462# 1069 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3463#elif defined(MFC_OpenMP)
3464# 1069 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3466# 1069 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3469 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
3471# 1072 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3472#if defined(MFC_OpenACC)
3473# 1072 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3475# 1072 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3476#elif defined(MFC_OpenMP)
3477# 1072 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3479# 1072 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3482 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j + 1)), kind=stp)
3485# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3486#if defined(MFC_OpenACC)
3487# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3489# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3490#elif defined(MFC_OpenMP)
3491# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3493# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3496 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j)), kind=stp)
3498# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3499#if defined(MFC_OpenACC)
3500# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3502# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3503#elif defined(MFC_OpenMP)
3504# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3506# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3509 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j)), kind=stp)
3512# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3513#if defined(MFC_OpenACC)
3514# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3516# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3517#elif defined(MFC_OpenMP)
3518# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3520# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3523 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dx(
j + 1)), kind=stp)
3525# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3526#if defined(MFC_OpenACC)
3527# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3529# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3530#elif defined(MFC_OpenMP)
3531# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3533# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3536 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dx(
j + 1)), kind=stp)
3539# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3540#if defined(MFC_OpenACC)
3541# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3543# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3544#elif defined(MFC_OpenMP)
3545# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3547# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3550 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dx(
j)), kind=stp)
3552# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3553#if defined(MFC_OpenACC)
3554# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3556# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3557#elif defined(MFC_OpenMP)
3558# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3560# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3563 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dx(
j)), kind=stp)
3566# 1097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3567#if defined(MFC_OpenACC)
3568# 1097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3570# 1097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3571#elif defined(MFC_OpenMP)
3572# 1097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3574# 1097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3577 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dx(
j + 1)), kind=stp)
3579# 1100 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3580#if defined(MFC_OpenACC)
3581# 1100 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3583# 1100 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3584#elif defined(MFC_OpenMP)
3585# 1100 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3587# 1100 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3590 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dx(
j + 1)), kind=stp)
3593# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3594#if defined(MFC_OpenACC)
3595# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3597# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3598#elif defined(MFC_OpenMP)
3599# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3601# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3604 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dx(
j)), kind=stp)
3606# 1107 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3607#if defined(MFC_OpenACC)
3608# 1107 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3610# 1107 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3611#elif defined(MFC_OpenMP)
3612# 1107 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3614# 1107 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3617 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dx(
j)), kind=stp)
3620# 1111 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3621#if defined(MFC_OpenACC)
3622# 1111 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3624# 1111 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3625#elif defined(MFC_OpenMP)
3626# 1111 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3628# 1111 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3631 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
3633# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3634#if defined(MFC_OpenACC)
3635# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3637# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3638#elif defined(MFC_OpenMP)
3639# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3641# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3644 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j + 1)), kind=stp)
3647# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3648#if defined(MFC_OpenACC)
3649# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3651# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3652#elif defined(MFC_OpenMP)
3653# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3655# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3658 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j)), kind=stp)
3660# 1121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3661#if defined(MFC_OpenACC)
3662# 1121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3664# 1121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3665#elif defined(MFC_OpenMP)
3666# 1121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3668# 1121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3671 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j)), kind=stp)
3674# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3675#if defined(MFC_OpenACC)
3676# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3678# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3679#elif defined(MFC_OpenMP)
3680# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3682# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3685 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
3687# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3688#if defined(MFC_OpenACC)
3689# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3691# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3692#elif defined(MFC_OpenMP)
3693# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3695# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3698 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j + 1)), kind=stp)
3701# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3702#if defined(MFC_OpenACC)
3703# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3705# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3706#elif defined(MFC_OpenMP)
3707# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3709# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3712 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j)), kind=stp)
3714# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3715#if defined(MFC_OpenACC)
3716# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3718# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3719#elif defined(MFC_OpenMP)
3720# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3722# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3725 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j)), kind=stp)
3728 e_l = 0._wp; e_r = 0._wp
3731# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3732#if defined(MFC_OpenACC)
3733# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3735# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3736#elif defined(MFC_OpenMP)
3737# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3739# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3746# 1147 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3747#if defined(MFC_OpenACC)
3748# 1147 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3750# 1147 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3751#elif defined(MFC_OpenMP)
3752# 1147 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3754# 1147 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3760 call s_get_derived_states(e_l, gamma_l, pi_inf_l, rho_l, vel_l, e_r, gamma_r, pi_inf_r, rho_r, &
3761 & vel_r, pres_l, pres_r, cfl)
3764# 1155 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3765#if defined(MFC_OpenACC)
3766# 1155 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3768# 1155 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3769#elif defined(MFC_OpenMP)
3770# 1155 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3772# 1155 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3774 do i = 1, num_fluids
3776# 1157 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3777#if defined(MFC_OpenACC)
3778# 1157 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3780# 1157 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3781#elif defined(MFC_OpenMP)
3782# 1157 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3784# 1157 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3786 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
3787 &
l) + real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
3788 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
3791# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3792#if defined(MFC_OpenACC)
3793# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3795# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3796#elif defined(MFC_OpenMP)
3797# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3799# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3801 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
3802 &
l) - real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
3803 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j))), kind=stp)
3806 if (num_fluids > 1)
then
3808# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3809#if defined(MFC_OpenACC)
3810# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3812# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3813#elif defined(MFC_OpenMP)
3814# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3816# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3818 do i = 1, num_fluids - 1
3820# 1171 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3821#if defined(MFC_OpenACC)
3822# 1171 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3824# 1171 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3825#elif defined(MFC_OpenMP)
3826# 1171 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3828# 1171 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3831 &
l) + real((0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
3832 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
3835# 1176 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3836#if defined(MFC_OpenACC)
3837# 1176 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3839# 1176 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3840#elif defined(MFC_OpenMP)
3841# 1176 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3843# 1176 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3847 &
l)*vel_l(1)*(1._wp/dx(
j + 1)), kind=stp)
3850# 1181 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3851#if defined(MFC_OpenACC)
3852# 1181 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3854# 1181 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3855#elif defined(MFC_OpenMP)
3856# 1181 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3858# 1181 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3861 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
3862 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j)), kind=stp)
3865# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3866#if defined(MFC_OpenACC)
3867# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3869# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3870#elif defined(MFC_OpenMP)
3871# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3873# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3877 &
l)*vel_l(1)*(1._wp/dx(
j)), kind=stp)
3882# 1193 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3883#if defined(MFC_OpenACC)
3884# 1193 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3886# 1193 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3887#elif defined(MFC_OpenMP)
3888# 1193 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3890# 1193 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3893 &
l) + real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j + 1)) &
3894 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j + 1))), kind=stp)
3897# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3898#if defined(MFC_OpenACC)
3899# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3901# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3902#elif defined(MFC_OpenMP)
3903# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3905# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3908 &
l) + real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j + 1)) &
3909 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j + 1))), kind=stp)
3912# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3913#if defined(MFC_OpenACC)
3914# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3916# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3917#elif defined(MFC_OpenMP)
3918# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3920# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3923 &
l) + real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(3)*(1._wp/dx(
j + 1)) &
3924 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dx(
j + 1))), kind=stp)
3927# 1208 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3928#if defined(MFC_OpenACC)
3929# 1208 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3931# 1208 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3932#elif defined(MFC_OpenMP)
3933# 1208 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3935# 1208 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3938 &
l) + real((0.5_wp*dt*(vel_l(1)*(e_l + pres_l))*(1._wp/dx(
j + 1)) - 0.5_wp*dt*cfl*(e_l) &
3939 & *(1._wp/dx(
j + 1))), kind=stp)
3942# 1213 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3943#if defined(MFC_OpenACC)
3944# 1213 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3946# 1213 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3947#elif defined(MFC_OpenMP)
3948# 1213 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3950# 1213 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3953 &
l) - real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j)) &
3954 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j))), kind=stp)
3957# 1218 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3958#if defined(MFC_OpenACC)
3959# 1218 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3961# 1218 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3962#elif defined(MFC_OpenMP)
3963# 1218 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3965# 1218 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3968 &
l) - real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j)) &
3969 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j))), kind=stp)
3972# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3973#if defined(MFC_OpenACC)
3974# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3976# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3977#elif defined(MFC_OpenMP)
3978# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3980# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3983 &
l) - real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(3)*(1._wp/dx(
j)) &
3984 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dx(
j))), kind=stp)
3987# 1228 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3988#if defined(MFC_OpenACC)
3989# 1228 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3991# 1228 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3992#elif defined(MFC_OpenMP)
3993# 1228 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3995# 1228 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3998 &
l) - real((0.5_wp*dt*(vel_l(1)*(e_l + pres_l))*(1._wp/dx(
j)) - 0.5_wp*dt*cfl*(e_l) &
3999 & *(1._wp/dx(
j))), kind=stp)
4002# 1233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4003#if defined(MFC_OpenACC)
4004# 1233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4006# 1233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4007#elif defined(MFC_OpenMP)
4008# 1233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4010# 1233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4012 do i = 1, num_fluids
4014# 1235 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4015#if defined(MFC_OpenACC)
4016# 1235 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4018# 1235 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4019#elif defined(MFC_OpenMP)
4020# 1235 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4022# 1235 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4024 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
4025 &
l) + real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
4026 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
4029# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4030#if defined(MFC_OpenACC)
4031# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4033# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4034#elif defined(MFC_OpenMP)
4035# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4037# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4039 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
4040 &
l) - real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
4041 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j))), kind=stp)
4044 if (num_fluids > 1)
then
4046# 1247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4047#if defined(MFC_OpenACC)
4048# 1247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4050# 1247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4051#elif defined(MFC_OpenMP)
4052# 1247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4054# 1247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4056 do i = 1, num_fluids - 1
4058# 1249 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4059#if defined(MFC_OpenACC)
4060# 1249 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4062# 1249 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4063#elif defined(MFC_OpenMP)
4064# 1249 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4066# 1249 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4069 &
l) + real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
4070 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
4073# 1254 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4074#if defined(MFC_OpenACC)
4075# 1254 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4077# 1254 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4078#elif defined(MFC_OpenMP)
4079# 1254 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4081# 1254 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4085 &
l)*vel_r(1)*(1._wp/dx(
j + 1))), kind=stp)
4088# 1259 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4089#if defined(MFC_OpenACC)
4090# 1259 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4092# 1259 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4093#elif defined(MFC_OpenMP)
4094# 1259 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4096# 1259 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4099 &
l) - real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
4100 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j))), kind=stp)
4103# 1264 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4104#if defined(MFC_OpenACC)
4105# 1264 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4107# 1264 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4108#elif defined(MFC_OpenMP)
4109# 1264 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4111# 1264 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4115 &
l)*vel_r(1)*(1._wp/dx(
j))), kind=stp)
4120# 1271 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4121#if defined(MFC_OpenACC)
4122# 1271 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4124# 1271 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4125#elif defined(MFC_OpenMP)
4126# 1271 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4128# 1271 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4131 &
l) + real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j + 1)) &
4132 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j + 1))), kind=stp)
4135# 1276 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4136#if defined(MFC_OpenACC)
4137# 1276 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4139# 1276 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4140#elif defined(MFC_OpenMP)
4141# 1276 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4143# 1276 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4146 &
l) + real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j + 1)) &
4147 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j + 1))), kind=stp)
4150# 1281 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4151#if defined(MFC_OpenACC)
4152# 1281 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4154# 1281 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4155#elif defined(MFC_OpenMP)
4156# 1281 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4158# 1281 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4161 &
l) + real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(3)*(1._wp/dx(
j + 1)) &
4162 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dx(
j + 1))), kind=stp)
4165# 1286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4166#if defined(MFC_OpenACC)
4167# 1286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4169# 1286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4170#elif defined(MFC_OpenMP)
4171# 1286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4173# 1286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4176 &
l) + real((0.5_wp*dt*(vel_r(1)*(e_r + pres_r))*(1._wp/dx(
j + 1)) + 0.5_wp*dt*cfl*(e_r) &
4177 & *(1._wp/dx(
j + 1))), kind=stp)
4180# 1291 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4181#if defined(MFC_OpenACC)
4182# 1291 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4184# 1291 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4185#elif defined(MFC_OpenMP)
4186# 1291 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4188# 1291 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4191 &
l) - real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j)) &
4192 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j))), kind=stp)
4195# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4196#if defined(MFC_OpenACC)
4197# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4199# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4200#elif defined(MFC_OpenMP)
4201# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4203# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4206 &
l) - real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j)) &
4207 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j))), kind=stp)
4210# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4211#if defined(MFC_OpenACC)
4212# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4214# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4215#elif defined(MFC_OpenMP)
4216# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4218# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4221 &
l) - real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(3)*(1._wp/dx(
j)) &
4222 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dx(
j))), kind=stp)
4225# 1306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4226#if defined(MFC_OpenACC)
4227# 1306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4229# 1306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4230#elif defined(MFC_OpenMP)
4231# 1306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4233# 1306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4236 &
l) - real((0.5_wp*dt*(vel_r(1)*(e_r + pres_r))*(1._wp/dx(
j)) + 0.5_wp*dt*cfl*(e_r) &
4237 & *(1._wp/dx(
j))), kind=stp)
4242# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4243#if defined(MFC_OpenACC)
4244# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4246# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4247#elif defined(MFC_OpenMP)
4248# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4250# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4252# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4254# 1315 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4256 else if (idir == 2)
then
4258# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4260# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4262# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4263#if defined(MFC_OpenACC)
4264# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4266# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4267#elif defined(MFC_OpenMP)
4268# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4270# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4272# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4274# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4276# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4278# 1322 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4286# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4288# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4289#if defined(MFC_OpenACC)
4290# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4292# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4293#elif defined(MFC_OpenMP)
4294# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4296# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4301# 1339 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4302#if defined(MFC_OpenACC)
4303# 1339 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4305# 1339 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4306#elif defined(MFC_OpenMP)
4307# 1339 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4309# 1339 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4314# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4315#if defined(MFC_OpenACC)
4316# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4318# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4319#elif defined(MFC_OpenMP)
4320# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4322# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4324 do r = 1, num_fluids
4327 rho_sf_small(
i) = rho_l
4334 &
l)/rho_sf_small(-1))
4337 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
4338 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(1))/3._wp
4341 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
4342 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(1))/3._wp
4347# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4348#if defined(MFC_OpenACC)
4349# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4351# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4352#elif defined(MFC_OpenMP)
4353# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4355# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4360# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4361#if defined(MFC_OpenACC)
4362# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4364# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4365#elif defined(MFC_OpenMP)
4366# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4368# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4370 do r = 1, num_fluids
4373 rho_sf_small(
i) = rho_l
4380 &
l)/rho_sf_small(-1))
4383 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
4384 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(2))/3._wp
4387 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
4388 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(2))/3._wp
4394# 1392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4395#if defined(MFC_OpenACC)
4396# 1392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4398# 1392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4399#elif defined(MFC_OpenMP)
4400# 1392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4402# 1392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4404 do i = 1, num_fluids
4405 alpha_rho_l(
i) = 0._wp
4406 alpha_rho_r(
i) = 0._wp
4411# 1399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4412#if defined(MFC_OpenACC)
4413# 1399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4415# 1399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4416#elif defined(MFC_OpenMP)
4417# 1399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4419# 1399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4427# 1405 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4428#if defined(MFC_OpenACC)
4429# 1405 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4431# 1405 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4432#elif defined(MFC_OpenMP)
4433# 1405 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4435# 1405 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4439# 1407 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4440#if defined(MFC_OpenACC)
4441# 1407 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4443# 1407 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4444#elif defined(MFC_OpenMP)
4445# 1407 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4447# 1407 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4449 do i = 1, num_fluids
4453 if (num_fluids > 1)
then
4455# 1413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4456#if defined(MFC_OpenACC)
4457# 1413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4459# 1413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4460#elif defined(MFC_OpenMP)
4461# 1413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4463# 1413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4465 do i = 1, num_fluids - 1
4473# 1421 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4474#if defined(MFC_OpenACC)
4475# 1421 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4477# 1421 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4478#elif defined(MFC_OpenMP)
4479# 1421 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4481# 1421 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4489# 1427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4490#if defined(MFC_OpenACC)
4491# 1427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4493# 1427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4494#elif defined(MFC_OpenMP)
4495# 1427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4497# 1427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4501# 1429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4502#if defined(MFC_OpenACC)
4503# 1429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4505# 1429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4506#elif defined(MFC_OpenMP)
4507# 1429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4509# 1429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4511 do i = 1, num_fluids
4515 if (num_fluids > 1)
then
4517# 1435 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4518#if defined(MFC_OpenACC)
4519# 1435 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4521# 1435 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4522#elif defined(MFC_OpenMP)
4523# 1435 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4525# 1435 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4527 do i = 1, num_fluids - 1
4535# 1443 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4536#if defined(MFC_OpenACC)
4537# 1443 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4539# 1443 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4540#elif defined(MFC_OpenMP)
4541# 1443 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4543# 1443 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4550 if (num_fluids > 1)
then
4551 alpha_l(num_fluids) = 1._wp
4552 alpha_r(num_fluids) = 1._wp
4555# 1453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4556#if defined(MFC_OpenACC)
4557# 1453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4559# 1453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4560#elif defined(MFC_OpenMP)
4561# 1453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4563# 1453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4565 do i = 1, num_fluids - 1
4566 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
4567 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
4571 rho_l = 0._wp; rho_r = 0._wp
4572 gamma_l = 0._wp; gamma_r = 0._wp
4573 pi_inf_l = 0._wp; pi_inf_r = 0._wp
4576# 1464 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4577#if defined(MFC_OpenACC)
4578# 1464 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4580# 1464 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4581#elif defined(MFC_OpenMP)
4582# 1464 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4584# 1464 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4586 do i = 1, num_fluids
4587 rho_l = rho_l + alpha_rho_l(
i)
4588 gamma_l = gamma_l + alpha_l(
i)*gammas(
i)
4589 pi_inf_l = pi_inf_l + alpha_l(
i)*pi_infs(
i)
4591 rho_r = rho_r + alpha_rho_r(
i)
4592 gamma_r = gamma_r + alpha_r(
i)*gammas(
i)
4593 pi_inf_r = pi_inf_r + alpha_r(
i)*pi_infs(
i)
4597# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4598#if defined(MFC_OpenACC)
4599# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4601# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4602#elif defined(MFC_OpenMP)
4603# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4605# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4608 vel_l(
i) = vel_l(
i)/rho_l
4609 vel_r(
i) = vel_r(
i)/rho_r
4616# 1484 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4617#if defined(MFC_OpenACC)
4618# 1484 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4620# 1484 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4621#elif defined(MFC_OpenMP)
4622# 1484 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4624# 1484 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4626 do i = 1, num_fluids
4632# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4633#if defined(MFC_OpenACC)
4634# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4636# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4637#elif defined(MFC_OpenMP)
4638# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4640# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4643 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
4645# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4646#if defined(MFC_OpenACC)
4647# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4649# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4650#elif defined(MFC_OpenMP)
4651# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4653# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4656 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k + 1)), kind=stp)
4659# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4660#if defined(MFC_OpenACC)
4661# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4663# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4664#elif defined(MFC_OpenMP)
4665# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4667# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4670 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k)), kind=stp)
4672# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4673#if defined(MFC_OpenACC)
4674# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4676# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4677#elif defined(MFC_OpenMP)
4678# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4680# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4683 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k)), kind=stp)
4686# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4687#if defined(MFC_OpenACC)
4688# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4690# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4691#elif defined(MFC_OpenMP)
4692# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4694# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4697 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
4699# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4700#if defined(MFC_OpenACC)
4701# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4703# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4704#elif defined(MFC_OpenMP)
4705# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4707# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4710 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k + 1)), kind=stp)
4713# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4714#if defined(MFC_OpenACC)
4715# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4717# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4718#elif defined(MFC_OpenMP)
4719# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4721# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4724 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k)), kind=stp)
4726# 1514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4727#if defined(MFC_OpenACC)
4728# 1514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4730# 1514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4731#elif defined(MFC_OpenMP)
4732# 1514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4734# 1514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4737 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k)), kind=stp)
4740# 1518 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4741#if defined(MFC_OpenACC)
4742# 1518 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4744# 1518 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4745#elif defined(MFC_OpenMP)
4746# 1518 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4748# 1518 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4751 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
4753# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4754#if defined(MFC_OpenACC)
4755# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4757# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4758#elif defined(MFC_OpenMP)
4759# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4761# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4764 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
4767# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4768#if defined(MFC_OpenACC)
4769# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4771# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4772#elif defined(MFC_OpenMP)
4773# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4775# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4778 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k)), kind=stp)
4780# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4781#if defined(MFC_OpenACC)
4782# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4784# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4785#elif defined(MFC_OpenMP)
4786# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4788# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4791 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
4794# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4795#if defined(MFC_OpenACC)
4796# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4798# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4799#elif defined(MFC_OpenMP)
4800# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4802# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4805 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
4807# 1535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4808#if defined(MFC_OpenACC)
4809# 1535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4811# 1535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4812#elif defined(MFC_OpenMP)
4813# 1535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4815# 1535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4818 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
4821# 1539 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4822#if defined(MFC_OpenACC)
4823# 1539 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4825# 1539 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4826#elif defined(MFC_OpenMP)
4827# 1539 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4829# 1539 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4832 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k)), kind=stp)
4834# 1542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4835#if defined(MFC_OpenACC)
4836# 1542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4838# 1542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4839#elif defined(MFC_OpenMP)
4840# 1542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4842# 1542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4845 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
4848 e_l = 0._wp; e_r = 0._wp
4849 f_l = 0._wp; f_r = 0._wp
4852# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4853#if defined(MFC_OpenACC)
4854# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4856# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4857#elif defined(MFC_OpenMP)
4858# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4860# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4868# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4869#if defined(MFC_OpenACC)
4870# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4872# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4873#elif defined(MFC_OpenMP)
4874# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4876# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4883 call s_get_derived_states(e_l, gamma_l, pi_inf_l, rho_l, vel_l, e_r, gamma_r, pi_inf_r, rho_r, &
4884 & vel_r, pres_l, pres_r, cfl)
4887# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4888#if defined(MFC_OpenACC)
4889# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4891# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4892#elif defined(MFC_OpenMP)
4893# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4895# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4897 do i = 1, num_fluids
4899# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4900#if defined(MFC_OpenACC)
4901# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4903# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4904#elif defined(MFC_OpenMP)
4905# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4907# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4909 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
4910 &
l) + real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
4911 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
4914# 1572 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4915#if defined(MFC_OpenACC)
4916# 1572 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4918# 1572 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4919#elif defined(MFC_OpenMP)
4920# 1572 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4922# 1572 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4924 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
4925 &
l) - real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
4926 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k)), kind=stp)
4929 if (num_fluids > 1)
then
4931# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4932#if defined(MFC_OpenACC)
4933# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4935# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4936#elif defined(MFC_OpenMP)
4937# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4939# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4941 do i = 1, num_fluids - 1
4943# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4944#if defined(MFC_OpenACC)
4945# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4947# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4948#elif defined(MFC_OpenMP)
4949# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4951# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4954 &
l) + real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
4955 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
4958# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4959#if defined(MFC_OpenACC)
4960# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4962# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4963#elif defined(MFC_OpenMP)
4964# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4966# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4970 &
l)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
4973# 1591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4974#if defined(MFC_OpenACC)
4975# 1591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4977# 1591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4978#elif defined(MFC_OpenMP)
4979# 1591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4981# 1591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4984 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
4985 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k)), kind=stp)
4988# 1596 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4989#if defined(MFC_OpenACC)
4990# 1596 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4992# 1596 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4993#elif defined(MFC_OpenMP)
4994# 1596 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4996# 1596 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5000 &
l)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
5005# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5006#if defined(MFC_OpenACC)
5007# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5009# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5010#elif defined(MFC_OpenMP)
5011# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5013# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5016 &
l) + real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k + 1)) &
5017 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k + 1)), kind=stp)
5020# 1608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5021#if defined(MFC_OpenACC)
5022# 1608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5024# 1608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5025#elif defined(MFC_OpenMP)
5026# 1608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5028# 1608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5031 &
l) + real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dy(
k + 1)) &
5032 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dy(
k + 1)), kind=stp)
5035# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5036#if defined(MFC_OpenACC)
5037# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5039# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5040#elif defined(MFC_OpenMP)
5041# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5043# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5046 &
l) + real(0.5_wp*dt*(vel_l(2)*(e_l + pres_l + f_l))*(1._wp/dy(
k + 1)) &
5047 & - 0.5_wp*dt*cfl*(e_l)*(1._wp/dy(
k + 1)), kind=stp)
5050# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5051#if defined(MFC_OpenACC)
5052# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5054# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5055#elif defined(MFC_OpenMP)
5056# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5058# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5061 &
l) - real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k)) &
5062 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k)), kind=stp)
5065# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5066#if defined(MFC_OpenACC)
5067# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5069# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5070#elif defined(MFC_OpenMP)
5071# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5073# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5076 &
l) - real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dy(
k)) - 0.5_wp*dt*cfl*(rho_l*vel_l(1) &
5077 & )*(1._wp/dy(
k)), kind=stp)
5080# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5081#if defined(MFC_OpenACC)
5082# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5084# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5085#elif defined(MFC_OpenMP)
5086# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5088# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5091 &
l) - real(0.5_wp*dt*(vel_l(2)*(e_l + pres_l + f_l))*(1._wp/dy(
k)) - 0.5_wp*dt*cfl*(e_l) &
5092 & *(1._wp/dy(
k)), kind=stp)
5095# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5096#if defined(MFC_OpenACC)
5097# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5099# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5100#elif defined(MFC_OpenMP)
5101# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5103# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5105 do i = 1, num_fluids
5107# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5108#if defined(MFC_OpenACC)
5109# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5111# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5112#elif defined(MFC_OpenMP)
5113# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5115# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5117 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
5118 &
l) + real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
5119 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
5121# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5122#if defined(MFC_OpenACC)
5123# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5125# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5126#elif defined(MFC_OpenMP)
5127# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5129# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5131 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
5132 &
l) - real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
5133 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k)), kind=stp)
5136 if (num_fluids > 1)
then
5138# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5139#if defined(MFC_OpenACC)
5140# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5142# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5143#elif defined(MFC_OpenMP)
5144# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5146# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5148 do i = 1, num_fluids - 1
5150# 1648 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5151#if defined(MFC_OpenACC)
5152# 1648 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5154# 1648 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5155#elif defined(MFC_OpenMP)
5156# 1648 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5158# 1648 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5161 &
l) + real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
5162 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
5165# 1653 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5166#if defined(MFC_OpenACC)
5167# 1653 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5169# 1653 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5170#elif defined(MFC_OpenMP)
5171# 1653 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5173# 1653 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5177 &
l)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
5180# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5181#if defined(MFC_OpenACC)
5182# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5184# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5185#elif defined(MFC_OpenMP)
5186# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5188# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5191 &
l) - real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
5192 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k)), kind=stp)
5195# 1663 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5196#if defined(MFC_OpenACC)
5197# 1663 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5199# 1663 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5200#elif defined(MFC_OpenMP)
5201# 1663 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5203# 1663 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5207 &
l)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
5211# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5212#if defined(MFC_OpenACC)
5213# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5215# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5216#elif defined(MFC_OpenMP)
5217# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5219# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5222 &
l) + real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k + 1)) &
5223 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k + 1)), kind=stp)
5225# 1673 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5226#if defined(MFC_OpenACC)
5227# 1673 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5229# 1673 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5230#elif defined(MFC_OpenMP)
5231# 1673 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5233# 1673 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5236 &
l) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(1)*(1._wp/dy(
k + 1)) &
5237 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dy(
k + 1)), kind=stp)
5239# 1677 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5240#if defined(MFC_OpenACC)
5241# 1677 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5243# 1677 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5244#elif defined(MFC_OpenMP)
5245# 1677 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5247# 1677 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5250 &
l) + real(0.5_wp*dt*(vel_r(2)*(e_r + pres_r + f_r))*(1._wp/dy(
k + 1)) &
5251 & + 0.5_wp*dt*cfl*(e_r)*(1._wp/dy(
k + 1)), kind=stp)
5253# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5254#if defined(MFC_OpenACC)
5255# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5257# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5258#elif defined(MFC_OpenMP)
5259# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5261# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5264 &
l) - real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k)) &
5265 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k)), kind=stp)
5267# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5268#if defined(MFC_OpenACC)
5269# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5271# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5272#elif defined(MFC_OpenMP)
5273# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5275# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5278 &
l) - real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(1)*(1._wp/dy(
k)) + 0.5_wp*dt*cfl*(rho_r*vel_r(1) &
5279 & )*(1._wp/dy(
k)), kind=stp)
5281# 1689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5282#if defined(MFC_OpenACC)
5283# 1689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5285# 1689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5286#elif defined(MFC_OpenMP)
5287# 1689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5289# 1689 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5292 &
l) - real(0.5_wp*dt*(vel_r(2)*(e_r + pres_r + f_r))*(1._wp/dy(
k)) + 0.5_wp*dt*cfl*(e_r) &
5293 & *(1._wp/dy(
k)), kind=stp)
5298# 1696 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5299#if defined(MFC_OpenACC)
5300# 1696 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5302# 1696 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5303#elif defined(MFC_OpenMP)
5304# 1696 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5306# 1696 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5308# 1696 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5310# 1698 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5312# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5314# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5316# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5317#if defined(MFC_OpenACC)
5318# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5320# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5321#elif defined(MFC_OpenMP)
5322# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5324# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5326# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5328# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5330# 1700 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5332# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5340# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5342# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5343#if defined(MFC_OpenACC)
5344# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5346# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5347#elif defined(MFC_OpenMP)
5348# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5350# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5355# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5356#if defined(MFC_OpenACC)
5357# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5359# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5360#elif defined(MFC_OpenMP)
5361# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5363# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5368# 1723 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5369#if defined(MFC_OpenACC)
5370# 1723 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5372# 1723 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5373#elif defined(MFC_OpenMP)
5374# 1723 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5376# 1723 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5378 do r = 1, num_fluids
5381 rho_sf_small(
i) = rho_l
5388 &
l)/rho_sf_small(-1))
5391 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
5392 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(1))/3._wp
5395 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
5396 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(1))/3._wp
5401# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5402#if defined(MFC_OpenACC)
5403# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5405# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5406#elif defined(MFC_OpenMP)
5407# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5409# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5414# 1749 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5415#if defined(MFC_OpenACC)
5416# 1749 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5418# 1749 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5419#elif defined(MFC_OpenMP)
5420# 1749 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5422# 1749 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5424 do r = 1, num_fluids
5427 rho_sf_small(
i) = rho_l
5434 &
l)/rho_sf_small(-1))
5437 &
l)/rho_sf_small(-1))
5440 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
5441 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(3))
5442 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(2))/3._wp
5445 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
5446 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(3))
5447 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(2))/3._wp
5452# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5453#if defined(MFC_OpenACC)
5454# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5456# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5457#elif defined(MFC_OpenMP)
5458# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5460# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5465# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5466#if defined(MFC_OpenACC)
5467# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5469# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5470#elif defined(MFC_OpenMP)
5471# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5473# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5475 do r = 1, num_fluids
5478 rho_sf_small(
i) = rho_l
5483 &
l - 1)/rho_sf_small(-1))
5486 &
l - 1)/rho_sf_small(-1))
5488 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(2))
5489 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(3))/3._wp
5492 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(2))
5493 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(3))/3._wp
5499# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5500#if defined(MFC_OpenACC)
5501# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5503# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5504#elif defined(MFC_OpenMP)
5505# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5507# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5509 do i = 1, num_fluids
5510 alpha_rho_l(
i) = 0._wp
5511 alpha_rho_r(
i) = 0._wp
5516# 1811 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5517#if defined(MFC_OpenACC)
5518# 1811 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5520# 1811 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5521#elif defined(MFC_OpenMP)
5522# 1811 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5524# 1811 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5532# 1817 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5533#if defined(MFC_OpenACC)
5534# 1817 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5536# 1817 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5537#elif defined(MFC_OpenMP)
5538# 1817 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5540# 1817 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5544# 1819 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5545#if defined(MFC_OpenACC)
5546# 1819 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5548# 1819 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5549#elif defined(MFC_OpenMP)
5550# 1819 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5552# 1819 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5554 do i = 1, num_fluids
5558 if (num_fluids > 1)
then
5560# 1825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5561#if defined(MFC_OpenACC)
5562# 1825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5564# 1825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5565#elif defined(MFC_OpenMP)
5566# 1825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5568# 1825 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5570 do i = 1, num_fluids - 1
5578# 1833 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5579#if defined(MFC_OpenACC)
5580# 1833 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5582# 1833 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5583#elif defined(MFC_OpenMP)
5584# 1833 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5586# 1833 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5594# 1839 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5595#if defined(MFC_OpenACC)
5596# 1839 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5598# 1839 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5599#elif defined(MFC_OpenMP)
5600# 1839 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5602# 1839 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5606# 1841 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5607#if defined(MFC_OpenACC)
5608# 1841 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5610# 1841 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5611#elif defined(MFC_OpenMP)
5612# 1841 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5614# 1841 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5616 do i = 1, num_fluids
5620 if (num_fluids > 1)
then
5622# 1847 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5623#if defined(MFC_OpenACC)
5624# 1847 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5626# 1847 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5627#elif defined(MFC_OpenMP)
5628# 1847 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5630# 1847 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5632 do i = 1, num_fluids - 1
5640# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5641#if defined(MFC_OpenACC)
5642# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5644# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5645#elif defined(MFC_OpenMP)
5646# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5648# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5655 if (num_fluids > 1)
then
5656 alpha_l(num_fluids) = 1._wp
5657 alpha_r(num_fluids) = 1._wp
5660# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5661#if defined(MFC_OpenACC)
5662# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5664# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5665#elif defined(MFC_OpenMP)
5666# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5668# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5670 do i = 1, num_fluids - 1
5671 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
5672 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
5676 rho_l = 0._wp; rho_r = 0._wp
5677 gamma_l = 0._wp; gamma_r = 0._wp
5678 pi_inf_l = 0._wp; pi_inf_r = 0._wp
5681# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5682#if defined(MFC_OpenACC)
5683# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5685# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5686#elif defined(MFC_OpenMP)
5687# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5689# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5691 do i = 1, num_fluids
5692 rho_l = rho_l + alpha_rho_l(
i)
5693 gamma_l = gamma_l + alpha_l(
i)*gammas(
i)
5694 pi_inf_l = pi_inf_l + alpha_l(
i)*pi_infs(
i)
5696 rho_r = rho_r + alpha_rho_r(
i)
5697 gamma_r = gamma_r + alpha_r(
i)*gammas(
i)
5698 pi_inf_r = pi_inf_r + alpha_r(
i)*pi_infs(
i)
5702# 1887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5703#if defined(MFC_OpenACC)
5704# 1887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5706# 1887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5707#elif defined(MFC_OpenMP)
5708# 1887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5710# 1887 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5713 vel_l(
i) = vel_l(
i)/rho_l
5714 vel_r(
i) = vel_r(
i)/rho_r
5721# 1896 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5722#if defined(MFC_OpenACC)
5723# 1896 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5725# 1896 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5726#elif defined(MFC_OpenMP)
5727# 1896 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5729# 1896 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5731 do i = 1, num_fluids
5737# 1902 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5738#if defined(MFC_OpenACC)
5739# 1902 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5741# 1902 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5742#elif defined(MFC_OpenMP)
5743# 1902 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5745# 1902 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5748 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
5750# 1905 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5751#if defined(MFC_OpenACC)
5752# 1905 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5754# 1905 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5755#elif defined(MFC_OpenMP)
5756# 1905 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5758# 1905 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5761 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k + 1)), kind=stp)
5764# 1909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5765#if defined(MFC_OpenACC)
5766# 1909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5768# 1909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5769#elif defined(MFC_OpenMP)
5770# 1909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5772# 1909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5775 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k)), kind=stp)
5777# 1912 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5778#if defined(MFC_OpenACC)
5779# 1912 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5781# 1912 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5782#elif defined(MFC_OpenMP)
5783# 1912 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5785# 1912 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5788 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k)), kind=stp)
5791# 1916 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5792#if defined(MFC_OpenACC)
5793# 1916 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5795# 1916 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5796#elif defined(MFC_OpenMP)
5797# 1916 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5799# 1916 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5802 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
5804# 1919 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5805#if defined(MFC_OpenACC)
5806# 1919 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5808# 1919 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5809#elif defined(MFC_OpenMP)
5810# 1919 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5812# 1919 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5815 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k + 1)), kind=stp)
5818# 1923 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5819#if defined(MFC_OpenACC)
5820# 1923 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5822# 1923 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5823#elif defined(MFC_OpenMP)
5824# 1923 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5826# 1923 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5829 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k)), kind=stp)
5831# 1926 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5832#if defined(MFC_OpenACC)
5833# 1926 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5835# 1926 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5836#elif defined(MFC_OpenMP)
5837# 1926 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5839# 1926 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5842 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k)), kind=stp)
5845# 1930 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5846#if defined(MFC_OpenACC)
5847# 1930 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5849# 1930 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5850#elif defined(MFC_OpenMP)
5851# 1930 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5853# 1930 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5856 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dy(
k + 1)), kind=stp)
5858# 1933 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5859#if defined(MFC_OpenACC)
5860# 1933 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5862# 1933 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5863#elif defined(MFC_OpenMP)
5864# 1933 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5866# 1933 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5869 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dy(
k + 1)), kind=stp)
5872# 1937 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5873#if defined(MFC_OpenACC)
5874# 1937 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5876# 1937 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5877#elif defined(MFC_OpenMP)
5878# 1937 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5880# 1937 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5883 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dy(
k)), kind=stp)
5885# 1940 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5886#if defined(MFC_OpenACC)
5887# 1940 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5889# 1940 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5890#elif defined(MFC_OpenMP)
5891# 1940 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5893# 1940 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5896 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dy(
k)), kind=stp)
5899# 1944 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5900#if defined(MFC_OpenACC)
5901# 1944 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5903# 1944 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5904#elif defined(MFC_OpenMP)
5905# 1944 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5907# 1944 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5910 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dy(
k + 1)), kind=stp)
5912# 1947 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5913#if defined(MFC_OpenACC)
5914# 1947 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5916# 1947 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5917#elif defined(MFC_OpenMP)
5918# 1947 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5920# 1947 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5923 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dy(
k + 1)), kind=stp)
5926# 1951 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5927#if defined(MFC_OpenACC)
5928# 1951 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5930# 1951 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5931#elif defined(MFC_OpenMP)
5932# 1951 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5934# 1951 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5937 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dy(
k)), kind=stp)
5939# 1954 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5940#if defined(MFC_OpenACC)
5941# 1954 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5943# 1954 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5944#elif defined(MFC_OpenMP)
5945# 1954 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5947# 1954 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5950 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dy(
k)), kind=stp)
5953# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5954#if defined(MFC_OpenACC)
5955# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5957# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5958#elif defined(MFC_OpenMP)
5959# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5961# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5964 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
5966# 1961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5967#if defined(MFC_OpenACC)
5968# 1961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5970# 1961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5971#elif defined(MFC_OpenMP)
5972# 1961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5974# 1961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5977 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
5980# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5981#if defined(MFC_OpenACC)
5982# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5984# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5985#elif defined(MFC_OpenMP)
5986# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5988# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5991 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k)), kind=stp)
5993# 1968 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5994#if defined(MFC_OpenACC)
5995# 1968 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5997# 1968 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5998#elif defined(MFC_OpenMP)
5999# 1968 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6001# 1968 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6004 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
6007# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6008#if defined(MFC_OpenACC)
6009# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6011# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6012#elif defined(MFC_OpenMP)
6013# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6015# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6018 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
6020# 1975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6021#if defined(MFC_OpenACC)
6022# 1975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6024# 1975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6025#elif defined(MFC_OpenMP)
6026# 1975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6028# 1975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6031 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
6034# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6035#if defined(MFC_OpenACC)
6036# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6038# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6039#elif defined(MFC_OpenMP)
6040# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6042# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6045 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k)), kind=stp)
6047# 1982 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6048#if defined(MFC_OpenACC)
6049# 1982 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6051# 1982 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6052#elif defined(MFC_OpenMP)
6053# 1982 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6055# 1982 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6058 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
6061 e_l = 0._wp; e_r = 0._wp
6062 f_l = 0._wp; f_r = 0._wp
6065# 1990 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6066#if defined(MFC_OpenACC)
6067# 1990 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6069# 1990 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6070#elif defined(MFC_OpenMP)
6071# 1990 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6073# 1990 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6081# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6082#if defined(MFC_OpenACC)
6083# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6085# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6086#elif defined(MFC_OpenMP)
6087# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6089# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6096 call s_get_derived_states(e_l, gamma_l, pi_inf_l, rho_l, vel_l, e_r, gamma_r, pi_inf_r, rho_r, &
6097 & vel_r, pres_l, pres_r, cfl)
6100# 2005 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6101#if defined(MFC_OpenACC)
6102# 2005 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6104# 2005 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6105#elif defined(MFC_OpenMP)
6106# 2005 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6108# 2005 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6110 do i = 1, num_fluids
6112# 2007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6113#if defined(MFC_OpenACC)
6114# 2007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6116# 2007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6117#elif defined(MFC_OpenMP)
6118# 2007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6120# 2007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6122 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
6123 &
l) + real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
6124 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
6127# 2012 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6128#if defined(MFC_OpenACC)
6129# 2012 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6131# 2012 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6132#elif defined(MFC_OpenMP)
6133# 2012 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6135# 2012 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6137 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
6138 &
l) - real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
6139 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k)), kind=stp)
6142 if (num_fluids > 1)
then
6144# 2019 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6145#if defined(MFC_OpenACC)
6146# 2019 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6148# 2019 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6149#elif defined(MFC_OpenMP)
6150# 2019 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6152# 2019 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6154 do i = 1, num_fluids - 1
6156# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6157#if defined(MFC_OpenACC)
6158# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6160# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6161#elif defined(MFC_OpenMP)
6162# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6164# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6167 &
l) + real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
6168 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
6171# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6172#if defined(MFC_OpenACC)
6173# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6175# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6176#elif defined(MFC_OpenMP)
6177# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6179# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6183 &
l)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
6186# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6187#if defined(MFC_OpenACC)
6188# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6190# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6191#elif defined(MFC_OpenMP)
6192# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6194# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6197 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
6198 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k)), kind=stp)
6201# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6202#if defined(MFC_OpenACC)
6203# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6205# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6206#elif defined(MFC_OpenMP)
6207# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6209# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6213 &
l)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
6218# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6219#if defined(MFC_OpenACC)
6220# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6222# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6223#elif defined(MFC_OpenMP)
6224# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6226# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6229 &
l) + real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k + 1)) &
6230 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k + 1)), kind=stp)
6233# 2048 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6234#if defined(MFC_OpenACC)
6235# 2048 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6237# 2048 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6238#elif defined(MFC_OpenMP)
6239# 2048 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6241# 2048 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6244 &
l) + real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dy(
k + 1)) &
6245 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dy(
k + 1)), kind=stp)
6248# 2053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6249#if defined(MFC_OpenACC)
6250# 2053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6252# 2053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6253#elif defined(MFC_OpenMP)
6254# 2053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6256# 2053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6259 &
l) + real(0.5_wp*dt*rho_l*vel_l(3)*vel_l(2)*(1._wp/dy(
k + 1)) &
6260 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dy(
k + 1)), kind=stp)
6263# 2058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6264#if defined(MFC_OpenACC)
6265# 2058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6267# 2058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6268#elif defined(MFC_OpenMP)
6269# 2058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6271# 2058 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6274 &
l) + real(0.5_wp*dt*(vel_l(2)*(e_l + pres_l + f_l))*(1._wp/dy(
k + 1)) &
6275 & - 0.5_wp*dt*cfl*(e_l)*(1._wp/dy(
k + 1)), kind=stp)
6278# 2063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6279#if defined(MFC_OpenACC)
6280# 2063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6282# 2063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6283#elif defined(MFC_OpenMP)
6284# 2063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6286# 2063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6289 &
l) - real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k)) &
6290 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k)), kind=stp)
6293# 2068 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6294#if defined(MFC_OpenACC)
6295# 2068 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6297# 2068 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6298#elif defined(MFC_OpenMP)
6299# 2068 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6301# 2068 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6304 &
l) - real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dy(
k)) - 0.5_wp*dt*cfl*(rho_l*vel_l(1) &
6305 & )*(1._wp/dy(
k)), kind=stp)
6308# 2073 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6309#if defined(MFC_OpenACC)
6310# 2073 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6312# 2073 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6313#elif defined(MFC_OpenMP)
6314# 2073 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6316# 2073 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6319 &
l) - real(0.5_wp*dt*rho_l*vel_l(3)*vel_l(2)*(1._wp/dy(
k)) - 0.5_wp*dt*cfl*(rho_l*vel_l(3) &
6320 & )*(1._wp/dy(
k)), kind=stp)
6323# 2078 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6324#if defined(MFC_OpenACC)
6325# 2078 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6327# 2078 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6328#elif defined(MFC_OpenMP)
6329# 2078 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6331# 2078 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6334 &
l) - real(0.5_wp*dt*(vel_l(2)*(e_l + pres_l + f_l))*(1._wp/dy(
k)) - 0.5_wp*dt*cfl*(e_l) &
6335 & *(1._wp/dy(
k)), kind=stp)
6338# 2083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6339#if defined(MFC_OpenACC)
6340# 2083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6342# 2083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6343#elif defined(MFC_OpenMP)
6344# 2083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6346# 2083 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6348 do i = 1, num_fluids
6350# 2085 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6351#if defined(MFC_OpenACC)
6352# 2085 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6354# 2085 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6355#elif defined(MFC_OpenMP)
6356# 2085 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6358# 2085 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6360 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
6361 &
l) + real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
6362 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
6365# 2090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6366#if defined(MFC_OpenACC)
6367# 2090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6369# 2090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6370#elif defined(MFC_OpenMP)
6371# 2090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6373# 2090 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6375 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
6376 &
l) - real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
6377 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k)), kind=stp)
6380 if (num_fluids > 1)
then
6382# 2097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6383#if defined(MFC_OpenACC)
6384# 2097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6386# 2097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6387#elif defined(MFC_OpenMP)
6388# 2097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6390# 2097 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6392 do i = 1, num_fluids - 1
6394# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6395#if defined(MFC_OpenACC)
6396# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6398# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6399#elif defined(MFC_OpenMP)
6400# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6402# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6405 &
l) + real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
6406 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
6409# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6410#if defined(MFC_OpenACC)
6411# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6413# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6414#elif defined(MFC_OpenMP)
6415# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6417# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6421 &
l)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
6424# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6425#if defined(MFC_OpenACC)
6426# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6428# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6429#elif defined(MFC_OpenMP)
6430# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6432# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6435 &
l) - real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
6436 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k)), kind=stp)
6439# 2114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6440#if defined(MFC_OpenACC)
6441# 2114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6443# 2114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6444#elif defined(MFC_OpenMP)
6445# 2114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6447# 2114 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6451 &
l)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
6456# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6457#if defined(MFC_OpenACC)
6458# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6460# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6461#elif defined(MFC_OpenMP)
6462# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6464# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6467 &
l) + real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k + 1)) &
6468 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k + 1)), kind=stp)
6471# 2126 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6472#if defined(MFC_OpenACC)
6473# 2126 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6475# 2126 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6476#elif defined(MFC_OpenMP)
6477# 2126 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6479# 2126 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6482 &
l) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(1)*(1._wp/dy(
k + 1)) &
6483 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dy(
k + 1)), kind=stp)
6486# 2131 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6487#if defined(MFC_OpenACC)
6488# 2131 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6490# 2131 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6491#elif defined(MFC_OpenMP)
6492# 2131 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6494# 2131 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6497 &
l) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(3)*(1._wp/dy(
k + 1)) &
6498 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dy(
k + 1)), kind=stp)
6501# 2136 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6502#if defined(MFC_OpenACC)
6503# 2136 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6505# 2136 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6506#elif defined(MFC_OpenMP)
6507# 2136 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6509# 2136 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6512 &
l) + real(0.5_wp*dt*(vel_r(2)*(e_r + pres_r + f_r))*(1._wp/dy(
k + 1)) &
6513 & + 0.5_wp*dt*cfl*(e_r)*(1._wp/dy(
k + 1)), kind=stp)
6516# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6517#if defined(MFC_OpenACC)
6518# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6520# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6521#elif defined(MFC_OpenMP)
6522# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6524# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6527 &
l) - real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k)) &
6528 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k)), kind=stp)
6531# 2146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6532#if defined(MFC_OpenACC)
6533# 2146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6535# 2146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6536#elif defined(MFC_OpenMP)
6537# 2146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6539# 2146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6542 &
l) - real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(1)*(1._wp/dy(
k)) + 0.5_wp*dt*cfl*(rho_r*vel_r(1) &
6543 & )*(1._wp/dy(
k)), kind=stp)
6546# 2151 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6547#if defined(MFC_OpenACC)
6548# 2151 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6550# 2151 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6551#elif defined(MFC_OpenMP)
6552# 2151 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6554# 2151 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6557 &
l) - real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(3)*(1._wp/dy(
k)) + 0.5_wp*dt*cfl*(rho_r*vel_r(3) &
6558 & )*(1._wp/dy(
k)), kind=stp)
6561# 2156 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6562#if defined(MFC_OpenACC)
6563# 2156 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6565# 2156 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6566#elif defined(MFC_OpenMP)
6567# 2156 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6569# 2156 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6572 &
l) - real(0.5_wp*dt*(vel_r(2)*(e_r + pres_r + f_r))*(1._wp/dy(
k)) + 0.5_wp*dt*cfl*(e_r) &
6573 & *(1._wp/dy(
k)), kind=stp)
6578# 2163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6579#if defined(MFC_OpenACC)
6580# 2163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6582# 2163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6583#elif defined(MFC_OpenMP)
6584# 2163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6586# 2163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6588# 2163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6590# 2165 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6592 else if (idir == 3)
then
6593# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6595# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6597# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6598#if defined(MFC_OpenACC)
6599# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6601# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6602#elif defined(MFC_OpenMP)
6603# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6605# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6607# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6609# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6611# 2168 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6613# 2171 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6621# 2185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6623# 2185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6624#if defined(MFC_OpenACC)
6625# 2185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6627# 2185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6628#elif defined(MFC_OpenMP)
6629# 2185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6631# 2185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6636# 2188 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6637#if defined(MFC_OpenACC)
6638# 2188 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6640# 2188 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6641#elif defined(MFC_OpenMP)
6642# 2188 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6644# 2188 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6649# 2191 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6650#if defined(MFC_OpenACC)
6651# 2191 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6653# 2191 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6654#elif defined(MFC_OpenMP)
6655# 2191 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6657# 2191 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6659 do r = 1, num_fluids
6662 rho_sf_small(
i) = rho_l
6669 &
l +
q)/rho_sf_small(-1))
6672 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(3))
6673 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(1))/3._wp
6676 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(3))
6677 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(1))/3._wp
6682# 2214 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6683#if defined(MFC_OpenACC)
6684# 2214 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6686# 2214 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6687#elif defined(MFC_OpenMP)
6688# 2214 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6690# 2214 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6695# 2217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6696#if defined(MFC_OpenACC)
6697# 2217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6699# 2217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6700#elif defined(MFC_OpenMP)
6701# 2217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6703# 2217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6705 do r = 1, num_fluids
6708 rho_sf_small(
i) = rho_l
6713 &
l +
q)/rho_sf_small(-1))
6716 &
l +
q)/rho_sf_small(-1))
6719 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(3))
6720 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(2))/3._wp
6723 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(3))
6724 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(2))/3._wp
6729# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6730#if defined(MFC_OpenACC)
6731# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6733# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6734#elif defined(MFC_OpenMP)
6735# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6737# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6742# 2244 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6743#if defined(MFC_OpenACC)
6744# 2244 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6746# 2244 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6747#elif defined(MFC_OpenMP)
6748# 2244 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6750# 2244 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6752 do r = 1, num_fluids
6755 rho_sf_small(
i) = rho_l
6759 &
l - 1 +
q)/rho_sf_small(-1))
6762 &
l - 1 +
q)/rho_sf_small(-1))
6765 &
l - 1 +
q)/rho_sf_small(-1))
6767 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
6768 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(2))
6769 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(3))/3._wp
6772 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
6773 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(2))
6774 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(3))/3._wp
6780# 2272 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6781#if defined(MFC_OpenACC)
6782# 2272 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6784# 2272 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6785#elif defined(MFC_OpenMP)
6786# 2272 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6788# 2272 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6790 do i = 1, num_fluids
6791 alpha_rho_l(
i) = 0._wp
6792 alpha_rho_r(
i) = 0._wp
6798# 2280 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6799#if defined(MFC_OpenACC)
6800# 2280 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6802# 2280 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6803#elif defined(MFC_OpenMP)
6804# 2280 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6806# 2280 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6814# 2286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6815#if defined(MFC_OpenACC)
6816# 2286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6818# 2286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6819#elif defined(MFC_OpenMP)
6820# 2286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6822# 2286 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6826# 2288 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6827#if defined(MFC_OpenACC)
6828# 2288 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6830# 2288 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6831#elif defined(MFC_OpenMP)
6832# 2288 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6834# 2288 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6836 do i = 1, num_fluids
6840 if (num_fluids > 1)
then
6842# 2294 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6843#if defined(MFC_OpenACC)
6844# 2294 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6846# 2294 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6847#elif defined(MFC_OpenMP)
6848# 2294 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6850# 2294 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6852 do i = 1, num_fluids - 1
6860# 2302 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6861#if defined(MFC_OpenACC)
6862# 2302 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6864# 2302 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6865#elif defined(MFC_OpenMP)
6866# 2302 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6868# 2302 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6876# 2308 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6877#if defined(MFC_OpenACC)
6878# 2308 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6880# 2308 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6881#elif defined(MFC_OpenMP)
6882# 2308 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6884# 2308 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6888# 2310 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6889#if defined(MFC_OpenACC)
6890# 2310 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6892# 2310 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6893#elif defined(MFC_OpenMP)
6894# 2310 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6896# 2310 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6898 do i = 1, num_fluids
6902 if (num_fluids > 1)
then
6904# 2316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6905#if defined(MFC_OpenACC)
6906# 2316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6908# 2316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6909#elif defined(MFC_OpenMP)
6910# 2316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6912# 2316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6914 do i = 1, num_fluids - 1
6922# 2324 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6923#if defined(MFC_OpenACC)
6924# 2324 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6926# 2324 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6927#elif defined(MFC_OpenMP)
6928# 2324 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6930# 2324 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6937 if (num_fluids > 1)
then
6938 alpha_l(num_fluids) = 1._wp
6939 alpha_r(num_fluids) = 1._wp
6942# 2334 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6943#if defined(MFC_OpenACC)
6944# 2334 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6946# 2334 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6947#elif defined(MFC_OpenMP)
6948# 2334 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6950# 2334 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6952 do i = 1, num_fluids - 1
6953 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
6954 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
6958 rho_l = 0._wp; rho_r = 0._wp
6959 gamma_l = 0._wp; gamma_r = 0._wp
6960 pi_inf_l = 0._wp; pi_inf_r = 0._wp
6963# 2345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6964#if defined(MFC_OpenACC)
6965# 2345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6967# 2345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6968#elif defined(MFC_OpenMP)
6969# 2345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6971# 2345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6973 do i = 1, num_fluids
6974 rho_l = rho_l + alpha_rho_l(
i)
6975 gamma_l = gamma_l + alpha_l(
i)*gammas(
i)
6976 pi_inf_l = pi_inf_l + alpha_l(
i)*pi_infs(
i)
6978 rho_r = rho_r + alpha_rho_r(
i)
6979 gamma_r = gamma_r + alpha_r(
i)*gammas(
i)
6980 pi_inf_r = pi_inf_r + alpha_r(
i)*pi_infs(
i)
6984# 2356 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6985#if defined(MFC_OpenACC)
6986# 2356 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6988# 2356 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6989#elif defined(MFC_OpenMP)
6990# 2356 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6992# 2356 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6995 vel_l(
i) = vel_l(
i)/rho_l
6996 vel_r(
i) = vel_r(
i)/rho_r
7003# 2365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7004#if defined(MFC_OpenACC)
7005# 2365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7007# 2365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7008#elif defined(MFC_OpenMP)
7009# 2365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7011# 2365 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7013 do i = 1, num_fluids
7019# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7020#if defined(MFC_OpenACC)
7021# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7023# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7024#elif defined(MFC_OpenMP)
7025# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7027# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7030 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dz(
l + 1)), kind=stp)
7032# 2374 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7033#if defined(MFC_OpenACC)
7034# 2374 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7036# 2374 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7037#elif defined(MFC_OpenMP)
7038# 2374 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7040# 2374 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7043 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dz(
l + 1)), kind=stp)
7046# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7047#if defined(MFC_OpenACC)
7048# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7050# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7051#elif defined(MFC_OpenMP)
7052# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7054# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7057 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dz(
l)), kind=stp)
7059# 2381 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7060#if defined(MFC_OpenACC)
7061# 2381 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7063# 2381 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7064#elif defined(MFC_OpenMP)
7065# 2381 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7067# 2381 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7070 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dz(
l)), kind=stp)
7073# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7074#if defined(MFC_OpenACC)
7075# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7077# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7078#elif defined(MFC_OpenMP)
7079# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7081# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7084 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dz(
l + 1)), kind=stp)
7086# 2388 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7087#if defined(MFC_OpenACC)
7088# 2388 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7090# 2388 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7091#elif defined(MFC_OpenMP)
7092# 2388 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7094# 2388 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7097 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dz(
l + 1)), kind=stp)
7100# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7101#if defined(MFC_OpenACC)
7102# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7104# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7105#elif defined(MFC_OpenMP)
7106# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7108# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7111 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dz(
l)), kind=stp)
7113# 2395 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7114#if defined(MFC_OpenACC)
7115# 2395 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7117# 2395 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7118#elif defined(MFC_OpenMP)
7119# 2395 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7121# 2395 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7124 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dz(
l)), kind=stp)
7127# 2399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7128#if defined(MFC_OpenACC)
7129# 2399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7131# 2399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7132#elif defined(MFC_OpenMP)
7133# 2399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7135# 2399 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7138 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dz(
l + 1)), kind=stp)
7140# 2402 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7141#if defined(MFC_OpenACC)
7142# 2402 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7144# 2402 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7145#elif defined(MFC_OpenMP)
7146# 2402 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7148# 2402 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7151 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(2)*(1._wp/dz(
l + 1)), kind=stp)
7154# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7155#if defined(MFC_OpenACC)
7156# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7158# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7159#elif defined(MFC_OpenMP)
7160# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7162# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7165 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dz(
l)), kind=stp)
7167# 2409 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7168#if defined(MFC_OpenACC)
7169# 2409 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7171# 2409 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7172#elif defined(MFC_OpenMP)
7173# 2409 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7175# 2409 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7178 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(2)*(1._wp/dz(
l)), kind=stp)
7181# 2413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7182#if defined(MFC_OpenACC)
7183# 2413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7185# 2413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7186#elif defined(MFC_OpenMP)
7187# 2413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7189# 2413 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7192 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dz(
l + 1)), kind=stp)
7194# 2416 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7195#if defined(MFC_OpenACC)
7196# 2416 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7198# 2416 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7199#elif defined(MFC_OpenMP)
7200# 2416 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7202# 2416 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7205 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(2)*(1._wp/dz(
l + 1)), kind=stp)
7208# 2420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7209#if defined(MFC_OpenACC)
7210# 2420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7212# 2420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7213#elif defined(MFC_OpenMP)
7214# 2420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7216# 2420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7219 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dz(
l)), kind=stp)
7221# 2423 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7222#if defined(MFC_OpenACC)
7223# 2423 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7225# 2423 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7226#elif defined(MFC_OpenMP)
7227# 2423 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7229# 2423 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7232 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(2)*(1._wp/dz(
l)), kind=stp)
7235# 2427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7236#if defined(MFC_OpenACC)
7237# 2427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7239# 2427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7240#elif defined(MFC_OpenMP)
7241# 2427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7243# 2427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7246 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dz(
l + 1)), kind=stp)
7248# 2430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7249#if defined(MFC_OpenACC)
7250# 2430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7252# 2430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7253#elif defined(MFC_OpenMP)
7254# 2430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7256# 2430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7259 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(3)*(1._wp/dz(
l + 1)), kind=stp)
7262# 2434 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7263#if defined(MFC_OpenACC)
7264# 2434 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7266# 2434 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7267#elif defined(MFC_OpenMP)
7268# 2434 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7270# 2434 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7273 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dz(
l)), kind=stp)
7275# 2437 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7276#if defined(MFC_OpenACC)
7277# 2437 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7279# 2437 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7280#elif defined(MFC_OpenMP)
7281# 2437 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7283# 2437 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7286 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(3)*(1._wp/dz(
l)), kind=stp)
7289# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7290#if defined(MFC_OpenACC)
7291# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7293# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7294#elif defined(MFC_OpenMP)
7295# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7297# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7300 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dz(
l + 1)), kind=stp)
7302# 2444 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7303#if defined(MFC_OpenACC)
7304# 2444 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7306# 2444 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7307#elif defined(MFC_OpenMP)
7308# 2444 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7310# 2444 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7313 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(3)*(1._wp/dz(
l + 1)), kind=stp)
7316# 2448 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7317#if defined(MFC_OpenACC)
7318# 2448 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7320# 2448 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7321#elif defined(MFC_OpenMP)
7322# 2448 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7324# 2448 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7327 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dz(
l)), kind=stp)
7329# 2451 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7330#if defined(MFC_OpenACC)
7331# 2451 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7333# 2451 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7334#elif defined(MFC_OpenMP)
7335# 2451 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7337# 2451 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7340 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(3)*(1._wp/dz(
l)), kind=stp)
7343 e_l = 0._wp; e_r = 0._wp
7344 f_l = 0._wp; f_r = 0._wp
7347# 2459 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7348#if defined(MFC_OpenACC)
7349# 2459 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7351# 2459 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7352#elif defined(MFC_OpenMP)
7353# 2459 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7355# 2459 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7363# 2465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7364#if defined(MFC_OpenACC)
7365# 2465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7367# 2465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7368#elif defined(MFC_OpenMP)
7369# 2465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7371# 2465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7378 call s_get_derived_states(e_l, gamma_l, pi_inf_l, rho_l, vel_l, e_r, gamma_r, pi_inf_r, rho_r, vel_r, &
7379 & pres_l, pres_r, cfl)
7382# 2474 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7383#if defined(MFC_OpenACC)
7384# 2474 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7386# 2474 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7387#elif defined(MFC_OpenMP)
7388# 2474 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7390# 2474 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7392 do i = 1, num_fluids
7394# 2476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7395#if defined(MFC_OpenACC)
7396# 2476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7398# 2476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7399#elif defined(MFC_OpenMP)
7400# 2476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7402# 2476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7404 rhs_vf(
i)%sf(
j,
k,
l + 1) = rhs_vf(
i)%sf(
j,
k, &
7405 &
l + 1) + real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(3))*(1._wp/dz(
l + 1)) &
7406 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dz(
l + 1)), kind=stp)
7409# 2481 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7410#if defined(MFC_OpenACC)
7411# 2481 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7413# 2481 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7414#elif defined(MFC_OpenMP)
7415# 2481 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7417# 2481 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7419 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
7420 &
l) - real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(3))*(1._wp/dz(
l)) &
7421 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dz(
l)), kind=stp)
7424 if (num_fluids > 1)
then
7426# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7427#if defined(MFC_OpenACC)
7428# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7430# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7431#elif defined(MFC_OpenMP)
7432# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7434# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7436 do i = 1, num_fluids - 1
7438# 2490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7439#if defined(MFC_OpenACC)
7440# 2490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7442# 2490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7443#elif defined(MFC_OpenMP)
7444# 2490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7446# 2490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7449 &
l + 1) + real(0.5_wp*dt*(alpha_l(
i)*vel_l(3))*(1._wp/dz(
l + 1)) &
7450 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dz(
l + 1)), kind=stp)
7453# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7454#if defined(MFC_OpenACC)
7455# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7457# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7458#elif defined(MFC_OpenMP)
7459# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7461# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7465 &
l + 1)*vel_l(3)*(1._wp/dz(
l + 1)), kind=stp)
7468# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7469#if defined(MFC_OpenACC)
7470# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7472# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7473#elif defined(MFC_OpenMP)
7474# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7476# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7479 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(3))*(1._wp/dz(
l)) - 0.5_wp*dt*cfl*(alpha_l(
i)) &
7480 & *(1._wp/dz(
l)), kind=stp)
7483# 2505 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7484#if defined(MFC_OpenACC)
7485# 2505 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7487# 2505 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7488#elif defined(MFC_OpenMP)
7489# 2505 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7491# 2505 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7500# 2512 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7501#if defined(MFC_OpenACC)
7502# 2512 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7504# 2512 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7505#elif defined(MFC_OpenMP)
7506# 2512 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7508# 2512 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7511 &
l + 1) + real(0.5_wp*dt*(rho_l*(vel_l(3))**2.0 + pres_l + f_l)*(1._wp/dz(
l + 1)) &
7512 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dz(
l + 1)), kind=stp)
7515# 2517 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7516#if defined(MFC_OpenACC)
7517# 2517 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7519# 2517 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7520#elif defined(MFC_OpenMP)
7521# 2517 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7523# 2517 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7526 &
l + 1) + real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(3)*(1._wp/dz(
l + 1)) &
7527 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dz(
l + 1)), kind=stp)
7530# 2522 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7531#if defined(MFC_OpenACC)
7532# 2522 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7534# 2522 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7535#elif defined(MFC_OpenMP)
7536# 2522 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7538# 2522 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7541 &
l + 1) + real(0.5_wp*dt*rho_l*vel_l(2)*vel_l(3)*(1._wp/dz(
l + 1)) &
7542 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dz(
l + 1)), kind=stp)
7545# 2527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7546#if defined(MFC_OpenACC)
7547# 2527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7549# 2527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7550#elif defined(MFC_OpenMP)
7551# 2527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7553# 2527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7556 &
l + 1) + real(0.5_wp*dt*(vel_l(3)*(e_l + pres_l + f_l))*(1._wp/dz(
l + 1)) &
7557 & - 0.5_wp*dt*cfl*(e_l)*(1._wp/dz(
l + 1)), kind=stp)
7560# 2532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7561#if defined(MFC_OpenACC)
7562# 2532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7564# 2532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7565#elif defined(MFC_OpenMP)
7566# 2532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7568# 2532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7571 &
l) - real(0.5_wp*dt*(rho_l*(vel_l(3))**2.0 + pres_l + f_l)*(1._wp/dz(
l)) &
7572 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dz(
l)), kind=stp)
7575# 2537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7576#if defined(MFC_OpenACC)
7577# 2537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7579# 2537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7580#elif defined(MFC_OpenMP)
7581# 2537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7583# 2537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7586 &
l) - real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(3)*(1._wp/dz(
l)) - 0.5_wp*dt*cfl*(rho_l*vel_l(1)) &
7587 & *(1._wp/dz(
l)), kind=stp)
7590# 2542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7591#if defined(MFC_OpenACC)
7592# 2542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7594# 2542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7595#elif defined(MFC_OpenMP)
7596# 2542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7598# 2542 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7601 &
l) - real(0.5_wp*dt*rho_l*vel_l(2)*vel_l(3)*(1._wp/dz(
l)) - 0.5_wp*dt*cfl*(rho_l*vel_l(2)) &
7602 & *(1._wp/dz(
l)), kind=stp)
7605# 2547 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7606#if defined(MFC_OpenACC)
7607# 2547 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7609# 2547 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7610#elif defined(MFC_OpenMP)
7611# 2547 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7613# 2547 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7616 &
l) - real(0.5_wp*dt*(vel_l(3)*(e_l + pres_l + f_l))*(1._wp/dz(
l)) - 0.5_wp*dt*cfl*(e_l) &
7617 & *(1._wp/dz(
l)), kind=stp)
7620# 2552 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7621#if defined(MFC_OpenACC)
7622# 2552 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7624# 2552 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7625#elif defined(MFC_OpenMP)
7626# 2552 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7628# 2552 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7630 do i = 1, num_fluids
7632# 2554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7633#if defined(MFC_OpenACC)
7634# 2554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7636# 2554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7637#elif defined(MFC_OpenMP)
7638# 2554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7640# 2554 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7642 rhs_vf(
i)%sf(
j,
k,
l + 1) = rhs_vf(
i)%sf(
j,
k, &
7643 &
l + 1) + real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(3))*(1._wp/dz(
l + 1)) &
7644 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dz(
l + 1)), kind=stp)
7647# 2559 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7648#if defined(MFC_OpenACC)
7649# 2559 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7651# 2559 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7652#elif defined(MFC_OpenMP)
7653# 2559 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7655# 2559 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7657 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
7658 &
l) - real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(3))*(1._wp/dz(
l)) &
7659 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dz(
l)), kind=stp)
7662 if (num_fluids > 1)
then
7664# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7665#if defined(MFC_OpenACC)
7666# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7668# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7669#elif defined(MFC_OpenMP)
7670# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7672# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7674 do i = 1, num_fluids - 1
7676# 2568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7677#if defined(MFC_OpenACC)
7678# 2568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7680# 2568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7681#elif defined(MFC_OpenMP)
7682# 2568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7684# 2568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7687 &
l + 1) + real(0.5_wp*dt*(alpha_r(
i)*vel_r(3))*(1._wp/dz(
l + 1)) &
7688 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dz(
l + 1)), kind=stp)
7691# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7692#if defined(MFC_OpenACC)
7693# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7695# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7696#elif defined(MFC_OpenMP)
7697# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7699# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7703 &
l + 1)*vel_r(3)*(1._wp/dz(
l + 1)), kind=stp)
7706# 2578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7707#if defined(MFC_OpenACC)
7708# 2578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7710# 2578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7711#elif defined(MFC_OpenMP)
7712# 2578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7714# 2578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7717 &
l) - real(0.5_wp*dt*(alpha_r(
i)*vel_r(3))*(1._wp/dz(
l)) + 0.5_wp*dt*cfl*(alpha_r(
i)) &
7718 & *(1._wp/dz(
l)), kind=stp)
7721# 2583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7722#if defined(MFC_OpenACC)
7723# 2583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7725# 2583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7726#elif defined(MFC_OpenMP)
7727# 2583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7729# 2583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7738# 2590 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7739#if defined(MFC_OpenACC)
7740# 2590 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7742# 2590 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7743#elif defined(MFC_OpenMP)
7744# 2590 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7746# 2590 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7749 &
l + 1) + real(0.5_wp*dt*(rho_r*(vel_r(3))**2.0 + pres_r + f_r)*(1._wp/dz(
l + 1)) &
7750 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dz(
l + 1)), kind=stp)
7753# 2595 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7754#if defined(MFC_OpenACC)
7755# 2595 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7757# 2595 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7758#elif defined(MFC_OpenMP)
7759# 2595 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7761# 2595 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7764 &
l + 1) + real(0.5_wp*dt*rho_r*vel_r(1)*vel_r(3)*(1._wp/dz(
l + 1)) &
7765 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dz(
l + 1)), kind=stp)
7768# 2600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7769#if defined(MFC_OpenACC)
7770# 2600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7772# 2600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7773#elif defined(MFC_OpenMP)
7774# 2600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7776# 2600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7779 &
l + 1) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(3)*(1._wp/dz(
l + 1)) &
7780 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dz(
l + 1)), kind=stp)
7783# 2605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7784#if defined(MFC_OpenACC)
7785# 2605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7787# 2605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7788#elif defined(MFC_OpenMP)
7789# 2605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7791# 2605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7794 &
l + 1) + real(0.5_wp*dt*(vel_r(3)*(e_r + pres_r + f_r))*(1._wp/dz(
l + 1)) &
7795 & + 0.5_wp*dt*cfl*(e_r)*(1._wp/dz(
l + 1)), kind=stp)
7798# 2610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7799#if defined(MFC_OpenACC)
7800# 2610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7802# 2610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7803#elif defined(MFC_OpenMP)
7804# 2610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7806# 2610 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7809 &
l) - real(0.5_wp*dt*(rho_r*(vel_r(3))**2.0 + pres_r + f_r)*(1._wp/dz(
l)) &
7810 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dz(
l)), kind=stp)
7813# 2615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7814#if defined(MFC_OpenACC)
7815# 2615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7817# 2615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7818#elif defined(MFC_OpenMP)
7819# 2615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7821# 2615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7824 &
l) - real(0.5_wp*dt*rho_r*vel_r(1)*vel_r(3)*(1._wp/dz(
l)) + 0.5_wp*dt*cfl*(rho_r*vel_r(1)) &
7825 & *(1._wp/dz(
l)), kind=stp)
7828# 2620 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7829#if defined(MFC_OpenACC)
7830# 2620 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7832# 2620 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7833#elif defined(MFC_OpenMP)
7834# 2620 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7836# 2620 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7839 &
l) - real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(3)*(1._wp/dz(
l)) + 0.5_wp*dt*cfl*(rho_r*vel_r(2)) &
7840 & *(1._wp/dz(
l)), kind=stp)
7843# 2625 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7844#if defined(MFC_OpenACC)
7845# 2625 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7847# 2625 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7848#elif defined(MFC_OpenMP)
7849# 2625 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7851# 2625 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7854 &
l) - real(0.5_wp*dt*(vel_r(3)*(e_r + pres_r + f_r))*(1._wp/dz(
l)) + 0.5_wp*dt*cfl*(e_r) &
7855 & *(1._wp/dz(
l)), kind=stp)
7860# 2632 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7861#if defined(MFC_OpenACC)
7862# 2632 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7864# 2632 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7865#elif defined(MFC_OpenMP)
7866# 2632 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7868# 2632 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7870# 2632 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7872# 2634 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"