470 & qR_prim_rsx_vf, dqR_prim_dx_vf, dqR_prim_dy_vf, dqR_prim_dz_vf, qR_prim_vf, &
471 & q_prim_vf, flux_vf, flux_src_vf, flux_gsrc_vf, norm_dir, ix, iy, iz)
473 real(wp),
dimension(idwbuff(1)%beg:,idwbuff(2)%beg:,idwbuff(3)%beg:,1:),
intent(inout) :: qL_prim_rsx_vf, qR_prim_rsx_vf
474 type(
scalar_field),
allocatable,
dimension(:),
intent(inout) :: dqL_prim_dx_vf, dqR_prim_dx_vf, dqL_prim_dy_vf, &
475 & dqR_prim_dy_vf, dqL_prim_dz_vf, dqR_prim_dz_vf
477 type(
scalar_field),
allocatable,
dimension(:),
intent(inout) :: qL_prim_vf, qR_prim_vf
478 type(
scalar_field),
dimension(sys_size),
intent(in) :: q_prim_vf
479 type(
scalar_field),
dimension(sys_size),
intent(inout) :: flux_vf, flux_src_vf, flux_gsrc_vf
480 integer,
intent(in) :: norm_dir
489# 89 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
490 real(wp),
dimension(num_fluids) :: alpha_L, alpha_R, alpha_rho_L, alpha_rho_R
491# 91 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
499 real(wp) :: S_L, S_R, s_M, S_Lstar, S_Rstar
500 real(wp) :: pTot_L, pTot_R, rhoL_star, rhoR_star
501 real(wp),
dimension(14) :: U_L, U_R
502 real(wp),
dimension(14) :: F_L, F_R, F_hlld
503 real(wp) :: us_c, uss_c
506 real(wp) :: rho_HLL, u_n_HLL_cons, tau_nn_HLL
507 real(wp) :: u_n_HLL_trace, u_t_HLL_trace
508 real(wp) :: p_face_HLL, tau_qq_face_HLL
513 real(wp) :: G_eff, G_eff_tol, C_NC, sqrtC_NC
514 real(wp) :: A_L, A_R, denomA, fac_L, fac_R
515 real(wp) :: u_n_L, u_t_L, u_n_R, u_t_R
516 real(wp) :: u_t2_L, u_t2_R
517 real(wp) :: tau_nn_L, tau_nt_L, tau_tt_L, tau_nn_R, tau_nt_R, tau_tt_R
518 real(wp) :: tau_nt2_L, tau_nt2_R, tau_t2t2_L, tau_t2t2_R, tau_t1t2_L, tau_t1t2_R
519 real(wp) :: tau_qq_L, tau_qq_R
521# 123 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
522 real(wp),
dimension(eqn_idx%stress%end - eqn_idx%stress%beg + 1) :: tau_e_L, tau_e_R
523# 125 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
525 real(wp) :: alpha1_L_star, alpha1_R_star, alpha2_L_star, alpha2_R_star
526 real(wp) :: u_t_star, tau_nt_star
527 real(wp) :: u_t2_star, tau_nt2_star
528 real(wp) :: tau_nn_L_star, tau_nn_R_star, tau_tt_L_star, tau_tt_R_star
529 real(wp) :: tau_tt_L_starstar, tau_tt_R_starstar
530 real(wp) :: tau_t2t2_L_star, tau_t2t2_R_star
531 real(wp) :: tau_t2t2_L_starstar, tau_t2t2_R_starstar
532 real(wp) :: tau_t1t2_L_star, tau_t1t2_R_star
533 real(wp) :: tau_t1t2_L_starstar, tau_t1t2_R_starstar
534 real(wp) :: tau_qq_L_star, tau_qq_R_star
535 real(wp) :: pTot_star
536 real(wp) :: E_L_star, E_R_star
537 real(wp) :: E_L_starstar, E_R_starstar
538 real(wp) :: p_face, tau_qq_face
539 real(wp) :: u_n_face, u_t_face
542 real(wp) :: tau_nn_hat, tau_nt_hat, tau_tt_hat, tau_qq_hat
543 real(wp) :: tau_nt2_hat, tau_t2t2_hat, tau_t1t2_hat
545# 153 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
546 real(wp),
dimension(num_fluids) :: alpha_hat, alpha_rho_hat
547 real(wp),
dimension(eqn_idx%stress%end - eqn_idx%stress%beg + 1) :: tau_e_hat
548# 156 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
550 real(wp) :: pres_hat, blkmod1_hat, blkmod2_hat, K_hat, alpha_hat_q, alpha_rho_hat_q
551 real(wp) :: C_hat_1, C_hat_2
552 real(wp) :: Sigma_L, Sigma_R, dSigma, Sigma_ref
553 real(wp) :: a_L_ref, a_R_ref, a_ref
554 real(wp) :: du_t, dtau_nt, du_t2, dtau_nt2
555 real(wp) :: sensor_ptot, sensor_vt, sensor_tnt, sensor_combined
557 real(wp),
parameter :: ADC_power = 1.0_wp
558 real(wp) :: alpha_L_sum, alpha_R_sum
559 logical :: degenerate, shear_degenerate, fan_fallback, shear_cond
560 integer :: i, j, k, l, ipass, zone
563 & qr_prim_rsx_vf, dqr_prim_dx_vf, dqr_prim_dy_vf, dqr_prim_dz_vf, norm_dir, ix, iy, iz)
567# 178 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
568# 179 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
569# 180 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
570 if (norm_dir == 1)
then
573# 184 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
579# 190 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
580# 191 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
581# 192 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
582# 193 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
585# 196 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
589# 213 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
590# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
592# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
594# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
595#if defined(MFC_OpenACC)
596# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
598# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
600# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
602# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
604# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
606# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
608# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
610# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
612# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
614# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
615#elif defined(MFC_OpenMP)
616# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
618# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
620# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
622# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
624# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
626# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
628# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
630# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
632# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
634# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
636# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
638# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
640# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
642# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
649 do i = 1, eqn_idx%cont%end
650 alpha_rho_l(i) = ql_prim_rsx_vf(j, k, l, i)
651 alpha_rho_r(i) = qr_prim_rsx_vf(j + 1, k, l, i)
660 vel%L(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%cont%end + i)
661 vel%R(i) = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%cont%end + i)
664 vel_rms%L = vel%L(1)**2 + vel%L(2)**2 + vel%L(3)**2
665 vel_rms%R = vel%R(1)**2 + vel%R(2)**2 + vel%R(3)**2
668 alpha_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%E + i)
669 alpha_r(i) = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%E + i)
677 alpha_rho_l(i) = max(0._wp, alpha_rho_l(i))
678 alpha_l(i) = min(max(0._wp, alpha_l(i)), 1._wp)
679 alpha_l_sum = alpha_l_sum + alpha_l(i)
681 alpha_l = alpha_l/max(alpha_l_sum, sgm_eps)
684 alpha_rho_r(i) = max(0._wp, alpha_rho_r(i))
685 alpha_r(i) = min(max(0._wp, alpha_r(i)), 1._wp)
686 alpha_r_sum = alpha_r_sum + alpha_r(i)
688 alpha_r = alpha_r/max(alpha_r_sum, sgm_eps)
691 pres%L = ql_prim_rsx_vf(j, k, l, eqn_idx%E)
692 pres%R = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%E)
696# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
697#if defined(MFC_OpenACC)
698# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
700# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
701#elif defined(MFC_OpenMP)
702# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
704# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
706 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
707 tau_e_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%stress%beg - 1 + i)
708 tau_e_r(i) = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%stress%beg - 1 + i)
711 u_t2_l = 0._wp; u_t2_r = 0._wp
712 tau_nt2_l = 0._wp; tau_nt2_r = 0._wp
713 tau_t2t2_l = 0._wp; tau_t2t2_r = 0._wp
714 tau_t1t2_l = 0._wp; tau_t1t2_r = 0._wp
721 else if (p == 0)
then
737 tau_qq_l = tau_e_l(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
738 tau_qq_r = tau_e_r(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
740 tau_qq_l = 0._wp; tau_qq_r = 0._wp
743 ptot_l = pres%L - tau_nn_l
744 ptot_r = pres%R - tau_nn_r
750 if (abs(ptot_r - ptot_l) < 1e-12_wp*max(abs(ptot_l), abs(ptot_r), 1._wp))
then
751 ptot_l = 5e-1_wp*(ptot_l + ptot_r)
758 g_l = 0._wp; g_r = 0._wp
760# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
761#if defined(MFC_OpenACC)
762# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
764# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
765#elif defined(MFC_OpenMP)
766# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
768# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
771 g_l = g_l + alpha_l(i)*
gs_rs(i)
772 g_r = g_r + alpha_r(i)*
gs_rs(i)
782# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
783#if defined(MFC_OpenACC)
784# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
786# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
787#elif defined(MFC_OpenMP)
788# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
790# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
792 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
793 shear_cond = (n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5))
803 s_l = min(u_n_l - sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)), &
804 & u_n_r - sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)))
805 s_r = max(u_n_r + sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)), &
806 & u_n_l + sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)))
808 if (p > 0 .and. .not. cyl_coord)
then
810 u_l(1) = alpha_rho_l(1); u_r(1) = alpha_rho_r(1)
811 u_l(2) = alpha_rho_l(2); u_r(2) = alpha_rho_r(2)
812 u_l(3) = rho%L*u_n_l; u_r(3) = rho%R*u_n_r
813 u_l(4) = rho%L*u_t_l; u_r(4) = rho%R*u_t_r
814 u_l(5) = rho%L*u_t2_l; u_r(5) = rho%R*u_t2_r
815 u_l(6) = e%L; u_r(6) = e%R
816 u_l(7) = alpha_l(1); u_r(7) = alpha_r(1)
817 u_l(8) = alpha_l(2); u_r(8) = alpha_r(2)
818 u_l(9) = rho%L*tau_nn_l; u_r(9) = rho%R*tau_nn_r
819 u_l(10) = rho%L*tau_nt_l; u_r(10) = rho%R*tau_nt_r
820 u_l(11) = rho%L*tau_nt2_l; u_r(11) = rho%R*tau_nt2_r
821 u_l(12) = rho%L*tau_tt_l; u_r(12) = rho%R*tau_tt_r
822 u_l(13) = rho%L*tau_t2t2_l; u_r(13) = rho%R*tau_t2t2_r
823 u_l(14) = rho%L*tau_t1t2_l; u_r(14) = rho%R*tau_t1t2_r
825 f_l(1) = u_l(1)*u_n_l; f_r(1) = u_r(1)*u_n_r
826 f_l(2) = u_l(2)*u_n_l; f_r(2) = u_r(2)*u_n_r
827 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
828 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
829 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
830 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
831 f_l(5) = rho%L*u_n_l*u_t2_l - tau_nt2_l
832 f_r(5) = rho%R*u_n_r*u_t2_r - tau_nt2_r
833 f_l(6) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l - u_t2_l*tau_nt2_l
834 f_r(6) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r - u_t2_r*tau_nt2_r
839 u_l(1) = alpha_rho_l(1)
840 u_l(2) = alpha_rho_l(2)
846 u_l(8) = rho%L*tau_nn_l
847 u_l(9) = rho%L*tau_nt_l
848 u_l(10) = rho%L*tau_tt_l
849 u_l(11) = rho%L*tau_qq_l
851 u_r(1) = alpha_rho_r(1)
852 u_r(2) = alpha_rho_r(2)
858 u_r(8) = rho%R*tau_nn_r
859 u_r(9) = rho%R*tau_nt_r
860 u_r(10) = rho%R*tau_tt_r
861 u_r(11) = rho%R*tau_qq_r
863 f_l(1) = u_l(1)*u_n_l
864 f_l(2) = u_l(2)*u_n_l
865 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
866 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
867 f_l(5) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l
869 f_r(1) = u_r(1)*u_n_r
870 f_r(2) = u_r(2)*u_n_r
871 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
872 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
873 f_r(5) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r
878 a_l = rho%L*(s_l - u_n_l)
879 a_r = rho%R*(s_r - u_n_r)
882 s_m = ((ptot_r - ptot_l) + a_l*u_n_l - a_r*u_n_r)/(a_l - a_r + verysmall)
886 degenerate = (abs(denoma) < verysmall .or. .not. (s_l - verysmall <= s_m .and. s_m <= s_r + verysmall))
888 if (.not. degenerate)
then
891 ptot_star = ptot_l + a_l*(s_m - u_n_l)
896 rhol_star = rho%L*(s_l - u_n_l)/min(s_l - s_m, -verysmall)
897 rhor_star = rho%R*(s_r - u_n_r)/max(s_r - s_m, verysmall)
898 fac_l = (s_l - u_n_l)/min(s_l - s_m, -verysmall)
899 fac_r = (s_r - u_n_r)/max(s_r - s_m, verysmall)
901 e_l_star = (e%L*(u_n_l - s_l) + u_n_l*ptot_l - s_m*ptot_star)/max(s_m - s_l, verysmall)
902 e_r_star = (e%R*(u_n_r - s_r) + u_n_r*ptot_r - s_m*ptot_star)/min(s_m - s_r, -verysmall)
904 if (riemann_hypo_adc)
then
908 dsigma = sigma_r - sigma_l
909 sigma_ref = max(max(abs(sigma_l), abs(sigma_r)), verysmall)
911 a_l_ref = sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L))
912 a_r_ref = sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R))
913 a_ref = max(max(a_l_ref, a_r_ref), verysmall)
916 dtau_nt = tau_nt_r - tau_nt_l
917 du_t2 = u_t2_r - u_t2_l
918 dtau_nt2 = tau_nt2_r - tau_nt2_l
920 sensor_ptot = (dsigma*dsigma)/((adc_kappa*sigma_ref)**2 + verysmall)
921 sensor_vt = (du_t*du_t + du_t2*du_t2)/((adc_kappa*a_ref)**2 + verysmall)
922 sensor_tnt = (dtau_nt*dtau_nt + dtau_nt2*dtau_nt2)/((adc_kappa*sigma_ref)**2 + verysmall)
924 sensor_combined = sensor_ptot + sensor_tnt + sensor_vt
926 phi = exp(-(sensor_combined**adc_power))
934# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
935#if defined(MFC_OpenACC)
936# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
938# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
939#elif defined(MFC_OpenMP)
940# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
942# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
945 do i = 1, eqn_idx%cont%end
946 alpha_rho_hat(i) = q_prim_vf(i)%sf(j + ipass - 1, k, l)
949 alpha_hat(i) = q_prim_vf(eqn_idx%E + i)%sf(j + ipass - 1, k, l)
952# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
953#if defined(MFC_OpenACC)
954# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
956# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
957#elif defined(MFC_OpenMP)
958# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
960# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
962 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
963 tau_e_hat(i) = q_prim_vf(eqn_idx%stress%beg - 1 + i)%sf(j + ipass - 1, k, l)
967 tau_nt2_hat = 0._wp; tau_t2t2_hat = 0._wp; tau_t1t2_hat = 0._wp
979 tau_qq_hat = tau_e_hat(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
986# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
987#if defined(MFC_OpenACC)
988# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
990# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
991#elif defined(MFC_OpenMP)
992# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
994# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
997 rho_hat = rho_hat + alpha_rho_hat(i)
1002# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1003#if defined(MFC_OpenACC)
1004# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1006# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1007#elif defined(MFC_OpenMP)
1008# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1010# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1012 do i = 1, num_fluids
1013 g_hat = g_hat + alpha_hat(i)*
gs_rs(i)
1018 if (alt_soundspeed)
then
1019 pres_hat = q_prim_vf(eqn_idx%E)%sf(j + ipass - 1, k, l)
1020 alpha_hat_q = alpha_hat(1)
1021 alpha_rho_hat_q = alpha_rho_hat(1)
1023 alpha_hat_q = alpha_hat(2)
1024 alpha_rho_hat_q = alpha_rho_hat(2)
1026 blkmod1_hat = blkmod1_hat + (4._wp/3._wp)*
gs_rs(1)
1027 blkmod2_hat = blkmod2_hat + (4._wp/3._wp)*
gs_rs(2)
1028 k_hat = alpha_hat(1)*alpha_hat(2)*(blkmod2_hat - blkmod1_hat)/(alpha_hat(1)*blkmod2_hat &
1029 & + alpha_hat(2)*blkmod1_hat + verysmall)
1031 c_hat_1 = alpha_hat(1) + k_hat
1032 c_hat_2 = alpha_hat(2) - k_hat
1034 if (p > 0 .and. .not. cyl_coord)
then
1036 f_l(7) = u_l(7)*u_n_l - c_hat_1*u_n_l
1037 f_r(7) = u_r(7)*u_n_r - c_hat_1*u_n_r
1038 f_l(8) = u_l(8)*u_n_l - c_hat_2*u_n_l
1039 f_r(8) = u_r(8)*u_n_r - c_hat_2*u_n_r
1040 f_l(9) = u_l(9)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
1041 f_r(9) = u_r(9)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
1042 f_l(10) = u_l(10)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
1043 f_r(10) = u_r(10)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
1044 f_l(11) = u_l(11)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t2_l
1045 f_r(11) = u_r(11)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t2_r
1046 f_l(12) = u_l(12)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
1047 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
1048 f_r(12) = u_r(12)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
1049 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
1050 f_l(13) = u_l(13)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_l &
1051 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_l
1052 f_r(13) = u_r(13)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_r &
1053 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_r
1054 f_l(14) = u_l(14)*u_n_l + rho_hat*tau_t1t2_hat*u_n_l - rho_hat*tau_nt2_hat*u_t_l &
1055 & - rho_hat*tau_nt_hat*u_t2_l
1056 f_r(14) = u_r(14)*u_n_r + rho_hat*tau_t1t2_hat*u_n_r - rho_hat*tau_nt2_hat*u_t_r &
1057 & - rho_hat*tau_nt_hat*u_t2_r
1060 f_l(6) = u_l(6)*u_n_l - c_hat_1*u_n_l
1061 f_r(6) = u_r(6)*u_n_r - c_hat_1*u_n_r
1062 f_l(7) = u_l(7)*u_n_l - c_hat_2*u_n_l
1063 f_r(7) = u_r(7)*u_n_r - c_hat_2*u_n_r
1064 f_l(8) = u_l(8)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
1065 f_r(8) = u_r(8)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
1066 f_l(9) = u_l(9)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
1067 f_r(9) = u_r(9)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
1068 f_l(10) = u_l(10)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
1069 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
1070 f_r(10) = u_r(10)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
1071 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
1072 f_l(11) = u_l(11)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_l
1073 f_r(11) = u_r(11)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_r
1080 fan_fallback = degenerate
1081 shear_degenerate = .false.
1082 if (.not. fan_fallback)
then
1083 g_eff = g_hat + tau_nn_hat
1086 if (.not. (g_eff >= -g_eff_tol))
then
1088 fan_fallback = .true.
1089 else if (g_eff <= g_eff_tol)
then
1090 shear_degenerate = .true.
1094 c_nc = rho_hat*g_eff
1095 if (.not. (c_nc > 0._wp))
then
1098 shear_degenerate = (rho_hat > 0._wp)
1099 fan_fallback = .not. shear_degenerate
1103 sqrtc_nc = sqrt(c_nc)
1104 s_lstar = s_m - sqrtc_nc/rhol_star
1105 s_rstar = s_m + sqrtc_nc/rhor_star
1107 fan_fallback = (s_lstar < s_l .or. s_rstar > s_r)
1112 if (fan_fallback)
then
1114 if (s_l < 0._wp .and. s_r > 0._wp)
then
1116 f_hlld(i) =
f_hll_flux(s_l, s_r, f_l(i), f_r(i), u_l(i), u_r(i))
1118 else if (s_l >= 0._wp)
then
1119 f_hlld(1:ncomp) = f_l(1:ncomp)
1121 f_hlld(1:ncomp) = f_r(1:ncomp)
1126 ptot_star = 5e-1_wp*(ptot_l + ptot_r)
1129 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
1130 tau_nn_l_star = tau_nn_l
1131 tau_nn_r_star = tau_nn_r
1132 tau_qq_l_star = tau_qq_l
1133 tau_qq_r_star = tau_qq_r
1137 ptot_star = ptot_l + a_l*(s_m - u_n_l)
1139 if (shear_degenerate)
then
1141 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
1142 tau_nt_star = 5e-1_wp*(tau_nt_l + tau_nt_r)
1143 u_t2_star = 5e-1_wp*(u_t2_l + u_t2_r)
1144 tau_nt2_star = 5e-1_wp*(tau_nt2_l + tau_nt2_r)
1146 u_t_star = 5e-1_wp*((tau_nt_r - tau_nt_l)/sqrtc_nc + (u_t_r + u_t_l))
1147 tau_nt_star = 5e-1_wp*((u_t_r - u_t_l)*sqrtc_nc + (tau_nt_r + tau_nt_l))
1148 u_t2_star = 5e-1_wp*((tau_nt2_r - tau_nt2_l)/sqrtc_nc + (u_t2_r + u_t2_l))
1149 tau_nt2_star = 5e-1_wp*((u_t2_r - u_t2_l)*sqrtc_nc + (tau_nt2_r + tau_nt2_l))
1155# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1156# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1157 tau_nn_l_star = tau_nn_l - (rho_hat*(g_hat*4._wp/3._wp + tau_nn_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
1158# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1159 tau_nn_r_star = tau_nn_r - (rho_hat*(g_hat*4._wp/3._wp + tau_nn_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
1160# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1161# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1162# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1163 tau_tt_l_star = tau_tt_l + (rho_hat*(g_hat*2._wp/3._wp + tau_tt_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
1164# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1165 tau_tt_r_star = tau_tt_r + (rho_hat*(g_hat*2._wp/3._wp + tau_tt_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
1166# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1167# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1168# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1169 tau_t2t2_l_star = tau_t2t2_l + (rho_hat*(g_hat*2._wp/3._wp + tau_t2t2_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
1170# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1171 tau_t2t2_r_star = tau_t2t2_r + (rho_hat*(g_hat*2._wp/3._wp + tau_t2t2_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
1172# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1173# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1174# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1175 tau_t1t2_l_star = tau_t1t2_l + (rho_hat*tau_t1t2_hat*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
1176# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1177 tau_t1t2_r_star = tau_t1t2_r + (rho_hat*tau_t1t2_hat*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
1178# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1179# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1180# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1181 tau_qq_l_star = tau_qq_l + (rho_hat*(g_hat*2._wp/3._wp + tau_qq_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
1182# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1183 tau_qq_r_star = tau_qq_r + (rho_hat*(g_hat*2._wp/3._wp + tau_qq_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
1184# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1185# 675 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1187 if (shear_degenerate)
then
1189 tau_tt_l_starstar = tau_tt_l_star
1190 tau_tt_r_starstar = tau_tt_r_star
1191 tau_t2t2_l_starstar = tau_t2t2_l_star
1192 tau_t2t2_r_starstar = tau_t2t2_r_star
1193 tau_t1t2_l_starstar = tau_t1t2_l_star
1194 tau_t1t2_r_starstar = tau_t1t2_r_star
1195 e_l_starstar = e_l_star
1196 e_r_starstar = e_r_star
1198 tau_tt_l_starstar = tau_tt_l_star + 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_l)
1199 tau_tt_r_starstar = tau_tt_r_star - 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_r)
1200 tau_t2t2_l_starstar = tau_t2t2_l_star + 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
1202 tau_t2t2_r_starstar = tau_t2t2_r_star - 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
1204 tau_t1t2_l_starstar = tau_t1t2_l_star + rho_hat*(tau_nt2_hat*(u_t_star - u_t_l) &
1205 & + tau_nt_hat*(u_t2_star - u_t2_l))/sqrtc_nc
1206 tau_t1t2_r_starstar = tau_t1t2_r_star - rho_hat*(tau_nt2_hat*(u_t_star - u_t_r) &
1207 & + tau_nt_hat*(u_t2_star - u_t2_r))/sqrtc_nc
1208 e_l_starstar = e_l_star + (rhol_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_l*tau_nt_l) &
1209 & + (u_t2_star*tau_nt2_star - u_t2_l*tau_nt2_l))
1210 e_r_starstar = e_r_star - (rhor_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_r*tau_nt_r) &
1211 & + (u_t2_star*tau_nt2_star - u_t2_r*tau_nt2_r))
1214 alpha1_l_star = (alpha_l(1)*(s_l - u_n_l) - c_hat_1*(s_m - u_n_l))/min(s_l - s_m, -verysmall)
1215 alpha1_r_star = (alpha_r(1)*(s_r - u_n_r) - c_hat_1*(s_m - u_n_r))/max(s_r - s_m, verysmall)
1217 alpha2_l_star = (alpha_l(2)*(s_l - u_n_l) - c_hat_2*(s_m - u_n_l))/min(s_l - s_m, -verysmall)
1218 alpha2_r_star = (alpha_r(2)*(s_r - u_n_r) - c_hat_2*(s_m - u_n_r))/max(s_r - s_m, verysmall)
1233 f_hlld(1:ncomp) = f_l(1:ncomp)
1234 else if (zone == 5)
then
1235 f_hlld(1:ncomp) = f_r(1:ncomp)
1236# 726 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1237 else if (zone == 1 .or. zone == 2)
then
1242 if (p > 0 .and. .not. cyl_coord)
then
1243# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1246 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
1247 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
1249# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1252 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
1253 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
1255# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1256 us_c = rhol_star*s_m
1258 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
1259 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
1261# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1262 us_c = rhol_star*u_t_l
1263 uss_c = rhol_star*u_t_star
1264 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
1265 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
1267# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1268 us_c = rhol_star*u_t2_l
1269 uss_c = rhol_star*u_t2_star
1270 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
1271 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
1273# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1275 uss_c = e_l_starstar
1276 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
1277 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
1279# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1280 us_c = alpha1_l_star
1282 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
1283 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
1285# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1286 us_c = alpha2_l_star
1288 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
1289 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
1291# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1292 us_c = rhol_star*tau_nn_l_star
1294 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
1295 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
1297# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1298 us_c = rhol_star*tau_nt_l
1299 uss_c = rhol_star*tau_nt_star
1300 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
1301 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
1303# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1304 us_c = rhol_star*tau_nt2_l
1305 uss_c = rhol_star*tau_nt2_star
1306 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
1307 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
1309# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1310 us_c = rhol_star*tau_tt_l_star
1311 uss_c = rhol_star*tau_tt_l_starstar
1312 f_hlld(12) = f_l(12) + s_l*(us_c - u_l(12))
1313 if (zone == 2) f_hlld(12) = f_hlld(12) + s_lstar*(uss_c &
1315# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1316 us_c = rhol_star*tau_t2t2_l_star
1317 uss_c = rhol_star*tau_t2t2_l_starstar
1318 f_hlld(13) = f_l(13) + s_l*(us_c - u_l(13))
1319 if (zone == 2) f_hlld(13) = f_hlld(13) + s_lstar*(uss_c &
1321# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1322 us_c = rhol_star*tau_t1t2_l_star
1323 uss_c = rhol_star*tau_t1t2_l_starstar
1324 f_hlld(14) = f_l(14) + s_l*(us_c - u_l(14))
1325 if (zone == 2) f_hlld(14) = f_hlld(14) + s_lstar*(uss_c &
1327# 739 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1329# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1332 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
1333 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
1335# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1338 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
1339 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
1341# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1342 us_c = rhol_star*s_m
1344 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
1345 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
1347# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1348 us_c = rhol_star*u_t_l
1349 uss_c = rhol_star*u_t_star
1350 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
1351 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
1353# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1355 uss_c = e_l_starstar
1356 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
1357 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
1359# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1360 us_c = alpha1_l_star
1362 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
1363 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
1365# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1366 us_c = alpha2_l_star
1368 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
1369 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
1371# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1372 us_c = rhol_star*tau_nn_l_star
1374 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
1375 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
1377# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1378 us_c = rhol_star*tau_nt_l
1379 uss_c = rhol_star*tau_nt_star
1380 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
1381 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
1383# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1384 us_c = rhol_star*tau_tt_l_star
1385 uss_c = rhol_star*tau_tt_l_starstar
1386 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
1387 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
1389# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1390 us_c = rhol_star*tau_qq_l_star
1392 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
1393 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
1395# 747 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1397# 726 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1398 else if (zone == 3 .or. zone == 4)
then
1403 if (p > 0 .and. .not. cyl_coord)
then
1404# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1407 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
1408 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
1410# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1413 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
1414 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
1416# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1417 us_c = rhor_star*s_m
1419 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
1420 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
1422# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1423 us_c = rhor_star*u_t_r
1424 uss_c = rhor_star*u_t_star
1425 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
1426 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
1428# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1429 us_c = rhor_star*u_t2_r
1430 uss_c = rhor_star*u_t2_star
1431 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
1432 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
1434# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1436 uss_c = e_r_starstar
1437 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
1438 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
1440# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1441 us_c = alpha1_r_star
1443 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
1444 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
1446# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1447 us_c = alpha2_r_star
1449 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
1450 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
1452# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1453 us_c = rhor_star*tau_nn_r_star
1455 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
1456 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
1458# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1459 us_c = rhor_star*tau_nt_r
1460 uss_c = rhor_star*tau_nt_star
1461 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
1462 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
1464# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1465 us_c = rhor_star*tau_nt2_r
1466 uss_c = rhor_star*tau_nt2_star
1467 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
1468 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
1470# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1471 us_c = rhor_star*tau_tt_r_star
1472 uss_c = rhor_star*tau_tt_r_starstar
1473 f_hlld(12) = f_r(12) + s_r*(us_c - u_r(12))
1474 if (zone == 3) f_hlld(12) = f_hlld(12) + s_rstar*(uss_c &
1476# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1477 us_c = rhor_star*tau_t2t2_r_star
1478 uss_c = rhor_star*tau_t2t2_r_starstar
1479 f_hlld(13) = f_r(13) + s_r*(us_c - u_r(13))
1480 if (zone == 3) f_hlld(13) = f_hlld(13) + s_rstar*(uss_c &
1482# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1483 us_c = rhor_star*tau_t1t2_r_star
1484 uss_c = rhor_star*tau_t1t2_r_starstar
1485 f_hlld(14) = f_r(14) + s_r*(us_c - u_r(14))
1486 if (zone == 3) f_hlld(14) = f_hlld(14) + s_rstar*(uss_c &
1488# 739 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1490# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1493 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
1494 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
1496# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1499 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
1500 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
1502# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1503 us_c = rhor_star*s_m
1505 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
1506 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
1508# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1509 us_c = rhor_star*u_t_r
1510 uss_c = rhor_star*u_t_star
1511 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
1512 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
1514# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1516 uss_c = e_r_starstar
1517 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
1518 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
1520# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1521 us_c = alpha1_r_star
1523 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
1524 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
1526# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1527 us_c = alpha2_r_star
1529 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
1530 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
1532# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1533 us_c = rhor_star*tau_nn_r_star
1535 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
1536 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
1538# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1539 us_c = rhor_star*tau_nt_r
1540 uss_c = rhor_star*tau_nt_star
1541 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
1542 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
1544# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1545 us_c = rhor_star*tau_tt_r_star
1546 uss_c = rhor_star*tau_tt_r_starstar
1547 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
1548 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
1550# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1551 us_c = rhor_star*tau_qq_r_star
1553 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
1554 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
1556# 747 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1558# 749 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1563 if (riemann_hypo_adc)
then
1569 if (s_l < 0._wp .and. s_r > 0._wp)
then
1571 f_hll_c =
f_hll_flux(s_l, s_r, f_l(i), f_r(i), u_l(i), u_r(i))
1572 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
1577 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
1587 if (cyl_coord .and. (riemann_hypo_adc .or. fan_fallback))
then
1588 if (0._wp <= s_l)
then
1589 u_n_hll_trace = u_n_l; u_t_hll_trace = u_t_l
1590 p_face_hll = pres%L; tau_qq_face_hll = tau_qq_l
1591 else if (s_r <= 0._wp)
then
1592 u_n_hll_trace = u_n_r; u_t_hll_trace = u_t_r
1593 p_face_hll = pres%R; tau_qq_face_hll = tau_qq_r
1595 u_n_hll_trace = (s_r*u_n_l - s_l*u_n_r)/(s_r - s_l + verysmall)
1596 u_t_hll_trace = (s_r*u_t_l - s_l*u_t_r)/(s_r - s_l + verysmall)
1598 u_hll_c = (s_r*u_r(1) - s_l*u_l(1) - (f_r(1) - f_l(1)))/(s_r - s_l + verysmall)
1600 u_hll_c = (s_r*u_r(2) - s_l*u_l(2) - (f_r(2) - f_l(2)))/(s_r - s_l + verysmall)
1601 rho_hll = rho_hll + u_hll_c
1602 u_hll_c = (s_r*u_r(3) - s_l*u_l(3) - (f_r(3) - f_l(3)))/(s_r - s_l + verysmall)
1603 u_n_hll_cons = u_hll_c/(rho_hll + verysmall)
1604 u_hll_c = (s_r*u_r(8) - s_l*u_l(8) - (f_r(8) - f_l(8)))/(s_r - s_l + verysmall)
1605 tau_nn_hll = u_hll_c/(rho_hll + verysmall)
1606 u_hll_c = (s_r*u_r(11) - s_l*u_l(11) - (f_r(11) - f_l(11)))/(s_r - s_l + verysmall)
1607 tau_qq_face_hll = u_hll_c/(rho_hll + verysmall)
1609 f_hll_c =
f_hll_flux(s_l, s_r, f_l(3), f_r(3), u_l(3), u_r(3))
1610 p_face_hll = f_hll_c - rho_hll*u_n_hll_cons*u_n_hll_cons + tau_nn_hll
1615# 806 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1616 if (ipass == 1)
then
1617 if (p > 0 .and. .not. cyl_coord)
then
1629 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(9)
1630 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
1631 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(12)
1632 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(11)
1633 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(14)
1634 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
1635 else if (
dir_idx(1) == 2)
then
1636 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
1637 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
1638 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(9)
1639 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(14)
1640 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
1641 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
1643 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
1644 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(14)
1645 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(13)
1646 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(10)
1647 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
1648 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(9)
1659 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(9)
1661 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(8)
1662 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(10)
1664 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(10)
1665 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(8)
1668 flux_rsx_vf(j, k, l, eqn_idx%stress%end) = f_hlld(11)
1672# 806 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1673 if (ipass == 2)
then
1674 if (p > 0 .and. .not. cyl_coord)
then
1692 else if (
dir_idx(1) == 2)
then
1729# 863 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1733 if (fan_fallback)
then
1736 u_n_face = u_n_hll_trace; u_t_face = u_t_hll_trace
1742 u_n_face = u_n_l; u_t_face = u_t_l
1743 else if (zone == 1)
then
1744 u_n_face = s_m; u_t_face = u_t_l
1745 else if (zone <= 3)
then
1746 u_n_face = s_m; u_t_face = u_t_star
1747 else if (zone == 4)
then
1748 u_n_face = s_m; u_t_face = u_t_r
1750 u_n_face = u_n_r; u_t_face = u_t_r
1754 if (riemann_hypo_adc)
then
1755 u_n_face = u_n_hll_trace + phi*(u_n_face - u_n_hll_trace)
1756 u_t_face = u_t_hll_trace + phi*(u_t_face - u_t_hll_trace)
1759# 893 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1760 if (ipass == 1)
then
1769# 893 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1770 if (ipass == 2)
then
1779# 903 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1788# 954 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1797# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1798#if defined(MFC_OpenACC)
1799# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1801# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1802#elif defined(MFC_OpenMP)
1803# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1805# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1807# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1810# 178 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1811# 179 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1812# 180 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1813 if (norm_dir == 2)
then
1816# 184 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1822# 190 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1823# 191 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1824# 192 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1825# 193 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1828# 196 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1832# 213 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1833# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1835# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1837# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1838#if defined(MFC_OpenACC)
1839# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1841# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1843# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1845# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1847# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1849# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1851# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1853# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1855# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1857# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1858#elif defined(MFC_OpenMP)
1859# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1861# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1863# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1865# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1867# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1869# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1871# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1873# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1875# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1877# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1879# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1881# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1883# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1885# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1892 do i = 1, eqn_idx%cont%end
1893 alpha_rho_l(i) = ql_prim_rsx_vf(j, k, l, i)
1894 alpha_rho_r(i) = qr_prim_rsx_vf(j, k + 1, l, i)
1903 vel%L(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%cont%end + i)
1904 vel%R(i) = qr_prim_rsx_vf(j, k + 1, l, eqn_idx%cont%end + i)
1907 vel_rms%L = vel%L(1)**2 + vel%L(2)**2 + vel%L(3)**2
1908 vel_rms%R = vel%R(1)**2 + vel%R(2)**2 + vel%R(3)**2
1910 do i = 1, num_fluids
1911 alpha_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%E + i)
1912 alpha_r(i) = qr_prim_rsx_vf(j, k + 1, l, eqn_idx%E + i)
1919 do i = 1, num_fluids
1920 alpha_rho_l(i) = max(0._wp, alpha_rho_l(i))
1921 alpha_l(i) = min(max(0._wp, alpha_l(i)), 1._wp)
1922 alpha_l_sum = alpha_l_sum + alpha_l(i)
1924 alpha_l = alpha_l/max(alpha_l_sum, sgm_eps)
1926 do i = 1, num_fluids
1927 alpha_rho_r(i) = max(0._wp, alpha_rho_r(i))
1928 alpha_r(i) = min(max(0._wp, alpha_r(i)), 1._wp)
1929 alpha_r_sum = alpha_r_sum + alpha_r(i)
1931 alpha_r = alpha_r/max(alpha_r_sum, sgm_eps)
1934 pres%L = ql_prim_rsx_vf(j, k, l, eqn_idx%E)
1935 pres%R = qr_prim_rsx_vf(j, k + 1, l, eqn_idx%E)
1939# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1940#if defined(MFC_OpenACC)
1941# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1943# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1944#elif defined(MFC_OpenMP)
1945# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1947# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1949 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
1950 tau_e_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%stress%beg - 1 + i)
1951 tau_e_r(i) = qr_prim_rsx_vf(j, k + 1, l, eqn_idx%stress%beg - 1 + i)
1954 u_t2_l = 0._wp; u_t2_r = 0._wp
1955 tau_nt2_l = 0._wp; tau_nt2_r = 0._wp
1956 tau_t2t2_l = 0._wp; tau_t2t2_r = 0._wp
1957 tau_t1t2_l = 0._wp; tau_t1t2_r = 0._wp
1964 else if (p == 0)
then
1980 tau_qq_l = tau_e_l(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
1981 tau_qq_r = tau_e_r(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
1983 tau_qq_l = 0._wp; tau_qq_r = 0._wp
1986 ptot_l = pres%L - tau_nn_l
1987 ptot_r = pres%R - tau_nn_r
1993 if (abs(ptot_r - ptot_l) < 1e-12_wp*max(abs(ptot_l), abs(ptot_r), 1._wp))
then
1994 ptot_l = 5e-1_wp*(ptot_l + ptot_r)
2001 g_l = 0._wp; g_r = 0._wp
2003# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2004#if defined(MFC_OpenACC)
2005# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2007# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2008#elif defined(MFC_OpenMP)
2009# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2011# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2013 do i = 1, num_fluids
2014 g_l = g_l + alpha_l(i)*
gs_rs(i)
2015 g_r = g_r + alpha_r(i)*
gs_rs(i)
2025# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2026#if defined(MFC_OpenACC)
2027# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2029# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2030#elif defined(MFC_OpenMP)
2031# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2033# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2035 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
2036 shear_cond = (n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5))
2046 s_l = min(u_n_l - sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)), &
2047 & u_n_r - sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)))
2048 s_r = max(u_n_r + sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)), &
2049 & u_n_l + sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)))
2051 if (p > 0 .and. .not. cyl_coord)
then
2053 u_l(1) = alpha_rho_l(1); u_r(1) = alpha_rho_r(1)
2054 u_l(2) = alpha_rho_l(2); u_r(2) = alpha_rho_r(2)
2055 u_l(3) = rho%L*u_n_l; u_r(3) = rho%R*u_n_r
2056 u_l(4) = rho%L*u_t_l; u_r(4) = rho%R*u_t_r
2057 u_l(5) = rho%L*u_t2_l; u_r(5) = rho%R*u_t2_r
2058 u_l(6) = e%L; u_r(6) = e%R
2059 u_l(7) = alpha_l(1); u_r(7) = alpha_r(1)
2060 u_l(8) = alpha_l(2); u_r(8) = alpha_r(2)
2061 u_l(9) = rho%L*tau_nn_l; u_r(9) = rho%R*tau_nn_r
2062 u_l(10) = rho%L*tau_nt_l; u_r(10) = rho%R*tau_nt_r
2063 u_l(11) = rho%L*tau_nt2_l; u_r(11) = rho%R*tau_nt2_r
2064 u_l(12) = rho%L*tau_tt_l; u_r(12) = rho%R*tau_tt_r
2065 u_l(13) = rho%L*tau_t2t2_l; u_r(13) = rho%R*tau_t2t2_r
2066 u_l(14) = rho%L*tau_t1t2_l; u_r(14) = rho%R*tau_t1t2_r
2068 f_l(1) = u_l(1)*u_n_l; f_r(1) = u_r(1)*u_n_r
2069 f_l(2) = u_l(2)*u_n_l; f_r(2) = u_r(2)*u_n_r
2070 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
2071 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
2072 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
2073 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
2074 f_l(5) = rho%L*u_n_l*u_t2_l - tau_nt2_l
2075 f_r(5) = rho%R*u_n_r*u_t2_r - tau_nt2_r
2076 f_l(6) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l - u_t2_l*tau_nt2_l
2077 f_r(6) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r - u_t2_r*tau_nt2_r
2082 u_l(1) = alpha_rho_l(1)
2083 u_l(2) = alpha_rho_l(2)
2084 u_l(3) = rho%L*u_n_l
2085 u_l(4) = rho%L*u_t_l
2089 u_l(8) = rho%L*tau_nn_l
2090 u_l(9) = rho%L*tau_nt_l
2091 u_l(10) = rho%L*tau_tt_l
2092 u_l(11) = rho%L*tau_qq_l
2094 u_r(1) = alpha_rho_r(1)
2095 u_r(2) = alpha_rho_r(2)
2096 u_r(3) = rho%R*u_n_r
2097 u_r(4) = rho%R*u_t_r
2101 u_r(8) = rho%R*tau_nn_r
2102 u_r(9) = rho%R*tau_nt_r
2103 u_r(10) = rho%R*tau_tt_r
2104 u_r(11) = rho%R*tau_qq_r
2106 f_l(1) = u_l(1)*u_n_l
2107 f_l(2) = u_l(2)*u_n_l
2108 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
2109 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
2110 f_l(5) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l
2112 f_r(1) = u_r(1)*u_n_r
2113 f_r(2) = u_r(2)*u_n_r
2114 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
2115 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
2116 f_r(5) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r
2121 a_l = rho%L*(s_l - u_n_l)
2122 a_r = rho%R*(s_r - u_n_r)
2123 denoma = (a_r - a_l)
2125 s_m = ((ptot_r - ptot_l) + a_l*u_n_l - a_r*u_n_r)/(a_l - a_r + verysmall)
2129 degenerate = (abs(denoma) < verysmall .or. .not. (s_l - verysmall <= s_m .and. s_m <= s_r + verysmall))
2131 if (.not. degenerate)
then
2134 ptot_star = ptot_l + a_l*(s_m - u_n_l)
2139 rhol_star = rho%L*(s_l - u_n_l)/min(s_l - s_m, -verysmall)
2140 rhor_star = rho%R*(s_r - u_n_r)/max(s_r - s_m, verysmall)
2141 fac_l = (s_l - u_n_l)/min(s_l - s_m, -verysmall)
2142 fac_r = (s_r - u_n_r)/max(s_r - s_m, verysmall)
2144 e_l_star = (e%L*(u_n_l - s_l) + u_n_l*ptot_l - s_m*ptot_star)/max(s_m - s_l, verysmall)
2145 e_r_star = (e%R*(u_n_r - s_r) + u_n_r*ptot_r - s_m*ptot_star)/min(s_m - s_r, -verysmall)
2147 if (riemann_hypo_adc)
then
2151 dsigma = sigma_r - sigma_l
2152 sigma_ref = max(max(abs(sigma_l), abs(sigma_r)), verysmall)
2154 a_l_ref = sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L))
2155 a_r_ref = sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R))
2156 a_ref = max(max(a_l_ref, a_r_ref), verysmall)
2158 du_t = u_t_r - u_t_l
2159 dtau_nt = tau_nt_r - tau_nt_l
2160 du_t2 = u_t2_r - u_t2_l
2161 dtau_nt2 = tau_nt2_r - tau_nt2_l
2163 sensor_ptot = (dsigma*dsigma)/((adc_kappa*sigma_ref)**2 + verysmall)
2164 sensor_vt = (du_t*du_t + du_t2*du_t2)/((adc_kappa*a_ref)**2 + verysmall)
2165 sensor_tnt = (dtau_nt*dtau_nt + dtau_nt2*dtau_nt2)/((adc_kappa*sigma_ref)**2 + verysmall)
2167 sensor_combined = sensor_ptot + sensor_tnt + sensor_vt
2169 phi = exp(-(sensor_combined**adc_power))
2177# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2178#if defined(MFC_OpenACC)
2179# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2181# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2182#elif defined(MFC_OpenMP)
2183# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2185# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2188 do i = 1, eqn_idx%cont%end
2189 alpha_rho_hat(i) = q_prim_vf(i)%sf(j, k + ipass - 1, l)
2191 do i = 1, num_fluids
2192 alpha_hat(i) = q_prim_vf(eqn_idx%E + i)%sf(j, k + ipass - 1, l)
2195# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2196#if defined(MFC_OpenACC)
2197# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2199# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2200#elif defined(MFC_OpenMP)
2201# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2203# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2205 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
2206 tau_e_hat(i) = q_prim_vf(eqn_idx%stress%beg - 1 + i)%sf(j, k + ipass - 1, l)
2210 tau_nt2_hat = 0._wp; tau_t2t2_hat = 0._wp; tau_t1t2_hat = 0._wp
2222 tau_qq_hat = tau_e_hat(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
2229# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2230#if defined(MFC_OpenACC)
2231# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2233# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2234#elif defined(MFC_OpenMP)
2235# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2237# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2239 do i = 1, num_fluids
2240 rho_hat = rho_hat + alpha_rho_hat(i)
2245# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2246#if defined(MFC_OpenACC)
2247# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2249# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2250#elif defined(MFC_OpenMP)
2251# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2253# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2255 do i = 1, num_fluids
2256 g_hat = g_hat + alpha_hat(i)*
gs_rs(i)
2261 if (alt_soundspeed)
then
2262 pres_hat = q_prim_vf(eqn_idx%E)%sf(j, k + ipass - 1, l)
2263 alpha_hat_q = alpha_hat(1)
2264 alpha_rho_hat_q = alpha_rho_hat(1)
2266 alpha_hat_q = alpha_hat(2)
2267 alpha_rho_hat_q = alpha_rho_hat(2)
2269 blkmod1_hat = blkmod1_hat + (4._wp/3._wp)*
gs_rs(1)
2270 blkmod2_hat = blkmod2_hat + (4._wp/3._wp)*
gs_rs(2)
2271 k_hat = alpha_hat(1)*alpha_hat(2)*(blkmod2_hat - blkmod1_hat)/(alpha_hat(1)*blkmod2_hat &
2272 & + alpha_hat(2)*blkmod1_hat + verysmall)
2274 c_hat_1 = alpha_hat(1) + k_hat
2275 c_hat_2 = alpha_hat(2) - k_hat
2277 if (p > 0 .and. .not. cyl_coord)
then
2279 f_l(7) = u_l(7)*u_n_l - c_hat_1*u_n_l
2280 f_r(7) = u_r(7)*u_n_r - c_hat_1*u_n_r
2281 f_l(8) = u_l(8)*u_n_l - c_hat_2*u_n_l
2282 f_r(8) = u_r(8)*u_n_r - c_hat_2*u_n_r
2283 f_l(9) = u_l(9)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
2284 f_r(9) = u_r(9)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
2285 f_l(10) = u_l(10)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
2286 f_r(10) = u_r(10)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
2287 f_l(11) = u_l(11)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t2_l
2288 f_r(11) = u_r(11)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t2_r
2289 f_l(12) = u_l(12)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
2290 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
2291 f_r(12) = u_r(12)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
2292 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
2293 f_l(13) = u_l(13)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_l &
2294 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_l
2295 f_r(13) = u_r(13)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_r &
2296 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_r
2297 f_l(14) = u_l(14)*u_n_l + rho_hat*tau_t1t2_hat*u_n_l - rho_hat*tau_nt2_hat*u_t_l &
2298 & - rho_hat*tau_nt_hat*u_t2_l
2299 f_r(14) = u_r(14)*u_n_r + rho_hat*tau_t1t2_hat*u_n_r - rho_hat*tau_nt2_hat*u_t_r &
2300 & - rho_hat*tau_nt_hat*u_t2_r
2303 f_l(6) = u_l(6)*u_n_l - c_hat_1*u_n_l
2304 f_r(6) = u_r(6)*u_n_r - c_hat_1*u_n_r
2305 f_l(7) = u_l(7)*u_n_l - c_hat_2*u_n_l
2306 f_r(7) = u_r(7)*u_n_r - c_hat_2*u_n_r
2307 f_l(8) = u_l(8)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
2308 f_r(8) = u_r(8)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
2309 f_l(9) = u_l(9)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
2310 f_r(9) = u_r(9)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
2311 f_l(10) = u_l(10)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
2312 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
2313 f_r(10) = u_r(10)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
2314 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
2315 f_l(11) = u_l(11)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_l
2316 f_r(11) = u_r(11)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_r
2323 fan_fallback = degenerate
2324 shear_degenerate = .false.
2325 if (.not. fan_fallback)
then
2326 g_eff = g_hat + tau_nn_hat
2329 if (.not. (g_eff >= -g_eff_tol))
then
2331 fan_fallback = .true.
2332 else if (g_eff <= g_eff_tol)
then
2333 shear_degenerate = .true.
2337 c_nc = rho_hat*g_eff
2338 if (.not. (c_nc > 0._wp))
then
2341 shear_degenerate = (rho_hat > 0._wp)
2342 fan_fallback = .not. shear_degenerate
2346 sqrtc_nc = sqrt(c_nc)
2347 s_lstar = s_m - sqrtc_nc/rhol_star
2348 s_rstar = s_m + sqrtc_nc/rhor_star
2350 fan_fallback = (s_lstar < s_l .or. s_rstar > s_r)
2355 if (fan_fallback)
then
2357 if (s_l < 0._wp .and. s_r > 0._wp)
then
2359 f_hlld(i) =
f_hll_flux(s_l, s_r, f_l(i), f_r(i), u_l(i), u_r(i))
2361 else if (s_l >= 0._wp)
then
2362 f_hlld(1:ncomp) = f_l(1:ncomp)
2364 f_hlld(1:ncomp) = f_r(1:ncomp)
2369 ptot_star = 5e-1_wp*(ptot_l + ptot_r)
2372 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
2373 tau_nn_l_star = tau_nn_l
2374 tau_nn_r_star = tau_nn_r
2375 tau_qq_l_star = tau_qq_l
2376 tau_qq_r_star = tau_qq_r
2380 ptot_star = ptot_l + a_l*(s_m - u_n_l)
2382 if (shear_degenerate)
then
2384 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
2385 tau_nt_star = 5e-1_wp*(tau_nt_l + tau_nt_r)
2386 u_t2_star = 5e-1_wp*(u_t2_l + u_t2_r)
2387 tau_nt2_star = 5e-1_wp*(tau_nt2_l + tau_nt2_r)
2389 u_t_star = 5e-1_wp*((tau_nt_r - tau_nt_l)/sqrtc_nc + (u_t_r + u_t_l))
2390 tau_nt_star = 5e-1_wp*((u_t_r - u_t_l)*sqrtc_nc + (tau_nt_r + tau_nt_l))
2391 u_t2_star = 5e-1_wp*((tau_nt2_r - tau_nt2_l)/sqrtc_nc + (u_t2_r + u_t2_l))
2392 tau_nt2_star = 5e-1_wp*((u_t2_r - u_t2_l)*sqrtc_nc + (tau_nt2_r + tau_nt2_l))
2398# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2399# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2400 tau_nn_l_star = tau_nn_l - (rho_hat*(g_hat*4._wp/3._wp + tau_nn_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
2401# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2402 tau_nn_r_star = tau_nn_r - (rho_hat*(g_hat*4._wp/3._wp + tau_nn_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
2403# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2404# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2405# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2406 tau_tt_l_star = tau_tt_l + (rho_hat*(g_hat*2._wp/3._wp + tau_tt_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
2407# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2408 tau_tt_r_star = tau_tt_r + (rho_hat*(g_hat*2._wp/3._wp + tau_tt_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
2409# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2410# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2411# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2412 tau_t2t2_l_star = tau_t2t2_l + (rho_hat*(g_hat*2._wp/3._wp + tau_t2t2_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
2413# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2414 tau_t2t2_r_star = tau_t2t2_r + (rho_hat*(g_hat*2._wp/3._wp + tau_t2t2_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
2415# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2416# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2417# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2418 tau_t1t2_l_star = tau_t1t2_l + (rho_hat*tau_t1t2_hat*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
2419# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2420 tau_t1t2_r_star = tau_t1t2_r + (rho_hat*tau_t1t2_hat*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
2421# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2422# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2423# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2424 tau_qq_l_star = tau_qq_l + (rho_hat*(g_hat*2._wp/3._wp + tau_qq_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
2425# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2426 tau_qq_r_star = tau_qq_r + (rho_hat*(g_hat*2._wp/3._wp + tau_qq_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
2427# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2428# 675 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2430 if (shear_degenerate)
then
2432 tau_tt_l_starstar = tau_tt_l_star
2433 tau_tt_r_starstar = tau_tt_r_star
2434 tau_t2t2_l_starstar = tau_t2t2_l_star
2435 tau_t2t2_r_starstar = tau_t2t2_r_star
2436 tau_t1t2_l_starstar = tau_t1t2_l_star
2437 tau_t1t2_r_starstar = tau_t1t2_r_star
2438 e_l_starstar = e_l_star
2439 e_r_starstar = e_r_star
2441 tau_tt_l_starstar = tau_tt_l_star + 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_l)
2442 tau_tt_r_starstar = tau_tt_r_star - 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_r)
2443 tau_t2t2_l_starstar = tau_t2t2_l_star + 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
2445 tau_t2t2_r_starstar = tau_t2t2_r_star - 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
2447 tau_t1t2_l_starstar = tau_t1t2_l_star + rho_hat*(tau_nt2_hat*(u_t_star - u_t_l) &
2448 & + tau_nt_hat*(u_t2_star - u_t2_l))/sqrtc_nc
2449 tau_t1t2_r_starstar = tau_t1t2_r_star - rho_hat*(tau_nt2_hat*(u_t_star - u_t_r) &
2450 & + tau_nt_hat*(u_t2_star - u_t2_r))/sqrtc_nc
2451 e_l_starstar = e_l_star + (rhol_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_l*tau_nt_l) &
2452 & + (u_t2_star*tau_nt2_star - u_t2_l*tau_nt2_l))
2453 e_r_starstar = e_r_star - (rhor_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_r*tau_nt_r) &
2454 & + (u_t2_star*tau_nt2_star - u_t2_r*tau_nt2_r))
2457 alpha1_l_star = (alpha_l(1)*(s_l - u_n_l) - c_hat_1*(s_m - u_n_l))/min(s_l - s_m, -verysmall)
2458 alpha1_r_star = (alpha_r(1)*(s_r - u_n_r) - c_hat_1*(s_m - u_n_r))/max(s_r - s_m, verysmall)
2460 alpha2_l_star = (alpha_l(2)*(s_l - u_n_l) - c_hat_2*(s_m - u_n_l))/min(s_l - s_m, -verysmall)
2461 alpha2_r_star = (alpha_r(2)*(s_r - u_n_r) - c_hat_2*(s_m - u_n_r))/max(s_r - s_m, verysmall)
2476 f_hlld(1:ncomp) = f_l(1:ncomp)
2477 else if (zone == 5)
then
2478 f_hlld(1:ncomp) = f_r(1:ncomp)
2479# 726 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2480 else if (zone == 1 .or. zone == 2)
then
2485 if (p > 0 .and. .not. cyl_coord)
then
2486# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2489 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
2490 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
2492# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2495 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
2496 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
2498# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2499 us_c = rhol_star*s_m
2501 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
2502 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
2504# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2505 us_c = rhol_star*u_t_l
2506 uss_c = rhol_star*u_t_star
2507 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
2508 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
2510# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2511 us_c = rhol_star*u_t2_l
2512 uss_c = rhol_star*u_t2_star
2513 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
2514 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
2516# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2518 uss_c = e_l_starstar
2519 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
2520 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
2522# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2523 us_c = alpha1_l_star
2525 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
2526 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
2528# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2529 us_c = alpha2_l_star
2531 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
2532 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
2534# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2535 us_c = rhol_star*tau_nn_l_star
2537 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
2538 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
2540# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2541 us_c = rhol_star*tau_nt_l
2542 uss_c = rhol_star*tau_nt_star
2543 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
2544 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
2546# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2547 us_c = rhol_star*tau_nt2_l
2548 uss_c = rhol_star*tau_nt2_star
2549 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
2550 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
2552# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2553 us_c = rhol_star*tau_tt_l_star
2554 uss_c = rhol_star*tau_tt_l_starstar
2555 f_hlld(12) = f_l(12) + s_l*(us_c - u_l(12))
2556 if (zone == 2) f_hlld(12) = f_hlld(12) + s_lstar*(uss_c &
2558# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2559 us_c = rhol_star*tau_t2t2_l_star
2560 uss_c = rhol_star*tau_t2t2_l_starstar
2561 f_hlld(13) = f_l(13) + s_l*(us_c - u_l(13))
2562 if (zone == 2) f_hlld(13) = f_hlld(13) + s_lstar*(uss_c &
2564# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2565 us_c = rhol_star*tau_t1t2_l_star
2566 uss_c = rhol_star*tau_t1t2_l_starstar
2567 f_hlld(14) = f_l(14) + s_l*(us_c - u_l(14))
2568 if (zone == 2) f_hlld(14) = f_hlld(14) + s_lstar*(uss_c &
2570# 739 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2572# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2575 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
2576 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
2578# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2581 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
2582 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
2584# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2585 us_c = rhol_star*s_m
2587 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
2588 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
2590# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2591 us_c = rhol_star*u_t_l
2592 uss_c = rhol_star*u_t_star
2593 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
2594 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
2596# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2598 uss_c = e_l_starstar
2599 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
2600 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
2602# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2603 us_c = alpha1_l_star
2605 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
2606 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
2608# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2609 us_c = alpha2_l_star
2611 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
2612 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
2614# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2615 us_c = rhol_star*tau_nn_l_star
2617 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
2618 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
2620# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2621 us_c = rhol_star*tau_nt_l
2622 uss_c = rhol_star*tau_nt_star
2623 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
2624 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
2626# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2627 us_c = rhol_star*tau_tt_l_star
2628 uss_c = rhol_star*tau_tt_l_starstar
2629 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
2630 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
2632# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2633 us_c = rhol_star*tau_qq_l_star
2635 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
2636 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
2638# 747 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2640# 726 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2641 else if (zone == 3 .or. zone == 4)
then
2646 if (p > 0 .and. .not. cyl_coord)
then
2647# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2650 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
2651 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
2653# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2656 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
2657 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
2659# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2660 us_c = rhor_star*s_m
2662 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
2663 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
2665# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2666 us_c = rhor_star*u_t_r
2667 uss_c = rhor_star*u_t_star
2668 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
2669 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
2671# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2672 us_c = rhor_star*u_t2_r
2673 uss_c = rhor_star*u_t2_star
2674 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
2675 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
2677# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2679 uss_c = e_r_starstar
2680 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
2681 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
2683# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2684 us_c = alpha1_r_star
2686 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
2687 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
2689# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2690 us_c = alpha2_r_star
2692 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
2693 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
2695# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2696 us_c = rhor_star*tau_nn_r_star
2698 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
2699 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
2701# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2702 us_c = rhor_star*tau_nt_r
2703 uss_c = rhor_star*tau_nt_star
2704 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
2705 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
2707# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2708 us_c = rhor_star*tau_nt2_r
2709 uss_c = rhor_star*tau_nt2_star
2710 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
2711 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
2713# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2714 us_c = rhor_star*tau_tt_r_star
2715 uss_c = rhor_star*tau_tt_r_starstar
2716 f_hlld(12) = f_r(12) + s_r*(us_c - u_r(12))
2717 if (zone == 3) f_hlld(12) = f_hlld(12) + s_rstar*(uss_c &
2719# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2720 us_c = rhor_star*tau_t2t2_r_star
2721 uss_c = rhor_star*tau_t2t2_r_starstar
2722 f_hlld(13) = f_r(13) + s_r*(us_c - u_r(13))
2723 if (zone == 3) f_hlld(13) = f_hlld(13) + s_rstar*(uss_c &
2725# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2726 us_c = rhor_star*tau_t1t2_r_star
2727 uss_c = rhor_star*tau_t1t2_r_starstar
2728 f_hlld(14) = f_r(14) + s_r*(us_c - u_r(14))
2729 if (zone == 3) f_hlld(14) = f_hlld(14) + s_rstar*(uss_c &
2731# 739 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2733# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2736 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
2737 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
2739# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2742 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
2743 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
2745# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2746 us_c = rhor_star*s_m
2748 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
2749 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
2751# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2752 us_c = rhor_star*u_t_r
2753 uss_c = rhor_star*u_t_star
2754 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
2755 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
2757# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2759 uss_c = e_r_starstar
2760 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
2761 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
2763# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2764 us_c = alpha1_r_star
2766 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
2767 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
2769# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2770 us_c = alpha2_r_star
2772 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
2773 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
2775# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2776 us_c = rhor_star*tau_nn_r_star
2778 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
2779 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
2781# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2782 us_c = rhor_star*tau_nt_r
2783 uss_c = rhor_star*tau_nt_star
2784 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
2785 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
2787# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2788 us_c = rhor_star*tau_tt_r_star
2789 uss_c = rhor_star*tau_tt_r_starstar
2790 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
2791 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
2793# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2794 us_c = rhor_star*tau_qq_r_star
2796 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
2797 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
2799# 747 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2801# 749 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2806 if (riemann_hypo_adc)
then
2812 if (s_l < 0._wp .and. s_r > 0._wp)
then
2814 f_hll_c =
f_hll_flux(s_l, s_r, f_l(i), f_r(i), u_l(i), u_r(i))
2815 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
2820 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
2830 if (cyl_coord .and. (riemann_hypo_adc .or. fan_fallback))
then
2831 if (0._wp <= s_l)
then
2832 u_n_hll_trace = u_n_l; u_t_hll_trace = u_t_l
2833 p_face_hll = pres%L; tau_qq_face_hll = tau_qq_l
2834 else if (s_r <= 0._wp)
then
2835 u_n_hll_trace = u_n_r; u_t_hll_trace = u_t_r
2836 p_face_hll = pres%R; tau_qq_face_hll = tau_qq_r
2838 u_n_hll_trace = (s_r*u_n_l - s_l*u_n_r)/(s_r - s_l + verysmall)
2839 u_t_hll_trace = (s_r*u_t_l - s_l*u_t_r)/(s_r - s_l + verysmall)
2841 u_hll_c = (s_r*u_r(1) - s_l*u_l(1) - (f_r(1) - f_l(1)))/(s_r - s_l + verysmall)
2843 u_hll_c = (s_r*u_r(2) - s_l*u_l(2) - (f_r(2) - f_l(2)))/(s_r - s_l + verysmall)
2844 rho_hll = rho_hll + u_hll_c
2845 u_hll_c = (s_r*u_r(3) - s_l*u_l(3) - (f_r(3) - f_l(3)))/(s_r - s_l + verysmall)
2846 u_n_hll_cons = u_hll_c/(rho_hll + verysmall)
2847 u_hll_c = (s_r*u_r(8) - s_l*u_l(8) - (f_r(8) - f_l(8)))/(s_r - s_l + verysmall)
2848 tau_nn_hll = u_hll_c/(rho_hll + verysmall)
2849 u_hll_c = (s_r*u_r(11) - s_l*u_l(11) - (f_r(11) - f_l(11)))/(s_r - s_l + verysmall)
2850 tau_qq_face_hll = u_hll_c/(rho_hll + verysmall)
2852 f_hll_c =
f_hll_flux(s_l, s_r, f_l(3), f_r(3), u_l(3), u_r(3))
2853 p_face_hll = f_hll_c - rho_hll*u_n_hll_cons*u_n_hll_cons + tau_nn_hll
2858# 806 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2859 if (ipass == 1)
then
2860 if (p > 0 .and. .not. cyl_coord)
then
2872 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(9)
2873 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
2874 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(12)
2875 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(11)
2876 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(14)
2877 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
2878 else if (
dir_idx(1) == 2)
then
2879 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
2880 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
2881 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(9)
2882 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(14)
2883 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
2884 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
2886 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
2887 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(14)
2888 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(13)
2889 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(10)
2890 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
2891 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(9)
2902 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(9)
2904 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(8)
2905 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(10)
2907 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(10)
2908 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(8)
2911 flux_rsx_vf(j, k, l, eqn_idx%stress%end) = f_hlld(11)
2915# 806 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2916 if (ipass == 2)
then
2917 if (p > 0 .and. .not. cyl_coord)
then
2935 else if (
dir_idx(1) == 2)
then
2972# 863 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2976 if (fan_fallback)
then
2979 u_n_face = u_n_hll_trace; u_t_face = u_t_hll_trace
2985 u_n_face = u_n_l; u_t_face = u_t_l
2986 else if (zone == 1)
then
2987 u_n_face = s_m; u_t_face = u_t_l
2988 else if (zone <= 3)
then
2989 u_n_face = s_m; u_t_face = u_t_star
2990 else if (zone == 4)
then
2991 u_n_face = s_m; u_t_face = u_t_r
2993 u_n_face = u_n_r; u_t_face = u_t_r
2997 if (riemann_hypo_adc)
then
2998 u_n_face = u_n_hll_trace + phi*(u_n_face - u_n_hll_trace)
2999 u_t_face = u_t_hll_trace + phi*(u_t_face - u_t_hll_trace)
3002# 893 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3003 if (ipass == 1)
then
3012# 893 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3013 if (ipass == 2)
then
3022# 903 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3031# 912 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3033# 915 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3034 if (ipass == 1)
then
3036# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3037#if defined(MFC_OpenACC)
3038# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3040# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3041#elif defined(MFC_OpenMP)
3042# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3044# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3049 if (fan_fallback)
then
3052 p_face = p_face_hll; tau_qq_face = tau_qq_face_hll
3056 p_face = pres%L; tau_qq_face = tau_qq_l
3057 else if (zone <= 2)
then
3058 p_face = ptot_star + tau_nn_l_star; tau_qq_face = tau_qq_l_star
3059 else if (zone <= 4)
then
3060 p_face = ptot_star + tau_nn_r_star; tau_qq_face = tau_qq_r_star
3062 p_face = pres%R; tau_qq_face = tau_qq_r
3066 if (riemann_hypo_adc)
then
3067 p_face = p_face_hll + phi*(p_face - p_face_hll)
3068 tau_qq_face = tau_qq_face_hll + phi*(tau_qq_face - tau_qq_face_hll)
3072 & eqn_idx%cont%end +
dir_idx(1)) - p_face + tau_qq_face
3076# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3077#if defined(MFC_OpenACC)
3078# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3080# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3081#elif defined(MFC_OpenMP)
3082# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3084# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3086 do i = eqn_idx%adv%beg, sys_size
3090# 915 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3091 if (ipass == 2)
then
3093# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3094#if defined(MFC_OpenACC)
3095# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3097# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3098#elif defined(MFC_OpenMP)
3099# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3101# 916 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3106 if (fan_fallback)
then
3109 p_face = p_face_hll; tau_qq_face = tau_qq_face_hll
3113 p_face = pres%L; tau_qq_face = tau_qq_l
3114 else if (zone <= 2)
then
3115 p_face = ptot_star + tau_nn_l_star; tau_qq_face = tau_qq_l_star
3116 else if (zone <= 4)
then
3117 p_face = ptot_star + tau_nn_r_star; tau_qq_face = tau_qq_r_star
3119 p_face = pres%R; tau_qq_face = tau_qq_r
3123 if (riemann_hypo_adc)
then
3124 p_face = p_face_hll + phi*(p_face - p_face_hll)
3125 tau_qq_face = tau_qq_face_hll + phi*(tau_qq_face - tau_qq_face_hll)
3129 & eqn_idx%cont%end +
dir_idx(1)) - p_face + tau_qq_face
3133# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3134#if defined(MFC_OpenACC)
3135# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3137# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3138#elif defined(MFC_OpenMP)
3139# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3141# 946 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3143 do i = eqn_idx%adv%beg, sys_size
3147# 952 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3149# 954 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3158# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3159#if defined(MFC_OpenACC)
3160# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3162# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3163#elif defined(MFC_OpenMP)
3164# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3166# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3168# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3171# 178 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3172# 179 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3173# 180 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3174 if (norm_dir == 3)
then
3177# 184 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3183# 190 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3184# 191 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3185# 192 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3186# 193 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3189# 196 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3193# 213 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3194# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3196# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3198# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3199#if defined(MFC_OpenACC)
3200# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3202# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3204# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3206# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3208# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3210# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3212# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3214# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3216# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3218# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3219#elif defined(MFC_OpenMP)
3220# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3222# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3224# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3226# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3228# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3230# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3232# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3234# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3236# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3238# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3240# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3242# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3244# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3246# 224 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3253 do i = 1, eqn_idx%cont%end
3254 alpha_rho_l(i) = ql_prim_rsx_vf(j, k, l, i)
3255 alpha_rho_r(i) = qr_prim_rsx_vf(j, k, l + 1, i)
3264 vel%L(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%cont%end + i)
3265 vel%R(i) = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%cont%end + i)
3268 vel_rms%L = vel%L(1)**2 + vel%L(2)**2 + vel%L(3)**2
3269 vel_rms%R = vel%R(1)**2 + vel%R(2)**2 + vel%R(3)**2
3271 do i = 1, num_fluids
3272 alpha_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%E + i)
3273 alpha_r(i) = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%E + i)
3280 do i = 1, num_fluids
3281 alpha_rho_l(i) = max(0._wp, alpha_rho_l(i))
3282 alpha_l(i) = min(max(0._wp, alpha_l(i)), 1._wp)
3283 alpha_l_sum = alpha_l_sum + alpha_l(i)
3285 alpha_l = alpha_l/max(alpha_l_sum, sgm_eps)
3287 do i = 1, num_fluids
3288 alpha_rho_r(i) = max(0._wp, alpha_rho_r(i))
3289 alpha_r(i) = min(max(0._wp, alpha_r(i)), 1._wp)
3290 alpha_r_sum = alpha_r_sum + alpha_r(i)
3292 alpha_r = alpha_r/max(alpha_r_sum, sgm_eps)
3295 pres%L = ql_prim_rsx_vf(j, k, l, eqn_idx%E)
3296 pres%R = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%E)
3300# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3301#if defined(MFC_OpenACC)
3302# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3304# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3305#elif defined(MFC_OpenMP)
3306# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3308# 276 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3310 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
3311 tau_e_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%stress%beg - 1 + i)
3312 tau_e_r(i) = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%stress%beg - 1 + i)
3315 u_t2_l = 0._wp; u_t2_r = 0._wp
3316 tau_nt2_l = 0._wp; tau_nt2_r = 0._wp
3317 tau_t2t2_l = 0._wp; tau_t2t2_r = 0._wp
3318 tau_t1t2_l = 0._wp; tau_t1t2_r = 0._wp
3325 else if (p == 0)
then
3341 tau_qq_l = tau_e_l(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
3342 tau_qq_r = tau_e_r(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
3344 tau_qq_l = 0._wp; tau_qq_r = 0._wp
3347 ptot_l = pres%L - tau_nn_l
3348 ptot_r = pres%R - tau_nn_r
3354 if (abs(ptot_r - ptot_l) < 1e-12_wp*max(abs(ptot_l), abs(ptot_r), 1._wp))
then
3355 ptot_l = 5e-1_wp*(ptot_l + ptot_r)
3362 g_l = 0._wp; g_r = 0._wp
3364# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3365#if defined(MFC_OpenACC)
3366# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3368# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3369#elif defined(MFC_OpenMP)
3370# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3372# 330 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3374 do i = 1, num_fluids
3375 g_l = g_l + alpha_l(i)*
gs_rs(i)
3376 g_r = g_r + alpha_r(i)*
gs_rs(i)
3386# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3387#if defined(MFC_OpenACC)
3388# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3390# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3391#elif defined(MFC_OpenMP)
3392# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3394# 342 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3396 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
3397 shear_cond = (n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5))
3407 s_l = min(u_n_l - sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)), &
3408 & u_n_r - sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)))
3409 s_r = max(u_n_r + sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)), &
3410 & u_n_l + sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)))
3412 if (p > 0 .and. .not. cyl_coord)
then
3414 u_l(1) = alpha_rho_l(1); u_r(1) = alpha_rho_r(1)
3415 u_l(2) = alpha_rho_l(2); u_r(2) = alpha_rho_r(2)
3416 u_l(3) = rho%L*u_n_l; u_r(3) = rho%R*u_n_r
3417 u_l(4) = rho%L*u_t_l; u_r(4) = rho%R*u_t_r
3418 u_l(5) = rho%L*u_t2_l; u_r(5) = rho%R*u_t2_r
3419 u_l(6) = e%L; u_r(6) = e%R
3420 u_l(7) = alpha_l(1); u_r(7) = alpha_r(1)
3421 u_l(8) = alpha_l(2); u_r(8) = alpha_r(2)
3422 u_l(9) = rho%L*tau_nn_l; u_r(9) = rho%R*tau_nn_r
3423 u_l(10) = rho%L*tau_nt_l; u_r(10) = rho%R*tau_nt_r
3424 u_l(11) = rho%L*tau_nt2_l; u_r(11) = rho%R*tau_nt2_r
3425 u_l(12) = rho%L*tau_tt_l; u_r(12) = rho%R*tau_tt_r
3426 u_l(13) = rho%L*tau_t2t2_l; u_r(13) = rho%R*tau_t2t2_r
3427 u_l(14) = rho%L*tau_t1t2_l; u_r(14) = rho%R*tau_t1t2_r
3429 f_l(1) = u_l(1)*u_n_l; f_r(1) = u_r(1)*u_n_r
3430 f_l(2) = u_l(2)*u_n_l; f_r(2) = u_r(2)*u_n_r
3431 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
3432 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
3433 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
3434 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
3435 f_l(5) = rho%L*u_n_l*u_t2_l - tau_nt2_l
3436 f_r(5) = rho%R*u_n_r*u_t2_r - tau_nt2_r
3437 f_l(6) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l - u_t2_l*tau_nt2_l
3438 f_r(6) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r - u_t2_r*tau_nt2_r
3443 u_l(1) = alpha_rho_l(1)
3444 u_l(2) = alpha_rho_l(2)
3445 u_l(3) = rho%L*u_n_l
3446 u_l(4) = rho%L*u_t_l
3450 u_l(8) = rho%L*tau_nn_l
3451 u_l(9) = rho%L*tau_nt_l
3452 u_l(10) = rho%L*tau_tt_l
3453 u_l(11) = rho%L*tau_qq_l
3455 u_r(1) = alpha_rho_r(1)
3456 u_r(2) = alpha_rho_r(2)
3457 u_r(3) = rho%R*u_n_r
3458 u_r(4) = rho%R*u_t_r
3462 u_r(8) = rho%R*tau_nn_r
3463 u_r(9) = rho%R*tau_nt_r
3464 u_r(10) = rho%R*tau_tt_r
3465 u_r(11) = rho%R*tau_qq_r
3467 f_l(1) = u_l(1)*u_n_l
3468 f_l(2) = u_l(2)*u_n_l
3469 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
3470 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
3471 f_l(5) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l
3473 f_r(1) = u_r(1)*u_n_r
3474 f_r(2) = u_r(2)*u_n_r
3475 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
3476 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
3477 f_r(5) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r
3482 a_l = rho%L*(s_l - u_n_l)
3483 a_r = rho%R*(s_r - u_n_r)
3484 denoma = (a_r - a_l)
3486 s_m = ((ptot_r - ptot_l) + a_l*u_n_l - a_r*u_n_r)/(a_l - a_r + verysmall)
3490 degenerate = (abs(denoma) < verysmall .or. .not. (s_l - verysmall <= s_m .and. s_m <= s_r + verysmall))
3492 if (.not. degenerate)
then
3495 ptot_star = ptot_l + a_l*(s_m - u_n_l)
3500 rhol_star = rho%L*(s_l - u_n_l)/min(s_l - s_m, -verysmall)
3501 rhor_star = rho%R*(s_r - u_n_r)/max(s_r - s_m, verysmall)
3502 fac_l = (s_l - u_n_l)/min(s_l - s_m, -verysmall)
3503 fac_r = (s_r - u_n_r)/max(s_r - s_m, verysmall)
3505 e_l_star = (e%L*(u_n_l - s_l) + u_n_l*ptot_l - s_m*ptot_star)/max(s_m - s_l, verysmall)
3506 e_r_star = (e%R*(u_n_r - s_r) + u_n_r*ptot_r - s_m*ptot_star)/min(s_m - s_r, -verysmall)
3508 if (riemann_hypo_adc)
then
3512 dsigma = sigma_r - sigma_l
3513 sigma_ref = max(max(abs(sigma_l), abs(sigma_r)), verysmall)
3515 a_l_ref = sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L))
3516 a_r_ref = sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R))
3517 a_ref = max(max(a_l_ref, a_r_ref), verysmall)
3519 du_t = u_t_r - u_t_l
3520 dtau_nt = tau_nt_r - tau_nt_l
3521 du_t2 = u_t2_r - u_t2_l
3522 dtau_nt2 = tau_nt2_r - tau_nt2_l
3524 sensor_ptot = (dsigma*dsigma)/((adc_kappa*sigma_ref)**2 + verysmall)
3525 sensor_vt = (du_t*du_t + du_t2*du_t2)/((adc_kappa*a_ref)**2 + verysmall)
3526 sensor_tnt = (dtau_nt*dtau_nt + dtau_nt2*dtau_nt2)/((adc_kappa*sigma_ref)**2 + verysmall)
3528 sensor_combined = sensor_ptot + sensor_tnt + sensor_vt
3530 phi = exp(-(sensor_combined**adc_power))
3538# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3539#if defined(MFC_OpenACC)
3540# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3542# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3543#elif defined(MFC_OpenMP)
3544# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3546# 484 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3549 do i = 1, eqn_idx%cont%end
3550 alpha_rho_hat(i) = q_prim_vf(i)%sf(j, k, l + ipass - 1)
3552 do i = 1, num_fluids
3553 alpha_hat(i) = q_prim_vf(eqn_idx%E + i)%sf(j, k, l + ipass - 1)
3556# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3557#if defined(MFC_OpenACC)
3558# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3560# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3561#elif defined(MFC_OpenMP)
3562# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3564# 492 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3566 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
3567 tau_e_hat(i) = q_prim_vf(eqn_idx%stress%beg - 1 + i)%sf(j, k, l + ipass - 1)
3571 tau_nt2_hat = 0._wp; tau_t2t2_hat = 0._wp; tau_t1t2_hat = 0._wp
3583 tau_qq_hat = tau_e_hat(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
3590# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3591#if defined(MFC_OpenACC)
3592# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3594# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3595#elif defined(MFC_OpenMP)
3596# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3598# 516 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3600 do i = 1, num_fluids
3601 rho_hat = rho_hat + alpha_rho_hat(i)
3606# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3607#if defined(MFC_OpenACC)
3608# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3610# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3611#elif defined(MFC_OpenMP)
3612# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3614# 522 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3616 do i = 1, num_fluids
3617 g_hat = g_hat + alpha_hat(i)*
gs_rs(i)
3622 if (alt_soundspeed)
then
3623 pres_hat = q_prim_vf(eqn_idx%E)%sf(j, k, l + ipass - 1)
3624 alpha_hat_q = alpha_hat(1)
3625 alpha_rho_hat_q = alpha_rho_hat(1)
3627 alpha_hat_q = alpha_hat(2)
3628 alpha_rho_hat_q = alpha_rho_hat(2)
3630 blkmod1_hat = blkmod1_hat + (4._wp/3._wp)*
gs_rs(1)
3631 blkmod2_hat = blkmod2_hat + (4._wp/3._wp)*
gs_rs(2)
3632 k_hat = alpha_hat(1)*alpha_hat(2)*(blkmod2_hat - blkmod1_hat)/(alpha_hat(1)*blkmod2_hat &
3633 & + alpha_hat(2)*blkmod1_hat + verysmall)
3635 c_hat_1 = alpha_hat(1) + k_hat
3636 c_hat_2 = alpha_hat(2) - k_hat
3638 if (p > 0 .and. .not. cyl_coord)
then
3640 f_l(7) = u_l(7)*u_n_l - c_hat_1*u_n_l
3641 f_r(7) = u_r(7)*u_n_r - c_hat_1*u_n_r
3642 f_l(8) = u_l(8)*u_n_l - c_hat_2*u_n_l
3643 f_r(8) = u_r(8)*u_n_r - c_hat_2*u_n_r
3644 f_l(9) = u_l(9)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
3645 f_r(9) = u_r(9)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
3646 f_l(10) = u_l(10)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
3647 f_r(10) = u_r(10)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
3648 f_l(11) = u_l(11)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t2_l
3649 f_r(11) = u_r(11)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t2_r
3650 f_l(12) = u_l(12)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
3651 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
3652 f_r(12) = u_r(12)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
3653 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
3654 f_l(13) = u_l(13)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_l &
3655 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_l
3656 f_r(13) = u_r(13)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_r &
3657 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_r
3658 f_l(14) = u_l(14)*u_n_l + rho_hat*tau_t1t2_hat*u_n_l - rho_hat*tau_nt2_hat*u_t_l &
3659 & - rho_hat*tau_nt_hat*u_t2_l
3660 f_r(14) = u_r(14)*u_n_r + rho_hat*tau_t1t2_hat*u_n_r - rho_hat*tau_nt2_hat*u_t_r &
3661 & - rho_hat*tau_nt_hat*u_t2_r
3664 f_l(6) = u_l(6)*u_n_l - c_hat_1*u_n_l
3665 f_r(6) = u_r(6)*u_n_r - c_hat_1*u_n_r
3666 f_l(7) = u_l(7)*u_n_l - c_hat_2*u_n_l
3667 f_r(7) = u_r(7)*u_n_r - c_hat_2*u_n_r
3668 f_l(8) = u_l(8)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
3669 f_r(8) = u_r(8)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
3670 f_l(9) = u_l(9)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
3671 f_r(9) = u_r(9)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
3672 f_l(10) = u_l(10)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
3673 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
3674 f_r(10) = u_r(10)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
3675 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
3676 f_l(11) = u_l(11)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_l
3677 f_r(11) = u_r(11)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_r
3684 fan_fallback = degenerate
3685 shear_degenerate = .false.
3686 if (.not. fan_fallback)
then
3687 g_eff = g_hat + tau_nn_hat
3690 if (.not. (g_eff >= -g_eff_tol))
then
3692 fan_fallback = .true.
3693 else if (g_eff <= g_eff_tol)
then
3694 shear_degenerate = .true.
3698 c_nc = rho_hat*g_eff
3699 if (.not. (c_nc > 0._wp))
then
3702 shear_degenerate = (rho_hat > 0._wp)
3703 fan_fallback = .not. shear_degenerate
3707 sqrtc_nc = sqrt(c_nc)
3708 s_lstar = s_m - sqrtc_nc/rhol_star
3709 s_rstar = s_m + sqrtc_nc/rhor_star
3711 fan_fallback = (s_lstar < s_l .or. s_rstar > s_r)
3716 if (fan_fallback)
then
3718 if (s_l < 0._wp .and. s_r > 0._wp)
then
3720 f_hlld(i) =
f_hll_flux(s_l, s_r, f_l(i), f_r(i), u_l(i), u_r(i))
3722 else if (s_l >= 0._wp)
then
3723 f_hlld(1:ncomp) = f_l(1:ncomp)
3725 f_hlld(1:ncomp) = f_r(1:ncomp)
3730 ptot_star = 5e-1_wp*(ptot_l + ptot_r)
3733 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
3734 tau_nn_l_star = tau_nn_l
3735 tau_nn_r_star = tau_nn_r
3736 tau_qq_l_star = tau_qq_l
3737 tau_qq_r_star = tau_qq_r
3741 ptot_star = ptot_l + a_l*(s_m - u_n_l)
3743 if (shear_degenerate)
then
3745 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
3746 tau_nt_star = 5e-1_wp*(tau_nt_l + tau_nt_r)
3747 u_t2_star = 5e-1_wp*(u_t2_l + u_t2_r)
3748 tau_nt2_star = 5e-1_wp*(tau_nt2_l + tau_nt2_r)
3750 u_t_star = 5e-1_wp*((tau_nt_r - tau_nt_l)/sqrtc_nc + (u_t_r + u_t_l))
3751 tau_nt_star = 5e-1_wp*((u_t_r - u_t_l)*sqrtc_nc + (tau_nt_r + tau_nt_l))
3752 u_t2_star = 5e-1_wp*((tau_nt2_r - tau_nt2_l)/sqrtc_nc + (u_t2_r + u_t2_l))
3753 tau_nt2_star = 5e-1_wp*((u_t2_r - u_t2_l)*sqrtc_nc + (tau_nt2_r + tau_nt2_l))
3759# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3760# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3761 tau_nn_l_star = tau_nn_l - (rho_hat*(g_hat*4._wp/3._wp + tau_nn_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
3762# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3763 tau_nn_r_star = tau_nn_r - (rho_hat*(g_hat*4._wp/3._wp + tau_nn_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
3764# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3765# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3766# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3767 tau_tt_l_star = tau_tt_l + (rho_hat*(g_hat*2._wp/3._wp + tau_tt_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
3768# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3769 tau_tt_r_star = tau_tt_r + (rho_hat*(g_hat*2._wp/3._wp + tau_tt_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
3770# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3771# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3772# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3773 tau_t2t2_l_star = tau_t2t2_l + (rho_hat*(g_hat*2._wp/3._wp + tau_t2t2_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
3774# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3775 tau_t2t2_r_star = tau_t2t2_r + (rho_hat*(g_hat*2._wp/3._wp + tau_t2t2_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
3776# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3777# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3778# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3779 tau_t1t2_l_star = tau_t1t2_l + (rho_hat*tau_t1t2_hat*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
3780# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3781 tau_t1t2_r_star = tau_t1t2_r + (rho_hat*tau_t1t2_hat*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
3782# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3783# 671 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3784# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3785 tau_qq_l_star = tau_qq_l + (rho_hat*(g_hat*2._wp/3._wp + tau_qq_hat)*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
3786# 672 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3787 tau_qq_r_star = tau_qq_r + (rho_hat*(g_hat*2._wp/3._wp + tau_qq_hat)*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
3788# 674 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3789# 675 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3791 if (shear_degenerate)
then
3793 tau_tt_l_starstar = tau_tt_l_star
3794 tau_tt_r_starstar = tau_tt_r_star
3795 tau_t2t2_l_starstar = tau_t2t2_l_star
3796 tau_t2t2_r_starstar = tau_t2t2_r_star
3797 tau_t1t2_l_starstar = tau_t1t2_l_star
3798 tau_t1t2_r_starstar = tau_t1t2_r_star
3799 e_l_starstar = e_l_star
3800 e_r_starstar = e_r_star
3802 tau_tt_l_starstar = tau_tt_l_star + 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_l)
3803 tau_tt_r_starstar = tau_tt_r_star - 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_r)
3804 tau_t2t2_l_starstar = tau_t2t2_l_star + 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
3806 tau_t2t2_r_starstar = tau_t2t2_r_star - 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
3808 tau_t1t2_l_starstar = tau_t1t2_l_star + rho_hat*(tau_nt2_hat*(u_t_star - u_t_l) &
3809 & + tau_nt_hat*(u_t2_star - u_t2_l))/sqrtc_nc
3810 tau_t1t2_r_starstar = tau_t1t2_r_star - rho_hat*(tau_nt2_hat*(u_t_star - u_t_r) &
3811 & + tau_nt_hat*(u_t2_star - u_t2_r))/sqrtc_nc
3812 e_l_starstar = e_l_star + (rhol_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_l*tau_nt_l) &
3813 & + (u_t2_star*tau_nt2_star - u_t2_l*tau_nt2_l))
3814 e_r_starstar = e_r_star - (rhor_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_r*tau_nt_r) &
3815 & + (u_t2_star*tau_nt2_star - u_t2_r*tau_nt2_r))
3818 alpha1_l_star = (alpha_l(1)*(s_l - u_n_l) - c_hat_1*(s_m - u_n_l))/min(s_l - s_m, -verysmall)
3819 alpha1_r_star = (alpha_r(1)*(s_r - u_n_r) - c_hat_1*(s_m - u_n_r))/max(s_r - s_m, verysmall)
3821 alpha2_l_star = (alpha_l(2)*(s_l - u_n_l) - c_hat_2*(s_m - u_n_l))/min(s_l - s_m, -verysmall)
3822 alpha2_r_star = (alpha_r(2)*(s_r - u_n_r) - c_hat_2*(s_m - u_n_r))/max(s_r - s_m, verysmall)
3837 f_hlld(1:ncomp) = f_l(1:ncomp)
3838 else if (zone == 5)
then
3839 f_hlld(1:ncomp) = f_r(1:ncomp)
3840# 726 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3841 else if (zone == 1 .or. zone == 2)
then
3846 if (p > 0 .and. .not. cyl_coord)
then
3847# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3850 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
3851 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
3853# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3856 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
3857 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
3859# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3860 us_c = rhol_star*s_m
3862 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
3863 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
3865# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3866 us_c = rhol_star*u_t_l
3867 uss_c = rhol_star*u_t_star
3868 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
3869 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
3871# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3872 us_c = rhol_star*u_t2_l
3873 uss_c = rhol_star*u_t2_star
3874 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
3875 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
3877# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3879 uss_c = e_l_starstar
3880 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
3881 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
3883# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3884 us_c = alpha1_l_star
3886 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
3887 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
3889# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3890 us_c = alpha2_l_star
3892 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
3893 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
3895# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3896 us_c = rhol_star*tau_nn_l_star
3898 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
3899 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
3901# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3902 us_c = rhol_star*tau_nt_l
3903 uss_c = rhol_star*tau_nt_star
3904 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
3905 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
3907# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3908 us_c = rhol_star*tau_nt2_l
3909 uss_c = rhol_star*tau_nt2_star
3910 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
3911 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
3913# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3914 us_c = rhol_star*tau_tt_l_star
3915 uss_c = rhol_star*tau_tt_l_starstar
3916 f_hlld(12) = f_l(12) + s_l*(us_c - u_l(12))
3917 if (zone == 2) f_hlld(12) = f_hlld(12) + s_lstar*(uss_c &
3919# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3920 us_c = rhol_star*tau_t2t2_l_star
3921 uss_c = rhol_star*tau_t2t2_l_starstar
3922 f_hlld(13) = f_l(13) + s_l*(us_c - u_l(13))
3923 if (zone == 2) f_hlld(13) = f_hlld(13) + s_lstar*(uss_c &
3925# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3926 us_c = rhol_star*tau_t1t2_l_star
3927 uss_c = rhol_star*tau_t1t2_l_starstar
3928 f_hlld(14) = f_l(14) + s_l*(us_c - u_l(14))
3929 if (zone == 2) f_hlld(14) = f_hlld(14) + s_lstar*(uss_c &
3931# 739 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3933# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3936 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
3937 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
3939# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3942 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
3943 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
3945# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3946 us_c = rhol_star*s_m
3948 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
3949 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
3951# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3952 us_c = rhol_star*u_t_l
3953 uss_c = rhol_star*u_t_star
3954 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
3955 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
3957# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3959 uss_c = e_l_starstar
3960 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
3961 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
3963# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3964 us_c = alpha1_l_star
3966 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
3967 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
3969# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3970 us_c = alpha2_l_star
3972 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
3973 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
3975# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3976 us_c = rhol_star*tau_nn_l_star
3978 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
3979 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
3981# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3982 us_c = rhol_star*tau_nt_l
3983 uss_c = rhol_star*tau_nt_star
3984 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
3985 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
3987# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3988 us_c = rhol_star*tau_tt_l_star
3989 uss_c = rhol_star*tau_tt_l_starstar
3990 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
3991 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
3993# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3994 us_c = rhol_star*tau_qq_l_star
3996 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
3997 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
3999# 747 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4001# 726 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4002 else if (zone == 3 .or. zone == 4)
then
4007 if (p > 0 .and. .not. cyl_coord)
then
4008# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4011 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
4012 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
4014# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4017 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
4018 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
4020# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4021 us_c = rhor_star*s_m
4023 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
4024 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
4026# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4027 us_c = rhor_star*u_t_r
4028 uss_c = rhor_star*u_t_star
4029 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
4030 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
4032# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4033 us_c = rhor_star*u_t2_r
4034 uss_c = rhor_star*u_t2_star
4035 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
4036 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
4038# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4040 uss_c = e_r_starstar
4041 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
4042 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
4044# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4045 us_c = alpha1_r_star
4047 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
4048 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
4050# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4051 us_c = alpha2_r_star
4053 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
4054 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
4056# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4057 us_c = rhor_star*tau_nn_r_star
4059 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
4060 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
4062# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4063 us_c = rhor_star*tau_nt_r
4064 uss_c = rhor_star*tau_nt_star
4065 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
4066 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
4068# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4069 us_c = rhor_star*tau_nt2_r
4070 uss_c = rhor_star*tau_nt2_star
4071 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
4072 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
4074# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4075 us_c = rhor_star*tau_tt_r_star
4076 uss_c = rhor_star*tau_tt_r_starstar
4077 f_hlld(12) = f_r(12) + s_r*(us_c - u_r(12))
4078 if (zone == 3) f_hlld(12) = f_hlld(12) + s_rstar*(uss_c &
4080# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4081 us_c = rhor_star*tau_t2t2_r_star
4082 uss_c = rhor_star*tau_t2t2_r_starstar
4083 f_hlld(13) = f_r(13) + s_r*(us_c - u_r(13))
4084 if (zone == 3) f_hlld(13) = f_hlld(13) + s_rstar*(uss_c &
4086# 733 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4087 us_c = rhor_star*tau_t1t2_r_star
4088 uss_c = rhor_star*tau_t1t2_r_starstar
4089 f_hlld(14) = f_r(14) + s_r*(us_c - u_r(14))
4090 if (zone == 3) f_hlld(14) = f_hlld(14) + s_rstar*(uss_c &
4092# 739 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4094# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4097 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
4098 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
4100# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4103 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
4104 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
4106# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4107 us_c = rhor_star*s_m
4109 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
4110 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
4112# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4113 us_c = rhor_star*u_t_r
4114 uss_c = rhor_star*u_t_star
4115 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
4116 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
4118# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4120 uss_c = e_r_starstar
4121 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
4122 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
4124# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4125 us_c = alpha1_r_star
4127 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
4128 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
4130# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4131 us_c = alpha2_r_star
4133 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
4134 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
4136# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4137 us_c = rhor_star*tau_nn_r_star
4139 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
4140 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
4142# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4143 us_c = rhor_star*tau_nt_r
4144 uss_c = rhor_star*tau_nt_star
4145 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
4146 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
4148# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4149 us_c = rhor_star*tau_tt_r_star
4150 uss_c = rhor_star*tau_tt_r_starstar
4151 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
4152 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
4154# 741 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4155 us_c = rhor_star*tau_qq_r_star
4157 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
4158 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
4160# 747 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4162# 749 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4167 if (riemann_hypo_adc)
then
4173 if (s_l < 0._wp .and. s_r > 0._wp)
then
4175 f_hll_c =
f_hll_flux(s_l, s_r, f_l(i), f_r(i), u_l(i), u_r(i))
4176 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
4181 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
4191 if (cyl_coord .and. (riemann_hypo_adc .or. fan_fallback))
then
4192 if (0._wp <= s_l)
then
4193 u_n_hll_trace = u_n_l; u_t_hll_trace = u_t_l
4194 p_face_hll = pres%L; tau_qq_face_hll = tau_qq_l
4195 else if (s_r <= 0._wp)
then
4196 u_n_hll_trace = u_n_r; u_t_hll_trace = u_t_r
4197 p_face_hll = pres%R; tau_qq_face_hll = tau_qq_r
4199 u_n_hll_trace = (s_r*u_n_l - s_l*u_n_r)/(s_r - s_l + verysmall)
4200 u_t_hll_trace = (s_r*u_t_l - s_l*u_t_r)/(s_r - s_l + verysmall)
4202 u_hll_c = (s_r*u_r(1) - s_l*u_l(1) - (f_r(1) - f_l(1)))/(s_r - s_l + verysmall)
4204 u_hll_c = (s_r*u_r(2) - s_l*u_l(2) - (f_r(2) - f_l(2)))/(s_r - s_l + verysmall)
4205 rho_hll = rho_hll + u_hll_c
4206 u_hll_c = (s_r*u_r(3) - s_l*u_l(3) - (f_r(3) - f_l(3)))/(s_r - s_l + verysmall)
4207 u_n_hll_cons = u_hll_c/(rho_hll + verysmall)
4208 u_hll_c = (s_r*u_r(8) - s_l*u_l(8) - (f_r(8) - f_l(8)))/(s_r - s_l + verysmall)
4209 tau_nn_hll = u_hll_c/(rho_hll + verysmall)
4210 u_hll_c = (s_r*u_r(11) - s_l*u_l(11) - (f_r(11) - f_l(11)))/(s_r - s_l + verysmall)
4211 tau_qq_face_hll = u_hll_c/(rho_hll + verysmall)
4213 f_hll_c =
f_hll_flux(s_l, s_r, f_l(3), f_r(3), u_l(3), u_r(3))
4214 p_face_hll = f_hll_c - rho_hll*u_n_hll_cons*u_n_hll_cons + tau_nn_hll
4219# 806 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4220 if (ipass == 1)
then
4221 if (p > 0 .and. .not. cyl_coord)
then
4233 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(9)
4234 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
4235 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(12)
4236 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(11)
4237 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(14)
4238 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
4239 else if (
dir_idx(1) == 2)
then
4240 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
4241 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
4242 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(9)
4243 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(14)
4244 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
4245 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
4247 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
4248 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(14)
4249 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(13)
4250 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(10)
4251 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
4252 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(9)
4263 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(9)
4265 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(8)
4266 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(10)
4268 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(10)
4269 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(8)
4272 flux_rsx_vf(j, k, l, eqn_idx%stress%end) = f_hlld(11)
4276# 806 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4277 if (ipass == 2)
then
4278 if (p > 0 .and. .not. cyl_coord)
then
4296 else if (
dir_idx(1) == 2)
then
4333# 863 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4337 if (fan_fallback)
then
4340 u_n_face = u_n_hll_trace; u_t_face = u_t_hll_trace
4346 u_n_face = u_n_l; u_t_face = u_t_l
4347 else if (zone == 1)
then
4348 u_n_face = s_m; u_t_face = u_t_l
4349 else if (zone <= 3)
then
4350 u_n_face = s_m; u_t_face = u_t_star
4351 else if (zone == 4)
then
4352 u_n_face = s_m; u_t_face = u_t_r
4354 u_n_face = u_n_r; u_t_face = u_t_r
4358 if (riemann_hypo_adc)
then
4359 u_n_face = u_n_hll_trace + phi*(u_n_face - u_n_hll_trace)
4360 u_t_face = u_t_hll_trace + phi*(u_t_face - u_t_hll_trace)
4363# 893 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4364 if (ipass == 1)
then
4373# 893 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4374 if (ipass == 2)
then
4383# 903 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4392# 954 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4401# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4402#if defined(MFC_OpenACC)
4403# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4405# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4406#elif defined(MFC_OpenMP)
4407# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4409# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4411# 961 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4414# 964 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"