2153 integer,
intent(in) :: cbc_dir_in, cbc_loc_in
2156 real(wp),
dimension(0:buff_size + 1) :: s_cb
2168 s_cb(i + 1) = s_cb(i) +
ds(i)
2172# 470 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2173 if (cbc_dir_in == 1 .and. recon_type == weno_type)
then
2174 if (weno_order == 1)
then
2181 elseif (weno_order == 3)
then
2188 pi_coef_x(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2194 fd_coef_x(0, cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) &
2195 - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2203 ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))* &
2204 (s_cb(1) - s_cb(3)))/((s_cb(1) - s_cb(4))* &
2205 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2207 ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))* &
2208 ((s_cb(1) - s_cb(3))*(s_cb(1) - s_cb(3)) - &
2209 (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + &
2210 (s_cb(4) - s_cb(1)))))/ &
2211 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2212 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2214 (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))* &
2215 (s_cb(1) - s_cb(3)) + ((s_cb(0) - s_cb(2)) + &
2216 (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/ &
2217 ((s_cb(2) - s_cb(0))*(s_cb(0) - s_cb(3))* &
2218 (s_cb(0) - s_cb(4)))
2220 ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))* &
2221 (s_cb(2) - s_cb(3)))/((s_cb(2) - s_cb(4))* &
2222 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2224 ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))* &
2225 ((s_cb(1) - s_cb(3))*(s_cb(2) - s_cb(3)) + &
2226 (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + &
2227 (s_cb(2) - s_cb(4)))))/ &
2228 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2229 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2231 ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))* &
2232 (s_cb(2) - s_cb(4)))/((s_cb(0) - s_cb(2))* &
2233 (s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2237# 470 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2238 if (cbc_dir_in == 2 .and. recon_type == weno_type)
then
2239 if (weno_order == 1)
then
2246 elseif (weno_order == 3)
then
2253 pi_coef_y(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2259 fd_coef_y(0, cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) &
2260 - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2268 ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))* &
2269 (s_cb(1) - s_cb(3)))/((s_cb(1) - s_cb(4))* &
2270 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2272 ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))* &
2273 ((s_cb(1) - s_cb(3))*(s_cb(1) - s_cb(3)) - &
2274 (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + &
2275 (s_cb(4) - s_cb(1)))))/ &
2276 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2277 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2279 (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))* &
2280 (s_cb(1) - s_cb(3)) + ((s_cb(0) - s_cb(2)) + &
2281 (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/ &
2282 ((s_cb(2) - s_cb(0))*(s_cb(0) - s_cb(3))* &
2283 (s_cb(0) - s_cb(4)))
2285 ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))* &
2286 (s_cb(2) - s_cb(3)))/((s_cb(2) - s_cb(4))* &
2287 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2289 ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))* &
2290 ((s_cb(1) - s_cb(3))*(s_cb(2) - s_cb(3)) + &
2291 (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + &
2292 (s_cb(2) - s_cb(4)))))/ &
2293 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2294 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2296 ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))* &
2297 (s_cb(2) - s_cb(4)))/((s_cb(0) - s_cb(2))* &
2298 (s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2302# 470 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2303 if (cbc_dir_in == 3 .and. recon_type == weno_type)
then
2304 if (weno_order == 1)
then
2311 elseif (weno_order == 3)
then
2318 pi_coef_z(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2324 fd_coef_z(0, cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) &
2325 - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2333 ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))* &
2334 (s_cb(1) - s_cb(3)))/((s_cb(1) - s_cb(4))* &
2335 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2337 ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))* &
2338 ((s_cb(1) - s_cb(3))*(s_cb(1) - s_cb(3)) - &
2339 (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + &
2340 (s_cb(4) - s_cb(1)))))/ &
2341 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2342 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2344 (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))* &
2345 (s_cb(1) - s_cb(3)) + ((s_cb(0) - s_cb(2)) + &
2346 (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/ &
2347 ((s_cb(2) - s_cb(0))*(s_cb(0) - s_cb(3))* &
2348 (s_cb(0) - s_cb(4)))
2350 ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))* &
2351 (s_cb(2) - s_cb(3)))/((s_cb(2) - s_cb(4))* &
2352 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2354 ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))* &
2355 ((s_cb(1) - s_cb(3))*(s_cb(2) - s_cb(3)) + &
2356 (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + &
2357 (s_cb(2) - s_cb(4)))))/ &
2358 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2359 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2361 ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))* &
2362 (s_cb(2) - s_cb(4)))/((s_cb(0) - s_cb(2))* &
2363 (s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2367# 535 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2462 subroutine s_cbc(q_prim_vf, flux_vf, flux_src_vf, &
2463 cbc_dir_norm, cbc_loc_norm, &
2466 type(scalar_field), &
2467 dimension(sys_size), &
2468 intent(in) :: q_prim_vf
2470 type(scalar_field), &
2471 dimension(sys_size), &
2472 intent(inout) :: flux_vf, flux_src_vf
2474 integer,
intent(in) :: cbc_dir_norm, cbc_loc_norm
2476 type(int_bounds_info),
intent(in) :: ix, iy, iz
2478 real(wp) :: dpres_dt
2479 real(wp) :: dgamma_dt
2480 real(wp) :: dpi_inf_dt
2482 real(wp) :: dpres_ds
2483# 643 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2484 real(wp),
dimension(sys_size) ::
l
2485# 645 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2486# 654 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2487 real(wp),
dimension(num_fluids) :: alpha_rho, dalpha_rho_ds, mf
2488 real(wp),
dimension(num_vels) :: vel, dvel_ds
2489 real(wp),
dimension(num_fluids) :: adv_local, dadv_ds
2490 real(wp),
dimension(num_fluids) :: dadv_dt
2491 real(wp),
dimension(num_dims) :: dvel_dt
2492 real(wp),
dimension(num_fluids) :: dalpha_rho_dt
2493 real(wp),
dimension(num_species) :: ys, h_k, dys_dt, dys_ds, xs, gamma_i, cp_i
2494# 662 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2495 real(wp),
dimension(2) :: re_cbc
2496 real(wp),
dimension(3) :: lambda
2507 real(wp) :: t, sum_enthalpies
2508 real(wp) :: cv, cp, e_mix, mw, r_gas
2510 real(wp) :: vel_k_sum, vel_dv_dt_sum
2512 integer :: i,
j,
k, r
2522# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2523#if defined(MFC_OpenACC)
2524# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2526# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2527#elif defined(MFC_OpenMP)
2528# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2530# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2538# 696 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2539 if (
cbc_dir == 1 .and. recon_type == weno_type)
then
2542 if (weno_order == 3 .or. dummy)
then
2547 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
2550# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2552# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2553#if defined(MFC_OpenACC)
2554# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2556# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2557#elif defined(MFC_OpenMP)
2558# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2560# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2562# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2564# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2566# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2568# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2581# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2583# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2584#if defined(MFC_OpenACC)
2585# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2587# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2588#elif defined(MFC_OpenMP)
2589# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2591# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2593# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2595# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2597# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2601# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2603# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2604#if defined(MFC_OpenACC)
2605# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2607# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2608#elif defined(MFC_OpenMP)
2609# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2611# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2613# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2615# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2617# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2619# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2632# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2634# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2635#if defined(MFC_OpenACC)
2636# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2638# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2639#elif defined(MFC_OpenMP)
2640# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2642# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2644# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2646# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2648# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2653 if (weno_order == 5 .or. dummy)
then
2657 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
2660# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2662# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2663#if defined(MFC_OpenACC)
2664# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2666# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2667#elif defined(MFC_OpenMP)
2668# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2670# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2672# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2674# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2676# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2678# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2699# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2701# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2702#if defined(MFC_OpenACC)
2703# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2705# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2706#elif defined(MFC_OpenMP)
2707# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2709# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2711# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2713# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2715# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2719# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2721# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2722#if defined(MFC_OpenACC)
2723# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2725# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2726#elif defined(MFC_OpenMP)
2727# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2729# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2731# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2733# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2735# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2737# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2758# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2760# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2761#if defined(MFC_OpenACC)
2762# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2764# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2765#elif defined(MFC_OpenMP)
2766# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2768# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2770# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2772# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2774# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2781# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2783# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2784#if defined(MFC_OpenACC)
2785# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2787# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2788#elif defined(MFC_OpenMP)
2789# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2791# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2793# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2795# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2797# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2799# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2806# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2807#if defined(MFC_OpenACC)
2808# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2810# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2811#elif defined(MFC_OpenMP)
2812# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2814# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2821# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2822#if defined(MFC_OpenACC)
2823# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2825# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2826#elif defined(MFC_OpenMP)
2827# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2829# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2837# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2838#if defined(MFC_OpenACC)
2839# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2841# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2842#elif defined(MFC_OpenMP)
2843# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2845# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2848 vel_k_sum = vel_k_sum + vel(i)**2._wp
2854# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2855#if defined(MFC_OpenACC)
2856# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2858# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2859#elif defined(MFC_OpenMP)
2860# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2862# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2864 do i = 1, advxe - e_idx
2868 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
2871# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2872#if defined(MFC_OpenACC)
2873# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2875# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2876#elif defined(MFC_OpenMP)
2877# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2879# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2882 mf(i) = alpha_rho(i)/rho
2887# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2888#if defined(MFC_OpenACC)
2889# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2891# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2892#elif defined(MFC_OpenMP)
2893# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2895# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2897 do i = chemxb, chemxe
2901 call get_mixture_molecular_weight(ys, mw)
2902 r_gas = gas_constant/mw
2904 call get_mixture_specific_heat_cp_mass(t, ys, cp)
2905 call get_mixture_energy_mass(t, ys, e_mix)
2906 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
2907 if (chem_params%gamma_method == 1)
then
2909 call get_mole_fractions(mw, ys, xs)
2910 call get_species_specific_heats_r(t, cp_i)
2911 gamma_i = cp_i/(cp_i - 1.0_wp)
2912 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
2913 else if (chem_params%gamma_method == 2)
then
2915 call get_mixture_specific_heat_cv_mass(t, ys, cv)
2916 gamma = 1.0_wp/(cp/cv - 1.0_wp)
2919 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
2925 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
2930# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2931#if defined(MFC_OpenACC)
2932# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2934# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2935#elif defined(MFC_OpenMP)
2936# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2938# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2941 dalpha_rho_ds(i) = 0._wp
2945# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2946#if defined(MFC_OpenACC)
2947# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2949# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2950#elif defined(MFC_OpenMP)
2951# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2953# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2961# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2962#if defined(MFC_OpenACC)
2963# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2965# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2966#elif defined(MFC_OpenMP)
2967# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2969# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2971 do i = 1, advxe - e_idx
2977# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2978#if defined(MFC_OpenACC)
2979# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2981# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2982#elif defined(MFC_OpenMP)
2983# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2985# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2987 do i = 1, num_species
2993# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2994#if defined(MFC_OpenACC)
2995# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2997# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2998#elif defined(MFC_OpenMP)
2999# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3001# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3006# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3007#if defined(MFC_OpenACC)
3008# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3010# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3011#elif defined(MFC_OpenMP)
3012# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3014# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3022# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3023#if defined(MFC_OpenACC)
3024# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3026# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3027#elif defined(MFC_OpenMP)
3028# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3030# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3042# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3043#if defined(MFC_OpenACC)
3044# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3046# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3047#elif defined(MFC_OpenMP)
3048# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3050# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3052 do i = 1, advxe - e_idx
3060# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3061#if defined(MFC_OpenACC)
3062# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3064# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3065#elif defined(MFC_OpenMP)
3066# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3068# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3070 do i = 1, num_species
3079 lambda(1) = vel(dir_idx(1)) - c
3080 lambda(2) = vel(dir_idx(1))
3081 lambda(3) = vel(dir_idx(1)) + c
3083 ma = vel(dir_idx(1))/c
3085 if ((
cbc_loc == -1 .and.
bcxb == bc_char_slip_wall) .or. &
3086 (
cbc_loc == 1 .and.
bcxe == bc_char_slip_wall))
then
3087 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
3088 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_buffer) .or. &
3089 (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_buffer))
then
3090 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
3091 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_inflow) .or. &
3092 (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_inflow))
then
3093 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
3095 if (bc_x%grcbc_in)
then
3097# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3098#if defined(MFC_OpenACC)
3099# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3101# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3102#elif defined(MFC_OpenMP)
3103# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3105# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3111 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(1, dir_idx(2)))/
del_in(1)
3113 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(1, dir_idx(3)))/
del_in(1)
3117# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3118#if defined(MFC_OpenACC)
3119# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3121# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3122#elif defined(MFC_OpenMP)
3123# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3125# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3127 do i = e_idx, advxe - 1
3128 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, 1))/
del_in(1)
3132 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_outflow) .or. &
3133 (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_outflow))
then
3134 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
3136 if (bc_x%grcbc_out)
then
3140 if (bc_x%grcbc_vel_out)
then
3141 l(advxe) =
l(advxe) + rho*c**2._wp*(1._wp - ma)*(vel(dir_idx(1)) +
vel_out(1, dir_idx(1))*sign(1,
cbc_loc))/
del_out(1)
3144 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_ff_sub_outflow) .or. &
3145 (
cbc_loc == 1 .and.
bcxe == bc_char_ff_sub_outflow))
then
3146 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
3147 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_cp_sub_outflow) .or. &
3148 (
cbc_loc == 1 .and.
bcxe == bc_char_cp_sub_outflow))
then
3149 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
3150 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_sup_inflow) .or. &
3151 (
cbc_loc == 1 .and.
bcxe == bc_char_sup_inflow))
then
3152 call s_compute_supersonic_inflow_l(
l)
3153 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_sup_outflow) .or. &
3154 (
cbc_loc == 1 .and.
bcxe == bc_char_sup_outflow))
then
3155 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
3160 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1)) &
3163 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
3167# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3168#if defined(MFC_OpenACC)
3169# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3171# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3172#elif defined(MFC_OpenMP)
3173# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3175# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3178 dalpha_rho_dt(i) = &
3179 -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
3183# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3184#if defined(MFC_OpenACC)
3185# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3187# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3188#elif defined(MFC_OpenMP)
3189# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3191# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3194 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))* &
3195 (
l(1) -
l(advxe))/(2._wp*rho*c) + &
3196 (dir_flg(dir_idx(i)) - 1._wp)* &
3200 vel_dv_dt_sum = 0._wp
3202# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3203#if defined(MFC_OpenACC)
3204# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3206# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3207#elif defined(MFC_OpenMP)
3208# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3210# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3213 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
3218# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3219#if defined(MFC_OpenACC)
3220# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3222# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3223#elif defined(MFC_OpenMP)
3224# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3226# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3228 do i = 1, num_species
3229 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
3236# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3237#if defined(MFC_OpenACC)
3238# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3240# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3241#elif defined(MFC_OpenMP)
3242# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3244# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3246 do i = 1, advxe - e_idx
3247 dadv_dt(i) = -
l(momxe + i)
3251# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3252#if defined(MFC_OpenACC)
3253# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3255# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3256#elif defined(MFC_OpenMP)
3257# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3259# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3261 do i = 1, advxe - e_idx
3262 dadv_dt(i) = -
l(momxe + i)
3266 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
3268 if (model_eqns == 1)
then
3269 drho_dt = dalpha_rho_dt(1)
3270 dgamma_dt = dadv_dt(1)
3271# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3272 dpi_inf_dt = dadv_dt(2)
3273# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3276# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3277#if defined(MFC_OpenACC)
3278# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3280# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3281#elif defined(MFC_OpenMP)
3282# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3284# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3286 do i = 1, num_fluids
3287 drho_dt = drho_dt + dalpha_rho_dt(i)
3288 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
3289 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
3290 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
3296# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3297#if defined(MFC_OpenACC)
3298# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3300# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3301#elif defined(MFC_OpenMP)
3302# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3304# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3308 +
ds(0)*dalpha_rho_dt(i)
3312# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3313#if defined(MFC_OpenACC)
3314# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3316# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3317#elif defined(MFC_OpenMP)
3318# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3320# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3324 +
ds(0)*(vel(i - contxe)*drho_dt &
3325 + rho*dvel_dt(i - contxe))
3330 call get_species_enthalpies_rt(t, h_k)
3331 sum_enthalpies = 0._wp
3333# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3334#if defined(MFC_OpenACC)
3335# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3337# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3338#elif defined(MFC_OpenMP)
3339# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3341# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3343 do i = 1, num_species
3345# 1064 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3346 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
3347 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
3348# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3351 +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) + sum_enthalpies)
3353# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3354#if defined(MFC_OpenACC)
3355# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3357# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3358#elif defined(MFC_OpenMP)
3359# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3361# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3363 do i = 1, num_species
3365 +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
3369 +
ds(0)*(pres*dgamma_dt &
3373 + rho*vel_dv_dt_sum &
3374 + 5.e-1_wp*drho_dt*vel_k_sum)
3377 if (riemann_solver == 1)
then
3379# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3380#if defined(MFC_OpenACC)
3381# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3383# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3384#elif defined(MFC_OpenMP)
3385# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3387# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3394# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3395#if defined(MFC_OpenACC)
3396# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3398# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3399#elif defined(MFC_OpenMP)
3400# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3402# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3406 1._wp/max(abs(vel(dir_idx(1))), sgm_eps) &
3407 *sign(1._wp, vel(dir_idx(1))) &
3411 +
ds(0)*dadv_dt(i - e_idx))
3417# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3418#if defined(MFC_OpenACC)
3419# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3421# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3422#elif defined(MFC_OpenMP)
3423# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3425# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3429 ds(0)*dadv_dt(i - e_idx)
3433# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3434#if defined(MFC_OpenACC)
3435# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3437# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3438#elif defined(MFC_OpenMP)
3439# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3441# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3453# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3455# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3456#if defined(MFC_OpenACC)
3457# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3459# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3460#elif defined(MFC_OpenMP)
3461# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3463# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3465# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3467# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3469# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3472# 696 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3473 if (
cbc_dir == 2 .and. recon_type == weno_type)
then
3476 if (weno_order == 3 .or. dummy)
then
3481 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
3484# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3486# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3487#if defined(MFC_OpenACC)
3488# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3490# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3491#elif defined(MFC_OpenMP)
3492# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3494# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3496# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3498# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3500# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3502# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3515# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3517# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3518#if defined(MFC_OpenACC)
3519# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3521# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3522#elif defined(MFC_OpenMP)
3523# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3525# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3527# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3529# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3531# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3535# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3537# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3538#if defined(MFC_OpenACC)
3539# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3541# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3542#elif defined(MFC_OpenMP)
3543# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3545# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3547# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3549# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3551# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3553# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3566# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3568# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3569#if defined(MFC_OpenACC)
3570# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3572# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3573#elif defined(MFC_OpenMP)
3574# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3576# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3578# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3580# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3582# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3587 if (weno_order == 5 .or. dummy)
then
3591 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
3594# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3596# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3597#if defined(MFC_OpenACC)
3598# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3600# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3601#elif defined(MFC_OpenMP)
3602# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3604# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3606# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3608# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3610# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3612# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3633# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3635# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3636#if defined(MFC_OpenACC)
3637# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3639# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3640#elif defined(MFC_OpenMP)
3641# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3643# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3645# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3647# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3649# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3653# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3655# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3656#if defined(MFC_OpenACC)
3657# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3659# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3660#elif defined(MFC_OpenMP)
3661# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3663# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3665# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3667# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3669# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3671# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3692# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3694# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3695#if defined(MFC_OpenACC)
3696# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3698# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3699#elif defined(MFC_OpenMP)
3700# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3702# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3704# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3706# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3708# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3715# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3717# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3718#if defined(MFC_OpenACC)
3719# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3721# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3722#elif defined(MFC_OpenMP)
3723# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3725# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3727# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3729# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3731# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3733# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3740# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3741#if defined(MFC_OpenACC)
3742# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3744# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3745#elif defined(MFC_OpenMP)
3746# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3748# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3755# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3756#if defined(MFC_OpenACC)
3757# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3759# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3760#elif defined(MFC_OpenMP)
3761# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3763# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3771# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3772#if defined(MFC_OpenACC)
3773# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3775# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3776#elif defined(MFC_OpenMP)
3777# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3779# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3782 vel_k_sum = vel_k_sum + vel(i)**2._wp
3788# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3789#if defined(MFC_OpenACC)
3790# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3792# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3793#elif defined(MFC_OpenMP)
3794# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3796# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3798 do i = 1, advxe - e_idx
3802 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
3805# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3806#if defined(MFC_OpenACC)
3807# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3809# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3810#elif defined(MFC_OpenMP)
3811# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3813# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3816 mf(i) = alpha_rho(i)/rho
3821# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3822#if defined(MFC_OpenACC)
3823# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3825# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3826#elif defined(MFC_OpenMP)
3827# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3829# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3831 do i = chemxb, chemxe
3835 call get_mixture_molecular_weight(ys, mw)
3836 r_gas = gas_constant/mw
3838 call get_mixture_specific_heat_cp_mass(t, ys, cp)
3839 call get_mixture_energy_mass(t, ys, e_mix)
3840 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
3841 if (chem_params%gamma_method == 1)
then
3843 call get_mole_fractions(mw, ys, xs)
3844 call get_species_specific_heats_r(t, cp_i)
3845 gamma_i = cp_i/(cp_i - 1.0_wp)
3846 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
3847 else if (chem_params%gamma_method == 2)
then
3849 call get_mixture_specific_heat_cv_mass(t, ys, cv)
3850 gamma = 1.0_wp/(cp/cv - 1.0_wp)
3853 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
3859 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
3864# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3865#if defined(MFC_OpenACC)
3866# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3868# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3869#elif defined(MFC_OpenMP)
3870# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3872# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3875 dalpha_rho_ds(i) = 0._wp
3879# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3880#if defined(MFC_OpenACC)
3881# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3883# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3884#elif defined(MFC_OpenMP)
3885# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3887# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3895# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3896#if defined(MFC_OpenACC)
3897# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3899# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3900#elif defined(MFC_OpenMP)
3901# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3903# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3905 do i = 1, advxe - e_idx
3911# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3912#if defined(MFC_OpenACC)
3913# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3915# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3916#elif defined(MFC_OpenMP)
3917# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3919# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3921 do i = 1, num_species
3927# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3928#if defined(MFC_OpenACC)
3929# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3931# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3932#elif defined(MFC_OpenMP)
3933# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3935# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3940# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3941#if defined(MFC_OpenACC)
3942# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3944# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3945#elif defined(MFC_OpenMP)
3946# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3948# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3956# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3957#if defined(MFC_OpenACC)
3958# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3960# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3961#elif defined(MFC_OpenMP)
3962# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3964# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3976# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3977#if defined(MFC_OpenACC)
3978# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3980# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3981#elif defined(MFC_OpenMP)
3982# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3984# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3986 do i = 1, advxe - e_idx
3994# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3995#if defined(MFC_OpenACC)
3996# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3998# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3999#elif defined(MFC_OpenMP)
4000# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4002# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4004 do i = 1, num_species
4013 lambda(1) = vel(dir_idx(1)) - c
4014 lambda(2) = vel(dir_idx(1))
4015 lambda(3) = vel(dir_idx(1)) + c
4017 ma = vel(dir_idx(1))/c
4019 if ((
cbc_loc == -1 .and.
bcyb == bc_char_slip_wall) .or. &
4020 (
cbc_loc == 1 .and.
bcye == bc_char_slip_wall))
then
4021 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4022 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_buffer) .or. &
4023 (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_buffer))
then
4024 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4025 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_inflow) .or. &
4026 (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_inflow))
then
4027 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4029 if (bc_y%grcbc_in)
then
4031# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4032#if defined(MFC_OpenACC)
4033# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4035# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4036#elif defined(MFC_OpenMP)
4037# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4039# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4045 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(2, dir_idx(2)))/
del_in(2)
4047 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(2, dir_idx(3)))/
del_in(2)
4051# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4052#if defined(MFC_OpenACC)
4053# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4055# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4056#elif defined(MFC_OpenMP)
4057# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4059# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4061 do i = e_idx, advxe - 1
4062 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, 2))/
del_in(2)
4066 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_outflow) .or. &
4067 (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_outflow))
then
4068 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4070 if (bc_y%grcbc_out)
then
4074 if (bc_y%grcbc_vel_out)
then
4075 l(advxe) =
l(advxe) + rho*c**2._wp*(1._wp - ma)*(vel(dir_idx(1)) +
vel_out(2, dir_idx(1))*sign(1,
cbc_loc))/
del_out(2)
4078 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_ff_sub_outflow) .or. &
4079 (
cbc_loc == 1 .and.
bcye == bc_char_ff_sub_outflow))
then
4080 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
4081 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_cp_sub_outflow) .or. &
4082 (
cbc_loc == 1 .and.
bcye == bc_char_cp_sub_outflow))
then
4083 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
4084 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_sup_inflow) .or. &
4085 (
cbc_loc == 1 .and.
bcye == bc_char_sup_inflow))
then
4086 call s_compute_supersonic_inflow_l(
l)
4087 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_sup_outflow) .or. &
4088 (
cbc_loc == 1 .and.
bcye == bc_char_sup_outflow))
then
4089 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4094 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1)) &
4097 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
4101# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4102#if defined(MFC_OpenACC)
4103# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4105# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4106#elif defined(MFC_OpenMP)
4107# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4109# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4112 dalpha_rho_dt(i) = &
4113 -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
4117# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4118#if defined(MFC_OpenACC)
4119# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4121# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4122#elif defined(MFC_OpenMP)
4123# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4125# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4128 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))* &
4129 (
l(1) -
l(advxe))/(2._wp*rho*c) + &
4130 (dir_flg(dir_idx(i)) - 1._wp)* &
4134 vel_dv_dt_sum = 0._wp
4136# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4137#if defined(MFC_OpenACC)
4138# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4140# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4141#elif defined(MFC_OpenMP)
4142# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4144# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4147 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
4152# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4153#if defined(MFC_OpenACC)
4154# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4156# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4157#elif defined(MFC_OpenMP)
4158# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4160# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4162 do i = 1, num_species
4163 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
4170# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4171#if defined(MFC_OpenACC)
4172# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4174# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4175#elif defined(MFC_OpenMP)
4176# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4178# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4180 do i = 1, advxe - e_idx
4181 dadv_dt(i) = -
l(momxe + i)
4185# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4186#if defined(MFC_OpenACC)
4187# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4189# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4190#elif defined(MFC_OpenMP)
4191# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4193# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4195 do i = 1, advxe - e_idx
4196 dadv_dt(i) = -
l(momxe + i)
4200 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
4202 if (model_eqns == 1)
then
4203 drho_dt = dalpha_rho_dt(1)
4204 dgamma_dt = dadv_dt(1)
4205# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4206 dpi_inf_dt = dadv_dt(2)
4207# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4210# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4211#if defined(MFC_OpenACC)
4212# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4214# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4215#elif defined(MFC_OpenMP)
4216# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4218# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4220 do i = 1, num_fluids
4221 drho_dt = drho_dt + dalpha_rho_dt(i)
4222 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
4223 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
4224 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
4230# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4231#if defined(MFC_OpenACC)
4232# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4234# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4235#elif defined(MFC_OpenMP)
4236# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4238# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4242 +
ds(0)*dalpha_rho_dt(i)
4246# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4247#if defined(MFC_OpenACC)
4248# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4250# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4251#elif defined(MFC_OpenMP)
4252# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4254# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4258 +
ds(0)*(vel(i - contxe)*drho_dt &
4259 + rho*dvel_dt(i - contxe))
4264 call get_species_enthalpies_rt(t, h_k)
4265 sum_enthalpies = 0._wp
4267# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4268#if defined(MFC_OpenACC)
4269# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4271# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4272#elif defined(MFC_OpenMP)
4273# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4275# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4277 do i = 1, num_species
4279# 1064 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4280 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
4281 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
4282# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4285 +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) + sum_enthalpies)
4287# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4288#if defined(MFC_OpenACC)
4289# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4291# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4292#elif defined(MFC_OpenMP)
4293# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4295# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4297 do i = 1, num_species
4299 +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
4303 +
ds(0)*(pres*dgamma_dt &
4307 + rho*vel_dv_dt_sum &
4308 + 5.e-1_wp*drho_dt*vel_k_sum)
4311 if (riemann_solver == 1)
then
4313# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4314#if defined(MFC_OpenACC)
4315# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4317# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4318#elif defined(MFC_OpenMP)
4319# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4321# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4328# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4329#if defined(MFC_OpenACC)
4330# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4332# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4333#elif defined(MFC_OpenMP)
4334# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4336# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4340 1._wp/max(abs(vel(dir_idx(1))), sgm_eps) &
4341 *sign(1._wp, vel(dir_idx(1))) &
4345 +
ds(0)*dadv_dt(i - e_idx))
4351# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4352#if defined(MFC_OpenACC)
4353# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4355# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4356#elif defined(MFC_OpenMP)
4357# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4359# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4363 ds(0)*dadv_dt(i - e_idx)
4367# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4368#if defined(MFC_OpenACC)
4369# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4371# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4372#elif defined(MFC_OpenMP)
4373# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4375# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4387# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4389# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4390#if defined(MFC_OpenACC)
4391# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4393# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4394#elif defined(MFC_OpenMP)
4395# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4397# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4399# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4401# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4403# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4406# 696 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4407 if (
cbc_dir == 3 .and. recon_type == weno_type)
then
4410 if (weno_order == 3 .or. dummy)
then
4415 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
4418# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4420# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4421#if defined(MFC_OpenACC)
4422# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4424# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4425#elif defined(MFC_OpenMP)
4426# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4428# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4430# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4432# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4434# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4436# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4449# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4451# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4452#if defined(MFC_OpenACC)
4453# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4455# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4456#elif defined(MFC_OpenMP)
4457# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4459# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4461# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4463# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4465# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4469# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4471# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4472#if defined(MFC_OpenACC)
4473# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4475# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4476#elif defined(MFC_OpenMP)
4477# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4479# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4481# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4483# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4485# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4487# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4500# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4502# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4503#if defined(MFC_OpenACC)
4504# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4506# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4507#elif defined(MFC_OpenMP)
4508# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4510# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4512# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4514# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4516# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4521 if (weno_order == 5 .or. dummy)
then
4525 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
4528# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4530# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4531#if defined(MFC_OpenACC)
4532# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4534# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4535#elif defined(MFC_OpenMP)
4536# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4538# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4540# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4542# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4544# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4546# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4567# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4569# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4570#if defined(MFC_OpenACC)
4571# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4573# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4574#elif defined(MFC_OpenMP)
4575# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4577# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4579# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4581# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4583# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4587# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4589# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4590#if defined(MFC_OpenACC)
4591# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4593# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4594#elif defined(MFC_OpenMP)
4595# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4597# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4599# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4601# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4603# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4605# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4626# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4628# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4629#if defined(MFC_OpenACC)
4630# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4632# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4633#elif defined(MFC_OpenMP)
4634# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4636# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4638# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4640# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4642# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4649# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4651# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4652#if defined(MFC_OpenACC)
4653# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4655# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4656#elif defined(MFC_OpenMP)
4657# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4659# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4661# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4663# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4665# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4667# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4674# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4675#if defined(MFC_OpenACC)
4676# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4678# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4679#elif defined(MFC_OpenMP)
4680# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4682# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4689# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4690#if defined(MFC_OpenACC)
4691# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4693# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4694#elif defined(MFC_OpenMP)
4695# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4697# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4705# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4706#if defined(MFC_OpenACC)
4707# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4709# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4710#elif defined(MFC_OpenMP)
4711# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4713# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4716 vel_k_sum = vel_k_sum + vel(i)**2._wp
4722# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4723#if defined(MFC_OpenACC)
4724# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4726# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4727#elif defined(MFC_OpenMP)
4728# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4730# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4732 do i = 1, advxe - e_idx
4736 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
4739# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4740#if defined(MFC_OpenACC)
4741# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4743# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4744#elif defined(MFC_OpenMP)
4745# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4747# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4750 mf(i) = alpha_rho(i)/rho
4755# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4756#if defined(MFC_OpenACC)
4757# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4759# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4760#elif defined(MFC_OpenMP)
4761# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4763# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4765 do i = chemxb, chemxe
4769 call get_mixture_molecular_weight(ys, mw)
4770 r_gas = gas_constant/mw
4772 call get_mixture_specific_heat_cp_mass(t, ys, cp)
4773 call get_mixture_energy_mass(t, ys, e_mix)
4774 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
4775 if (chem_params%gamma_method == 1)
then
4777 call get_mole_fractions(mw, ys, xs)
4778 call get_species_specific_heats_r(t, cp_i)
4779 gamma_i = cp_i/(cp_i - 1.0_wp)
4780 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
4781 else if (chem_params%gamma_method == 2)
then
4783 call get_mixture_specific_heat_cv_mass(t, ys, cv)
4784 gamma = 1.0_wp/(cp/cv - 1.0_wp)
4787 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
4793 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
4798# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4799#if defined(MFC_OpenACC)
4800# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4802# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4803#elif defined(MFC_OpenMP)
4804# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4806# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4809 dalpha_rho_ds(i) = 0._wp
4813# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4814#if defined(MFC_OpenACC)
4815# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4817# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4818#elif defined(MFC_OpenMP)
4819# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4821# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4829# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4830#if defined(MFC_OpenACC)
4831# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4833# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4834#elif defined(MFC_OpenMP)
4835# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4837# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4839 do i = 1, advxe - e_idx
4845# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4846#if defined(MFC_OpenACC)
4847# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4849# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4850#elif defined(MFC_OpenMP)
4851# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4853# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4855 do i = 1, num_species
4861# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4862#if defined(MFC_OpenACC)
4863# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4865# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4866#elif defined(MFC_OpenMP)
4867# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4869# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4874# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4875#if defined(MFC_OpenACC)
4876# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4878# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4879#elif defined(MFC_OpenMP)
4880# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4882# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4890# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4891#if defined(MFC_OpenACC)
4892# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4894# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4895#elif defined(MFC_OpenMP)
4896# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4898# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4910# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4911#if defined(MFC_OpenACC)
4912# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4914# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4915#elif defined(MFC_OpenMP)
4916# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4918# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4920 do i = 1, advxe - e_idx
4928# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4929#if defined(MFC_OpenACC)
4930# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4932# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4933#elif defined(MFC_OpenMP)
4934# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4936# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4938 do i = 1, num_species
4947 lambda(1) = vel(dir_idx(1)) - c
4948 lambda(2) = vel(dir_idx(1))
4949 lambda(3) = vel(dir_idx(1)) + c
4951 ma = vel(dir_idx(1))/c
4953 if ((
cbc_loc == -1 .and.
bczb == bc_char_slip_wall) .or. &
4954 (
cbc_loc == 1 .and.
bcze == bc_char_slip_wall))
then
4955 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4956 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_buffer) .or. &
4957 (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_buffer))
then
4958 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4959 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_inflow) .or. &
4960 (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_inflow))
then
4961 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4963 if (bc_z%grcbc_in)
then
4965# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4966#if defined(MFC_OpenACC)
4967# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4969# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4970#elif defined(MFC_OpenMP)
4971# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4973# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4979 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(3, dir_idx(2)))/
del_in(3)
4981 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(3, dir_idx(3)))/
del_in(3)
4985# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4986#if defined(MFC_OpenACC)
4987# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4989# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4990#elif defined(MFC_OpenMP)
4991# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4993# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4995 do i = e_idx, advxe - 1
4996 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, 3))/
del_in(3)
5000 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_outflow) .or. &
5001 (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_outflow))
then
5002 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
5004 if (bc_z%grcbc_out)
then
5008 if (bc_z%grcbc_vel_out)
then
5009 l(advxe) =
l(advxe) + rho*c**2._wp*(1._wp - ma)*(vel(dir_idx(1)) +
vel_out(3, dir_idx(1))*sign(1,
cbc_loc))/
del_out(3)
5012 else if ((
cbc_loc == -1 .and.
bczb == bc_char_ff_sub_outflow) .or. &
5013 (
cbc_loc == 1 .and.
bcze == bc_char_ff_sub_outflow))
then
5014 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
5015 else if ((
cbc_loc == -1 .and.
bczb == bc_char_cp_sub_outflow) .or. &
5016 (
cbc_loc == 1 .and.
bcze == bc_char_cp_sub_outflow))
then
5017 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
5018 else if ((
cbc_loc == -1 .and.
bczb == bc_char_sup_inflow) .or. &
5019 (
cbc_loc == 1 .and.
bcze == bc_char_sup_inflow))
then
5020 call s_compute_supersonic_inflow_l(
l)
5021 else if ((
cbc_loc == -1 .and.
bczb == bc_char_sup_outflow) .or. &
5022 (
cbc_loc == 1 .and.
bcze == bc_char_sup_outflow))
then
5023 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
5028 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1)) &
5031 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
5035# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5036#if defined(MFC_OpenACC)
5037# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5039# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5040#elif defined(MFC_OpenMP)
5041# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5043# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5046 dalpha_rho_dt(i) = &
5047 -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
5051# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5052#if defined(MFC_OpenACC)
5053# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5055# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5056#elif defined(MFC_OpenMP)
5057# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5059# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5062 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))* &
5063 (
l(1) -
l(advxe))/(2._wp*rho*c) + &
5064 (dir_flg(dir_idx(i)) - 1._wp)* &
5068 vel_dv_dt_sum = 0._wp
5070# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5071#if defined(MFC_OpenACC)
5072# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5074# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5075#elif defined(MFC_OpenMP)
5076# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5078# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5081 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
5086# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5087#if defined(MFC_OpenACC)
5088# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5090# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5091#elif defined(MFC_OpenMP)
5092# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5094# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5096 do i = 1, num_species
5097 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
5104# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5105#if defined(MFC_OpenACC)
5106# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5108# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5109#elif defined(MFC_OpenMP)
5110# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5112# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5114 do i = 1, advxe - e_idx
5115 dadv_dt(i) = -
l(momxe + i)
5119# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5120#if defined(MFC_OpenACC)
5121# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5123# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5124#elif defined(MFC_OpenMP)
5125# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5127# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5129 do i = 1, advxe - e_idx
5130 dadv_dt(i) = -
l(momxe + i)
5134 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
5136 if (model_eqns == 1)
then
5137 drho_dt = dalpha_rho_dt(1)
5138 dgamma_dt = dadv_dt(1)
5139# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5140 dpi_inf_dt = dadv_dt(2)
5141# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5144# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5145#if defined(MFC_OpenACC)
5146# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5148# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5149#elif defined(MFC_OpenMP)
5150# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5152# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5154 do i = 1, num_fluids
5155 drho_dt = drho_dt + dalpha_rho_dt(i)
5156 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
5157 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
5158 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
5164# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5165#if defined(MFC_OpenACC)
5166# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5168# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5169#elif defined(MFC_OpenMP)
5170# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5172# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5176 +
ds(0)*dalpha_rho_dt(i)
5180# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5181#if defined(MFC_OpenACC)
5182# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5184# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5185#elif defined(MFC_OpenMP)
5186# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5188# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5192 +
ds(0)*(vel(i - contxe)*drho_dt &
5193 + rho*dvel_dt(i - contxe))
5198 call get_species_enthalpies_rt(t, h_k)
5199 sum_enthalpies = 0._wp
5201# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5202#if defined(MFC_OpenACC)
5203# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5205# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5206#elif defined(MFC_OpenMP)
5207# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5209# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5211 do i = 1, num_species
5213# 1064 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5214 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
5215 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
5216# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5219 +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) + sum_enthalpies)
5221# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5222#if defined(MFC_OpenACC)
5223# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5225# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5226#elif defined(MFC_OpenMP)
5227# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5229# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5231 do i = 1, num_species
5233 +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
5237 +
ds(0)*(pres*dgamma_dt &
5241 + rho*vel_dv_dt_sum &
5242 + 5.e-1_wp*drho_dt*vel_k_sum)
5245 if (riemann_solver == 1)
then
5247# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5248#if defined(MFC_OpenACC)
5249# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5251# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5252#elif defined(MFC_OpenMP)
5253# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5255# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5262# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5263#if defined(MFC_OpenACC)
5264# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5266# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5267#elif defined(MFC_OpenMP)
5268# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5270# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5274 1._wp/max(abs(vel(dir_idx(1))), sgm_eps) &
5275 *sign(1._wp, vel(dir_idx(1))) &
5279 +
ds(0)*dadv_dt(i - e_idx))
5285# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5286#if defined(MFC_OpenACC)
5287# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5289# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5290#elif defined(MFC_OpenMP)
5291# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5293# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5297 ds(0)*dadv_dt(i - e_idx)
5301# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5302#if defined(MFC_OpenACC)
5303# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5305# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5306#elif defined(MFC_OpenMP)
5307# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5309# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5321# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5323# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5324#if defined(MFC_OpenACC)
5325# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5327# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5328#elif defined(MFC_OpenMP)
5329# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5331# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5333# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5335# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5337# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5340# 1123 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"