2155 integer,
intent(in) :: cbc_dir_in, cbc_loc_in
2158 real(wp),
dimension(0:buff_size + 1) :: s_cb
2170 s_cb(i + 1) = s_cb(i) +
ds(i)
2174# 470 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2175 if (cbc_dir_in == 1 .and. recon_type == weno_type)
then
2176 if (weno_order == 1)
then
2183 elseif (weno_order == 3)
then
2190 pi_coef_x(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2196 fd_coef_x(0, cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) &
2197 - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2205 ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))* &
2206 (s_cb(1) - s_cb(3)))/((s_cb(1) - s_cb(4))* &
2207 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2209 ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))* &
2210 ((s_cb(1) - s_cb(3))*(s_cb(1) - s_cb(3)) - &
2211 (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + &
2212 (s_cb(4) - s_cb(1)))))/ &
2213 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2214 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2216 (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))* &
2217 (s_cb(1) - s_cb(3)) + ((s_cb(0) - s_cb(2)) + &
2218 (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/ &
2219 ((s_cb(2) - s_cb(0))*(s_cb(0) - s_cb(3))* &
2220 (s_cb(0) - s_cb(4)))
2222 ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))* &
2223 (s_cb(2) - s_cb(3)))/((s_cb(2) - s_cb(4))* &
2224 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2226 ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))* &
2227 ((s_cb(1) - s_cb(3))*(s_cb(2) - s_cb(3)) + &
2228 (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + &
2229 (s_cb(2) - s_cb(4)))))/ &
2230 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2231 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2233 ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))* &
2234 (s_cb(2) - s_cb(4)))/((s_cb(0) - s_cb(2))* &
2235 (s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2239# 470 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2240 if (cbc_dir_in == 2 .and. recon_type == weno_type)
then
2241 if (weno_order == 1)
then
2248 elseif (weno_order == 3)
then
2255 pi_coef_y(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2261 fd_coef_y(0, cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) &
2262 - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2270 ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))* &
2271 (s_cb(1) - s_cb(3)))/((s_cb(1) - s_cb(4))* &
2272 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2274 ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))* &
2275 ((s_cb(1) - s_cb(3))*(s_cb(1) - s_cb(3)) - &
2276 (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + &
2277 (s_cb(4) - s_cb(1)))))/ &
2278 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2279 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2281 (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))* &
2282 (s_cb(1) - s_cb(3)) + ((s_cb(0) - s_cb(2)) + &
2283 (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/ &
2284 ((s_cb(2) - s_cb(0))*(s_cb(0) - s_cb(3))* &
2285 (s_cb(0) - s_cb(4)))
2287 ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))* &
2288 (s_cb(2) - s_cb(3)))/((s_cb(2) - s_cb(4))* &
2289 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2291 ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))* &
2292 ((s_cb(1) - s_cb(3))*(s_cb(2) - s_cb(3)) + &
2293 (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + &
2294 (s_cb(2) - s_cb(4)))))/ &
2295 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2296 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2298 ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))* &
2299 (s_cb(2) - s_cb(4)))/((s_cb(0) - s_cb(2))* &
2300 (s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2304# 470 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2305 if (cbc_dir_in == 3 .and. recon_type == weno_type)
then
2306 if (weno_order == 1)
then
2313 elseif (weno_order == 3)
then
2320 pi_coef_z(0, 0, cbc_loc_in) = (s_cb(0) - s_cb(1))/(s_cb(0) - s_cb(2))
2326 fd_coef_z(0, cbc_loc_in) = -50._wp/(25._wp*
ds(0) + 2._wp*
ds(1) &
2327 - 1.e1_wp*
ds(2) + 1.e1_wp*
ds(3) &
2335 ((s_cb(0) - s_cb(1))*(s_cb(1) - s_cb(2))* &
2336 (s_cb(1) - s_cb(3)))/((s_cb(1) - s_cb(4))* &
2337 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(2)))
2339 ((s_cb(1) - s_cb(0))*(s_cb(1) - s_cb(2))* &
2340 ((s_cb(1) - s_cb(3))*(s_cb(1) - s_cb(3)) - &
2341 (s_cb(0) - s_cb(4))*((s_cb(3) - s_cb(1)) + &
2342 (s_cb(4) - s_cb(1)))))/ &
2343 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2344 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2346 (s_cb(1) - s_cb(0))*((s_cb(1) - s_cb(2))* &
2347 (s_cb(1) - s_cb(3)) + ((s_cb(0) - s_cb(2)) + &
2348 (s_cb(1) - s_cb(3)))*(s_cb(0) - s_cb(4)))/ &
2349 ((s_cb(2) - s_cb(0))*(s_cb(0) - s_cb(3))* &
2350 (s_cb(0) - s_cb(4)))
2352 ((s_cb(0) - s_cb(2))*(s_cb(2) - s_cb(1))* &
2353 (s_cb(2) - s_cb(3)))/((s_cb(2) - s_cb(4))* &
2354 (s_cb(4) - s_cb(0))*(s_cb(4) - s_cb(1)))
2356 ((s_cb(0) - s_cb(2))*(s_cb(1) - s_cb(2))* &
2357 ((s_cb(1) - s_cb(3))*(s_cb(2) - s_cb(3)) + &
2358 (s_cb(0) - s_cb(4))*((s_cb(1) - s_cb(3)) + &
2359 (s_cb(2) - s_cb(4)))))/ &
2360 ((s_cb(0) - s_cb(3))*(s_cb(1) - s_cb(3))* &
2361 (s_cb(0) - s_cb(4))*(s_cb(1) - s_cb(4)))
2363 ((s_cb(1) - s_cb(2))*(s_cb(2) - s_cb(3))* &
2364 (s_cb(2) - s_cb(4)))/((s_cb(0) - s_cb(2))* &
2365 (s_cb(0) - s_cb(3))*(s_cb(0) - s_cb(4)))
2369# 535 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2464 subroutine s_cbc(q_prim_vf, flux_vf, flux_src_vf, &
2465 cbc_dir_norm, cbc_loc_norm, &
2468 type(scalar_field), &
2469 dimension(sys_size), &
2470 intent(in) :: q_prim_vf
2472 type(scalar_field), &
2473 dimension(sys_size), &
2474 intent(inout) :: flux_vf, flux_src_vf
2476 integer,
intent(in) :: cbc_dir_norm, cbc_loc_norm
2478 type(int_bounds_info),
intent(in) :: ix, iy, iz
2480 real(wp) :: dpres_dt
2481 real(wp) :: dgamma_dt
2482 real(wp) :: dpi_inf_dt
2484 real(wp) :: dpres_ds
2485# 643 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2486 real(wp),
dimension(sys_size) ::
l
2487# 645 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2488# 654 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2489 real(wp),
dimension(num_fluids) :: alpha_rho, dalpha_rho_ds, mf
2490 real(wp),
dimension(num_vels) :: vel, dvel_ds
2491 real(wp),
dimension(num_fluids) :: adv_local, dadv_ds
2492 real(wp),
dimension(num_fluids) :: dadv_dt
2493 real(wp),
dimension(num_dims) :: dvel_dt
2494 real(wp),
dimension(num_fluids) :: dalpha_rho_dt
2495 real(wp),
dimension(num_species) :: ys, h_k, dys_dt, dys_ds, xs, gamma_i, cp_i
2496# 662 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2497 real(wp),
dimension(2) :: re_cbc
2498 real(wp),
dimension(3) :: lambda
2509 real(wp) :: t, sum_enthalpies
2510 real(wp) :: cv, cp, e_mix, mw, r_gas
2512 real(wp) :: vel_k_sum, vel_dv_dt_sum
2514 integer :: i,
j,
k, r
2524# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2525#if defined(MFC_OpenACC)
2526# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2528# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2529#elif defined(MFC_OpenMP)
2530# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2532# 688 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2540# 696 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2541 if (
cbc_dir == 1 .and. recon_type == weno_type)
then
2544 if (weno_order == 3 .or. dummy)
then
2549 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
2552# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2554# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2555#if defined(MFC_OpenACC)
2556# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2558# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2559#elif defined(MFC_OpenMP)
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"
2570# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2583# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2585# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2586#if defined(MFC_OpenACC)
2587# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2589# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2590#elif defined(MFC_OpenMP)
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"
2599# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2603# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2605# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2606#if defined(MFC_OpenACC)
2607# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2609# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2610#elif defined(MFC_OpenMP)
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"
2621# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2634# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2636# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2637#if defined(MFC_OpenACC)
2638# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2640# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2641#elif defined(MFC_OpenMP)
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"
2650# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2655 if (weno_order == 5 .or. dummy)
then
2659 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
2662# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2664# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2665#if defined(MFC_OpenACC)
2666# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2668# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2669#elif defined(MFC_OpenMP)
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"
2680# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2701# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2703# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2704#if defined(MFC_OpenACC)
2705# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2707# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2708#elif defined(MFC_OpenMP)
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"
2717# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2721# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2723# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2724#if defined(MFC_OpenACC)
2725# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2727# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2728#elif defined(MFC_OpenMP)
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"
2739# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2760# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2762# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2763#if defined(MFC_OpenACC)
2764# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2766# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2767#elif defined(MFC_OpenMP)
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"
2776# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2783# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2785# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2786#if defined(MFC_OpenACC)
2787# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2789# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2790#elif defined(MFC_OpenMP)
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"
2801# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2808# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2809#if defined(MFC_OpenACC)
2810# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2812# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2813#elif defined(MFC_OpenMP)
2814# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2816# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2823# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2824#if defined(MFC_OpenACC)
2825# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2827# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2828#elif defined(MFC_OpenMP)
2829# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2831# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2839# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2840#if defined(MFC_OpenACC)
2841# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2843# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2844#elif defined(MFC_OpenMP)
2845# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2847# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2850 vel_k_sum = vel_k_sum + vel(i)**2._wp
2856# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2857#if defined(MFC_OpenACC)
2858# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2860# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2861#elif defined(MFC_OpenMP)
2862# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2864# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2866 do i = 1, advxe - e_idx
2870 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
2873# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2874#if defined(MFC_OpenACC)
2875# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2877# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2878#elif defined(MFC_OpenMP)
2879# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2881# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2884 mf(i) = alpha_rho(i)/rho
2889# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2890#if defined(MFC_OpenACC)
2891# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2893# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2894#elif defined(MFC_OpenMP)
2895# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2897# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2899 do i = chemxb, chemxe
2903 call get_mixture_molecular_weight(ys, mw)
2904 r_gas = gas_constant/mw
2906 call get_mixture_specific_heat_cp_mass(t, ys, cp)
2907 call get_mixture_energy_mass(t, ys, e_mix)
2908 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
2909 if (chem_params%gamma_method == 1)
then
2911 call get_mole_fractions(mw, ys, xs)
2912 call get_species_specific_heats_r(t, cp_i)
2913 gamma_i = cp_i/(cp_i - 1.0_wp)
2914 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
2915 else if (chem_params%gamma_method == 2)
then
2917 call get_mixture_specific_heat_cv_mass(t, ys, cv)
2918 gamma = 1.0_wp/(cp/cv - 1.0_wp)
2921 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
2927 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
2932# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2933#if defined(MFC_OpenACC)
2934# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2936# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2937#elif defined(MFC_OpenMP)
2938# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2940# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2943 dalpha_rho_ds(i) = 0._wp
2947# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2948#if defined(MFC_OpenACC)
2949# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2951# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2952#elif defined(MFC_OpenMP)
2953# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2955# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2963# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2964#if defined(MFC_OpenACC)
2965# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2967# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2968#elif defined(MFC_OpenMP)
2969# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2971# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2973 do i = 1, advxe - e_idx
2979# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2980#if defined(MFC_OpenACC)
2981# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2983# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2984#elif defined(MFC_OpenMP)
2985# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2987# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2989 do i = 1, num_species
2995# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2996#if defined(MFC_OpenACC)
2997# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
2999# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3000#elif defined(MFC_OpenMP)
3001# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3003# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3008# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3009#if defined(MFC_OpenACC)
3010# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3012# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3013#elif defined(MFC_OpenMP)
3014# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3016# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3024# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3025#if defined(MFC_OpenACC)
3026# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3028# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3029#elif defined(MFC_OpenMP)
3030# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3032# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3044# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3045#if defined(MFC_OpenACC)
3046# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3048# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3049#elif defined(MFC_OpenMP)
3050# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3052# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3054 do i = 1, advxe - e_idx
3062# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3063#if defined(MFC_OpenACC)
3064# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3066# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3067#elif defined(MFC_OpenMP)
3068# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3070# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3072 do i = 1, num_species
3081 lambda(1) = vel(dir_idx(1)) - c
3082 lambda(2) = vel(dir_idx(1))
3083 lambda(3) = vel(dir_idx(1)) + c
3085 ma = vel(dir_idx(1))/c
3087 if ((
cbc_loc == -1 .and.
bcxb == bc_char_slip_wall) .or. &
3088 (
cbc_loc == 1 .and.
bcxe == bc_char_slip_wall))
then
3089 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
3090 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_buffer) .or. &
3091 (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_buffer))
then
3092 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
3093 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_inflow) .or. &
3094 (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_inflow))
then
3095 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
3097 if (bc_x%grcbc_in)
then
3099# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3100#if defined(MFC_OpenACC)
3101# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3103# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3104#elif defined(MFC_OpenMP)
3105# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3107# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3113 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(1, dir_idx(2)))/
del_in(1)
3115 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(1, dir_idx(3)))/
del_in(1)
3119# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3120#if defined(MFC_OpenACC)
3121# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3123# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3124#elif defined(MFC_OpenMP)
3125# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3127# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3129 do i = e_idx, advxe - 1
3130 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, 1))/
del_in(1)
3134 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_nr_sub_outflow) .or. &
3135 (
cbc_loc == 1 .and.
bcxe == bc_char_nr_sub_outflow))
then
3136 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
3138 if (bc_x%grcbc_out)
then
3142 if (bc_x%grcbc_vel_out)
then
3143 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)
3146 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_ff_sub_outflow) .or. &
3147 (
cbc_loc == 1 .and.
bcxe == bc_char_ff_sub_outflow))
then
3148 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
3149 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_cp_sub_outflow) .or. &
3150 (
cbc_loc == 1 .and.
bcxe == bc_char_cp_sub_outflow))
then
3151 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
3152 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_sup_inflow) .or. &
3153 (
cbc_loc == 1 .and.
bcxe == bc_char_sup_inflow))
then
3154 call s_compute_supersonic_inflow_l(
l)
3155 else if ((
cbc_loc == -1 .and.
bcxb == bc_char_sup_outflow) .or. &
3156 (
cbc_loc == 1 .and.
bcxe == bc_char_sup_outflow))
then
3157 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
3162 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1)) &
3165 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
3169# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3170#if defined(MFC_OpenACC)
3171# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3173# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3174#elif defined(MFC_OpenMP)
3175# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3177# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3180 dalpha_rho_dt(i) = &
3181 -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
3185# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3186#if defined(MFC_OpenACC)
3187# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3189# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3190#elif defined(MFC_OpenMP)
3191# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3193# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3196 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))* &
3197 (
l(1) -
l(advxe))/(2._wp*rho*c) + &
3198 (dir_flg(dir_idx(i)) - 1._wp)* &
3202 vel_dv_dt_sum = 0._wp
3204# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3205#if defined(MFC_OpenACC)
3206# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3208# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3209#elif defined(MFC_OpenMP)
3210# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3212# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3215 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
3220# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3221#if defined(MFC_OpenACC)
3222# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3224# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3225#elif defined(MFC_OpenMP)
3226# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3228# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3230 do i = 1, num_species
3231 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
3238# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3239#if defined(MFC_OpenACC)
3240# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3242# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3243#elif defined(MFC_OpenMP)
3244# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3246# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3248 do i = 1, advxe - e_idx
3249 dadv_dt(i) = -
l(momxe + i)
3253# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3254#if defined(MFC_OpenACC)
3255# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3257# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3258#elif defined(MFC_OpenMP)
3259# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3261# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3263 do i = 1, advxe - e_idx
3264 dadv_dt(i) = -
l(momxe + i)
3268 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
3270 if (model_eqns == 1)
then
3271 drho_dt = dalpha_rho_dt(1)
3272 dgamma_dt = dadv_dt(1)
3273# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3274 dpi_inf_dt = dadv_dt(2)
3275# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3278# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3279#if defined(MFC_OpenACC)
3280# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3282# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3283#elif defined(MFC_OpenMP)
3284# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3286# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3288 do i = 1, num_fluids
3289 drho_dt = drho_dt + dalpha_rho_dt(i)
3290 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
3291 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
3292 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
3298# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3299#if defined(MFC_OpenACC)
3300# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3302# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3303#elif defined(MFC_OpenMP)
3304# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3306# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3310 +
ds(0)*dalpha_rho_dt(i)
3314# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3315#if defined(MFC_OpenACC)
3316# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3318# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3319#elif defined(MFC_OpenMP)
3320# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3322# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3326 +
ds(0)*(vel(i - contxe)*drho_dt &
3327 + rho*dvel_dt(i - contxe))
3332 call get_species_enthalpies_rt(t, h_k)
3333 sum_enthalpies = 0._wp
3335# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3336#if defined(MFC_OpenACC)
3337# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3339# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3340#elif defined(MFC_OpenMP)
3341# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3343# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3345 do i = 1, num_species
3347# 1064 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3348 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
3349 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
3350# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3353 +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) + sum_enthalpies)
3355# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3356#if defined(MFC_OpenACC)
3357# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3359# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3360#elif defined(MFC_OpenMP)
3361# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3363# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3365 do i = 1, num_species
3367 +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
3371 +
ds(0)*(pres*dgamma_dt &
3375 + rho*vel_dv_dt_sum &
3376 + 5.e-1_wp*drho_dt*vel_k_sum)
3379 if (riemann_solver == 1)
then
3381# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3382#if defined(MFC_OpenACC)
3383# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3385# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3386#elif defined(MFC_OpenMP)
3387# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3389# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3396# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3397#if defined(MFC_OpenACC)
3398# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3400# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3401#elif defined(MFC_OpenMP)
3402# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3404# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3408 1._wp/max(abs(vel(dir_idx(1))), sgm_eps) &
3409 *sign(1._wp, vel(dir_idx(1))) &
3413 +
ds(0)*dadv_dt(i - e_idx))
3419# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3420#if defined(MFC_OpenACC)
3421# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3423# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3424#elif defined(MFC_OpenMP)
3425# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3427# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3431 ds(0)*dadv_dt(i - e_idx)
3435# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3436#if defined(MFC_OpenACC)
3437# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3439# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3440#elif defined(MFC_OpenMP)
3441# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3443# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3455# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3457# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3458#if defined(MFC_OpenACC)
3459# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3461# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3462#elif defined(MFC_OpenMP)
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"
3471# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3474# 696 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3475 if (
cbc_dir == 2 .and. recon_type == weno_type)
then
3478 if (weno_order == 3 .or. dummy)
then
3483 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
3486# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3488# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3489#if defined(MFC_OpenACC)
3490# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3492# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3493#elif defined(MFC_OpenMP)
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"
3504# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3517# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3519# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3520#if defined(MFC_OpenACC)
3521# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3523# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3524#elif defined(MFC_OpenMP)
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"
3533# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3537# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3539# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3540#if defined(MFC_OpenACC)
3541# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3543# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3544#elif defined(MFC_OpenMP)
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"
3555# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3568# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3570# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3571#if defined(MFC_OpenACC)
3572# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3574# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3575#elif defined(MFC_OpenMP)
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"
3584# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3589 if (weno_order == 5 .or. dummy)
then
3593 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
3596# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3598# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3599#if defined(MFC_OpenACC)
3600# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3602# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3603#elif defined(MFC_OpenMP)
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"
3614# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3635# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3637# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3638#if defined(MFC_OpenACC)
3639# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3641# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3642#elif defined(MFC_OpenMP)
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"
3651# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3655# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3657# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3658#if defined(MFC_OpenACC)
3659# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3661# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3662#elif defined(MFC_OpenMP)
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"
3673# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3694# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3696# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3697#if defined(MFC_OpenACC)
3698# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3700# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3701#elif defined(MFC_OpenMP)
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"
3710# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3717# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3719# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3720#if defined(MFC_OpenACC)
3721# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3723# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3724#elif defined(MFC_OpenMP)
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"
3735# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3742# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3743#if defined(MFC_OpenACC)
3744# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3746# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3747#elif defined(MFC_OpenMP)
3748# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3750# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3757# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3758#if defined(MFC_OpenACC)
3759# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3761# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3762#elif defined(MFC_OpenMP)
3763# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3765# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3773# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3774#if defined(MFC_OpenACC)
3775# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3777# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3778#elif defined(MFC_OpenMP)
3779# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3781# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3784 vel_k_sum = vel_k_sum + vel(i)**2._wp
3790# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3791#if defined(MFC_OpenACC)
3792# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3794# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3795#elif defined(MFC_OpenMP)
3796# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3798# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3800 do i = 1, advxe - e_idx
3804 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
3807# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3808#if defined(MFC_OpenACC)
3809# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3811# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3812#elif defined(MFC_OpenMP)
3813# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3815# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3818 mf(i) = alpha_rho(i)/rho
3823# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3824#if defined(MFC_OpenACC)
3825# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3827# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3828#elif defined(MFC_OpenMP)
3829# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3831# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3833 do i = chemxb, chemxe
3837 call get_mixture_molecular_weight(ys, mw)
3838 r_gas = gas_constant/mw
3840 call get_mixture_specific_heat_cp_mass(t, ys, cp)
3841 call get_mixture_energy_mass(t, ys, e_mix)
3842 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
3843 if (chem_params%gamma_method == 1)
then
3845 call get_mole_fractions(mw, ys, xs)
3846 call get_species_specific_heats_r(t, cp_i)
3847 gamma_i = cp_i/(cp_i - 1.0_wp)
3848 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
3849 else if (chem_params%gamma_method == 2)
then
3851 call get_mixture_specific_heat_cv_mass(t, ys, cv)
3852 gamma = 1.0_wp/(cp/cv - 1.0_wp)
3855 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
3861 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
3866# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3867#if defined(MFC_OpenACC)
3868# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3870# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3871#elif defined(MFC_OpenMP)
3872# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3874# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3877 dalpha_rho_ds(i) = 0._wp
3881# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3882#if defined(MFC_OpenACC)
3883# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3885# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3886#elif defined(MFC_OpenMP)
3887# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3889# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3897# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3898#if defined(MFC_OpenACC)
3899# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3901# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3902#elif defined(MFC_OpenMP)
3903# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3905# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3907 do i = 1, advxe - e_idx
3913# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3914#if defined(MFC_OpenACC)
3915# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3917# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3918#elif defined(MFC_OpenMP)
3919# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3921# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3923 do i = 1, num_species
3929# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3930#if defined(MFC_OpenACC)
3931# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3933# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3934#elif defined(MFC_OpenMP)
3935# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3937# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3942# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3943#if defined(MFC_OpenACC)
3944# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3946# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3947#elif defined(MFC_OpenMP)
3948# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3950# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3958# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3959#if defined(MFC_OpenACC)
3960# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3962# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3963#elif defined(MFC_OpenMP)
3964# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3966# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3978# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3979#if defined(MFC_OpenACC)
3980# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3982# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3983#elif defined(MFC_OpenMP)
3984# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3986# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3988 do i = 1, advxe - e_idx
3996# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
3997#if defined(MFC_OpenACC)
3998# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4000# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4001#elif defined(MFC_OpenMP)
4002# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4004# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4006 do i = 1, num_species
4015 lambda(1) = vel(dir_idx(1)) - c
4016 lambda(2) = vel(dir_idx(1))
4017 lambda(3) = vel(dir_idx(1)) + c
4019 ma = vel(dir_idx(1))/c
4021 if ((
cbc_loc == -1 .and.
bcyb == bc_char_slip_wall) .or. &
4022 (
cbc_loc == 1 .and.
bcye == bc_char_slip_wall))
then
4023 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4024 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_buffer) .or. &
4025 (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_buffer))
then
4026 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4027 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_inflow) .or. &
4028 (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_inflow))
then
4029 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4031 if (bc_y%grcbc_in)
then
4033# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4034#if defined(MFC_OpenACC)
4035# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4037# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4038#elif defined(MFC_OpenMP)
4039# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4041# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4047 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(2, dir_idx(2)))/
del_in(2)
4049 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(2, dir_idx(3)))/
del_in(2)
4053# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4054#if defined(MFC_OpenACC)
4055# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4057# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4058#elif defined(MFC_OpenMP)
4059# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4061# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4063 do i = e_idx, advxe - 1
4064 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, 2))/
del_in(2)
4068 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_nr_sub_outflow) .or. &
4069 (
cbc_loc == 1 .and.
bcye == bc_char_nr_sub_outflow))
then
4070 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4072 if (bc_y%grcbc_out)
then
4076 if (bc_y%grcbc_vel_out)
then
4077 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)
4080 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_ff_sub_outflow) .or. &
4081 (
cbc_loc == 1 .and.
bcye == bc_char_ff_sub_outflow))
then
4082 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
4083 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_cp_sub_outflow) .or. &
4084 (
cbc_loc == 1 .and.
bcye == bc_char_cp_sub_outflow))
then
4085 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
4086 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_sup_inflow) .or. &
4087 (
cbc_loc == 1 .and.
bcye == bc_char_sup_inflow))
then
4088 call s_compute_supersonic_inflow_l(
l)
4089 else if ((
cbc_loc == -1 .and.
bcyb == bc_char_sup_outflow) .or. &
4090 (
cbc_loc == 1 .and.
bcye == bc_char_sup_outflow))
then
4091 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4096 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1)) &
4099 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
4103# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4104#if defined(MFC_OpenACC)
4105# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4107# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4108#elif defined(MFC_OpenMP)
4109# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4111# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4114 dalpha_rho_dt(i) = &
4115 -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
4119# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4120#if defined(MFC_OpenACC)
4121# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4123# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4124#elif defined(MFC_OpenMP)
4125# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4127# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4130 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))* &
4131 (
l(1) -
l(advxe))/(2._wp*rho*c) + &
4132 (dir_flg(dir_idx(i)) - 1._wp)* &
4136 vel_dv_dt_sum = 0._wp
4138# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4139#if defined(MFC_OpenACC)
4140# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4142# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4143#elif defined(MFC_OpenMP)
4144# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4146# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4149 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
4154# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4155#if defined(MFC_OpenACC)
4156# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4158# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4159#elif defined(MFC_OpenMP)
4160# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4162# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4164 do i = 1, num_species
4165 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
4172# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4173#if defined(MFC_OpenACC)
4174# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4176# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4177#elif defined(MFC_OpenMP)
4178# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4180# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4182 do i = 1, advxe - e_idx
4183 dadv_dt(i) = -
l(momxe + i)
4187# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4188#if defined(MFC_OpenACC)
4189# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4191# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4192#elif defined(MFC_OpenMP)
4193# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4195# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4197 do i = 1, advxe - e_idx
4198 dadv_dt(i) = -
l(momxe + i)
4202 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
4204 if (model_eqns == 1)
then
4205 drho_dt = dalpha_rho_dt(1)
4206 dgamma_dt = dadv_dt(1)
4207# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4208 dpi_inf_dt = dadv_dt(2)
4209# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4212# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4213#if defined(MFC_OpenACC)
4214# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4216# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4217#elif defined(MFC_OpenMP)
4218# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4220# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4222 do i = 1, num_fluids
4223 drho_dt = drho_dt + dalpha_rho_dt(i)
4224 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
4225 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
4226 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
4232# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4233#if defined(MFC_OpenACC)
4234# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4236# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4237#elif defined(MFC_OpenMP)
4238# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4240# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4244 +
ds(0)*dalpha_rho_dt(i)
4248# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4249#if defined(MFC_OpenACC)
4250# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4252# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4253#elif defined(MFC_OpenMP)
4254# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4256# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4260 +
ds(0)*(vel(i - contxe)*drho_dt &
4261 + rho*dvel_dt(i - contxe))
4266 call get_species_enthalpies_rt(t, h_k)
4267 sum_enthalpies = 0._wp
4269# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4270#if defined(MFC_OpenACC)
4271# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4273# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4274#elif defined(MFC_OpenMP)
4275# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4277# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4279 do i = 1, num_species
4281# 1064 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4282 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
4283 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
4284# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4287 +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) + sum_enthalpies)
4289# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4290#if defined(MFC_OpenACC)
4291# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4293# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4294#elif defined(MFC_OpenMP)
4295# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4297# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4299 do i = 1, num_species
4301 +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
4305 +
ds(0)*(pres*dgamma_dt &
4309 + rho*vel_dv_dt_sum &
4310 + 5.e-1_wp*drho_dt*vel_k_sum)
4313 if (riemann_solver == 1)
then
4315# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4316#if defined(MFC_OpenACC)
4317# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4319# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4320#elif defined(MFC_OpenMP)
4321# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4323# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4330# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4331#if defined(MFC_OpenACC)
4332# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4334# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4335#elif defined(MFC_OpenMP)
4336# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4338# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4342 1._wp/max(abs(vel(dir_idx(1))), sgm_eps) &
4343 *sign(1._wp, vel(dir_idx(1))) &
4347 +
ds(0)*dadv_dt(i - e_idx))
4353# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4354#if defined(MFC_OpenACC)
4355# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4357# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4358#elif defined(MFC_OpenMP)
4359# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4361# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4365 ds(0)*dadv_dt(i - e_idx)
4369# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4370#if defined(MFC_OpenACC)
4371# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4373# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4374#elif defined(MFC_OpenMP)
4375# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4377# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4389# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4391# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4392#if defined(MFC_OpenACC)
4393# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4395# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4396#elif defined(MFC_OpenMP)
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"
4405# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4408# 696 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4409 if (
cbc_dir == 3 .and. recon_type == weno_type)
then
4412 if (weno_order == 3 .or. dummy)
then
4417 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
4420# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4422# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4423#if defined(MFC_OpenACC)
4424# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4426# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4427#elif defined(MFC_OpenMP)
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"
4438# 706 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4451# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4453# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4454#if defined(MFC_OpenACC)
4455# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4457# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4458#elif defined(MFC_OpenMP)
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"
4467# 717 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4471# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4473# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4474#if defined(MFC_OpenACC)
4475# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4477# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4478#elif defined(MFC_OpenMP)
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"
4489# 719 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4502# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4504# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4505#if defined(MFC_OpenACC)
4506# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4508# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4509#elif defined(MFC_OpenMP)
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"
4518# 730 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4523 if (weno_order == 5 .or. dummy)
then
4527 is1,
is2,
is3, idwbuff(2)%beg, idwbuff(3)%beg)
4530# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4532# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4533#if defined(MFC_OpenACC)
4534# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4536# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4537#elif defined(MFC_OpenMP)
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"
4548# 740 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4569# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4571# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4572#if defined(MFC_OpenACC)
4573# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4575# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4576#elif defined(MFC_OpenMP)
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"
4585# 759 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4589# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4591# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4592#if defined(MFC_OpenACC)
4593# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4595# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4596#elif defined(MFC_OpenMP)
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"
4607# 761 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4628# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4630# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4631#if defined(MFC_OpenACC)
4632# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4634# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4635#elif defined(MFC_OpenMP)
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"
4644# 780 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4651# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4653# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4654#if defined(MFC_OpenACC)
4655# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4657# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4658#elif defined(MFC_OpenMP)
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"
4669# 785 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4676# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4677#if defined(MFC_OpenACC)
4678# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4680# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4681#elif defined(MFC_OpenMP)
4682# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4684# 790 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4691# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4692#if defined(MFC_OpenACC)
4693# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4695# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4696#elif defined(MFC_OpenMP)
4697# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4699# 795 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4707# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4708#if defined(MFC_OpenACC)
4709# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4711# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4712#elif defined(MFC_OpenMP)
4713# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4715# 801 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4718 vel_k_sum = vel_k_sum + vel(i)**2._wp
4724# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4725#if defined(MFC_OpenACC)
4726# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4728# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4729#elif defined(MFC_OpenMP)
4730# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4732# 808 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4734 do i = 1, advxe - e_idx
4738 call s_convert_species_to_mixture_variables_acc(rho, gamma, pi_inf, qv, adv_local, alpha_rho, re_cbc)
4741# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4742#if defined(MFC_OpenACC)
4743# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4745# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4746#elif defined(MFC_OpenMP)
4747# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4749# 815 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4752 mf(i) = alpha_rho(i)/rho
4757# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4758#if defined(MFC_OpenACC)
4759# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4761# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4762#elif defined(MFC_OpenMP)
4763# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4765# 821 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4767 do i = chemxb, chemxe
4771 call get_mixture_molecular_weight(ys, mw)
4772 r_gas = gas_constant/mw
4774 call get_mixture_specific_heat_cp_mass(t, ys, cp)
4775 call get_mixture_energy_mass(t, ys, e_mix)
4776 e = rho*e_mix + 5.e-1_wp*rho*vel_k_sum
4777 if (chem_params%gamma_method == 1)
then
4779 call get_mole_fractions(mw, ys, xs)
4780 call get_species_specific_heats_r(t, cp_i)
4781 gamma_i = cp_i/(cp_i - 1.0_wp)
4782 gamma = sum(xs(:)/(gamma_i(:) - 1.0_wp))
4783 else if (chem_params%gamma_method == 2)
then
4785 call get_mixture_specific_heat_cv_mass(t, ys, cv)
4786 gamma = 1.0_wp/(cp/cv - 1.0_wp)
4789 e = gamma*pres + pi_inf + 5.e-1_wp*rho*vel_k_sum
4795 call s_compute_speed_of_sound(pres, rho, gamma, pi_inf, h, adv_local, vel_k_sum, 0._wp, c, qv)
4800# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4801#if defined(MFC_OpenACC)
4802# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4804# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4805#elif defined(MFC_OpenMP)
4806# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4808# 854 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4811 dalpha_rho_ds(i) = 0._wp
4815# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4816#if defined(MFC_OpenACC)
4817# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4819# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4820#elif defined(MFC_OpenMP)
4821# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4823# 859 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4831# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4832#if defined(MFC_OpenACC)
4833# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4835# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4836#elif defined(MFC_OpenMP)
4837# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4839# 865 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4841 do i = 1, advxe - e_idx
4847# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4848#if defined(MFC_OpenACC)
4849# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4851# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4852#elif defined(MFC_OpenMP)
4853# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4855# 871 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4857 do i = 1, num_species
4863# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4864#if defined(MFC_OpenACC)
4865# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4867# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4868#elif defined(MFC_OpenMP)
4869# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4871# 877 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4876# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4877#if defined(MFC_OpenACC)
4878# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4880# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4881#elif defined(MFC_OpenMP)
4882# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4884# 880 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4892# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4893#if defined(MFC_OpenACC)
4894# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4896# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4897#elif defined(MFC_OpenMP)
4898# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4900# 886 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4912# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4913#if defined(MFC_OpenACC)
4914# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4916# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4917#elif defined(MFC_OpenMP)
4918# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4920# 896 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4922 do i = 1, advxe - e_idx
4930# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4931#if defined(MFC_OpenACC)
4932# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4934# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4935#elif defined(MFC_OpenMP)
4936# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4938# 904 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4940 do i = 1, num_species
4949 lambda(1) = vel(dir_idx(1)) - c
4950 lambda(2) = vel(dir_idx(1))
4951 lambda(3) = vel(dir_idx(1)) + c
4953 ma = vel(dir_idx(1))/c
4955 if ((
cbc_loc == -1 .and.
bczb == bc_char_slip_wall) .or. &
4956 (
cbc_loc == 1 .and.
bcze == bc_char_slip_wall))
then
4957 call s_compute_slip_wall_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4958 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_buffer) .or. &
4959 (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_buffer))
then
4960 call s_compute_nonreflecting_subsonic_buffer_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
4961 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_inflow) .or. &
4962 (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_inflow))
then
4963 call s_compute_nonreflecting_subsonic_inflow_l(lambda,
l, rho, c, dpres_ds, dvel_ds)
4965 if (bc_z%grcbc_in)
then
4967# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4968#if defined(MFC_OpenACC)
4969# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4971# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4972#elif defined(MFC_OpenMP)
4973# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4975# 931 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4981 l(momxb + 1) = c*ma*(vel(dir_idx(2)) -
vel_in(3, dir_idx(2)))/
del_in(3)
4983 l(momxb + 2) = c*ma*(vel(dir_idx(3)) -
vel_in(3, dir_idx(3)))/
del_in(3)
4987# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4988#if defined(MFC_OpenACC)
4989# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4991# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4992#elif defined(MFC_OpenMP)
4993# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4995# 941 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
4997 do i = e_idx, advxe - 1
4998 l(i) = c*ma*(adv_local(i + 1 - e_idx) -
alpha_in(i + 1 - e_idx, 3))/
del_in(3)
5002 else if ((
cbc_loc == -1 .and.
bczb == bc_char_nr_sub_outflow) .or. &
5003 (
cbc_loc == 1 .and.
bcze == bc_char_nr_sub_outflow))
then
5004 call s_compute_nonreflecting_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
5006 if (bc_z%grcbc_out)
then
5010 if (bc_z%grcbc_vel_out)
then
5011 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)
5014 else if ((
cbc_loc == -1 .and.
bczb == bc_char_ff_sub_outflow) .or. &
5015 (
cbc_loc == 1 .and.
bcze == bc_char_ff_sub_outflow))
then
5016 call s_compute_force_free_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
5017 else if ((
cbc_loc == -1 .and.
bczb == bc_char_cp_sub_outflow) .or. &
5018 (
cbc_loc == 1 .and.
bcze == bc_char_cp_sub_outflow))
then
5019 call s_compute_constant_pressure_subsonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
5020 else if ((
cbc_loc == -1 .and.
bczb == bc_char_sup_inflow) .or. &
5021 (
cbc_loc == 1 .and.
bcze == bc_char_sup_inflow))
then
5022 call s_compute_supersonic_inflow_l(
l)
5023 else if ((
cbc_loc == -1 .and.
bczb == bc_char_sup_outflow) .or. &
5024 (
cbc_loc == 1 .and.
bcze == bc_char_sup_outflow))
then
5025 call s_compute_supersonic_outflow_l(lambda,
l, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dys_ds)
5030 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1)) + rho*c*c*vel(dir_idx(1)) &
5033 dpres_dt = -5.e-1_wp*(
l(advxe) +
l(1))
5037# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5038#if defined(MFC_OpenACC)
5039# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5041# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5042#elif defined(MFC_OpenMP)
5043# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5045# 981 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5048 dalpha_rho_dt(i) = &
5049 -(
l(i + 1) - mf(i)*dpres_dt)/(c*c)
5053# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5054#if defined(MFC_OpenACC)
5055# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5057# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5058#elif defined(MFC_OpenMP)
5059# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5061# 987 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5064 dvel_dt(dir_idx(i)) = dir_flg(dir_idx(i))* &
5065 (
l(1) -
l(advxe))/(2._wp*rho*c) + &
5066 (dir_flg(dir_idx(i)) - 1._wp)* &
5070 vel_dv_dt_sum = 0._wp
5072# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5073#if defined(MFC_OpenACC)
5074# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5076# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5077#elif defined(MFC_OpenMP)
5078# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5080# 996 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5083 vel_dv_dt_sum = vel_dv_dt_sum + vel(i)*dvel_dt(i)
5088# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5089#if defined(MFC_OpenACC)
5090# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5092# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5093#elif defined(MFC_OpenMP)
5094# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5096# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5098 do i = 1, num_species
5099 dys_dt(i) = -1._wp*
l(chemxb + i - 1)
5106# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5107#if defined(MFC_OpenACC)
5108# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5110# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5111#elif defined(MFC_OpenMP)
5112# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5114# 1010 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5116 do i = 1, advxe - e_idx
5117 dadv_dt(i) = -
l(momxe + i)
5121# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5122#if defined(MFC_OpenACC)
5123# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5125# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5126#elif defined(MFC_OpenMP)
5127# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5129# 1015 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5131 do i = 1, advxe - e_idx
5132 dadv_dt(i) = -
l(momxe + i)
5136 drho_dt = 0._wp; dgamma_dt = 0._wp; dpi_inf_dt = 0._wp; dqv_dt = 0._wp
5138 if (model_eqns == 1)
then
5139 drho_dt = dalpha_rho_dt(1)
5140 dgamma_dt = dadv_dt(1)
5141# 1027 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5142 dpi_inf_dt = dadv_dt(2)
5143# 1029 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5146# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5147#if defined(MFC_OpenACC)
5148# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5150# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5151#elif defined(MFC_OpenMP)
5152# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5154# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5156 do i = 1, num_fluids
5157 drho_dt = drho_dt + dalpha_rho_dt(i)
5158 dgamma_dt = dgamma_dt + dadv_dt(i)*gammas(i)
5159 dpi_inf_dt = dpi_inf_dt + dadv_dt(i)*pi_infs(i)
5160 dqv_dt = dqv_dt + dalpha_rho_dt(i)*qvs(i)
5166# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5167#if defined(MFC_OpenACC)
5168# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5170# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5171#elif defined(MFC_OpenMP)
5172# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5174# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5178 +
ds(0)*dalpha_rho_dt(i)
5182# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5183#if defined(MFC_OpenACC)
5184# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5186# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5187#elif defined(MFC_OpenMP)
5188# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5190# 1046 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5194 +
ds(0)*(vel(i - contxe)*drho_dt &
5195 + rho*dvel_dt(i - contxe))
5200 call get_species_enthalpies_rt(t, h_k)
5201 sum_enthalpies = 0._wp
5203# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5204#if defined(MFC_OpenACC)
5205# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5207# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5208#elif defined(MFC_OpenMP)
5209# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5211# 1057 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5213 do i = 1, num_species
5215# 1064 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5216 h_k(i) = h_k(i)*gas_constant/molecular_weights(i)*t
5217 sum_enthalpies = sum_enthalpies + (rho*h_k(i) - pres*mw/molecular_weights(i)*cp/r_gas)*dys_dt(i)
5218# 1067 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5221 +
ds(0)*((e/rho + pres/rho)*drho_dt + rho*vel_dv_dt_sum + cp*t*
l(2)/(c*c) + sum_enthalpies)
5223# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5224#if defined(MFC_OpenACC)
5225# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5227# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5228#elif defined(MFC_OpenMP)
5229# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5231# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5233 do i = 1, num_species
5235 +
ds(0)*(drho_dt*ys(i) + rho*dys_dt(i))
5239 +
ds(0)*(pres*dgamma_dt &
5243 + rho*vel_dv_dt_sum &
5244 + 5.e-1_wp*drho_dt*vel_k_sum)
5247 if (riemann_solver == 1)
then
5249# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5250#if defined(MFC_OpenACC)
5251# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5253# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5254#elif defined(MFC_OpenMP)
5255# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5257# 1086 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5264# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5265#if defined(MFC_OpenACC)
5266# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5268# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5269#elif defined(MFC_OpenMP)
5270# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5272# 1091 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5276 1._wp/max(abs(vel(dir_idx(1))), sgm_eps) &
5277 *sign(1._wp, vel(dir_idx(1))) &
5281 +
ds(0)*dadv_dt(i - e_idx))
5287# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5288#if defined(MFC_OpenACC)
5289# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5291# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5292#elif defined(MFC_OpenMP)
5293# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5295# 1104 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5299 ds(0)*dadv_dt(i - e_idx)
5303# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5304#if defined(MFC_OpenACC)
5305# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5307# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5308#elif defined(MFC_OpenMP)
5309# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5311# 1110 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5323# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5325# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5326#if defined(MFC_OpenACC)
5327# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5329# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5330#elif defined(MFC_OpenMP)
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"
5339# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"
5342# 1123 "/home/runner/work/MFC/MFC/src/simulation/m_cbc.fpp"