1798 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf
1799 type(scalar_field),
dimension(sys_size),
intent(inout) ::
q_cons_vf
1800 integer,
intent(in) :: idir
1802 real(wp) :: rho_l, gamma_l, pi_inf_l, e_l, mu_l, f_l, pres_l
1804 real(wp) :: rho_r, gamma_r, pi_inf_r, e_r, mu_r, f_r, pres_r
1805 real(wp),
dimension(3) :: vflux_l_arr, vflux_r_arr
1806 real(wp),
dimension(-1:1) :: rho_sf_small
1807# 417 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1808 real(wp),
dimension(num_fluids_max) :: alpha_rho_l, alpha_l, alpha_r, alpha_rho_r
1809 real(wp),
dimension(3) :: vel_l, vel_r
1810 real(wp),
dimension(3, 3) :: dvel
1811 real(wp),
dimension(3) :: dvel_small
1812# 427 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1816# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1818# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1820# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1821#if defined(MFC_OpenACC)
1822# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1824# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1826# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1827#elif defined(MFC_OpenMP)
1828# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1830# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1832# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1834# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1836# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1838# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1840# 431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1842# 434 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1849# 447 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1851# 447 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1852#if defined(MFC_OpenACC)
1853# 447 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1855# 447 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1856#elif defined(MFC_OpenMP)
1857# 447 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1859# 447 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1864# 450 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1865#if defined(MFC_OpenACC)
1866# 450 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1868# 450 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1869#elif defined(MFC_OpenMP)
1870# 450 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1872# 450 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1877# 453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1878#if defined(MFC_OpenACC)
1879# 453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1881# 453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1882#elif defined(MFC_OpenMP)
1883# 453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1885# 453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1887 do r = 1, num_fluids
1890 rho_sf_small(
i) = rho_l
1895 &
l)/rho_sf_small(-1))
1901# 467 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1902#if defined(MFC_OpenACC)
1903# 467 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1905# 467 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1906#elif defined(MFC_OpenMP)
1907# 467 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1909# 467 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1912 dvel(
i, 1) = dvel_small(
i)
1917 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
1918 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(1))/3._wp
1921 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
1922 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(1))/3._wp
1927# 483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1928#if defined(MFC_OpenACC)
1929# 483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1931# 483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1932#elif defined(MFC_OpenMP)
1933# 483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1935# 483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1940# 486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1941#if defined(MFC_OpenACC)
1942# 486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1944# 486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1945#elif defined(MFC_OpenMP)
1946# 486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1948# 486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1950 do r = 1, num_fluids
1953 rho_sf_small(
i) = rho_l
1963# 499 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1964#if defined(MFC_OpenACC)
1965# 499 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1967# 499 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1968#elif defined(MFC_OpenMP)
1969# 499 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1971# 499 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1974 dvel(
i, 2) = dvel_small(
i)
1979 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
1980 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(2))/3._wp
1983 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
1984 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(2))/3._wp
1989 & 1)**2._wp + dvel(2, 2)**2._wp + (dvel(1, 1) + dvel(2, 2))**2._wp), kind=stp)
1994# 520 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1995#if defined(MFC_OpenACC)
1996# 520 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1998# 520 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
1999#elif defined(MFC_OpenMP)
2000# 520 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2002# 520 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2004 do i = 1, num_fluids
2005 alpha_rho_l(
i) = 0._wp
2006 alpha_rho_r(
i) = 0._wp
2011# 527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2012#if defined(MFC_OpenACC)
2013# 527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2015# 527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2016#elif defined(MFC_OpenMP)
2017# 527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2019# 527 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2027# 533 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2028#if defined(MFC_OpenACC)
2029# 533 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2031# 533 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2032#elif defined(MFC_OpenMP)
2033# 533 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2035# 533 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2039# 535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2040#if defined(MFC_OpenACC)
2041# 535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2043# 535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2044#elif defined(MFC_OpenMP)
2045# 535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2047# 535 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2049 do i = 1, num_fluids
2053 if (num_fluids > 1)
then
2055# 541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2056#if defined(MFC_OpenACC)
2057# 541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2059# 541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2060#elif defined(MFC_OpenMP)
2061# 541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2063# 541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2065 do i = 1, num_fluids - 1
2073# 549 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2074#if defined(MFC_OpenACC)
2075# 549 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2077# 549 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2078#elif defined(MFC_OpenMP)
2079# 549 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2081# 549 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2089# 555 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2090#if defined(MFC_OpenACC)
2091# 555 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2093# 555 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2094#elif defined(MFC_OpenMP)
2095# 555 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2097# 555 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2101# 557 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2102#if defined(MFC_OpenACC)
2103# 557 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2105# 557 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2106#elif defined(MFC_OpenMP)
2107# 557 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2109# 557 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2111 do i = 1, num_fluids
2115 if (num_fluids > 1)
then
2117# 563 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2118#if defined(MFC_OpenACC)
2119# 563 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2121# 563 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2122#elif defined(MFC_OpenMP)
2123# 563 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2125# 563 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2127 do i = 1, num_fluids - 1
2135# 571 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2136#if defined(MFC_OpenACC)
2137# 571 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2139# 571 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2140#elif defined(MFC_OpenMP)
2141# 571 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2143# 571 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2150 if (num_fluids > 1)
then
2151 alpha_l(num_fluids) = 1._wp
2152 alpha_r(num_fluids) = 1._wp
2155# 581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2156#if defined(MFC_OpenACC)
2157# 581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2159# 581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2160#elif defined(MFC_OpenMP)
2161# 581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2163# 581 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2165 do i = 1, num_fluids - 1
2166 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
2167 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
2171 call s_compute_mixture_coefficients(alpha_rho_l, alpha_l, rho_l, gamma_l, pi_inf_l, qv_igr)
2172 call s_compute_mixture_coefficients(alpha_rho_r, alpha_r, rho_r, gamma_r, pi_inf_r, qv_igr)
2175# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2176#if defined(MFC_OpenACC)
2177# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2179# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2180#elif defined(MFC_OpenMP)
2181# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2183# 591 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2186 vel_l(
i) = vel_l(
i)/rho_l
2187 vel_r(
i) = vel_r(
i)/rho_r
2191 mu_l = 0._wp; mu_r = 0._wp
2193# 599 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2194#if defined(MFC_OpenACC)
2195# 599 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2197# 599 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2198#elif defined(MFC_OpenMP)
2199# 599 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2201# 599 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2203 do i = 1, num_fluids
2209# 605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2210#if defined(MFC_OpenACC)
2211# 605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2213# 605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2214#elif defined(MFC_OpenMP)
2215# 605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2217# 605 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2220 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
2222# 608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2223#if defined(MFC_OpenACC)
2224# 608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2226# 608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2227#elif defined(MFC_OpenMP)
2228# 608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2230# 608 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2233 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j + 1)), kind=stp)
2236# 612 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2237#if defined(MFC_OpenACC)
2238# 612 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2240# 612 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2241#elif defined(MFC_OpenMP)
2242# 612 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2244# 612 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2247 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j)), kind=stp)
2249# 615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2250#if defined(MFC_OpenACC)
2251# 615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2253# 615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2254#elif defined(MFC_OpenMP)
2255# 615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2257# 615 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2260 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j)), kind=stp)
2263# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2264#if defined(MFC_OpenACC)
2265# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2267# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2268#elif defined(MFC_OpenMP)
2269# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2271# 619 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2274 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
2276# 622 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2277#if defined(MFC_OpenACC)
2278# 622 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2280# 622 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2281#elif defined(MFC_OpenMP)
2282# 622 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2284# 622 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2287 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j + 1)), kind=stp)
2290# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2291#if defined(MFC_OpenACC)
2292# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2294# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2295#elif defined(MFC_OpenMP)
2296# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2298# 626 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2301 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j)), kind=stp)
2303# 629 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2304#if defined(MFC_OpenACC)
2305# 629 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2307# 629 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2308#elif defined(MFC_OpenMP)
2309# 629 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2311# 629 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2314 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j)), kind=stp)
2317# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2318#if defined(MFC_OpenACC)
2319# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2321# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2322#elif defined(MFC_OpenMP)
2323# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2325# 633 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2328 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
2330# 636 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2331#if defined(MFC_OpenACC)
2332# 636 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2334# 636 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2335#elif defined(MFC_OpenMP)
2336# 636 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2338# 636 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2341 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j + 1)), kind=stp)
2344# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2345#if defined(MFC_OpenACC)
2346# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2348# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2349#elif defined(MFC_OpenMP)
2350# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2352# 640 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2355 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j)), kind=stp)
2357# 643 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2358#if defined(MFC_OpenACC)
2359# 643 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2361# 643 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2362#elif defined(MFC_OpenMP)
2363# 643 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2365# 643 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2368 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j)), kind=stp)
2371# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2372#if defined(MFC_OpenACC)
2373# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2375# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2376#elif defined(MFC_OpenMP)
2377# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2379# 647 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2382 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
2384# 650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2385#if defined(MFC_OpenACC)
2386# 650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2388# 650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2389#elif defined(MFC_OpenMP)
2390# 650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2392# 650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2395 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j + 1)), kind=stp)
2398# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2399#if defined(MFC_OpenACC)
2400# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2402# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2403#elif defined(MFC_OpenMP)
2404# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2406# 654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2409 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j)), kind=stp)
2411# 657 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2412#if defined(MFC_OpenACC)
2413# 657 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2415# 657 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2416#elif defined(MFC_OpenMP)
2417# 657 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2419# 657 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2422 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j)), kind=stp)
2425 e_l = 0._wp; e_r = 0._wp
2428# 664 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2429#if defined(MFC_OpenACC)
2430# 664 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2432# 664 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2433#elif defined(MFC_OpenMP)
2434# 664 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2436# 664 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2443# 669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2444#if defined(MFC_OpenACC)
2445# 669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2447# 669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2448#elif defined(MFC_OpenMP)
2449# 669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2451# 669 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2457 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, &
2458 & vel_r, pres_l, pres_r, cfl)
2460 do i = 1, num_fluids
2462# 678 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2463#if defined(MFC_OpenACC)
2464# 678 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2466# 678 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2467#elif defined(MFC_OpenMP)
2468# 678 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2470# 678 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2472 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
2473 &
l) + real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
2474 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
2477# 683 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2478#if defined(MFC_OpenACC)
2479# 683 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2481# 683 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2482#elif defined(MFC_OpenMP)
2483# 683 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2485# 683 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2487 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
2488 &
l) - real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
2489 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j))), kind=stp)
2492 if (num_fluids > 1)
then
2494# 690 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2495#if defined(MFC_OpenACC)
2496# 690 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2498# 690 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2499#elif defined(MFC_OpenMP)
2500# 690 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2502# 690 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2504 do i = 1, num_fluids - 1
2506# 692 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2507#if defined(MFC_OpenACC)
2508# 692 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2510# 692 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2511#elif defined(MFC_OpenMP)
2512# 692 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2514# 692 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2517 &
l) + real((0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
2518 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
2521# 697 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2522#if defined(MFC_OpenACC)
2523# 697 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2525# 697 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2526#elif defined(MFC_OpenMP)
2527# 697 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2529# 697 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2533 &
l)*vel_l(1)*(1._wp/dx(
j + 1))), kind=stp)
2536# 702 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2537#if defined(MFC_OpenACC)
2538# 702 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2540# 702 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2541#elif defined(MFC_OpenMP)
2542# 702 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2544# 702 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2547 &
l) - real((0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
2548 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j))), kind=stp)
2551# 707 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2552#if defined(MFC_OpenACC)
2553# 707 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2555# 707 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2556#elif defined(MFC_OpenMP)
2557# 707 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2559# 707 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2563 &
l)*vel_l(1)*(1._wp/dx(
j))), kind=stp)
2568# 714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2569#if defined(MFC_OpenACC)
2570# 714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2572# 714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2573#elif defined(MFC_OpenMP)
2574# 714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2576# 714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2579 &
l) + real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j + 1)) &
2580 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j + 1))), kind=stp)
2583# 719 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2584#if defined(MFC_OpenACC)
2585# 719 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2587# 719 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2588#elif defined(MFC_OpenMP)
2589# 719 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2591# 719 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2594 &
l) + real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j + 1)) &
2595 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j + 1))), kind=stp)
2598# 724 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2599#if defined(MFC_OpenACC)
2600# 724 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2602# 724 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2603#elif defined(MFC_OpenMP)
2604# 724 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2606# 724 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2609 &
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) &
2610 & *(1._wp/dx(
j + 1))), kind=stp)
2613# 729 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2614#if defined(MFC_OpenACC)
2615# 729 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2617# 729 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2618#elif defined(MFC_OpenMP)
2619# 729 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2621# 729 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2624 &
l) - real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j)) &
2625 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j))), kind=stp)
2628# 734 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2629#if defined(MFC_OpenACC)
2630# 734 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2632# 734 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2633#elif defined(MFC_OpenMP)
2634# 734 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2636# 734 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2639 &
l) - real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j)) &
2640 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j))), kind=stp)
2643# 739 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2644#if defined(MFC_OpenACC)
2645# 739 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2647# 739 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2648#elif defined(MFC_OpenMP)
2649# 739 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2651# 739 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2654 &
l) - real((0.5_wp*dt*(vel_l(1)*(e_l + pres_l))*(1._wp/dx(
j)) - 0.5_wp*dt*cfl*(e_l) &
2655 & *(1._wp/dx(
j))), kind=stp)
2658# 744 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2659#if defined(MFC_OpenACC)
2660# 744 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2662# 744 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2663#elif defined(MFC_OpenMP)
2664# 744 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2666# 744 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2668 do i = 1, num_fluids
2670# 746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2671#if defined(MFC_OpenACC)
2672# 746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2674# 746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2675#elif defined(MFC_OpenMP)
2676# 746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2678# 746 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2680 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
2681 &
l) + real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
2682 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
2685# 751 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2686#if defined(MFC_OpenACC)
2687# 751 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2689# 751 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2690#elif defined(MFC_OpenMP)
2691# 751 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2693# 751 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2695 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
2696 &
l) - real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
2697 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j))), kind=stp)
2700 if (num_fluids > 1)
then
2702# 758 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2703#if defined(MFC_OpenACC)
2704# 758 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2706# 758 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2707#elif defined(MFC_OpenMP)
2708# 758 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2710# 758 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2712 do i = 1, num_fluids - 1
2714# 760 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2715#if defined(MFC_OpenACC)
2716# 760 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2718# 760 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2719#elif defined(MFC_OpenMP)
2720# 760 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2722# 760 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2725 &
l) + real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
2726 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
2729# 765 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2730#if defined(MFC_OpenACC)
2731# 765 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2733# 765 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2734#elif defined(MFC_OpenMP)
2735# 765 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2737# 765 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2741 &
l)*vel_r(1)*(1._wp/dx(
j + 1))), kind=stp)
2744# 770 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2745#if defined(MFC_OpenACC)
2746# 770 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2748# 770 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2749#elif defined(MFC_OpenMP)
2750# 770 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2752# 770 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2755 &
l) - real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
2756 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j))), kind=stp)
2759# 775 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2760#if defined(MFC_OpenACC)
2761# 775 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2763# 775 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2764#elif defined(MFC_OpenMP)
2765# 775 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2767# 775 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2771 &
l)*vel_r(1)*(1._wp/dx(
j))), kind=stp)
2776# 782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2777#if defined(MFC_OpenACC)
2778# 782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2780# 782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2781#elif defined(MFC_OpenMP)
2782# 782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2784# 782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2787 &
l) + real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j + 1)) &
2788 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j + 1))), kind=stp)
2791# 787 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2792#if defined(MFC_OpenACC)
2793# 787 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2795# 787 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2796#elif defined(MFC_OpenMP)
2797# 787 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2799# 787 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2802 &
l) + real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j + 1)) &
2803 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j + 1))), kind=stp)
2806# 792 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2807#if defined(MFC_OpenACC)
2808# 792 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2810# 792 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2811#elif defined(MFC_OpenMP)
2812# 792 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2814# 792 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2817 &
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) &
2818 & *(1._wp/dx(
j + 1))), kind=stp)
2821# 797 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2822#if defined(MFC_OpenACC)
2823# 797 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2825# 797 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2826#elif defined(MFC_OpenMP)
2827# 797 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2829# 797 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2832 &
l) - real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j)) &
2833 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j))), kind=stp)
2836# 802 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2837#if defined(MFC_OpenACC)
2838# 802 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2840# 802 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2841#elif defined(MFC_OpenMP)
2842# 802 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2844# 802 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2847 &
l) - real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j)) &
2848 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j))), kind=stp)
2851# 807 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2852#if defined(MFC_OpenACC)
2853# 807 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2855# 807 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2856#elif defined(MFC_OpenMP)
2857# 807 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2859# 807 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2862 &
l) - real((0.5_wp*dt*(vel_r(1)*(e_r + pres_r))*(1._wp/dx(
j)) + 0.5_wp*dt*cfl*(e_r) &
2863 & *(1._wp/dx(
j))), kind=stp)
2868# 814 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2869#if defined(MFC_OpenACC)
2870# 814 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2872# 814 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2873#elif defined(MFC_OpenMP)
2874# 814 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2876# 814 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2878# 814 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2880# 816 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2882# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2884# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2886# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2887#if defined(MFC_OpenACC)
2888# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2890# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2892# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2893#elif defined(MFC_OpenMP)
2894# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2896# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2898# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2900# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2902# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2904# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2906# 818 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2908# 821 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2915# 834 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2917# 834 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2918#if defined(MFC_OpenACC)
2919# 834 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2921# 834 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2922#elif defined(MFC_OpenMP)
2923# 834 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2925# 834 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2930# 837 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2931#if defined(MFC_OpenACC)
2932# 837 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2934# 837 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2935#elif defined(MFC_OpenMP)
2936# 837 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2938# 837 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2943# 840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2944#if defined(MFC_OpenACC)
2945# 840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2947# 840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2948#elif defined(MFC_OpenMP)
2949# 840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2951# 840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2953 do r = 1, num_fluids
2956 rho_sf_small(
i) = rho_l
2968# 855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2969#if defined(MFC_OpenACC)
2970# 855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2972# 855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2973#elif defined(MFC_OpenMP)
2974# 855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2976# 855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2979 dvel(
i, 1) = dvel_small(
i)
2984 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
2985 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(3))
2986 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(1))/3._wp
2989 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
2990 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(3))
2991 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(1))/3._wp
2996# 873 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
2997#if defined(MFC_OpenACC)
2998# 873 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3000# 873 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3001#elif defined(MFC_OpenMP)
3002# 873 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3004# 873 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3009# 876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3010#if defined(MFC_OpenACC)
3011# 876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3013# 876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3014#elif defined(MFC_OpenMP)
3015# 876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3017# 876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3019 do r = 1, num_fluids
3022 rho_sf_small(
i) = rho_l
3033# 890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3034#if defined(MFC_OpenACC)
3035# 890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3037# 890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3038#elif defined(MFC_OpenMP)
3039# 890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3041# 890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3044 dvel(
i, 2) = dvel_small(
i)
3049 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
3050 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(2))/3._wp
3053 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
3054 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(2))/3._wp
3059# 906 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3060#if defined(MFC_OpenACC)
3061# 906 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3063# 906 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3064#elif defined(MFC_OpenMP)
3065# 906 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3067# 906 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3072# 909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3073#if defined(MFC_OpenACC)
3074# 909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3076# 909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3077#elif defined(MFC_OpenMP)
3078# 909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3080# 909 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3082 do r = 1, num_fluids
3085 rho_sf_small(
i) = rho_l
3094 &
l - 1)/rho_sf_small(-1))
3097# 924 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3098#if defined(MFC_OpenACC)
3099# 924 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3101# 924 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3102#elif defined(MFC_OpenMP)
3103# 924 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3105# 924 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3108 dvel(
i, 3) = dvel_small(
i)
3113 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(1))
3114 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(3))/3._wp
3117 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(1))
3118 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(3))/3._wp
3122 jac_rhs(
j,
k,
l) = real(
alf_igr*(2._wp*(dvel(1, 2)*dvel(2, 1) + dvel(1, 3)*dvel(3, &
3123 & 1) + dvel(2, 3)*dvel(3, 2)) + dvel(1, 1)**2._wp + dvel(2, 2)**2._wp + dvel(3, &
3124 & 3)**2._wp + (dvel(1, 1) + dvel(2, 2) + dvel(3, 3))**2._wp), kind=stp)
3129# 946 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3130#if defined(MFC_OpenACC)
3131# 946 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3133# 946 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3134#elif defined(MFC_OpenMP)
3135# 946 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3137# 946 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3139 do i = 1, num_fluids
3140 alpha_rho_l(
i) = 0._wp
3141 alpha_rho_r(
i) = 0._wp
3146# 953 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3147#if defined(MFC_OpenACC)
3148# 953 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3150# 953 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3151#elif defined(MFC_OpenMP)
3152# 953 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3154# 953 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3162# 959 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3163#if defined(MFC_OpenACC)
3164# 959 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3166# 959 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3167#elif defined(MFC_OpenMP)
3168# 959 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3170# 959 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3174# 961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3175#if defined(MFC_OpenACC)
3176# 961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3178# 961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3179#elif defined(MFC_OpenMP)
3180# 961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3182# 961 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3184 do i = 1, num_fluids
3188 if (num_fluids > 1)
then
3190# 967 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3191#if defined(MFC_OpenACC)
3192# 967 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3194# 967 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3195#elif defined(MFC_OpenMP)
3196# 967 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3198# 967 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3200 do i = 1, num_fluids - 1
3208# 975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3209#if defined(MFC_OpenACC)
3210# 975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3212# 975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3213#elif defined(MFC_OpenMP)
3214# 975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3216# 975 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3224# 981 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3225#if defined(MFC_OpenACC)
3226# 981 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3228# 981 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3229#elif defined(MFC_OpenMP)
3230# 981 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3232# 981 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3236# 983 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3237#if defined(MFC_OpenACC)
3238# 983 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3240# 983 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3241#elif defined(MFC_OpenMP)
3242# 983 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3244# 983 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3246 do i = 1, num_fluids
3250 if (num_fluids > 1)
then
3252# 989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3253#if defined(MFC_OpenACC)
3254# 989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3256# 989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3257#elif defined(MFC_OpenMP)
3258# 989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3260# 989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3262 do i = 1, num_fluids - 1
3270# 997 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3271#if defined(MFC_OpenACC)
3272# 997 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3274# 997 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3275#elif defined(MFC_OpenMP)
3276# 997 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3278# 997 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3285 if (num_fluids > 1)
then
3286 alpha_l(num_fluids) = 1._wp
3287 alpha_r(num_fluids) = 1._wp
3290# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3291#if defined(MFC_OpenACC)
3292# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3294# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3295#elif defined(MFC_OpenMP)
3296# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3298# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3300 do i = 1, num_fluids - 1
3301 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
3302 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
3306 call s_compute_mixture_coefficients(alpha_rho_l, alpha_l, rho_l, gamma_l, pi_inf_l, qv_igr)
3307 call s_compute_mixture_coefficients(alpha_rho_r, alpha_r, rho_r, gamma_r, pi_inf_r, qv_igr)
3310# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3311#if defined(MFC_OpenACC)
3312# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3314# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3315#elif defined(MFC_OpenMP)
3316# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3318# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3321 vel_l(
i) = vel_l(
i)/rho_l
3322 vel_r(
i) = vel_r(
i)/rho_r
3329# 1026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3330#if defined(MFC_OpenACC)
3331# 1026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3333# 1026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3334#elif defined(MFC_OpenMP)
3335# 1026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3337# 1026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3339 do i = 1, num_fluids
3345# 1032 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3346#if defined(MFC_OpenACC)
3347# 1032 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3349# 1032 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3350#elif defined(MFC_OpenMP)
3351# 1032 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3353# 1032 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3356 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
3358# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3359#if defined(MFC_OpenACC)
3360# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3362# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3363#elif defined(MFC_OpenMP)
3364# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3366# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3369 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j + 1)), kind=stp)
3372# 1039 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3373#if defined(MFC_OpenACC)
3374# 1039 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3376# 1039 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3377#elif defined(MFC_OpenMP)
3378# 1039 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3380# 1039 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3383 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dx(
j)), kind=stp)
3385# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3386#if defined(MFC_OpenACC)
3387# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3389# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3390#elif defined(MFC_OpenMP)
3391# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3393# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3396 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(2)*(1._wp/dx(
j)), kind=stp)
3399# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3400#if defined(MFC_OpenACC)
3401# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3403# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3404#elif defined(MFC_OpenMP)
3405# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3407# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3410 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j + 1)), kind=stp)
3412# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3413#if defined(MFC_OpenACC)
3414# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3416# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3417#elif defined(MFC_OpenMP)
3418# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3420# 1049 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3423 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j + 1)), kind=stp)
3426# 1053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3427#if defined(MFC_OpenACC)
3428# 1053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3430# 1053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3431#elif defined(MFC_OpenMP)
3432# 1053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3434# 1053 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3437 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dx(
j)), kind=stp)
3439# 1056 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3440#if defined(MFC_OpenACC)
3441# 1056 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3443# 1056 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3444#elif defined(MFC_OpenMP)
3445# 1056 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3447# 1056 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3450 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(2)*(1._wp/dx(
j)), kind=stp)
3453# 1060 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3454#if defined(MFC_OpenACC)
3455# 1060 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3457# 1060 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3458#elif defined(MFC_OpenMP)
3459# 1060 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3461# 1060 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3464 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dx(
j + 1)), kind=stp)
3466# 1063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3467#if defined(MFC_OpenACC)
3468# 1063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3470# 1063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3471#elif defined(MFC_OpenMP)
3472# 1063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3474# 1063 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3477 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dx(
j + 1)), kind=stp)
3480# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3481#if defined(MFC_OpenACC)
3482# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3484# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3485#elif defined(MFC_OpenMP)
3486# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3488# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3491 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dx(
j)), kind=stp)
3493# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3494#if defined(MFC_OpenACC)
3495# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3497# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3498#elif defined(MFC_OpenMP)
3499# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3501# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3504 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dx(
j)), kind=stp)
3507# 1074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3508#if defined(MFC_OpenACC)
3509# 1074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3511# 1074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3512#elif defined(MFC_OpenMP)
3513# 1074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3515# 1074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3518 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dx(
j + 1)), kind=stp)
3520# 1077 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3521#if defined(MFC_OpenACC)
3522# 1077 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3524# 1077 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3525#elif defined(MFC_OpenMP)
3526# 1077 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3528# 1077 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3531 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dx(
j + 1)), kind=stp)
3534# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3535#if defined(MFC_OpenACC)
3536# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3538# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3539#elif defined(MFC_OpenMP)
3540# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3542# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3545 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dx(
j)), kind=stp)
3547# 1084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3548#if defined(MFC_OpenACC)
3549# 1084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3551# 1084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3552#elif defined(MFC_OpenMP)
3553# 1084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3555# 1084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3558 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dx(
j)), kind=stp)
3561# 1088 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3562#if defined(MFC_OpenACC)
3563# 1088 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3565# 1088 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3566#elif defined(MFC_OpenMP)
3567# 1088 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3569# 1088 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3572 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
3574# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3575#if defined(MFC_OpenACC)
3576# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3578# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3579#elif defined(MFC_OpenMP)
3580# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3582# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3585 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j + 1)), kind=stp)
3588# 1095 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3589#if defined(MFC_OpenACC)
3590# 1095 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3592# 1095 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3593#elif defined(MFC_OpenMP)
3594# 1095 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3596# 1095 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3599 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dx(
j)), kind=stp)
3601# 1098 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3602#if defined(MFC_OpenACC)
3603# 1098 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3605# 1098 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3606#elif defined(MFC_OpenMP)
3607# 1098 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3609# 1098 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3612 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(1)*(1._wp/dx(
j)), kind=stp)
3615# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3616#if defined(MFC_OpenACC)
3617# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3619# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3620#elif defined(MFC_OpenMP)
3621# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3623# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3626 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j + 1)), kind=stp)
3628# 1105 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3629#if defined(MFC_OpenACC)
3630# 1105 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3632# 1105 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3633#elif defined(MFC_OpenMP)
3634# 1105 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3636# 1105 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3639 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j + 1)), kind=stp)
3642# 1109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3643#if defined(MFC_OpenACC)
3644# 1109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3646# 1109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3647#elif defined(MFC_OpenMP)
3648# 1109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3650# 1109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3653 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dx(
j)), kind=stp)
3655# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3656#if defined(MFC_OpenACC)
3657# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3659# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3660#elif defined(MFC_OpenMP)
3661# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3663# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3666 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(1)*(1._wp/dx(
j)), kind=stp)
3669 e_l = 0._wp; e_r = 0._wp
3672# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3673#if defined(MFC_OpenACC)
3674# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3676# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3677#elif defined(MFC_OpenMP)
3678# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3680# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3687# 1124 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3688#if defined(MFC_OpenACC)
3689# 1124 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3691# 1124 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3692#elif defined(MFC_OpenMP)
3693# 1124 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3695# 1124 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3701 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, &
3702 & vel_r, pres_l, pres_r, cfl)
3705# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3706#if defined(MFC_OpenACC)
3707# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3709# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3710#elif defined(MFC_OpenMP)
3711# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3713# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3715 do i = 1, num_fluids
3717# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3718#if defined(MFC_OpenACC)
3719# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3721# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3722#elif defined(MFC_OpenMP)
3723# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3725# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3727 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
3728 &
l) + real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
3729 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
3732# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3733#if defined(MFC_OpenACC)
3734# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3736# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3737#elif defined(MFC_OpenMP)
3738# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3740# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3742 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
3743 &
l) - real((0.5_wp*dt*(alpha_rho_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
3744 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dx(
j))), kind=stp)
3747 if (num_fluids > 1)
then
3749# 1146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3750#if defined(MFC_OpenACC)
3751# 1146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3753# 1146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3754#elif defined(MFC_OpenMP)
3755# 1146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3757# 1146 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3759 do i = 1, num_fluids - 1
3761# 1148 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3762#if defined(MFC_OpenACC)
3763# 1148 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3765# 1148 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3766#elif defined(MFC_OpenMP)
3767# 1148 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3769# 1148 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3772 &
l) + real((0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j + 1)) &
3773 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j + 1))), kind=stp)
3776# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3777#if defined(MFC_OpenACC)
3778# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3780# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3781#elif defined(MFC_OpenMP)
3782# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3784# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3788 &
l)*vel_l(1)*(1._wp/dx(
j + 1)), kind=stp)
3791# 1158 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3792#if defined(MFC_OpenACC)
3793# 1158 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3795# 1158 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3796#elif defined(MFC_OpenMP)
3797# 1158 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3799# 1158 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3802 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(1))*(1._wp/dx(
j)) &
3803 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dx(
j)), kind=stp)
3806# 1163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3807#if defined(MFC_OpenACC)
3808# 1163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3810# 1163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3811#elif defined(MFC_OpenMP)
3812# 1163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3814# 1163 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3818 &
l)*vel_l(1)*(1._wp/dx(
j)), kind=stp)
3823# 1170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3824#if defined(MFC_OpenACC)
3825# 1170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3827# 1170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3828#elif defined(MFC_OpenMP)
3829# 1170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3831# 1170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3834 &
l) + real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j + 1)) &
3835 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j + 1))), kind=stp)
3838# 1175 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3839#if defined(MFC_OpenACC)
3840# 1175 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3842# 1175 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3843#elif defined(MFC_OpenMP)
3844# 1175 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3846# 1175 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3849 &
l) + real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j + 1)) &
3850 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j + 1))), kind=stp)
3853# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3854#if defined(MFC_OpenACC)
3855# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3857# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3858#elif defined(MFC_OpenMP)
3859# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3861# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3864 &
l) + real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(3)*(1._wp/dx(
j + 1)) &
3865 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dx(
j + 1))), kind=stp)
3868# 1185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3869#if defined(MFC_OpenACC)
3870# 1185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3872# 1185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3873#elif defined(MFC_OpenMP)
3874# 1185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3876# 1185 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3879 &
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) &
3880 & *(1._wp/dx(
j + 1))), kind=stp)
3883# 1190 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3884#if defined(MFC_OpenACC)
3885# 1190 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3887# 1190 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3888#elif defined(MFC_OpenMP)
3889# 1190 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3891# 1190 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3894 &
l) - real((0.5_wp*dt*(rho_l*(vel_l(1))**2.0 + pres_l)*(1._wp/dx(
j)) &
3895 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dx(
j))), kind=stp)
3898# 1195 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3899#if defined(MFC_OpenACC)
3900# 1195 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3902# 1195 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3903#elif defined(MFC_OpenMP)
3904# 1195 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3906# 1195 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3909 &
l) - real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dx(
j)) &
3910 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dx(
j))), kind=stp)
3913# 1200 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3914#if defined(MFC_OpenACC)
3915# 1200 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3917# 1200 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3918#elif defined(MFC_OpenMP)
3919# 1200 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3921# 1200 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3924 &
l) - real((0.5_wp*dt*rho_l*vel_l(1)*vel_l(3)*(1._wp/dx(
j)) &
3925 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dx(
j))), kind=stp)
3928# 1205 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3929#if defined(MFC_OpenACC)
3930# 1205 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3932# 1205 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3933#elif defined(MFC_OpenMP)
3934# 1205 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3936# 1205 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3939 &
l) - real((0.5_wp*dt*(vel_l(1)*(e_l + pres_l))*(1._wp/dx(
j)) - 0.5_wp*dt*cfl*(e_l) &
3940 & *(1._wp/dx(
j))), kind=stp)
3943# 1210 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3944#if defined(MFC_OpenACC)
3945# 1210 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3947# 1210 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3948#elif defined(MFC_OpenMP)
3949# 1210 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3951# 1210 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3953 do i = 1, num_fluids
3955# 1212 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3956#if defined(MFC_OpenACC)
3957# 1212 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3959# 1212 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3960#elif defined(MFC_OpenMP)
3961# 1212 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3963# 1212 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3965 rhs_vf(
i)%sf(
j + 1,
k,
l) = rhs_vf(
i)%sf(
j + 1,
k, &
3966 &
l) + real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
3967 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
3970# 1217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3971#if defined(MFC_OpenACC)
3972# 1217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3974# 1217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3975#elif defined(MFC_OpenMP)
3976# 1217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3978# 1217 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3980 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
3981 &
l) - real((0.5_wp*dt*(alpha_rho_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
3982 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dx(
j))), kind=stp)
3985 if (num_fluids > 1)
then
3987# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3988#if defined(MFC_OpenACC)
3989# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3991# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3992#elif defined(MFC_OpenMP)
3993# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3995# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
3997 do i = 1, num_fluids - 1
3999# 1226 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4000#if defined(MFC_OpenACC)
4001# 1226 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4003# 1226 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4004#elif defined(MFC_OpenMP)
4005# 1226 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4007# 1226 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4010 &
l) + real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j + 1)) &
4011 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j + 1))), kind=stp)
4014# 1231 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4015#if defined(MFC_OpenACC)
4016# 1231 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4018# 1231 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4019#elif defined(MFC_OpenMP)
4020# 1231 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4022# 1231 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4026 &
l)*vel_r(1)*(1._wp/dx(
j + 1))), kind=stp)
4029# 1236 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4030#if defined(MFC_OpenACC)
4031# 1236 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4033# 1236 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4034#elif defined(MFC_OpenMP)
4035# 1236 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4037# 1236 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4040 &
l) - real((0.5_wp*dt*(alpha_r(
i)*vel_r(1))*(1._wp/dx(
j)) &
4041 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dx(
j))), kind=stp)
4044# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4045#if defined(MFC_OpenACC)
4046# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4048# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4049#elif defined(MFC_OpenMP)
4050# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4052# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4056 &
l)*vel_r(1)*(1._wp/dx(
j))), kind=stp)
4061# 1248 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4062#if defined(MFC_OpenACC)
4063# 1248 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4065# 1248 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4066#elif defined(MFC_OpenMP)
4067# 1248 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4069# 1248 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4072 &
l) + real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j + 1)) &
4073 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j + 1))), kind=stp)
4076# 1253 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4077#if defined(MFC_OpenACC)
4078# 1253 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4080# 1253 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4081#elif defined(MFC_OpenMP)
4082# 1253 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4084# 1253 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4087 &
l) + real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j + 1)) &
4088 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j + 1))), kind=stp)
4091# 1258 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4092#if defined(MFC_OpenACC)
4093# 1258 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4095# 1258 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4096#elif defined(MFC_OpenMP)
4097# 1258 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4099# 1258 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4102 &
l) + real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(3)*(1._wp/dx(
j + 1)) &
4103 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dx(
j + 1))), kind=stp)
4106# 1263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4107#if defined(MFC_OpenACC)
4108# 1263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4110# 1263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4111#elif defined(MFC_OpenMP)
4112# 1263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4114# 1263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4117 &
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) &
4118 & *(1._wp/dx(
j + 1))), kind=stp)
4121# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4122#if defined(MFC_OpenACC)
4123# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4125# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4126#elif defined(MFC_OpenMP)
4127# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4129# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4132 &
l) - real((0.5_wp*dt*(rho_r*(vel_r(1))**2.0 + pres_r)*(1._wp/dx(
j)) &
4133 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dx(
j))), kind=stp)
4136# 1273 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4137#if defined(MFC_OpenACC)
4138# 1273 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4140# 1273 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4141#elif defined(MFC_OpenMP)
4142# 1273 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4144# 1273 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4147 &
l) - real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(2)*(1._wp/dx(
j)) &
4148 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dx(
j))), kind=stp)
4151# 1278 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4152#if defined(MFC_OpenACC)
4153# 1278 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4155# 1278 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4156#elif defined(MFC_OpenMP)
4157# 1278 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4159# 1278 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4162 &
l) - real((0.5_wp*dt*rho_r*vel_r(1)*vel_r(3)*(1._wp/dx(
j)) &
4163 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dx(
j))), kind=stp)
4166# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4167#if defined(MFC_OpenACC)
4168# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4170# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4171#elif defined(MFC_OpenMP)
4172# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4174# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4177 &
l) - real((0.5_wp*dt*(vel_r(1)*(e_r + pres_r))*(1._wp/dx(
j)) + 0.5_wp*dt*cfl*(e_r) &
4178 & *(1._wp/dx(
j))), kind=stp)
4183# 1290 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4184#if defined(MFC_OpenACC)
4185# 1290 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4187# 1290 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4188#elif defined(MFC_OpenMP)
4189# 1290 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4191# 1290 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4193# 1290 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4195# 1292 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4197 else if (idir == 2)
then
4199# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4201# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4203# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4204#if defined(MFC_OpenACC)
4205# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4207# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4209# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4210#elif defined(MFC_OpenMP)
4211# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4213# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4215# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4217# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4219# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4221# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4223# 1296 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4225# 1299 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4233# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4235# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4236#if defined(MFC_OpenACC)
4237# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4239# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4240#elif defined(MFC_OpenMP)
4241# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4243# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4248# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4249#if defined(MFC_OpenACC)
4250# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4252# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4253#elif defined(MFC_OpenMP)
4254# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4256# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4261# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4262#if defined(MFC_OpenACC)
4263# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4265# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4266#elif defined(MFC_OpenMP)
4267# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4269# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4271 do r = 1, num_fluids
4274 rho_sf_small(
i) = rho_l
4281 &
l)/rho_sf_small(-1))
4284 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
4285 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(1))/3._wp
4288 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
4289 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(1))/3._wp
4294# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4295#if defined(MFC_OpenACC)
4296# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4298# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4299#elif defined(MFC_OpenMP)
4300# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4302# 1342 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4307# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4308#if defined(MFC_OpenACC)
4309# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4311# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4312#elif defined(MFC_OpenMP)
4313# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4315# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4317 do r = 1, num_fluids
4320 rho_sf_small(
i) = rho_l
4327 &
l)/rho_sf_small(-1))
4330 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
4331 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(2))/3._wp
4334 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
4335 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(2))/3._wp
4341# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4342#if defined(MFC_OpenACC)
4343# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4345# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4346#elif defined(MFC_OpenMP)
4347# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4349# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4351 do i = 1, num_fluids
4352 alpha_rho_l(
i) = 0._wp
4353 alpha_rho_r(
i) = 0._wp
4358# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4359#if defined(MFC_OpenACC)
4360# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4362# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4363#elif defined(MFC_OpenMP)
4364# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4366# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4374# 1382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4375#if defined(MFC_OpenACC)
4376# 1382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4378# 1382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4379#elif defined(MFC_OpenMP)
4380# 1382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4382# 1382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4386# 1384 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4387#if defined(MFC_OpenACC)
4388# 1384 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4390# 1384 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4391#elif defined(MFC_OpenMP)
4392# 1384 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4394# 1384 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4396 do i = 1, num_fluids
4400 if (num_fluids > 1)
then
4402# 1390 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4403#if defined(MFC_OpenACC)
4404# 1390 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4406# 1390 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4407#elif defined(MFC_OpenMP)
4408# 1390 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4410# 1390 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4412 do i = 1, num_fluids - 1
4420# 1398 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4421#if defined(MFC_OpenACC)
4422# 1398 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4424# 1398 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4425#elif defined(MFC_OpenMP)
4426# 1398 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4428# 1398 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4436# 1404 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4437#if defined(MFC_OpenACC)
4438# 1404 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4440# 1404 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4441#elif defined(MFC_OpenMP)
4442# 1404 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4444# 1404 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4448# 1406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4449#if defined(MFC_OpenACC)
4450# 1406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4452# 1406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4453#elif defined(MFC_OpenMP)
4454# 1406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4456# 1406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4458 do i = 1, num_fluids
4462 if (num_fluids > 1)
then
4464# 1412 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4465#if defined(MFC_OpenACC)
4466# 1412 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4468# 1412 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4469#elif defined(MFC_OpenMP)
4470# 1412 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4472# 1412 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4474 do i = 1, num_fluids - 1
4482# 1420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4483#if defined(MFC_OpenACC)
4484# 1420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4486# 1420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4487#elif defined(MFC_OpenMP)
4488# 1420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4490# 1420 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4497 if (num_fluids > 1)
then
4498 alpha_l(num_fluids) = 1._wp
4499 alpha_r(num_fluids) = 1._wp
4502# 1430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4503#if defined(MFC_OpenACC)
4504# 1430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4506# 1430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4507#elif defined(MFC_OpenMP)
4508# 1430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4510# 1430 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4512 do i = 1, num_fluids - 1
4513 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
4514 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
4518 call s_compute_mixture_coefficients(alpha_rho_l, alpha_l, rho_l, gamma_l, pi_inf_l, qv_igr)
4519 call s_compute_mixture_coefficients(alpha_rho_r, alpha_r, rho_r, gamma_r, pi_inf_r, qv_igr)
4522# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4523#if defined(MFC_OpenACC)
4524# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4526# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4527#elif defined(MFC_OpenMP)
4528# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4530# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4533 vel_l(
i) = vel_l(
i)/rho_l
4534 vel_r(
i) = vel_r(
i)/rho_r
4541# 1449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4542#if defined(MFC_OpenACC)
4543# 1449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4545# 1449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4546#elif defined(MFC_OpenMP)
4547# 1449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4549# 1449 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4551 do i = 1, num_fluids
4557# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4558#if defined(MFC_OpenACC)
4559# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4561# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4562#elif defined(MFC_OpenMP)
4563# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4565# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4568 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
4570# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4571#if defined(MFC_OpenACC)
4572# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4574# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4575#elif defined(MFC_OpenMP)
4576# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4578# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4581 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k + 1)), kind=stp)
4584# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4585#if defined(MFC_OpenACC)
4586# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4588# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4589#elif defined(MFC_OpenMP)
4590# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4592# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4595 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k)), kind=stp)
4597# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4598#if defined(MFC_OpenACC)
4599# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4601# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4602#elif defined(MFC_OpenMP)
4603# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4605# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4608 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k)), kind=stp)
4611# 1469 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4612#if defined(MFC_OpenACC)
4613# 1469 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4615# 1469 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4616#elif defined(MFC_OpenMP)
4617# 1469 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4619# 1469 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4622 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
4624# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4625#if defined(MFC_OpenACC)
4626# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4628# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4629#elif defined(MFC_OpenMP)
4630# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4632# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4635 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k + 1)), kind=stp)
4638# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4639#if defined(MFC_OpenACC)
4640# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4642# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4643#elif defined(MFC_OpenMP)
4644# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4646# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4649 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k)), kind=stp)
4651# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4652#if defined(MFC_OpenACC)
4653# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4655# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4656#elif defined(MFC_OpenMP)
4657# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4659# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4662 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k)), kind=stp)
4665# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4666#if defined(MFC_OpenACC)
4667# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4669# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4670#elif defined(MFC_OpenMP)
4671# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4673# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4676 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
4678# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4679#if defined(MFC_OpenACC)
4680# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4682# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4683#elif defined(MFC_OpenMP)
4684# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4686# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4689 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
4692# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4693#if defined(MFC_OpenACC)
4694# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4696# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4697#elif defined(MFC_OpenMP)
4698# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4700# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4703 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k)), kind=stp)
4705# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4706#if defined(MFC_OpenACC)
4707# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4709# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4710#elif defined(MFC_OpenMP)
4711# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4713# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4716 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
4719# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4720#if defined(MFC_OpenACC)
4721# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4723# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4724#elif defined(MFC_OpenMP)
4725# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4727# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4730 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
4732# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4733#if defined(MFC_OpenACC)
4734# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4736# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4737#elif defined(MFC_OpenMP)
4738# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4740# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4743 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
4746# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4747#if defined(MFC_OpenACC)
4748# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4750# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4751#elif defined(MFC_OpenMP)
4752# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4754# 1504 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4757 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k)), kind=stp)
4759# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4760#if defined(MFC_OpenACC)
4761# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4763# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4764#elif defined(MFC_OpenMP)
4765# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4767# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4770 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
4773 e_l = 0._wp; e_r = 0._wp
4774 f_l = 0._wp; f_r = 0._wp
4777# 1515 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4778#if defined(MFC_OpenACC)
4779# 1515 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4781# 1515 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4782#elif defined(MFC_OpenMP)
4783# 1515 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4785# 1515 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4793# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4794#if defined(MFC_OpenACC)
4795# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4797# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4798#elif defined(MFC_OpenMP)
4799# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4801# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4808 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, &
4809 & vel_r, pres_l, pres_r, cfl)
4812# 1530 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4813#if defined(MFC_OpenACC)
4814# 1530 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4816# 1530 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4817#elif defined(MFC_OpenMP)
4818# 1530 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4820# 1530 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4822 do i = 1, num_fluids
4824# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4825#if defined(MFC_OpenACC)
4826# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4828# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4829#elif defined(MFC_OpenMP)
4830# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4832# 1532 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4834 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
4835 &
l) + real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
4836 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
4839# 1537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4840#if defined(MFC_OpenACC)
4841# 1537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4843# 1537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4844#elif defined(MFC_OpenMP)
4845# 1537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4847# 1537 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4849 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
4850 &
l) - real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
4851 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k)), kind=stp)
4854 if (num_fluids > 1)
then
4856# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4857#if defined(MFC_OpenACC)
4858# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4860# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4861#elif defined(MFC_OpenMP)
4862# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4864# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4866 do i = 1, num_fluids - 1
4868# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4869#if defined(MFC_OpenACC)
4870# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4872# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4873#elif defined(MFC_OpenMP)
4874# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4876# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4879 &
l) + real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
4880 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
4883# 1551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4884#if defined(MFC_OpenACC)
4885# 1551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4887# 1551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4888#elif defined(MFC_OpenMP)
4889# 1551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4891# 1551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4895 &
l)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
4898# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4899#if defined(MFC_OpenACC)
4900# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4902# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4903#elif defined(MFC_OpenMP)
4904# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4906# 1556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4909 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
4910 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k)), kind=stp)
4913# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4914#if defined(MFC_OpenACC)
4915# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4917# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4918#elif defined(MFC_OpenMP)
4919# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4921# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4925 &
l)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
4930# 1568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4931#if defined(MFC_OpenACC)
4932# 1568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4934# 1568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4935#elif defined(MFC_OpenMP)
4936# 1568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4938# 1568 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4941 &
l) + real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k + 1)) &
4942 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k + 1)), kind=stp)
4945# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4946#if defined(MFC_OpenACC)
4947# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4949# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4950#elif defined(MFC_OpenMP)
4951# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4953# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4956 &
l) + real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dy(
k + 1)) &
4957 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dy(
k + 1)), kind=stp)
4960# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4961#if defined(MFC_OpenACC)
4962# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4964# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4965#elif defined(MFC_OpenMP)
4966# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4968# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4971 &
l) + real(0.5_wp*dt*(vel_l(2)*(e_l + pres_l + f_l))*(1._wp/dy(
k + 1)) &
4972 & - 0.5_wp*dt*cfl*(e_l)*(1._wp/dy(
k + 1)), kind=stp)
4975# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4976#if defined(MFC_OpenACC)
4977# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4979# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4980#elif defined(MFC_OpenMP)
4981# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4983# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4986 &
l) - real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k)) &
4987 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k)), kind=stp)
4990# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4991#if defined(MFC_OpenACC)
4992# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4994# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4995#elif defined(MFC_OpenMP)
4996# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
4998# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5001 &
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) &
5002 & )*(1._wp/dy(
k)), kind=stp)
5005# 1593 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5006#if defined(MFC_OpenACC)
5007# 1593 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5009# 1593 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5010#elif defined(MFC_OpenMP)
5011# 1593 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5013# 1593 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5016 &
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) &
5017 & *(1._wp/dy(
k)), kind=stp)
5020# 1598 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5021#if defined(MFC_OpenACC)
5022# 1598 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5024# 1598 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5025#elif defined(MFC_OpenMP)
5026# 1598 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5028# 1598 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5030 do i = 1, num_fluids
5032# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5033#if defined(MFC_OpenACC)
5034# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5036# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5037#elif defined(MFC_OpenMP)
5038# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5040# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5042 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
5043 &
l) + real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
5044 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
5046# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5047#if defined(MFC_OpenACC)
5048# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5050# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5051#elif defined(MFC_OpenMP)
5052# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5054# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5056 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
5057 &
l) - real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
5058 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k)), kind=stp)
5061 if (num_fluids > 1)
then
5063# 1611 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5064#if defined(MFC_OpenACC)
5065# 1611 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5067# 1611 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5068#elif defined(MFC_OpenMP)
5069# 1611 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5071# 1611 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5073 do i = 1, num_fluids - 1
5075# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5076#if defined(MFC_OpenACC)
5077# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5079# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5080#elif defined(MFC_OpenMP)
5081# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5083# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5086 &
l) + real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
5087 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
5090# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5091#if defined(MFC_OpenACC)
5092# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5094# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5095#elif defined(MFC_OpenMP)
5096# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5098# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5102 &
l)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
5105# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5106#if defined(MFC_OpenACC)
5107# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5109# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5110#elif defined(MFC_OpenMP)
5111# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5113# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5116 &
l) - real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
5117 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k)), kind=stp)
5120# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5121#if defined(MFC_OpenACC)
5122# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5124# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5125#elif defined(MFC_OpenMP)
5126# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5128# 1628 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5132 &
l)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
5136# 1634 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5137#if defined(MFC_OpenACC)
5138# 1634 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5140# 1634 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5141#elif defined(MFC_OpenMP)
5142# 1634 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5144# 1634 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5147 &
l) + real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k + 1)) &
5148 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k + 1)), kind=stp)
5150# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5151#if defined(MFC_OpenACC)
5152# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5154# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5155#elif defined(MFC_OpenMP)
5156# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5158# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5161 &
l) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(1)*(1._wp/dy(
k + 1)) &
5162 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dy(
k + 1)), kind=stp)
5164# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5165#if defined(MFC_OpenACC)
5166# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5168# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5169#elif defined(MFC_OpenMP)
5170# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5172# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5175 &
l) + real(0.5_wp*dt*(vel_r(2)*(e_r + pres_r + f_r))*(1._wp/dy(
k + 1)) &
5176 & + 0.5_wp*dt*cfl*(e_r)*(1._wp/dy(
k + 1)), kind=stp)
5178# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5179#if defined(MFC_OpenACC)
5180# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5182# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5183#elif defined(MFC_OpenMP)
5184# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5186# 1646 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5189 &
l) - real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k)) &
5190 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k)), kind=stp)
5192# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5193#if defined(MFC_OpenACC)
5194# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5196# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5197#elif defined(MFC_OpenMP)
5198# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5200# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5203 &
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) &
5204 & )*(1._wp/dy(
k)), kind=stp)
5206# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5207#if defined(MFC_OpenACC)
5208# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5210# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5211#elif defined(MFC_OpenMP)
5212# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5214# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5217 &
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) &
5218 & *(1._wp/dy(
k)), kind=stp)
5223# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5224#if defined(MFC_OpenACC)
5225# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5227# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5228#elif defined(MFC_OpenMP)
5229# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5231# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5233# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5235# 1663 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5237# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5239# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5241# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5242#if defined(MFC_OpenACC)
5243# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5245# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5247# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5248#elif defined(MFC_OpenMP)
5249# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5251# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5253# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5255# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5257# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5259# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5261# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5263# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5271# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5273# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5274#if defined(MFC_OpenACC)
5275# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5277# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5278#elif defined(MFC_OpenMP)
5279# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5281# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5286# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5287#if defined(MFC_OpenACC)
5288# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5290# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5291#elif defined(MFC_OpenMP)
5292# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5294# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5299# 1688 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5300#if defined(MFC_OpenACC)
5301# 1688 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5303# 1688 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5304#elif defined(MFC_OpenMP)
5305# 1688 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5307# 1688 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5309 do r = 1, num_fluids
5312 rho_sf_small(
i) = rho_l
5319 &
l)/rho_sf_small(-1))
5322 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(2))
5323 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(1))/3._wp
5326 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(2))
5327 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(1))/3._wp
5332# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5333#if defined(MFC_OpenACC)
5334# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5336# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5337#elif defined(MFC_OpenMP)
5338# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5340# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5345# 1714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5346#if defined(MFC_OpenACC)
5347# 1714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5349# 1714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5350#elif defined(MFC_OpenMP)
5351# 1714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5353# 1714 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5355 do r = 1, num_fluids
5358 rho_sf_small(
i) = rho_l
5365 &
l)/rho_sf_small(-1))
5368 &
l)/rho_sf_small(-1))
5371 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
5372 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(3))
5373 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(2))/3._wp
5376 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
5377 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(3))
5378 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(2))/3._wp
5383# 1742 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5384#if defined(MFC_OpenACC)
5385# 1742 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5387# 1742 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5388#elif defined(MFC_OpenMP)
5389# 1742 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5391# 1742 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5396# 1745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5397#if defined(MFC_OpenACC)
5398# 1745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5400# 1745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5401#elif defined(MFC_OpenMP)
5402# 1745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5404# 1745 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5406 do r = 1, num_fluids
5409 rho_sf_small(
i) = rho_l
5414 &
l - 1)/rho_sf_small(-1))
5417 &
l - 1)/rho_sf_small(-1))
5419 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(2))
5420 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(3))/3._wp
5423 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(2))
5424 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(3))/3._wp
5430# 1769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5431#if defined(MFC_OpenACC)
5432# 1769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5434# 1769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5435#elif defined(MFC_OpenMP)
5436# 1769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5438# 1769 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5440 do i = 1, num_fluids
5441 alpha_rho_l(
i) = 0._wp
5442 alpha_rho_r(
i) = 0._wp
5447# 1776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5448#if defined(MFC_OpenACC)
5449# 1776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5451# 1776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5452#elif defined(MFC_OpenMP)
5453# 1776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5455# 1776 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5463# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5464#if defined(MFC_OpenACC)
5465# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5467# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5468#elif defined(MFC_OpenMP)
5469# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5471# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5475# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5476#if defined(MFC_OpenACC)
5477# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5479# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5480#elif defined(MFC_OpenMP)
5481# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5483# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5485 do i = 1, num_fluids
5489 if (num_fluids > 1)
then
5491# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5492#if defined(MFC_OpenACC)
5493# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5495# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5496#elif defined(MFC_OpenMP)
5497# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5499# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5501 do i = 1, num_fluids - 1
5509# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5510#if defined(MFC_OpenACC)
5511# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5513# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5514#elif defined(MFC_OpenMP)
5515# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5517# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5525# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5526#if defined(MFC_OpenACC)
5527# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5529# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5530#elif defined(MFC_OpenMP)
5531# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5533# 1804 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5537# 1806 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5538#if defined(MFC_OpenACC)
5539# 1806 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5541# 1806 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5542#elif defined(MFC_OpenMP)
5543# 1806 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5545# 1806 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5547 do i = 1, num_fluids
5551 if (num_fluids > 1)
then
5553# 1812 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5554#if defined(MFC_OpenACC)
5555# 1812 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5557# 1812 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5558#elif defined(MFC_OpenMP)
5559# 1812 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5561# 1812 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5563 do i = 1, num_fluids - 1
5571# 1820 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5572#if defined(MFC_OpenACC)
5573# 1820 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5575# 1820 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5576#elif defined(MFC_OpenMP)
5577# 1820 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5579# 1820 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5586 if (num_fluids > 1)
then
5587 alpha_l(num_fluids) = 1._wp
5588 alpha_r(num_fluids) = 1._wp
5591# 1830 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5592#if defined(MFC_OpenACC)
5593# 1830 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5595# 1830 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5596#elif defined(MFC_OpenMP)
5597# 1830 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5599# 1830 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5601 do i = 1, num_fluids - 1
5602 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
5603 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
5607 call s_compute_mixture_coefficients(alpha_rho_l, alpha_l, rho_l, gamma_l, pi_inf_l, qv_igr)
5608 call s_compute_mixture_coefficients(alpha_rho_r, alpha_r, rho_r, gamma_r, pi_inf_r, qv_igr)
5611# 1840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5612#if defined(MFC_OpenACC)
5613# 1840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5615# 1840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5616#elif defined(MFC_OpenMP)
5617# 1840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5619# 1840 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5622 vel_l(
i) = vel_l(
i)/rho_l
5623 vel_r(
i) = vel_r(
i)/rho_r
5630# 1849 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5631#if defined(MFC_OpenACC)
5632# 1849 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5634# 1849 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5635#elif defined(MFC_OpenMP)
5636# 1849 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5638# 1849 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5640 do i = 1, num_fluids
5646# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5647#if defined(MFC_OpenACC)
5648# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5650# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5651#elif defined(MFC_OpenMP)
5652# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5654# 1855 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5657 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
5659# 1858 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5660#if defined(MFC_OpenACC)
5661# 1858 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5663# 1858 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5664#elif defined(MFC_OpenMP)
5665# 1858 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5667# 1858 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5670 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k + 1)), kind=stp)
5673# 1862 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5674#if defined(MFC_OpenACC)
5675# 1862 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5677# 1862 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5678#elif defined(MFC_OpenMP)
5679# 1862 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5681# 1862 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5684 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dy(
k)), kind=stp)
5686# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5687#if defined(MFC_OpenACC)
5688# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5690# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5691#elif defined(MFC_OpenMP)
5692# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5694# 1865 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5697 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dy(
k)), kind=stp)
5700# 1869 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5701#if defined(MFC_OpenACC)
5702# 1869 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5704# 1869 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5705#elif defined(MFC_OpenMP)
5706# 1869 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5708# 1869 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5711 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k + 1)), kind=stp)
5713# 1872 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5714#if defined(MFC_OpenACC)
5715# 1872 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5717# 1872 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5718#elif defined(MFC_OpenMP)
5719# 1872 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5721# 1872 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5724 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k + 1)), kind=stp)
5727# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5728#if defined(MFC_OpenACC)
5729# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5731# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5732#elif defined(MFC_OpenMP)
5733# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5735# 1876 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5738 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dy(
k)), kind=stp)
5740# 1879 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5741#if defined(MFC_OpenACC)
5742# 1879 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5744# 1879 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5745#elif defined(MFC_OpenMP)
5746# 1879 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5748# 1879 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5751 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dy(
k)), kind=stp)
5754# 1883 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5755#if defined(MFC_OpenACC)
5756# 1883 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5758# 1883 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5759#elif defined(MFC_OpenMP)
5760# 1883 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5762# 1883 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5765 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dy(
k + 1)), kind=stp)
5767# 1886 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5768#if defined(MFC_OpenACC)
5769# 1886 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5771# 1886 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5772#elif defined(MFC_OpenMP)
5773# 1886 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5775# 1886 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5778 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dy(
k + 1)), kind=stp)
5781# 1890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5782#if defined(MFC_OpenACC)
5783# 1890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5785# 1890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5786#elif defined(MFC_OpenMP)
5787# 1890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5789# 1890 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5792 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dy(
k)), kind=stp)
5794# 1893 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5795#if defined(MFC_OpenACC)
5796# 1893 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5798# 1893 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5799#elif defined(MFC_OpenMP)
5800# 1893 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5802# 1893 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5805 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(3)*(1._wp/dy(
k)), kind=stp)
5808# 1897 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5809#if defined(MFC_OpenACC)
5810# 1897 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5812# 1897 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5813#elif defined(MFC_OpenMP)
5814# 1897 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5816# 1897 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5819 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dy(
k + 1)), kind=stp)
5821# 1900 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5822#if defined(MFC_OpenACC)
5823# 1900 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5825# 1900 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5826#elif defined(MFC_OpenMP)
5827# 1900 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5829# 1900 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5832 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dy(
k + 1)), kind=stp)
5835# 1904 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5836#if defined(MFC_OpenACC)
5837# 1904 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5839# 1904 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5840#elif defined(MFC_OpenMP)
5841# 1904 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5843# 1904 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5846 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dy(
k)), kind=stp)
5848# 1907 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5849#if defined(MFC_OpenACC)
5850# 1907 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5852# 1907 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5853#elif defined(MFC_OpenMP)
5854# 1907 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5856# 1907 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5859 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(3)*(1._wp/dy(
k)), kind=stp)
5862# 1911 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5863#if defined(MFC_OpenACC)
5864# 1911 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5866# 1911 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5867#elif defined(MFC_OpenMP)
5868# 1911 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5870# 1911 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5873 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
5875# 1914 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5876#if defined(MFC_OpenACC)
5877# 1914 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5879# 1914 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5880#elif defined(MFC_OpenMP)
5881# 1914 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5883# 1914 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5886 &
l) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
5889# 1918 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5890#if defined(MFC_OpenACC)
5891# 1918 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5893# 1918 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5894#elif defined(MFC_OpenMP)
5895# 1918 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5897# 1918 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5900 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dy(
k)), kind=stp)
5902# 1921 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5903#if defined(MFC_OpenACC)
5904# 1921 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5906# 1921 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5907#elif defined(MFC_OpenMP)
5908# 1921 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5910# 1921 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5913 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
5916# 1925 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5917#if defined(MFC_OpenACC)
5918# 1925 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5920# 1925 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5921#elif defined(MFC_OpenMP)
5922# 1925 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5924# 1925 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5927 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k + 1)), kind=stp)
5929# 1928 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5930#if defined(MFC_OpenACC)
5931# 1928 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5933# 1928 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5934#elif defined(MFC_OpenMP)
5935# 1928 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5937# 1928 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5940 &
l) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
5943# 1932 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5944#if defined(MFC_OpenACC)
5945# 1932 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5947# 1932 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5948#elif defined(MFC_OpenMP)
5949# 1932 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5951# 1932 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5954 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dy(
k)), kind=stp)
5956# 1935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5957#if defined(MFC_OpenACC)
5958# 1935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5960# 1935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5961#elif defined(MFC_OpenMP)
5962# 1935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5964# 1935 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5967 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
5970 e_l = 0._wp; e_r = 0._wp
5971 f_l = 0._wp; f_r = 0._wp
5974# 1943 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5975#if defined(MFC_OpenACC)
5976# 1943 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5978# 1943 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5979#elif defined(MFC_OpenMP)
5980# 1943 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5982# 1943 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5990# 1949 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5991#if defined(MFC_OpenACC)
5992# 1949 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5994# 1949 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5995#elif defined(MFC_OpenMP)
5996# 1949 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
5998# 1949 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6005 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, &
6006 & vel_r, pres_l, pres_r, cfl)
6009# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6010#if defined(MFC_OpenACC)
6011# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6013# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6014#elif defined(MFC_OpenMP)
6015# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6017# 1958 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6019 do i = 1, num_fluids
6021# 1960 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6022#if defined(MFC_OpenACC)
6023# 1960 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6025# 1960 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6026#elif defined(MFC_OpenMP)
6027# 1960 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6029# 1960 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6031 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
6032 &
l) + real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
6033 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
6036# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6037#if defined(MFC_OpenACC)
6038# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6040# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6041#elif defined(MFC_OpenMP)
6042# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6044# 1965 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6046 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
6047 &
l) - real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
6048 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dy(
k)), kind=stp)
6051 if (num_fluids > 1)
then
6053# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6054#if defined(MFC_OpenACC)
6055# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6057# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6058#elif defined(MFC_OpenMP)
6059# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6061# 1972 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6063 do i = 1, num_fluids - 1
6065# 1974 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6066#if defined(MFC_OpenACC)
6067# 1974 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6069# 1974 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6070#elif defined(MFC_OpenMP)
6071# 1974 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6073# 1974 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6076 &
l) + real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k + 1)) &
6077 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k + 1)), kind=stp)
6080# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6081#if defined(MFC_OpenACC)
6082# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6084# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6085#elif defined(MFC_OpenMP)
6086# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6088# 1979 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6092 &
l)*vel_l(2)*(1._wp/dy(
k + 1)), kind=stp)
6095# 1984 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6096#if defined(MFC_OpenACC)
6097# 1984 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6099# 1984 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6100#elif defined(MFC_OpenMP)
6101# 1984 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6103# 1984 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6106 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(2))*(1._wp/dy(
k)) &
6107 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dy(
k)), kind=stp)
6110# 1989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6111#if defined(MFC_OpenACC)
6112# 1989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6114# 1989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6115#elif defined(MFC_OpenMP)
6116# 1989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6118# 1989 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6122 &
l)*vel_l(2)*(1._wp/dy(
k)), kind=stp)
6127# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6128#if defined(MFC_OpenACC)
6129# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6131# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6132#elif defined(MFC_OpenMP)
6133# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6135# 1996 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6138 &
l) + real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k + 1)) &
6139 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k + 1)), kind=stp)
6142# 2001 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6143#if defined(MFC_OpenACC)
6144# 2001 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6146# 2001 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6147#elif defined(MFC_OpenMP)
6148# 2001 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6150# 2001 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6153 &
l) + real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(2)*(1._wp/dy(
k + 1)) &
6154 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dy(
k + 1)), kind=stp)
6157# 2006 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6158#if defined(MFC_OpenACC)
6159# 2006 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6161# 2006 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6162#elif defined(MFC_OpenMP)
6163# 2006 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6165# 2006 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6168 &
l) + real(0.5_wp*dt*rho_l*vel_l(3)*vel_l(2)*(1._wp/dy(
k + 1)) &
6169 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dy(
k + 1)), kind=stp)
6172# 2011 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6173#if defined(MFC_OpenACC)
6174# 2011 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6176# 2011 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6177#elif defined(MFC_OpenMP)
6178# 2011 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6180# 2011 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6183 &
l) + real(0.5_wp*dt*(vel_l(2)*(e_l + pres_l + f_l))*(1._wp/dy(
k + 1)) &
6184 & - 0.5_wp*dt*cfl*(e_l)*(1._wp/dy(
k + 1)), kind=stp)
6187# 2016 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6188#if defined(MFC_OpenACC)
6189# 2016 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6191# 2016 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6192#elif defined(MFC_OpenMP)
6193# 2016 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6195# 2016 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6198 &
l) - real(0.5_wp*dt*(rho_l*(vel_l(2))**2.0 + pres_l + f_l)*(1._wp/dy(
k)) &
6199 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dy(
k)), kind=stp)
6202# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6203#if defined(MFC_OpenACC)
6204# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6206# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6207#elif defined(MFC_OpenMP)
6208# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6210# 2021 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6213 &
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) &
6214 & )*(1._wp/dy(
k)), kind=stp)
6217# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6218#if defined(MFC_OpenACC)
6219# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6221# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6222#elif defined(MFC_OpenMP)
6223# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6225# 2026 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6228 &
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) &
6229 & )*(1._wp/dy(
k)), kind=stp)
6232# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6233#if defined(MFC_OpenACC)
6234# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6236# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6237#elif defined(MFC_OpenMP)
6238# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6240# 2031 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6243 &
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) &
6244 & *(1._wp/dy(
k)), kind=stp)
6247# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6248#if defined(MFC_OpenACC)
6249# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6251# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6252#elif defined(MFC_OpenMP)
6253# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6255# 2036 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6257 do i = 1, num_fluids
6259# 2038 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6260#if defined(MFC_OpenACC)
6261# 2038 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6263# 2038 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6264#elif defined(MFC_OpenMP)
6265# 2038 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6267# 2038 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6269 rhs_vf(
i)%sf(
j,
k + 1,
l) = rhs_vf(
i)%sf(
j,
k + 1, &
6270 &
l) + real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
6271 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
6274# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6275#if defined(MFC_OpenACC)
6276# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6278# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6279#elif defined(MFC_OpenMP)
6280# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6282# 2043 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6284 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
6285 &
l) - real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
6286 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dy(
k)), kind=stp)
6289 if (num_fluids > 1)
then
6291# 2050 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6292#if defined(MFC_OpenACC)
6293# 2050 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6295# 2050 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6296#elif defined(MFC_OpenMP)
6297# 2050 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6299# 2050 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6301 do i = 1, num_fluids - 1
6303# 2052 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6304#if defined(MFC_OpenACC)
6305# 2052 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6307# 2052 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6308#elif defined(MFC_OpenMP)
6309# 2052 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6311# 2052 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6314 &
l) + real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k + 1)) &
6315 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k + 1)), kind=stp)
6318# 2057 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6319#if defined(MFC_OpenACC)
6320# 2057 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6322# 2057 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6323#elif defined(MFC_OpenMP)
6324# 2057 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6326# 2057 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6330 &
l)*vel_r(2)*(1._wp/dy(
k + 1)), kind=stp)
6333# 2062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6334#if defined(MFC_OpenACC)
6335# 2062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6337# 2062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6338#elif defined(MFC_OpenMP)
6339# 2062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6341# 2062 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6344 &
l) - real(0.5_wp*dt*(alpha_r(
i)*vel_r(2))*(1._wp/dy(
k)) &
6345 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dy(
k)), kind=stp)
6348# 2067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6349#if defined(MFC_OpenACC)
6350# 2067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6352# 2067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6353#elif defined(MFC_OpenMP)
6354# 2067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6356# 2067 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6360 &
l)*vel_r(2)*(1._wp/dy(
k)), kind=stp)
6365# 2074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6366#if defined(MFC_OpenACC)
6367# 2074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6369# 2074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6370#elif defined(MFC_OpenMP)
6371# 2074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6373# 2074 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6376 &
l) + real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k + 1)) &
6377 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k + 1)), kind=stp)
6380# 2079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6381#if defined(MFC_OpenACC)
6382# 2079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6384# 2079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6385#elif defined(MFC_OpenMP)
6386# 2079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6388# 2079 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6391 &
l) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(1)*(1._wp/dy(
k + 1)) &
6392 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dy(
k + 1)), kind=stp)
6395# 2084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6396#if defined(MFC_OpenACC)
6397# 2084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6399# 2084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6400#elif defined(MFC_OpenMP)
6401# 2084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6403# 2084 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6406 &
l) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(3)*(1._wp/dy(
k + 1)) &
6407 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dy(
k + 1)), kind=stp)
6410# 2089 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6411#if defined(MFC_OpenACC)
6412# 2089 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6414# 2089 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6415#elif defined(MFC_OpenMP)
6416# 2089 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6418# 2089 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6421 &
l) + real(0.5_wp*dt*(vel_r(2)*(e_r + pres_r + f_r))*(1._wp/dy(
k + 1)) &
6422 & + 0.5_wp*dt*cfl*(e_r)*(1._wp/dy(
k + 1)), kind=stp)
6425# 2094 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6426#if defined(MFC_OpenACC)
6427# 2094 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6429# 2094 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6430#elif defined(MFC_OpenMP)
6431# 2094 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6433# 2094 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6436 &
l) - real(0.5_wp*dt*(rho_r*(vel_r(2))**2.0 + pres_r + f_r)*(1._wp/dy(
k)) &
6437 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dy(
k)), kind=stp)
6440# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6441#if defined(MFC_OpenACC)
6442# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6444# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6445#elif defined(MFC_OpenMP)
6446# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6448# 2099 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6451 &
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) &
6452 & )*(1._wp/dy(
k)), kind=stp)
6455# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6456#if defined(MFC_OpenACC)
6457# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6459# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6460#elif defined(MFC_OpenMP)
6461# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6463# 2104 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6466 &
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) &
6467 & )*(1._wp/dy(
k)), kind=stp)
6470# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6471#if defined(MFC_OpenACC)
6472# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6474# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6475#elif defined(MFC_OpenMP)
6476# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6478# 2109 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6481 &
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) &
6482 & *(1._wp/dy(
k)), kind=stp)
6487# 2116 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6488#if defined(MFC_OpenACC)
6489# 2116 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6491# 2116 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6492#elif defined(MFC_OpenMP)
6493# 2116 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6495# 2116 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6497# 2116 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6499# 2118 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6501 else if (idir == 3)
then
6502# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6504# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6506# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6507#if defined(MFC_OpenACC)
6508# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6510# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6512# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6513#elif defined(MFC_OpenMP)
6514# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6516# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6518# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6520# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6522# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6524# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6526# 2121 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6528# 2124 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6536# 2138 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6538# 2138 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6539#if defined(MFC_OpenACC)
6540# 2138 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6542# 2138 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6543#elif defined(MFC_OpenMP)
6544# 2138 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6546# 2138 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6551# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6552#if defined(MFC_OpenACC)
6553# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6555# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6556#elif defined(MFC_OpenMP)
6557# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6559# 2141 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6564# 2144 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6565#if defined(MFC_OpenACC)
6566# 2144 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6568# 2144 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6569#elif defined(MFC_OpenMP)
6570# 2144 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6572# 2144 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6574 do r = 1, num_fluids
6577 rho_sf_small(
i) = rho_l
6584 &
l +
q)/rho_sf_small(-1))
6587 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(3))
6588 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(1))/3._wp
6591 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(3))
6592 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(1))/3._wp
6597# 2167 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6598#if defined(MFC_OpenACC)
6599# 2167 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6601# 2167 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6602#elif defined(MFC_OpenMP)
6603# 2167 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6605# 2167 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6610# 2170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6611#if defined(MFC_OpenACC)
6612# 2170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6614# 2170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6615#elif defined(MFC_OpenMP)
6616# 2170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6618# 2170 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6620 do r = 1, num_fluids
6623 rho_sf_small(
i) = rho_l
6628 &
l +
q)/rho_sf_small(-1))
6631 &
l +
q)/rho_sf_small(-1))
6634 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(3))
6635 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(-2._wp*dvel_small(2))/3._wp
6638 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(3))
6639 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(-2._wp*dvel_small(2))/3._wp
6644# 2194 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6645#if defined(MFC_OpenACC)
6646# 2194 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6648# 2194 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6649#elif defined(MFC_OpenMP)
6650# 2194 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6652# 2194 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6657# 2197 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6658#if defined(MFC_OpenACC)
6659# 2197 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6661# 2197 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6662#elif defined(MFC_OpenMP)
6663# 2197 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6665# 2197 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6667 do r = 1, num_fluids
6670 rho_sf_small(
i) = rho_l
6674 &
l - 1 +
q)/rho_sf_small(-1))
6677 &
l - 1 +
q)/rho_sf_small(-1))
6680 &
l - 1 +
q)/rho_sf_small(-1))
6682 vflux_l_arr(1) = vflux_l_arr(1) +
coeff_l(
q)*(dvel_small(1))
6683 vflux_l_arr(2) = vflux_l_arr(2) +
coeff_l(
q)*(dvel_small(2))
6684 vflux_l_arr(3) = vflux_l_arr(3) +
coeff_l(
q)*(4._wp*dvel_small(3))/3._wp
6687 vflux_r_arr(1) = vflux_r_arr(1) +
coeff_r(
q)*(dvel_small(1))
6688 vflux_r_arr(2) = vflux_r_arr(2) +
coeff_r(
q)*(dvel_small(2))
6689 vflux_r_arr(3) = vflux_r_arr(3) +
coeff_r(
q)*(4._wp*dvel_small(3))/3._wp
6695# 2225 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6696#if defined(MFC_OpenACC)
6697# 2225 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6699# 2225 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6700#elif defined(MFC_OpenMP)
6701# 2225 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6703# 2225 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6705 do i = 1, num_fluids
6706 alpha_rho_l(
i) = 0._wp
6707 alpha_rho_r(
i) = 0._wp
6713# 2233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6714#if defined(MFC_OpenACC)
6715# 2233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6717# 2233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6718#elif defined(MFC_OpenMP)
6719# 2233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6721# 2233 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6729# 2239 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6730#if defined(MFC_OpenACC)
6731# 2239 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6733# 2239 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6734#elif defined(MFC_OpenMP)
6735# 2239 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6737# 2239 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6741# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6742#if defined(MFC_OpenACC)
6743# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6745# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6746#elif defined(MFC_OpenMP)
6747# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6749# 2241 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6751 do i = 1, num_fluids
6755 if (num_fluids > 1)
then
6757# 2247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6758#if defined(MFC_OpenACC)
6759# 2247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6761# 2247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6762#elif defined(MFC_OpenMP)
6763# 2247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6765# 2247 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6767 do i = 1, num_fluids - 1
6775# 2255 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6776#if defined(MFC_OpenACC)
6777# 2255 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6779# 2255 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6780#elif defined(MFC_OpenMP)
6781# 2255 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6783# 2255 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6791# 2261 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6792#if defined(MFC_OpenACC)
6793# 2261 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6795# 2261 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6796#elif defined(MFC_OpenMP)
6797# 2261 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6799# 2261 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6803# 2263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6804#if defined(MFC_OpenACC)
6805# 2263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6807# 2263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6808#elif defined(MFC_OpenMP)
6809# 2263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6811# 2263 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6813 do i = 1, num_fluids
6817 if (num_fluids > 1)
then
6819# 2269 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6820#if defined(MFC_OpenACC)
6821# 2269 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6823# 2269 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6824#elif defined(MFC_OpenMP)
6825# 2269 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6827# 2269 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6829 do i = 1, num_fluids - 1
6837# 2277 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6838#if defined(MFC_OpenACC)
6839# 2277 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6841# 2277 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6842#elif defined(MFC_OpenMP)
6843# 2277 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6845# 2277 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6852 if (num_fluids > 1)
then
6853 alpha_l(num_fluids) = 1._wp
6854 alpha_r(num_fluids) = 1._wp
6857# 2287 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6858#if defined(MFC_OpenACC)
6859# 2287 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6861# 2287 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6862#elif defined(MFC_OpenMP)
6863# 2287 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6865# 2287 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6867 do i = 1, num_fluids - 1
6868 alpha_l(num_fluids) = alpha_l(num_fluids) - alpha_l(
i)
6869 alpha_r(num_fluids) = alpha_r(num_fluids) - alpha_r(
i)
6873 call s_compute_mixture_coefficients(alpha_rho_l, alpha_l, rho_l, gamma_l, pi_inf_l, qv_igr)
6874 call s_compute_mixture_coefficients(alpha_rho_r, alpha_r, rho_r, gamma_r, pi_inf_r, qv_igr)
6877# 2297 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6878#if defined(MFC_OpenACC)
6879# 2297 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6881# 2297 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6882#elif defined(MFC_OpenMP)
6883# 2297 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6885# 2297 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6888 vel_l(
i) = vel_l(
i)/rho_l
6889 vel_r(
i) = vel_r(
i)/rho_r
6896# 2306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6897#if defined(MFC_OpenACC)
6898# 2306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6900# 2306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6901#elif defined(MFC_OpenMP)
6902# 2306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6904# 2306 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6906 do i = 1, num_fluids
6912# 2312 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6913#if defined(MFC_OpenACC)
6914# 2312 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6916# 2312 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6917#elif defined(MFC_OpenMP)
6918# 2312 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6920# 2312 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6923 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dz(
l + 1)), kind=stp)
6925# 2315 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6926#if defined(MFC_OpenACC)
6927# 2315 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6929# 2315 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6930#elif defined(MFC_OpenMP)
6931# 2315 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6933# 2315 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6936 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dz(
l + 1)), kind=stp)
6939# 2319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6940#if defined(MFC_OpenACC)
6941# 2319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6943# 2319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6944#elif defined(MFC_OpenMP)
6945# 2319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6947# 2319 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6950 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*(1._wp/dz(
l)), kind=stp)
6952# 2322 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6953#if defined(MFC_OpenACC)
6954# 2322 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6956# 2322 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6957#elif defined(MFC_OpenMP)
6958# 2322 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6960# 2322 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6963 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(1)*vel_l(1)*(1._wp/dz(
l)), kind=stp)
6966# 2326 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6967#if defined(MFC_OpenACC)
6968# 2326 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6970# 2326 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6971#elif defined(MFC_OpenMP)
6972# 2326 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6974# 2326 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6977 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dz(
l + 1)), kind=stp)
6979# 2329 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6980#if defined(MFC_OpenACC)
6981# 2329 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6983# 2329 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6984#elif defined(MFC_OpenMP)
6985# 2329 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6987# 2329 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6990 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dz(
l + 1)), kind=stp)
6993# 2333 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6994#if defined(MFC_OpenACC)
6995# 2333 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6997# 2333 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
6998#elif defined(MFC_OpenMP)
6999# 2333 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7001# 2333 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7004 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*(1._wp/dz(
l)), kind=stp)
7006# 2336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7007#if defined(MFC_OpenACC)
7008# 2336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7010# 2336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7011#elif defined(MFC_OpenMP)
7012# 2336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7014# 2336 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7017 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(1)*vel_r(1)*(1._wp/dz(
l)), kind=stp)
7020# 2340 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7021#if defined(MFC_OpenACC)
7022# 2340 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7024# 2340 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7025#elif defined(MFC_OpenMP)
7026# 2340 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7028# 2340 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7031 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dz(
l + 1)), kind=stp)
7033# 2343 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7034#if defined(MFC_OpenACC)
7035# 2343 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7037# 2343 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7038#elif defined(MFC_OpenMP)
7039# 2343 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7041# 2343 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7044 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(2)*(1._wp/dz(
l + 1)), kind=stp)
7047# 2347 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7048#if defined(MFC_OpenACC)
7049# 2347 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7051# 2347 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7052#elif defined(MFC_OpenMP)
7053# 2347 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7055# 2347 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7058 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*(1._wp/dz(
l)), kind=stp)
7060# 2350 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7061#if defined(MFC_OpenACC)
7062# 2350 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7064# 2350 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7065#elif defined(MFC_OpenMP)
7066# 2350 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7068# 2350 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7071 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(2)*vel_l(2)*(1._wp/dz(
l)), kind=stp)
7074# 2354 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7075#if defined(MFC_OpenACC)
7076# 2354 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7078# 2354 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7079#elif defined(MFC_OpenMP)
7080# 2354 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7082# 2354 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7085 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dz(
l + 1)), kind=stp)
7087# 2357 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7088#if defined(MFC_OpenACC)
7089# 2357 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7091# 2357 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7092#elif defined(MFC_OpenMP)
7093# 2357 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7095# 2357 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7098 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(2)*(1._wp/dz(
l + 1)), kind=stp)
7101# 2361 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7102#if defined(MFC_OpenACC)
7103# 2361 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7105# 2361 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7106#elif defined(MFC_OpenMP)
7107# 2361 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7109# 2361 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7112 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*(1._wp/dz(
l)), kind=stp)
7114# 2364 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7115#if defined(MFC_OpenACC)
7116# 2364 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7118# 2364 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7119#elif defined(MFC_OpenMP)
7120# 2364 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7122# 2364 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7125 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(2)*vel_r(2)*(1._wp/dz(
l)), kind=stp)
7128# 2368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7129#if defined(MFC_OpenACC)
7130# 2368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7132# 2368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7133#elif defined(MFC_OpenMP)
7134# 2368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7136# 2368 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7139 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dz(
l + 1)), kind=stp)
7141# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7142#if defined(MFC_OpenACC)
7143# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7145# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7146#elif defined(MFC_OpenMP)
7147# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7149# 2371 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7152 &
l + 1) - real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(3)*(1._wp/dz(
l + 1)), kind=stp)
7155# 2375 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7156#if defined(MFC_OpenACC)
7157# 2375 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7159# 2375 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7160#elif defined(MFC_OpenMP)
7161# 2375 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7163# 2375 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7166 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*(1._wp/dz(
l)), kind=stp)
7168# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7169#if defined(MFC_OpenACC)
7170# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7172# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7173#elif defined(MFC_OpenMP)
7174# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7176# 2378 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7179 &
l) + real(0.5_wp*dt*mu_l*vflux_l_arr(3)*vel_l(3)*(1._wp/dz(
l)), kind=stp)
7182# 2382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7183#if defined(MFC_OpenACC)
7184# 2382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7186# 2382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7187#elif defined(MFC_OpenMP)
7188# 2382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7190# 2382 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7193 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dz(
l + 1)), kind=stp)
7195# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7196#if defined(MFC_OpenACC)
7197# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7199# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7200#elif defined(MFC_OpenMP)
7201# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7203# 2385 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7206 &
l + 1) - real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(3)*(1._wp/dz(
l + 1)), kind=stp)
7209# 2389 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7210#if defined(MFC_OpenACC)
7211# 2389 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7213# 2389 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7214#elif defined(MFC_OpenMP)
7215# 2389 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7217# 2389 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7220 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*(1._wp/dz(
l)), kind=stp)
7222# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7223#if defined(MFC_OpenACC)
7224# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7226# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7227#elif defined(MFC_OpenMP)
7228# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7230# 2392 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7233 &
l) + real(0.5_wp*dt*mu_r*vflux_r_arr(3)*vel_r(3)*(1._wp/dz(
l)), kind=stp)
7236 e_l = 0._wp; e_r = 0._wp
7237 f_l = 0._wp; f_r = 0._wp
7240# 2400 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7241#if defined(MFC_OpenACC)
7242# 2400 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7244# 2400 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7245#elif defined(MFC_OpenMP)
7246# 2400 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7248# 2400 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7256# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7257#if defined(MFC_OpenACC)
7258# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7260# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7261#elif defined(MFC_OpenMP)
7262# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7264# 2406 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7271 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, &
7272 & pres_l, pres_r, cfl)
7275# 2415 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7276#if defined(MFC_OpenACC)
7277# 2415 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7279# 2415 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7280#elif defined(MFC_OpenMP)
7281# 2415 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7283# 2415 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7285 do i = 1, num_fluids
7287# 2417 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7288#if defined(MFC_OpenACC)
7289# 2417 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7291# 2417 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7292#elif defined(MFC_OpenMP)
7293# 2417 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7295# 2417 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7297 rhs_vf(
i)%sf(
j,
k,
l + 1) = rhs_vf(
i)%sf(
j,
k, &
7298 &
l + 1) + real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(3))*(1._wp/dz(
l + 1)) &
7299 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dz(
l + 1)), kind=stp)
7302# 2422 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7303#if defined(MFC_OpenACC)
7304# 2422 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7306# 2422 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7307#elif defined(MFC_OpenMP)
7308# 2422 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7310# 2422 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7312 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
7313 &
l) - real(0.5_wp*dt*(alpha_rho_l(
i)*vel_l(3))*(1._wp/dz(
l)) &
7314 & - 0.5_wp*dt*cfl*(alpha_rho_l(
i))*(1._wp/dz(
l)), kind=stp)
7317 if (num_fluids > 1)
then
7319# 2429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7320#if defined(MFC_OpenACC)
7321# 2429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7323# 2429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7324#elif defined(MFC_OpenMP)
7325# 2429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7327# 2429 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7329 do i = 1, num_fluids - 1
7331# 2431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7332#if defined(MFC_OpenACC)
7333# 2431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7335# 2431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7336#elif defined(MFC_OpenMP)
7337# 2431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7339# 2431 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7342 &
l + 1) + real(0.5_wp*dt*(alpha_l(
i)*vel_l(3))*(1._wp/dz(
l + 1)) &
7343 & - 0.5_wp*dt*cfl*(alpha_l(
i))*(1._wp/dz(
l + 1)), kind=stp)
7346# 2436 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7347#if defined(MFC_OpenACC)
7348# 2436 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7350# 2436 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7351#elif defined(MFC_OpenMP)
7352# 2436 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7354# 2436 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7358 &
l + 1)*vel_l(3)*(1._wp/dz(
l + 1)), kind=stp)
7361# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7362#if defined(MFC_OpenACC)
7363# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7365# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7366#elif defined(MFC_OpenMP)
7367# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7369# 2441 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7372 &
l) - real(0.5_wp*dt*(alpha_l(
i)*vel_l(3))*(1._wp/dz(
l)) - 0.5_wp*dt*cfl*(alpha_l(
i)) &
7373 & *(1._wp/dz(
l)), kind=stp)
7376# 2446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7377#if defined(MFC_OpenACC)
7378# 2446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7380# 2446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7381#elif defined(MFC_OpenMP)
7382# 2446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7384# 2446 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7393# 2453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7394#if defined(MFC_OpenACC)
7395# 2453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7397# 2453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7398#elif defined(MFC_OpenMP)
7399# 2453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7401# 2453 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7404 &
l + 1) + real(0.5_wp*dt*(rho_l*(vel_l(3))**2.0 + pres_l + f_l)*(1._wp/dz(
l + 1)) &
7405 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dz(
l + 1)), kind=stp)
7408# 2458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7409#if defined(MFC_OpenACC)
7410# 2458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7412# 2458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7413#elif defined(MFC_OpenMP)
7414# 2458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7416# 2458 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7419 &
l + 1) + real(0.5_wp*dt*rho_l*vel_l(1)*vel_l(3)*(1._wp/dz(
l + 1)) &
7420 & - 0.5_wp*dt*cfl*(rho_l*vel_l(1))*(1._wp/dz(
l + 1)), kind=stp)
7423# 2463 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7424#if defined(MFC_OpenACC)
7425# 2463 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7427# 2463 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7428#elif defined(MFC_OpenMP)
7429# 2463 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7431# 2463 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7434 &
l + 1) + real(0.5_wp*dt*rho_l*vel_l(2)*vel_l(3)*(1._wp/dz(
l + 1)) &
7435 & - 0.5_wp*dt*cfl*(rho_l*vel_l(2))*(1._wp/dz(
l + 1)), kind=stp)
7438# 2468 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7439#if defined(MFC_OpenACC)
7440# 2468 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7442# 2468 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7443#elif defined(MFC_OpenMP)
7444# 2468 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7446# 2468 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7449 &
l + 1) + real(0.5_wp*dt*(vel_l(3)*(e_l + pres_l + f_l))*(1._wp/dz(
l + 1)) &
7450 & - 0.5_wp*dt*cfl*(e_l)*(1._wp/dz(
l + 1)), kind=stp)
7453# 2473 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7454#if defined(MFC_OpenACC)
7455# 2473 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7457# 2473 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7458#elif defined(MFC_OpenMP)
7459# 2473 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7461# 2473 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7464 &
l) - real(0.5_wp*dt*(rho_l*(vel_l(3))**2.0 + pres_l + f_l)*(1._wp/dz(
l)) &
7465 & - 0.5_wp*dt*cfl*(rho_l*vel_l(3))*(1._wp/dz(
l)), kind=stp)
7468# 2478 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7469#if defined(MFC_OpenACC)
7470# 2478 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7472# 2478 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7473#elif defined(MFC_OpenMP)
7474# 2478 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7476# 2478 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7479 &
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)) &
7480 & *(1._wp/dz(
l)), kind=stp)
7483# 2483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7484#if defined(MFC_OpenACC)
7485# 2483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7487# 2483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7488#elif defined(MFC_OpenMP)
7489# 2483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7491# 2483 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7494 &
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)) &
7495 & *(1._wp/dz(
l)), kind=stp)
7498# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7499#if defined(MFC_OpenACC)
7500# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7502# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7503#elif defined(MFC_OpenMP)
7504# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7506# 2488 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7509 &
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) &
7510 & *(1._wp/dz(
l)), kind=stp)
7513# 2493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7514#if defined(MFC_OpenACC)
7515# 2493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7517# 2493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7518#elif defined(MFC_OpenMP)
7519# 2493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7521# 2493 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7523 do i = 1, num_fluids
7525# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7526#if defined(MFC_OpenACC)
7527# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7529# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7530#elif defined(MFC_OpenMP)
7531# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7533# 2495 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7535 rhs_vf(
i)%sf(
j,
k,
l + 1) = rhs_vf(
i)%sf(
j,
k, &
7536 &
l + 1) + real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(3))*(1._wp/dz(
l + 1)) &
7537 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dz(
l + 1)), kind=stp)
7540# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7541#if defined(MFC_OpenACC)
7542# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7544# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7545#elif defined(MFC_OpenMP)
7546# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7548# 2500 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7550 rhs_vf(
i)%sf(
j,
k,
l) = rhs_vf(
i)%sf(
j,
k, &
7551 &
l) - real(0.5_wp*dt*(alpha_rho_r(
i)*vel_r(3))*(1._wp/dz(
l)) &
7552 & + 0.5_wp*dt*cfl*(alpha_rho_r(
i))*(1._wp/dz(
l)), kind=stp)
7555 if (num_fluids > 1)
then
7557# 2507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7558#if defined(MFC_OpenACC)
7559# 2507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7561# 2507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7562#elif defined(MFC_OpenMP)
7563# 2507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7565# 2507 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7567 do i = 1, num_fluids - 1
7569# 2509 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7570#if defined(MFC_OpenACC)
7571# 2509 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7573# 2509 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7574#elif defined(MFC_OpenMP)
7575# 2509 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7577# 2509 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7580 &
l + 1) + real(0.5_wp*dt*(alpha_r(
i)*vel_r(3))*(1._wp/dz(
l + 1)) &
7581 & + 0.5_wp*dt*cfl*(alpha_r(
i))*(1._wp/dz(
l + 1)), kind=stp)
7584# 2514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7585#if defined(MFC_OpenACC)
7586# 2514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7588# 2514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7589#elif defined(MFC_OpenMP)
7590# 2514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7592# 2514 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7596 &
l + 1)*vel_r(3)*(1._wp/dz(
l + 1)), kind=stp)
7599# 2519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7600#if defined(MFC_OpenACC)
7601# 2519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7603# 2519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7604#elif defined(MFC_OpenMP)
7605# 2519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7607# 2519 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7610 &
l) - real(0.5_wp*dt*(alpha_r(
i)*vel_r(3))*(1._wp/dz(
l)) + 0.5_wp*dt*cfl*(alpha_r(
i)) &
7611 & *(1._wp/dz(
l)), kind=stp)
7614# 2524 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7615#if defined(MFC_OpenACC)
7616# 2524 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7618# 2524 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7619#elif defined(MFC_OpenMP)
7620# 2524 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7622# 2524 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7631# 2531 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7632#if defined(MFC_OpenACC)
7633# 2531 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7635# 2531 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7636#elif defined(MFC_OpenMP)
7637# 2531 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7639# 2531 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7642 &
l + 1) + real(0.5_wp*dt*(rho_r*(vel_r(3))**2.0 + pres_r + f_r)*(1._wp/dz(
l + 1)) &
7643 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dz(
l + 1)), kind=stp)
7646# 2536 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7647#if defined(MFC_OpenACC)
7648# 2536 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7650# 2536 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7651#elif defined(MFC_OpenMP)
7652# 2536 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7654# 2536 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7657 &
l + 1) + real(0.5_wp*dt*rho_r*vel_r(1)*vel_r(3)*(1._wp/dz(
l + 1)) &
7658 & + 0.5_wp*dt*cfl*(rho_r*vel_r(1))*(1._wp/dz(
l + 1)), kind=stp)
7661# 2541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7662#if defined(MFC_OpenACC)
7663# 2541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7665# 2541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7666#elif defined(MFC_OpenMP)
7667# 2541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7669# 2541 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7672 &
l + 1) + real(0.5_wp*dt*rho_r*vel_r(2)*vel_r(3)*(1._wp/dz(
l + 1)) &
7673 & + 0.5_wp*dt*cfl*(rho_r*vel_r(2))*(1._wp/dz(
l + 1)), kind=stp)
7676# 2546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7677#if defined(MFC_OpenACC)
7678# 2546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7680# 2546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7681#elif defined(MFC_OpenMP)
7682# 2546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7684# 2546 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7687 &
l + 1) + real(0.5_wp*dt*(vel_r(3)*(e_r + pres_r + f_r))*(1._wp/dz(
l + 1)) &
7688 & + 0.5_wp*dt*cfl*(e_r)*(1._wp/dz(
l + 1)), kind=stp)
7691# 2551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7692#if defined(MFC_OpenACC)
7693# 2551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7695# 2551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7696#elif defined(MFC_OpenMP)
7697# 2551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7699# 2551 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7702 &
l) - real(0.5_wp*dt*(rho_r*(vel_r(3))**2.0 + pres_r + f_r)*(1._wp/dz(
l)) &
7703 & + 0.5_wp*dt*cfl*(rho_r*vel_r(3))*(1._wp/dz(
l)), kind=stp)
7706# 2556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7707#if defined(MFC_OpenACC)
7708# 2556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7710# 2556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7711#elif defined(MFC_OpenMP)
7712# 2556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7714# 2556 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7717 &
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)) &
7718 & *(1._wp/dz(
l)), kind=stp)
7721# 2561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7722#if defined(MFC_OpenACC)
7723# 2561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7725# 2561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7726#elif defined(MFC_OpenMP)
7727# 2561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7729# 2561 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7732 &
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)) &
7733 & *(1._wp/dz(
l)), kind=stp)
7736# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7737#if defined(MFC_OpenACC)
7738# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7740# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7741#elif defined(MFC_OpenMP)
7742# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7744# 2566 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7747 &
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) &
7748 & *(1._wp/dz(
l)), kind=stp)
7753# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7754#if defined(MFC_OpenACC)
7755# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7757# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7758#elif defined(MFC_OpenMP)
7759# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7761# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7763# 2573 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"
7765# 2575 "/home/runner/work/MFC/MFC/src/simulation/m_igr.fpp"