2101 integer,
intent(in) :: cbc_dir_in, cbc_loc_in
2104 real(wp),
dimension(0:buff_size + 1) :: s_cb
2117 s_cb(i + 1) = s_cb(i) +
ds(i)
2121# 379 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2122 if (cbc_dir_in == 1 .and. recon_type == weno_type)
then
2123 if (weno_order == 1)
then
2129 else if (weno_order == 3)
then
2135 pi_coef_x(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2141 & cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2149 & cbc_loc_in) = ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))*(s_cb(1) - s_cb(3)))/((s_cb(1) &
2150 & - s_cb(4))*(s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2152 & cbc_loc_in) = ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))*((s_cb(1) - s_cb(3))*(s_cb(1) &
2153 & - s_cb(3)) - (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + (s_cb(4) - s_cb(1)))))/((s_cb(0) &
2154 & - s_cb(3))*(s_cb(1) - s_cb(3))*(s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2156 & cbc_loc_in) = (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))*(s_cb(1) - s_cb(3)) + ((s_cb(0) &
2157 & - s_cb(2)) + (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/((s_cb(2) - s_cb(0))*(s_cb(0) &
2158 & - s_cb(3))*(s_cb(0) - s_cb(4)))
2160 & cbc_loc_in) = ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))*(s_cb(2) - s_cb(3)))/((s_cb(2) &
2161 & - s_cb(4))*(s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2163 & cbc_loc_in) = ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))*((s_cb(1) - s_cb(3))*(s_cb(2) &
2164 & - s_cb(3)) + (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + (s_cb(2) - s_cb(4)))))/((s_cb(0) &
2165 & - s_cb(3))*(s_cb(1) - s_cb(3))*(s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2167 & cbc_loc_in) = ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))*(s_cb(2) - s_cb(4)))/((s_cb(0) &
2168 & - s_cb(2))*(s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2171# 379 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2172 if (cbc_dir_in == 2 .and. recon_type == weno_type)
then
2173 if (weno_order == 1)
then
2179 else if (weno_order == 3)
then
2185 pi_coef_y(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2191 & cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2199 & cbc_loc_in) = ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))*(s_cb(1) - s_cb(3)))/((s_cb(1) &
2200 & - s_cb(4))*(s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2202 & cbc_loc_in) = ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))*((s_cb(1) - s_cb(3))*(s_cb(1) &
2203 & - s_cb(3)) - (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + (s_cb(4) - s_cb(1)))))/((s_cb(0) &
2204 & - s_cb(3))*(s_cb(1) - s_cb(3))*(s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2206 & cbc_loc_in) = (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))*(s_cb(1) - s_cb(3)) + ((s_cb(0) &
2207 & - s_cb(2)) + (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/((s_cb(2) - s_cb(0))*(s_cb(0) &
2208 & - s_cb(3))*(s_cb(0) - s_cb(4)))
2210 & cbc_loc_in) = ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))*(s_cb(2) - s_cb(3)))/((s_cb(2) &
2211 & - s_cb(4))*(s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2213 & cbc_loc_in) = ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))*((s_cb(1) - s_cb(3))*(s_cb(2) &
2214 & - s_cb(3)) + (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + (s_cb(2) - s_cb(4)))))/((s_cb(0) &
2215 & - s_cb(3))*(s_cb(1) - s_cb(3))*(s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2217 & cbc_loc_in) = ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))*(s_cb(2) - s_cb(4)))/((s_cb(0) &
2218 & - s_cb(2))*(s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2221# 379 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2222 if (cbc_dir_in == 3 .and. recon_type == weno_type)
then
2223 if (weno_order == 1)
then
2229 else if (weno_order == 3)
then
2235 pi_coef_z(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2241 & cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2249 & cbc_loc_in) = ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))*(s_cb(1) - s_cb(3)))/((s_cb(1) &
2250 & - s_cb(4))*(s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2252 & cbc_loc_in) = ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))*((s_cb(1) - s_cb(3))*(s_cb(1) &
2253 & - s_cb(3)) - (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + (s_cb(4) - s_cb(1)))))/((s_cb(0) &
2254 & - s_cb(3))*(s_cb(1) - s_cb(3))*(s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2256 & cbc_loc_in) = (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))*(s_cb(1) - s_cb(3)) + ((s_cb(0) &
2257 & - s_cb(2)) + (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/((s_cb(2) - s_cb(0))*(s_cb(0) &
2258 & - s_cb(3))*(s_cb(0) - s_cb(4)))
2260 & cbc_loc_in) = ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))*(s_cb(2) - s_cb(3)))/((s_cb(2) &
2261 & - s_cb(4))*(s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2263 & cbc_loc_in) = ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))*((s_cb(1) - s_cb(3))*(s_cb(2) &
2264 & - s_cb(3)) + (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + (s_cb(2) - s_cb(4)))))/((s_cb(0) &
2265 & - s_cb(3))*(s_cb(1) - s_cb(3))*(s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2267 & cbc_loc_in) = ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))*(s_cb(2) - s_cb(4)))/((s_cb(0) &
2268 & - s_cb(2))*(s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2271# 429 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2343 subroutine s_cbc(q_prim_vf, flux_vf, flux_src_vf, cbc_dir_norm, cbc_loc_norm, ix, iy, iz)
2345 type(scalar_field),
dimension(sys_size),
intent(in) :: q_prim_vf
2346 type(scalar_field),
dimension(sys_size),
intent(inout) :: flux_vf, flux_src_vf
2347 integer,
intent(in) :: cbc_dir_norm, cbc_loc_norm
2348 type(int_bounds_info),
intent(in) :: ix, iy, iz
2350 real(wp) :: dpres_dt
2351 real(wp) :: dgamma_dt
2352 real(wp) :: dpi_inf_dt
2354 real(wp) :: dpres_ds
2356# 506 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2357 real(wp),
dimension(sys_size) ::
l
2358# 508 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2359# 517 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2360 real(wp),
dimension(num_fluids) :: alpha_rho, dalpha_rho_ds, mf
2361 real(wp),
dimension(num_vels) :: vel, dvel_ds
2362 real(wp),
dimension(num_fluids) :: adv_local, dadv_ds
2363 real(wp),
dimension(num_fluids) :: dadv_dt
2364 real(wp),
dimension(num_dims) :: dvel_dt
2365 real(wp),
dimension(num_fluids) :: dalpha_rho_dt
2366 real(wp),
dimension(num_species) :: ys, h_k, dys_dt, dys_ds, xs, gamma_i, cp_i
2367# 525 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2368 real(wp),
dimension(2) :: re_cbc
2369 real(wp),
dimension(3) :: lambda
2379 real(wp) :: t, sum_enthalpies
2380 real(wp) :: cv, cp, e_mix, mw, r_gas
2381 real(wp) :: vel_k_sum, vel_dv_dt_sum
2382 integer :: i,
j,
k, r
2390# 546 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2391#if defined(MFC_OpenACC)
2392# 546 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2394# 546 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2395#elif defined(MFC_OpenMP)
2396# 546 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2398# 546 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2405# 553 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2406 if (
cbc_dir == 1 .and. recon_type == weno_type)
then
2408 if (weno_order == 3 .or. dummy)
then
2410 &
is3, idwbuff(2)%beg, idwbuff(3)%beg)
2413# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2415# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2416#if defined(MFC_OpenACC)
2417# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2419# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2420#elif defined(MFC_OpenMP)
2421# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2423# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2425# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2427# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2429# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2440# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2441#if defined(MFC_OpenACC)
2442# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2444# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2445#elif defined(MFC_OpenMP)
2446# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2448# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2450# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2454# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2456# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2457#if defined(MFC_OpenACC)
2458# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2460# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2461#elif defined(MFC_OpenMP)
2462# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2464# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2466# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2468# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2470# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2481# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2482#if defined(MFC_OpenACC)
2483# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2485# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2486#elif defined(MFC_OpenMP)
2487# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2489# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2491# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2496 if (weno_order == 5 .or. dummy)
then
2498 &
is3, idwbuff(2)%beg, idwbuff(3)%beg)
2501# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2503# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2504#if defined(MFC_OpenACC)
2505# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2507# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2508#elif defined(MFC_OpenMP)
2509# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2511# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2513# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2515# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2517# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2533# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2534#if defined(MFC_OpenACC)
2535# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2537# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2538#elif defined(MFC_OpenMP)
2539# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2541# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2543# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2547# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2549# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2550#if defined(MFC_OpenACC)
2551# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2553# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2554#elif defined(MFC_OpenMP)
2555# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2557# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2559# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2561# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2563# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2580# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2581#if defined(MFC_OpenACC)
2582# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2584# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2585#elif defined(MFC_OpenMP)
2586# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2588# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2590# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2596# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2598# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2599#if defined(MFC_OpenACC)
2600# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2602# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2603#elif defined(MFC_OpenMP)
2604# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2606# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2608# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2610# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2612# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2614# 627 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2619# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2620#if defined(MFC_OpenACC)
2621# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2623# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2624#elif defined(MFC_OpenMP)
2625# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2627# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2634# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2635#if defined(MFC_OpenACC)
2636# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2638# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2639#elif defined(MFC_OpenMP)
2640# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2642# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2650# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2651#if defined(MFC_OpenACC)
2652# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2654# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2655#elif defined(MFC_OpenMP)
2656# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2658# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2661 vel_k_sum = vel_k_sum + vel(i)**2._wp
2667# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2668#if defined(MFC_OpenACC)
2669# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2671# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2672#elif defined(MFC_OpenMP)
2673# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2675# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2677 do i = 1, advxe - e_idx
2681 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
2684# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2685#if defined(MFC_OpenACC)
2686# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2688# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2689#elif defined(MFC_OpenMP)
2690# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2692# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2695 mf(i) = alpha_rho(i)/rho
2700# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2701#if defined(MFC_OpenACC)
2702# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2704# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2705#elif defined(MFC_OpenMP)
2706# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2708# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2710 do i = chemxb, chemxe
2714 call get_mixture_molecular_weight(ys, mw)
2715 r_gas = gas_constant/mw
2717 call get_mixture_specific_heat_cp_mass(t, ys, cp)
2718 call get_mixture_energy_mass(t, ys, e_mix)
2719 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
2720 if (chem_params%gamma_method == 1)
then
2722 call get_mole_fractions(mw, ys, xs)
2723 call get_species_specific_heats_r(t, cp_i)
2724 gamma_i = cp_i/(cp_i - 1.0_wp)
2725 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
2726 else if (chem_params%gamma_method == 2)
then
2728 call get_mixture_specific_heat_cv_mass(t, ys, cv)
2729 gamma = 1.0_wp/(cp/cv - 1.0_wp)
2732 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
2738 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
2743# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2744#if defined(MFC_OpenACC)
2745# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2747# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2748#elif defined(MFC_OpenMP)
2749# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2751# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2754 dalpha_rho_ds(i) = 0._wp
2758# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2759#if defined(MFC_OpenACC)
2760# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2762# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2763#elif defined(MFC_OpenMP)
2764# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2766# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2774# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2775#if defined(MFC_OpenACC)
2776# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2778# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2779#elif defined(MFC_OpenMP)
2780# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2782# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2784 do i = 1, advxe - e_idx
2790# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2791#if defined(MFC_OpenACC)
2792# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2794# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2795#elif defined(MFC_OpenMP)
2796# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2798# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2800 do i = 1, num_species
2806# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2807#if defined(MFC_OpenACC)
2808# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2810# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2811#elif defined(MFC_OpenMP)
2812# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2814# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2818# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2819#if defined(MFC_OpenACC)
2820# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2822# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2823#elif defined(MFC_OpenMP)
2824# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2826# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2832# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2833#if defined(MFC_OpenACC)
2834# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2836# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2837#elif defined(MFC_OpenMP)
2838# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2840# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2848# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2849#if defined(MFC_OpenACC)
2850# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2852# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2853#elif defined(MFC_OpenMP)
2854# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2856# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2858 do i = 1, advxe - e_idx
2864# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2865#if defined(MFC_OpenACC)
2866# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2868# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2869#elif defined(MFC_OpenMP)
2870# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2872# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2874 do i = 1, num_species
2882 lambda(1) = vel(dir_idx(1)) - c
2883 lambda(2) = vel(dir_idx(1))
2884 lambda(3) = vel(dir_idx(1)) + c
2886 ma = vel(dir_idx(1))/c
2888 if ((
cbc_loc == -1 .and.
bcxb == bc_char_slip_wall) &
2889 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_slip_wall))
then
2890 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
2891 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_buffer) &
2892 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_buffer))
then
2893 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
2894 & dvel_ds, dadv_ds, dys_ds)
2895 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_inflow) &
2896 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_inflow))
then
2897 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
2899 if (bc_x%grcbc_in)
then
2901# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2902#if defined(MFC_OpenACC)
2903# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2905# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2906#elif defined(MFC_OpenMP)
2907# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2909# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2912 l(i) = c**3._wp*ma*(alpha_rho(i - 1) -
alpha_rho_in(i - 1, &
2916 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(1, dir_idx(2)))/
del_in(1)
2918 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(1, dir_idx(3)))/
del_in(1)
2922# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2923#if defined(MFC_OpenACC)
2924# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2926# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2927#elif defined(MFC_OpenMP)
2928# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2930# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2932 do i = e_idx, advxe - 1
2933 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, &
2936 l(advxe) = rho*c**2._wp*(1._wp + ma)*(vel(dir_idx(1)) +
vel_in(1, dir_idx(1))*sign(1, &
2939 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_outflow) &
2940 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_outflow))
then
2941 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
2942 & dvel_ds, dadv_ds, dys_ds)
2944 if (bc_x%grcbc_out)
then
2948 if (bc_x%grcbc_vel_out)
then
2949 l(advxe) =
l(advxe) + rho*c**2._wp*(1._wp - ma)*(vel(dir_idx(1)) +
vel_out(1, &
2953 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_ff_sub_outflow) &
2954 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_ff_sub_outflow))
then
2955 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, &
2957 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_cp_sub_outflow) &
2958 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_cp_sub_outflow))
then
2959 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
2961 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_sup_inflow) &
2962 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_sup_inflow))
then
2963 call s_compute_supersonic_inflow_l(
l)
2964 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_sup_outflow) &
2965 & .or. (
cbc_loc == 1 .and.
bcxe == bc_char_sup_outflow))
then
2966 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, &
2972 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1))/y_cc(n)
2974 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
2978# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2979#if defined(MFC_OpenACC)
2980# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2982# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2983#elif defined(MFC_OpenMP)
2984# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2986# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2989 dalpha_rho_dt(i) = -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
2993# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2994#if defined(MFC_OpenACC)
2995# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2997# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2998#elif defined(MFC_OpenMP)
2999# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3001# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3004 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))*(
l(1) -
l(advxe))/(2._wp*rho*c) + (dir_flg(dir_idx(i)) &
3005 & - 1._wp)*
l(momxb + i - 1)
3008 vel_dv_dt_sum = 0._wp
3010# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3011#if defined(MFC_OpenACC)
3012# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3014# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3015#elif defined(MFC_OpenMP)
3016# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3018# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3021 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
3026# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3027#if defined(MFC_OpenACC)
3028# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3030# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3031#elif defined(MFC_OpenMP)
3032# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3034# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3036 do i = 1, num_species
3037 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
3044# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3045#if defined(MFC_OpenACC)
3046# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3048# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3049#elif defined(MFC_OpenMP)
3050# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3052# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3054 do i = 1, advxe - e_idx
3055 dadv_dt(i) = -
l(momxe + i)
3059# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3060#if defined(MFC_OpenACC)
3061# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3063# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3064#elif defined(MFC_OpenMP)
3065# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3067# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3069 do i = 1, advxe - e_idx
3070 dadv_dt(i) = -
l(momxe + i)
3074 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
3076 if (model_eqns == 1)
then
3077 drho_dt = dalpha_rho_dt(1)
3078 dgamma_dt = dadv_dt(1)
3079# 862 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3080 dpi_inf_dt = dadv_dt(2)
3081# 864 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3084# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3085#if defined(MFC_OpenACC)
3086# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3088# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3089#elif defined(MFC_OpenMP)
3090# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3092# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3094 do i = 1, num_fluids
3095 drho_dt = drho_dt + dalpha_rho_dt(i)
3096 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
3097 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
3098 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
3104# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3105#if defined(MFC_OpenACC)
3106# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3108# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3109#elif defined(MFC_OpenMP)
3110# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3112# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3119# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3120#if defined(MFC_OpenACC)
3121# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3123# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3124#elif defined(MFC_OpenMP)
3125# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3127# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3131 & i) +
ds(0)*(vel(i - contxe)*drho_dt + rho*dvel_dt(i - contxe))
3137 call get_species_enthalpies_rt(t, h_k)
3138 sum_enthalpies = 0._wp
3140# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3141#if defined(MFC_OpenACC)
3142# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3144# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3145#elif defined(MFC_OpenMP)
3146# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3148# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3150 do i = 1, num_species
3151# 898 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3152 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
3153 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
3154# 901 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3157 & e_idx) +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) &
3160# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3161#if defined(MFC_OpenACC)
3162# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3164# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3165#elif defined(MFC_OpenMP)
3166# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3168# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3170 do i = 1, num_species
3172 & chemxb + i - 1) +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
3176 & e_idx) +
ds(0)*(pres*dgamma_dt + gamma*dpres_dt + dpi_inf_dt + dqv_dt &
3177 & + rho*vel_dv_dt_sum + 5.e-1_wp*drho_dt*vel_k_sum)
3180 if (riemann_solver == 1)
then
3182# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3183#if defined(MFC_OpenACC)
3184# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3186# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3187#elif defined(MFC_OpenMP)
3188# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3190# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3197# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3198#if defined(MFC_OpenACC)
3199# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3201# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3202#elif defined(MFC_OpenMP)
3203# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3205# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3208 flux_src_rsx_vf_l(-1,
k, r, i) = 1._wp/max(abs(vel(dir_idx(1))), sgm_eps)*sign(1._wp, &
3211 & i) +
ds(0)*dadv_dt(i - e_idx))
3215# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3216#if defined(MFC_OpenACC)
3217# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3219# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3220#elif defined(MFC_OpenMP)
3221# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3223# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3230# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3231#if defined(MFC_OpenACC)
3232# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3234# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3235#elif defined(MFC_OpenMP)
3236# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3238# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3248# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3249#if defined(MFC_OpenACC)
3250# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3252# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3253#elif defined(MFC_OpenMP)
3254# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3256# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3258# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3261# 553 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3262 if (
cbc_dir == 2 .and. recon_type == weno_type)
then
3264 if (weno_order == 3 .or. dummy)
then
3266 &
is3, idwbuff(2)%beg, idwbuff(3)%beg)
3269# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3271# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3272#if defined(MFC_OpenACC)
3273# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3275# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3276#elif defined(MFC_OpenMP)
3277# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3279# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3281# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3283# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3285# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3296# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3297#if defined(MFC_OpenACC)
3298# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3300# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3301#elif defined(MFC_OpenMP)
3302# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3304# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3306# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3310# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3312# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3313#if defined(MFC_OpenACC)
3314# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3316# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3317#elif defined(MFC_OpenMP)
3318# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3320# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3322# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3324# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3326# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3337# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3338#if defined(MFC_OpenACC)
3339# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3341# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3342#elif defined(MFC_OpenMP)
3343# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3345# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3347# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3352 if (weno_order == 5 .or. dummy)
then
3354 &
is3, idwbuff(2)%beg, idwbuff(3)%beg)
3357# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3359# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3360#if defined(MFC_OpenACC)
3361# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3363# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3364#elif defined(MFC_OpenMP)
3365# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3367# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3369# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3371# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3373# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3389# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3390#if defined(MFC_OpenACC)
3391# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3393# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3394#elif defined(MFC_OpenMP)
3395# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3397# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3399# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3403# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3405# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3406#if defined(MFC_OpenACC)
3407# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3409# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3410#elif defined(MFC_OpenMP)
3411# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3413# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3415# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3417# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3419# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3436# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3437#if defined(MFC_OpenACC)
3438# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3440# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3441#elif defined(MFC_OpenMP)
3442# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3444# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3446# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3452# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3454# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3455#if defined(MFC_OpenACC)
3456# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3458# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3459#elif defined(MFC_OpenMP)
3460# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3462# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3464# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3466# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3468# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3470# 627 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3475# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3476#if defined(MFC_OpenACC)
3477# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3479# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3480#elif defined(MFC_OpenMP)
3481# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3483# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3490# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3491#if defined(MFC_OpenACC)
3492# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3494# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3495#elif defined(MFC_OpenMP)
3496# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3498# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3506# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3507#if defined(MFC_OpenACC)
3508# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3510# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3511#elif defined(MFC_OpenMP)
3512# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3514# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3517 vel_k_sum = vel_k_sum + vel(i)**2._wp
3523# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3524#if defined(MFC_OpenACC)
3525# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3527# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3528#elif defined(MFC_OpenMP)
3529# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3531# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3533 do i = 1, advxe - e_idx
3537 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
3540# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3541#if defined(MFC_OpenACC)
3542# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3544# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3545#elif defined(MFC_OpenMP)
3546# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3548# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3551 mf(i) = alpha_rho(i)/rho
3556# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3557#if defined(MFC_OpenACC)
3558# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3560# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3561#elif defined(MFC_OpenMP)
3562# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3564# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3566 do i = chemxb, chemxe
3570 call get_mixture_molecular_weight(ys, mw)
3571 r_gas = gas_constant/mw
3573 call get_mixture_specific_heat_cp_mass(t, ys, cp)
3574 call get_mixture_energy_mass(t, ys, e_mix)
3575 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
3576 if (chem_params%gamma_method == 1)
then
3578 call get_mole_fractions(mw, ys, xs)
3579 call get_species_specific_heats_r(t, cp_i)
3580 gamma_i = cp_i/(cp_i - 1.0_wp)
3581 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
3582 else if (chem_params%gamma_method == 2)
then
3584 call get_mixture_specific_heat_cv_mass(t, ys, cv)
3585 gamma = 1.0_wp/(cp/cv - 1.0_wp)
3588 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
3594 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
3599# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3600#if defined(MFC_OpenACC)
3601# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3603# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3604#elif defined(MFC_OpenMP)
3605# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3607# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3610 dalpha_rho_ds(i) = 0._wp
3614# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3615#if defined(MFC_OpenACC)
3616# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3618# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3619#elif defined(MFC_OpenMP)
3620# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3622# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3630# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3631#if defined(MFC_OpenACC)
3632# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3634# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3635#elif defined(MFC_OpenMP)
3636# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3638# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3640 do i = 1, advxe - e_idx
3646# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3647#if defined(MFC_OpenACC)
3648# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3650# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3651#elif defined(MFC_OpenMP)
3652# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3654# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3656 do i = 1, num_species
3662# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3663#if defined(MFC_OpenACC)
3664# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3666# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3667#elif defined(MFC_OpenMP)
3668# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3670# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3674# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3675#if defined(MFC_OpenACC)
3676# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3678# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3679#elif defined(MFC_OpenMP)
3680# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3682# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3688# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3689#if defined(MFC_OpenACC)
3690# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3692# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3693#elif defined(MFC_OpenMP)
3694# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3696# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3704# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3705#if defined(MFC_OpenACC)
3706# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3708# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3709#elif defined(MFC_OpenMP)
3710# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3712# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3714 do i = 1, advxe - e_idx
3720# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3721#if defined(MFC_OpenACC)
3722# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3724# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3725#elif defined(MFC_OpenMP)
3726# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3728# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3730 do i = 1, num_species
3738 lambda(1) = vel(dir_idx(1)) - c
3739 lambda(2) = vel(dir_idx(1))
3740 lambda(3) = vel(dir_idx(1)) + c
3742 ma = vel(dir_idx(1))/c
3744 if ((
cbc_loc == -1 .and.
bcyb == bc_char_slip_wall) &
3745 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_slip_wall))
then
3746 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
3747 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_buffer) &
3748 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_buffer))
then
3749 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
3750 & dvel_ds, dadv_ds, dys_ds)
3751 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_inflow) &
3752 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_inflow))
then
3753 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
3755 if (bc_y%grcbc_in)
then
3757# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3758#if defined(MFC_OpenACC)
3759# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3761# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3762#elif defined(MFC_OpenMP)
3763# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3765# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3768 l(i) = c**3._wp*ma*(alpha_rho(i - 1) -
alpha_rho_in(i - 1, &
3772 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(2, dir_idx(2)))/
del_in(2)
3774 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(2, dir_idx(3)))/
del_in(2)
3778# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3779#if defined(MFC_OpenACC)
3780# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3782# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3783#elif defined(MFC_OpenMP)
3784# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3786# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3788 do i = e_idx, advxe - 1
3789 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, &
3792 l(advxe) = rho*c**2._wp*(1._wp + ma)*(vel(dir_idx(1)) +
vel_in(2, dir_idx(1))*sign(1, &
3795 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_outflow) &
3796 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_outflow))
then
3797 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
3798 & dvel_ds, dadv_ds, dys_ds)
3800 if (bc_y%grcbc_out)
then
3804 if (bc_y%grcbc_vel_out)
then
3805 l(advxe) =
l(advxe) + rho*c**2._wp*(1._wp - ma)*(vel(dir_idx(1)) +
vel_out(2, &
3809 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_ff_sub_outflow) &
3810 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_ff_sub_outflow))
then
3811 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, &
3813 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_cp_sub_outflow) &
3814 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_cp_sub_outflow))
then
3815 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
3817 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_sup_inflow) &
3818 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_sup_inflow))
then
3819 call s_compute_supersonic_inflow_l(
l)
3820 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_sup_outflow) &
3821 & .or. (
cbc_loc == 1 .and.
bcye == bc_char_sup_outflow))
then
3822 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, &
3828 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1))/y_cc(n)
3830 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
3834# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3835#if defined(MFC_OpenACC)
3836# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3838# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3839#elif defined(MFC_OpenMP)
3840# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3842# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3845 dalpha_rho_dt(i) = -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
3849# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3850#if defined(MFC_OpenACC)
3851# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3853# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3854#elif defined(MFC_OpenMP)
3855# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3857# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3860 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))*(
l(1) -
l(advxe))/(2._wp*rho*c) + (dir_flg(dir_idx(i)) &
3861 & - 1._wp)*
l(momxb + i - 1)
3864 vel_dv_dt_sum = 0._wp
3866# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3867#if defined(MFC_OpenACC)
3868# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3870# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3871#elif defined(MFC_OpenMP)
3872# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3874# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3877 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
3882# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3883#if defined(MFC_OpenACC)
3884# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3886# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3887#elif defined(MFC_OpenMP)
3888# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3890# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3892 do i = 1, num_species
3893 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
3900# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3901#if defined(MFC_OpenACC)
3902# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3904# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3905#elif defined(MFC_OpenMP)
3906# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3908# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3910 do i = 1, advxe - e_idx
3911 dadv_dt(i) = -
l(momxe + i)
3915# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3916#if defined(MFC_OpenACC)
3917# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3919# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3920#elif defined(MFC_OpenMP)
3921# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3923# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3925 do i = 1, advxe - e_idx
3926 dadv_dt(i) = -
l(momxe + i)
3930 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
3932 if (model_eqns == 1)
then
3933 drho_dt = dalpha_rho_dt(1)
3934 dgamma_dt = dadv_dt(1)
3935# 862 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3936 dpi_inf_dt = dadv_dt(2)
3937# 864 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3940# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3941#if defined(MFC_OpenACC)
3942# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3944# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3945#elif defined(MFC_OpenMP)
3946# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3948# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3950 do i = 1, num_fluids
3951 drho_dt = drho_dt + dalpha_rho_dt(i)
3952 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
3953 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
3954 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
3960# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3961#if defined(MFC_OpenACC)
3962# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3964# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3965#elif defined(MFC_OpenMP)
3966# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3968# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3975# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3976#if defined(MFC_OpenACC)
3977# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3979# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3980#elif defined(MFC_OpenMP)
3981# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3983# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3987 & i) +
ds(0)*(vel(i - contxe)*drho_dt + rho*dvel_dt(i - contxe))
3993 call get_species_enthalpies_rt(t, h_k)
3994 sum_enthalpies = 0._wp
3996# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3997#if defined(MFC_OpenACC)
3998# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4000# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4001#elif defined(MFC_OpenMP)
4002# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4004# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4006 do i = 1, num_species
4007# 898 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4008 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
4009 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
4010# 901 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4013 & e_idx) +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) &
4016# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4017#if defined(MFC_OpenACC)
4018# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4020# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4021#elif defined(MFC_OpenMP)
4022# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4024# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4026 do i = 1, num_species
4028 & chemxb + i - 1) +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
4032 & e_idx) +
ds(0)*(pres*dgamma_dt + gamma*dpres_dt + dpi_inf_dt + dqv_dt &
4033 & + rho*vel_dv_dt_sum + 5.e-1_wp*drho_dt*vel_k_sum)
4036 if (riemann_solver == 1)
then
4038# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4039#if defined(MFC_OpenACC)
4040# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4042# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4043#elif defined(MFC_OpenMP)
4044# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4046# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4053# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4054#if defined(MFC_OpenACC)
4055# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4057# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4058#elif defined(MFC_OpenMP)
4059# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4061# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4064 flux_src_rsy_vf_l(-1,
k, r, i) = 1._wp/max(abs(vel(dir_idx(1))), sgm_eps)*sign(1._wp, &
4067 & i) +
ds(0)*dadv_dt(i - e_idx))
4071# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4072#if defined(MFC_OpenACC)
4073# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4075# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4076#elif defined(MFC_OpenMP)
4077# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4079# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4086# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4087#if defined(MFC_OpenACC)
4088# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4090# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4091#elif defined(MFC_OpenMP)
4092# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4094# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4104# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4105#if defined(MFC_OpenACC)
4106# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4108# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4109#elif defined(MFC_OpenMP)
4110# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4112# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4114# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4117# 553 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4118 if (
cbc_dir == 3 .and. recon_type == weno_type)
then
4120 if (weno_order == 3 .or. dummy)
then
4122 &
is3, idwbuff(2)%beg, idwbuff(3)%beg)
4125# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4127# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4128#if defined(MFC_OpenACC)
4129# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4131# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4132#elif defined(MFC_OpenMP)
4133# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4135# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4137# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4139# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4141# 559 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4152# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4153#if defined(MFC_OpenACC)
4154# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4156# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4157#elif defined(MFC_OpenMP)
4158# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4160# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4162# 568 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4166# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4168# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4169#if defined(MFC_OpenACC)
4170# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4172# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4173#elif defined(MFC_OpenMP)
4174# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4176# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4178# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4180# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4182# 570 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4193# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4194#if defined(MFC_OpenACC)
4195# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4197# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4198#elif defined(MFC_OpenMP)
4199# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4201# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4203# 579 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4208 if (weno_order == 5 .or. dummy)
then
4210 &
is3, idwbuff(2)%beg, idwbuff(3)%beg)
4213# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4215# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4216#if defined(MFC_OpenACC)
4217# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4219# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4220#elif defined(MFC_OpenMP)
4221# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4223# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4225# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4227# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4229# 587 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4245# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4246#if defined(MFC_OpenACC)
4247# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4249# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4250#elif defined(MFC_OpenMP)
4251# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4253# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4255# 601 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4259# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4261# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4262#if defined(MFC_OpenACC)
4263# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4265# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4266#elif defined(MFC_OpenMP)
4267# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4269# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4271# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4273# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4275# 603 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4292# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4293#if defined(MFC_OpenACC)
4294# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4296# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4297#elif defined(MFC_OpenMP)
4298# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4300# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4302# 618 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4308# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4310# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4311#if defined(MFC_OpenACC)
4312# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4314# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4315#elif defined(MFC_OpenMP)
4316# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4318# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4320# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4322# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4324# 622 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4326# 627 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4331# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4332#if defined(MFC_OpenACC)
4333# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4335# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4336#elif defined(MFC_OpenMP)
4337# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4339# 630 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4346# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4347#if defined(MFC_OpenACC)
4348# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4350# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4351#elif defined(MFC_OpenMP)
4352# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4354# 635 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4362# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4363#if defined(MFC_OpenACC)
4364# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4366# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4367#elif defined(MFC_OpenMP)
4368# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4370# 641 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4373 vel_k_sum = vel_k_sum + vel(i)**2._wp
4379# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4380#if defined(MFC_OpenACC)
4381# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4383# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4384#elif defined(MFC_OpenMP)
4385# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4387# 648 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4389 do i = 1, advxe - e_idx
4393 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
4396# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4397#if defined(MFC_OpenACC)
4398# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4400# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4401#elif defined(MFC_OpenMP)
4402# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4404# 655 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4407 mf(i) = alpha_rho(i)/rho
4412# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4413#if defined(MFC_OpenACC)
4414# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4416# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4417#elif defined(MFC_OpenMP)
4418# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4420# 661 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4422 do i = chemxb, chemxe
4426 call get_mixture_molecular_weight(ys, mw)
4427 r_gas = gas_constant/mw
4429 call get_mixture_specific_heat_cp_mass(t, ys, cp)
4430 call get_mixture_energy_mass(t, ys, e_mix)
4431 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
4432 if (chem_params%gamma_method == 1)
then
4434 call get_mole_fractions(mw, ys, xs)
4435 call get_species_specific_heats_r(t, cp_i)
4436 gamma_i = cp_i/(cp_i - 1.0_wp)
4437 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
4438 else if (chem_params%gamma_method == 2)
then
4440 call get_mixture_specific_heat_cv_mass(t, ys, cv)
4441 gamma = 1.0_wp/(cp/cv - 1.0_wp)
4444 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
4450 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
4455# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4456#if defined(MFC_OpenACC)
4457# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4459# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4460#elif defined(MFC_OpenMP)
4461# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4463# 694 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4466 dalpha_rho_ds(i) = 0._wp
4470# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4471#if defined(MFC_OpenACC)
4472# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4474# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4475#elif defined(MFC_OpenMP)
4476# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4478# 699 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4486# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4487#if defined(MFC_OpenACC)
4488# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4490# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4491#elif defined(MFC_OpenMP)
4492# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4494# 705 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4496 do i = 1, advxe - e_idx
4502# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4503#if defined(MFC_OpenACC)
4504# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4506# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4507#elif defined(MFC_OpenMP)
4508# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4510# 711 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4512 do i = 1, num_species
4518# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4519#if defined(MFC_OpenACC)
4520# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4522# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4523#elif defined(MFC_OpenMP)
4524# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4526# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4530# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4531#if defined(MFC_OpenACC)
4532# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4534# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4535#elif defined(MFC_OpenMP)
4536# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4538# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4544# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4545#if defined(MFC_OpenACC)
4546# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4548# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4549#elif defined(MFC_OpenMP)
4550# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4552# 723 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4560# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4561#if defined(MFC_OpenACC)
4562# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4564# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4565#elif defined(MFC_OpenMP)
4566# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4568# 729 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4570 do i = 1, advxe - e_idx
4576# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4577#if defined(MFC_OpenACC)
4578# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4580# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4581#elif defined(MFC_OpenMP)
4582# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4584# 735 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4586 do i = 1, num_species
4594 lambda(1) = vel(dir_idx(1)) - c
4595 lambda(2) = vel(dir_idx(1))
4596 lambda(3) = vel(dir_idx(1)) + c
4598 ma = vel(dir_idx(1))/c
4600 if ((
cbc_loc == -1 .and.
bczb == bc_char_slip_wall) &
4601 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_slip_wall))
then
4602 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4603 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_buffer) &
4604 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_buffer))
then
4605 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
4606 & dvel_ds, dadv_ds, dys_ds)
4607 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_inflow) &
4608 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_inflow))
then
4609 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4611 if (bc_z%grcbc_in)
then
4613# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4614#if defined(MFC_OpenACC)
4615# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4617# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4618#elif defined(MFC_OpenMP)
4619# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4621# 762 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4624 l(i) = c**3._wp*ma*(alpha_rho(i - 1) -
alpha_rho_in(i - 1, &
4628 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(3, dir_idx(2)))/
del_in(3)
4630 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(3, dir_idx(3)))/
del_in(3)
4634# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4635#if defined(MFC_OpenACC)
4636# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4638# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4639#elif defined(MFC_OpenMP)
4640# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4642# 773 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4644 do i = e_idx, advxe - 1
4645 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, &
4648 l(advxe) = rho*c**2._wp*(1._wp + ma)*(vel(dir_idx(1)) +
vel_in(3, dir_idx(1))*sign(1, &
4651 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_outflow) &
4652 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_outflow))
then
4653 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
4654 & dvel_ds, dadv_ds, dys_ds)
4656 if (bc_z%grcbc_out)
then
4660 if (bc_z%grcbc_vel_out)
then
4661 l(advxe) =
l(advxe) + rho*c**2._wp*(1._wp - ma)*(vel(dir_idx(1)) +
vel_out(3, &
4665 else if ((
cbc_loc == -1 .and.
bczb == bc_char_ff_sub_outflow) &
4666 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_ff_sub_outflow))
then
4667 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, &
4669 else if ((
cbc_loc == -1 .and.
bczb == bc_char_cp_sub_outflow) &
4670 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_cp_sub_outflow))
then
4671 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, &
4673 else if ((
cbc_loc == -1 .and.
bczb == bc_char_sup_inflow) &
4674 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_sup_inflow))
then
4675 call s_compute_supersonic_inflow_l(
l)
4676 else if ((
cbc_loc == -1 .and.
bczb == bc_char_sup_outflow) &
4677 & .or. (
cbc_loc == 1 .and.
bcze == bc_char_sup_outflow))
then
4678 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, &
4684 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1))/y_cc(n)
4686 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
4690# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4691#if defined(MFC_OpenACC)
4692# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4694# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4695#elif defined(MFC_OpenMP)
4696# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4698# 819 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4701 dalpha_rho_dt(i) = -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
4705# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4706#if defined(MFC_OpenACC)
4707# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4709# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4710#elif defined(MFC_OpenMP)
4711# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4713# 824 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4716 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))*(
l(1) -
l(advxe))/(2._wp*rho*c) + (dir_flg(dir_idx(i)) &
4717 & - 1._wp)*
l(momxb + i - 1)
4720 vel_dv_dt_sum = 0._wp
4722# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4723#if defined(MFC_OpenACC)
4724# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4726# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4727#elif defined(MFC_OpenMP)
4728# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4730# 831 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4733 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
4738# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4739#if defined(MFC_OpenACC)
4740# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4742# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4743#elif defined(MFC_OpenMP)
4744# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4746# 837 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4748 do i = 1, num_species
4749 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
4756# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4757#if defined(MFC_OpenACC)
4758# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4760# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4761#elif defined(MFC_OpenMP)
4762# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4764# 845 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4766 do i = 1, advxe - e_idx
4767 dadv_dt(i) = -
l(momxe + i)
4771# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4772#if defined(MFC_OpenACC)
4773# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4775# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4776#elif defined(MFC_OpenMP)
4777# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4779# 850 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4781 do i = 1, advxe - e_idx
4782 dadv_dt(i) = -
l(momxe + i)
4786 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
4788 if (model_eqns == 1)
then
4789 drho_dt = dalpha_rho_dt(1)
4790 dgamma_dt = dadv_dt(1)
4791# 862 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4792 dpi_inf_dt = dadv_dt(2)
4793# 864 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4796# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4797#if defined(MFC_OpenACC)
4798# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4800# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4801#elif defined(MFC_OpenMP)
4802# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4804# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4806 do i = 1, num_fluids
4807 drho_dt = drho_dt + dalpha_rho_dt(i)
4808 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
4809 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
4810 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
4816# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4817#if defined(MFC_OpenACC)
4818# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4820# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4821#elif defined(MFC_OpenMP)
4822# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4824# 875 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4831# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4832#if defined(MFC_OpenACC)
4833# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4835# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4836#elif defined(MFC_OpenMP)
4837# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4839# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4843 & i) +
ds(0)*(vel(i - contxe)*drho_dt + rho*dvel_dt(i - contxe))
4849 call get_species_enthalpies_rt(t, h_k)
4850 sum_enthalpies = 0._wp
4852# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4853#if defined(MFC_OpenACC)
4854# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4856# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4857#elif defined(MFC_OpenMP)
4858# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4860# 891 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4862 do i = 1, num_species
4863# 898 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4864 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
4865 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
4866# 901 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4869 & e_idx) +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) &
4872# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4873#if defined(MFC_OpenACC)
4874# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4876# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4877#elif defined(MFC_OpenMP)
4878# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4880# 905 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4882 do i = 1, num_species
4884 & chemxb + i - 1) +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
4888 & e_idx) +
ds(0)*(pres*dgamma_dt + gamma*dpres_dt + dpi_inf_dt + dqv_dt &
4889 & + rho*vel_dv_dt_sum + 5.e-1_wp*drho_dt*vel_k_sum)
4892 if (riemann_solver == 1)
then
4894# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4895#if defined(MFC_OpenACC)
4896# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4898# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4899#elif defined(MFC_OpenMP)
4900# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4902# 917 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4909# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4910#if defined(MFC_OpenACC)
4911# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4913# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4914#elif defined(MFC_OpenMP)
4915# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4917# 922 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4920 flux_src_rsz_vf_l(-1,
k, r, i) = 1._wp/max(abs(vel(dir_idx(1))), sgm_eps)*sign(1._wp, &
4923 & i) +
ds(0)*dadv_dt(i - e_idx))
4927# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4928#if defined(MFC_OpenACC)
4929# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4931# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4932#elif defined(MFC_OpenMP)
4933# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4935# 930 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4942# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4943#if defined(MFC_OpenACC)
4944# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4946# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4947#elif defined(MFC_OpenMP)
4948# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4950# 935 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4960# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4961#if defined(MFC_OpenACC)
4962# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4964# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4965#elif defined(MFC_OpenMP)
4966# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4968# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4970# 943 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4973# 946 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4986 type(scalar_field),
dimension(sys_size),
intent(in) :: q_prim_vf
4987 type(scalar_field),
dimension(sys_size),
intent(in) :: flux_vf, flux_src_vf
4988 type(int_bounds_info),
intent(in) :: ix, iy, iz
4989 integer :: i, j, k, r
4998 dir_idx = (/1, 2, 3/); dir_flg = (/1._wp, 0._wp, 0._wp/)
5001 dir_idx = (/2, 1, 3/); dir_flg = (/0._wp, 1._wp, 0._wp/)
5004 dir_idx = (/3, 1, 2/); dir_flg = (/0._wp, 0._wp, 1._wp/)
5009# 980 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5010#if defined(MFC_OpenACC)
5011# 980 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5013# 980 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5014#elif defined(MFC_OpenMP)
5015# 980 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5017# 980 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5020# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5021#if defined(MFC_OpenACC)
5022# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5024# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5025#elif defined(MFC_OpenMP)
5026# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5028# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5034# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5036# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5037#if defined(MFC_OpenACC)
5038# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5040# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5041#elif defined(MFC_OpenMP)
5042# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5044# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5046# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5048# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5050# 985 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5062# 995 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5063#if defined(MFC_OpenACC)
5064# 995 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5066# 995 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5067#elif defined(MFC_OpenMP)
5068# 995 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5070# 995 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5072# 995 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5076# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5078# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5079#if defined(MFC_OpenACC)
5080# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5082# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5083#elif defined(MFC_OpenMP)
5084# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5086# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5088# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5090# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5092# 997 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5102# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5103#if defined(MFC_OpenACC)
5104# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5106# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5107#elif defined(MFC_OpenMP)
5108# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5110# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5112# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5116# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5118# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5119#if defined(MFC_OpenACC)
5120# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5122# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5123#elif defined(MFC_OpenMP)
5124# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5126# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5128# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5130# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5132# 1007 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5137 do j = -1, buff_size
5144# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5145#if defined(MFC_OpenACC)
5146# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5148# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5149#elif defined(MFC_OpenMP)
5150# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5152# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5154# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5158# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5160# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5161#if defined(MFC_OpenACC)
5162# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5164# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5165#elif defined(MFC_OpenMP)
5166# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5168# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5170# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5172# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5174# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5178 do j = -1, buff_size
5179 flux_rsx_vf_l(j, k, r, momxb) = flux_vf(momxb)%sf(
dj*((m - 1) - 2*j) + j, k, r)
5184# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5185#if defined(MFC_OpenACC)
5186# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5188# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5189#elif defined(MFC_OpenMP)
5190# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5192# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5194# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5197 if (riemann_solver == 1)
then
5199# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5201# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5202#if defined(MFC_OpenACC)
5203# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5205# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5206#elif defined(MFC_OpenMP)
5207# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5209# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5211# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5213# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5215# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5220 do j = -1, buff_size
5227# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5228#if defined(MFC_OpenACC)
5229# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5231# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5232#elif defined(MFC_OpenMP)
5233# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5235# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5237# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5241# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5243# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5244#if defined(MFC_OpenACC)
5245# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5247# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5248#elif defined(MFC_OpenMP)
5249# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5251# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5253# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5255# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5257# 1042 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5261 do j = -1, buff_size
5262 flux_src_rsx_vf_l(j, k, r, advxb) = flux_src_vf(advxb)%sf(
dj*((m - 1) - 2*j) + j, k, r)*sign(1._wp, &
5268# 1051 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5269#if defined(MFC_OpenACC)
5270# 1051 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5272# 1051 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5273#elif defined(MFC_OpenMP)
5274# 1051 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5276# 1051 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5278# 1051 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5287# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5289# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5290#if defined(MFC_OpenACC)
5291# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5293# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5294#elif defined(MFC_OpenMP)
5295# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5297# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5299# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5301# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5303# 1058 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5315# 1068 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5316#if defined(MFC_OpenACC)
5317# 1068 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5319# 1068 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5320#elif defined(MFC_OpenMP)
5321# 1068 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5323# 1068 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5325# 1068 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5329# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5331# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5332#if defined(MFC_OpenACC)
5333# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5335# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5336#elif defined(MFC_OpenMP)
5337# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5339# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5341# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5343# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5345# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5350 q_prim_rsy_vf(j, k, r, momxb + 1) = q_prim_vf(momxb + 1)%sf(k,
dj*(n - 2*j) + j, r)*sign(1._wp, &
5356# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5357#if defined(MFC_OpenACC)
5358# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5360# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5361#elif defined(MFC_OpenMP)
5362# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5364# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5366# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5370# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5372# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5373#if defined(MFC_OpenACC)
5374# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5376# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5377#elif defined(MFC_OpenMP)
5378# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5380# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5382# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5384# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5386# 1081 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5391 do j = -1, buff_size
5398# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5399#if defined(MFC_OpenACC)
5400# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5402# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5403#elif defined(MFC_OpenMP)
5404# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5406# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5408# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5412# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5414# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5415#if defined(MFC_OpenACC)
5416# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5418# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5419#elif defined(MFC_OpenMP)
5420# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5422# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5424# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5426# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5428# 1093 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5432 do j = -1, buff_size
5433 flux_rsy_vf_l(j, k, r, momxb + 1) = flux_vf(momxb + 1)%sf(k,
dj*((n - 1) - 2*j) + j, r)
5438# 1101 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5439#if defined(MFC_OpenACC)
5440# 1101 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5442# 1101 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5443#elif defined(MFC_OpenMP)
5444# 1101 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5446# 1101 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5448# 1101 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5451 if (riemann_solver == 1)
then
5453# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5455# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5456#if defined(MFC_OpenACC)
5457# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5459# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5460#elif defined(MFC_OpenMP)
5461# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5463# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5465# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5467# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5469# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5474 do j = -1, buff_size
5481# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5482#if defined(MFC_OpenACC)
5483# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5485# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5486#elif defined(MFC_OpenMP)
5487# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5489# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5491# 1114 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5495# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5497# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5498#if defined(MFC_OpenACC)
5499# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5501# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5502#elif defined(MFC_OpenMP)
5503# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5505# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5507# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5509# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5511# 1116 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5515 do j = -1, buff_size
5516 flux_src_rsy_vf_l(j, k, r, advxb) = flux_src_vf(advxb)%sf(k,
dj*((n - 1) - 2*j) + j, r)*sign(1._wp, &
5522# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5523#if defined(MFC_OpenACC)
5524# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5526# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5527#elif defined(MFC_OpenMP)
5528# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5530# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5532# 1125 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5541# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5543# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5544#if defined(MFC_OpenACC)
5545# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5547# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5548#elif defined(MFC_OpenMP)
5549# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5551# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5553# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5555# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5557# 1132 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5569# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5570#if defined(MFC_OpenACC)
5571# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5573# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5574#elif defined(MFC_OpenMP)
5575# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5577# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5579# 1142 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5583# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5585# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5586#if defined(MFC_OpenACC)
5587# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5589# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5590#elif defined(MFC_OpenMP)
5591# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5593# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5595# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5597# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5599# 1144 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5609# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5610#if defined(MFC_OpenACC)
5611# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5613# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5614#elif defined(MFC_OpenMP)
5615# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5617# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5619# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5623# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5625# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5626#if defined(MFC_OpenACC)
5627# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5629# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5630#elif defined(MFC_OpenMP)
5631# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5633# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5635# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5637# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5639# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5644 do j = -1, buff_size
5651# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5652#if defined(MFC_OpenACC)
5653# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5655# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5656#elif defined(MFC_OpenMP)
5657# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5659# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5661# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5665# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5667# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5668#if defined(MFC_OpenACC)
5669# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5671# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5672#elif defined(MFC_OpenMP)
5673# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5675# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5677# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5679# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5681# 1166 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5685 do j = -1, buff_size
5686 flux_rsz_vf_l(j, k, r, momxe) = flux_vf(momxe)%sf(r, k,
dj*((p - 1) - 2*j) + j)
5691# 1174 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5692#if defined(MFC_OpenACC)
5693# 1174 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5695# 1174 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5696#elif defined(MFC_OpenMP)
5697# 1174 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5699# 1174 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5701# 1174 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5704 if (riemann_solver == 1)
then
5706# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5708# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5709#if defined(MFC_OpenACC)
5710# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5712# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5713#elif defined(MFC_OpenMP)
5714# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5716# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5718# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5720# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5722# 1177 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5727 do j = -1, buff_size
5734# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5735#if defined(MFC_OpenACC)
5736# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5738# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5739#elif defined(MFC_OpenMP)
5740# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5742# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5744# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5748# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5750# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5751#if defined(MFC_OpenACC)
5752# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5754# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5755#elif defined(MFC_OpenMP)
5756# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5758# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5760# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5762# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5764# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5768 do j = -1, buff_size
5769 flux_src_rsz_vf_l(j, k, r, advxb) = flux_src_vf(advxb)%sf(r, k,
dj*((p - 1) - 2*j) + j)*sign(1._wp, &
5775# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5776#if defined(MFC_OpenACC)
5777# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5779# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5780#elif defined(MFC_OpenMP)
5781# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5783# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5785# 1198 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"