428 & qR_prim_rsx_vf, dqR_prim_dx_vf, dqR_prim_dy_vf, dqR_prim_dz_vf, qR_prim_vf, &
429 & q_prim_vf, flux_vf, flux_src_vf, flux_gsrc_vf, norm_dir, ix, iy, iz)
431 real(wp),
dimension(idwbuff(1)%beg:,idwbuff(2)%beg:,idwbuff(3)%beg:,1:),
intent(inout) :: qL_prim_rsx_vf, qR_prim_rsx_vf
432 type(
scalar_field),
allocatable,
dimension(:),
intent(inout) :: dqL_prim_dx_vf, dqR_prim_dx_vf, dqL_prim_dy_vf, &
433 & dqR_prim_dy_vf, dqL_prim_dz_vf, dqR_prim_dz_vf
435 type(
scalar_field),
allocatable,
dimension(:),
intent(inout) :: qL_prim_vf, qR_prim_vf
436 type(
scalar_field),
dimension(sys_size),
intent(in) :: q_prim_vf
437 type(
scalar_field),
dimension(sys_size),
intent(inout) :: flux_vf, flux_src_vf, flux_gsrc_vf
438 integer,
intent(in) :: norm_dir
447# 83 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
448 real(wp),
dimension(num_fluids) :: alpha_L, alpha_R, alpha_rho_L, alpha_rho_R
449# 85 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
457 real(wp) :: S_L, S_R, s_M, S_Lstar, S_Rstar
458 real(wp) :: pTot_L, pTot_R, rhoL_star, rhoR_star
459 real(wp),
dimension(14) :: U_L, U_R
460 real(wp),
dimension(14) :: F_L, F_R, F_hlld
461 real(wp) :: us_c, uss_c
464 real(wp) :: rho_HLL, u_n_HLL_cons, tau_nn_HLL
465 real(wp) :: u_n_HLL_trace, u_t_HLL_trace
466 real(wp) :: p_face_HLL, tau_qq_face_HLL
471 real(wp) :: G_eff, G_eff_tol, C_NC, sqrtC_NC
472 real(wp) :: A_L, A_R, denomA, fac_L, fac_R
473 real(wp) :: u_n_L, u_t_L, u_n_R, u_t_R
474 real(wp) :: u_t2_L, u_t2_R
475 real(wp) :: tau_nn_L, tau_nt_L, tau_tt_L, tau_nn_R, tau_nt_R, tau_tt_R
476 real(wp) :: tau_nt2_L, tau_nt2_R, tau_t2t2_L, tau_t2t2_R, tau_t1t2_L, tau_t1t2_R
477 real(wp) :: tau_qq_L, tau_qq_R
479# 117 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
480 real(wp),
dimension(eqn_idx%stress%end - eqn_idx%stress%beg + 1) :: tau_e_L, tau_e_R
481# 119 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
483 real(wp) :: alpha1_L_star, alpha1_R_star, alpha2_L_star, alpha2_R_star
484 real(wp) :: u_t_star, tau_nt_star
485 real(wp) :: u_t2_star, tau_nt2_star
486 real(wp) :: tau_nn_L_star, tau_nn_R_star, tau_tt_L_star, tau_tt_R_star
487 real(wp) :: tau_tt_L_starstar, tau_tt_R_starstar
488 real(wp) :: tau_t2t2_L_star, tau_t2t2_R_star
489 real(wp) :: tau_t2t2_L_starstar, tau_t2t2_R_starstar
490 real(wp) :: tau_t1t2_L_star, tau_t1t2_R_star
491 real(wp) :: tau_t1t2_L_starstar, tau_t1t2_R_starstar
492 real(wp) :: tau_qq_L_star, tau_qq_R_star
493 real(wp) :: pTot_star
494 real(wp) :: E_L_star, E_R_star
495 real(wp) :: E_L_starstar, E_R_starstar
496 real(wp) :: p_face, tau_qq_face
497 real(wp) :: u_n_face, u_t_face
500 real(wp) :: tau_nn_hat, tau_nt_hat, tau_tt_hat, tau_qq_hat
501 real(wp) :: tau_nt2_hat, tau_t2t2_hat, tau_t1t2_hat
503# 147 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
504 real(wp),
dimension(num_fluids) :: alpha_hat, alpha_rho_hat
505 real(wp),
dimension(eqn_idx%stress%end - eqn_idx%stress%beg + 1) :: tau_e_hat
506# 150 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
508 real(wp) :: pres_hat, blkmod1_hat, blkmod2_hat, K_hat
509 real(wp) :: C_hat_1, C_hat_2
510 real(wp) :: Sigma_L, Sigma_R, dSigma, Sigma_ref
511 real(wp) :: a_L_ref, a_R_ref, a_ref
512 real(wp) :: du_t, dtau_nt, du_t2, dtau_nt2
513 real(wp) :: sensor_ptot, sensor_vt, sensor_tnt, sensor_combined
515 real(wp),
parameter :: ADC_power = 1.0_wp
516 real(wp) :: alpha_L_sum, alpha_R_sum
517 logical :: degenerate, shear_degenerate, fan_fallback
518 integer :: i, j, k, l, ipass, zone
521 & qr_prim_rsx_vf, dqr_prim_dx_vf, dqr_prim_dy_vf, dqr_prim_dz_vf, norm_dir, ix, iy, iz)
525# 172 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
526# 173 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
527# 174 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
528 if (norm_dir == 1)
then
531# 178 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
537# 184 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
538# 185 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
539# 186 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
540# 187 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
543# 190 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
547# 207 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
548# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
550# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
552# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
553#if defined(MFC_OpenACC)
554# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
556# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
558# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
560# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
562# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
564# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
566# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
568# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
570# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
572# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
573#elif defined(MFC_OpenMP)
574# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
576# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
578# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
580# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
582# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
584# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
586# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
588# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
590# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
592# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
594# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
596# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
598# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
600# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
607 do i = 1, eqn_idx%cont%end
608 alpha_rho_l(i) = ql_prim_rsx_vf(j, k, l, i)
609 alpha_rho_r(i) = qr_prim_rsx_vf(j + 1, k, l, i)
618 vel%L(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%cont%end + i)
619 vel%R(i) = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%cont%end + i)
622 vel_rms%L = vel%L(1)**2 + vel%L(2)**2 + vel%L(3)**2
623 vel_rms%R = vel%R(1)**2 + vel%R(2)**2 + vel%R(3)**2
626 alpha_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%E + i)
627 alpha_r(i) = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%E + i)
635 alpha_rho_l(i) = max(0._wp, alpha_rho_l(i))
636 alpha_l(i) = min(max(0._wp, alpha_l(i)), 1._wp)
637 alpha_l_sum = alpha_l_sum + alpha_l(i)
639 alpha_l = alpha_l/max(alpha_l_sum, sgm_eps)
642 alpha_rho_r(i) = max(0._wp, alpha_rho_r(i))
643 alpha_r(i) = min(max(0._wp, alpha_r(i)), 1._wp)
644 alpha_r_sum = alpha_r_sum + alpha_r(i)
646 alpha_r = alpha_r/max(alpha_r_sum, sgm_eps)
649 pres%L = ql_prim_rsx_vf(j, k, l, eqn_idx%E)
650 pres%R = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%E)
654# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
655#if defined(MFC_OpenACC)
656# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
658# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
659#elif defined(MFC_OpenMP)
660# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
662# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
664 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
665 tau_e_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%stress%beg - 1 + i)
666 tau_e_r(i) = qr_prim_rsx_vf(j + 1, k, l, eqn_idx%stress%beg - 1 + i)
669 u_t2_l = 0._wp; u_t2_r = 0._wp
670 tau_nt2_l = 0._wp; tau_nt2_r = 0._wp
671 tau_t2t2_l = 0._wp; tau_t2t2_r = 0._wp
672 tau_t1t2_l = 0._wp; tau_t1t2_r = 0._wp
679 else if (p == 0)
then
695 tau_qq_l = tau_e_l(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
696 tau_qq_r = tau_e_r(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
698 tau_qq_l = 0._wp; tau_qq_r = 0._wp
701 ptot_l = pres%L - tau_nn_l
702 ptot_r = pres%R - tau_nn_r
708 if (abs(ptot_r - ptot_l) < 1e-12_wp*max(abs(ptot_l), abs(ptot_r), 1._wp))
then
709 ptot_l = 5e-1_wp*(ptot_l + ptot_r)
714 rho%L = 0._wp; gamma%L = 0._wp; pi_inf%L = 0._wp; qv%L = 0._wp
715 rho%R = 0._wp; gamma%R = 0._wp; pi_inf%R = 0._wp; qv%R = 0._wp
717# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
718#if defined(MFC_OpenACC)
719# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
721# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
722#elif defined(MFC_OpenMP)
723# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
725# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
728 rho%L = rho%L + alpha_rho_l(i)
729 gamma%L = gamma%L + alpha_l(i)*gammas(i)
730 pi_inf%L = pi_inf%L + alpha_l(i)*pi_infs(i)
731 qv%L = qv%L + alpha_rho_l(i)*qvs(i)
733 rho%R = rho%R + alpha_rho_r(i)
734 gamma%R = gamma%R + alpha_r(i)*gammas(i)
735 pi_inf%R = pi_inf%R + alpha_r(i)*pi_infs(i)
736 qv%R = qv%R + alpha_rho_r(i)*qvs(i)
739 g_l = 0._wp; g_r = 0._wp
741# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
742#if defined(MFC_OpenACC)
743# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
745# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
746#elif defined(MFC_OpenMP)
747# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
749# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
752 g_l = g_l + alpha_l(i)*
gs_rs(i)
753 g_r = g_r + alpha_r(i)*
gs_rs(i)
756 e%L = gamma%L*pres%L + pi_inf%L + 5e-1_wp*rho%L*vel_rms%L + qv%L
757 e%R = gamma%R*pres%R + pi_inf%R + 5e-1_wp*rho%R*vel_rms%R + qv%R
761 h%L = (e%L + pres%L)/rho%L
762 h%R = (e%R + pres%R)/rho%R
765# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
766#if defined(MFC_OpenACC)
767# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
769# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
770#elif defined(MFC_OpenMP)
771# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
773# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
775 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
776 if (g_l > verysmall)
then
777# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
778 e%L = e%L + (tau_e_l(i)*tau_e_l(i))/max(4._wp*g_l, verysmall)
779# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
780 if ((n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5)))
then
781# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
782 e%L = e%L + (tau_e_l(i)*tau_e_l(i))/max(4._wp*g_l, verysmall)
783# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
785# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
787# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
788 if (g_r > verysmall)
then
789# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
790 e%R = e%R + (tau_e_r(i)*tau_e_r(i))/max(4._wp*g_r, verysmall)
791# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
792 if ((n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5)))
then
793# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
794 e%R = e%R + (tau_e_r(i)*tau_e_r(i))/max(4._wp*g_r, verysmall)
795# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
797# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
799# 356 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
809 s_l = min(u_n_l - sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)), &
810 & u_n_r - sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)))
811 s_r = max(u_n_r + sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)), &
812 & u_n_l + sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)))
814 if (p > 0 .and. .not. cyl_coord)
then
816 u_l(1) = alpha_rho_l(1); u_r(1) = alpha_rho_r(1)
817 u_l(2) = alpha_rho_l(2); u_r(2) = alpha_rho_r(2)
818 u_l(3) = rho%L*u_n_l; u_r(3) = rho%R*u_n_r
819 u_l(4) = rho%L*u_t_l; u_r(4) = rho%R*u_t_r
820 u_l(5) = rho%L*u_t2_l; u_r(5) = rho%R*u_t2_r
821 u_l(6) = e%L; u_r(6) = e%R
822 u_l(7) = alpha_l(1); u_r(7) = alpha_r(1)
823 u_l(8) = alpha_l(2); u_r(8) = alpha_r(2)
824 u_l(9) = rho%L*tau_nn_l; u_r(9) = rho%R*tau_nn_r
825 u_l(10) = rho%L*tau_nt_l; u_r(10) = rho%R*tau_nt_r
826 u_l(11) = rho%L*tau_nt2_l; u_r(11) = rho%R*tau_nt2_r
827 u_l(12) = rho%L*tau_tt_l; u_r(12) = rho%R*tau_tt_r
828 u_l(13) = rho%L*tau_t2t2_l; u_r(13) = rho%R*tau_t2t2_r
829 u_l(14) = rho%L*tau_t1t2_l; u_r(14) = rho%R*tau_t1t2_r
831 f_l(1) = u_l(1)*u_n_l; f_r(1) = u_r(1)*u_n_r
832 f_l(2) = u_l(2)*u_n_l; f_r(2) = u_r(2)*u_n_r
833 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
834 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
835 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
836 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
837 f_l(5) = rho%L*u_n_l*u_t2_l - tau_nt2_l
838 f_r(5) = rho%R*u_n_r*u_t2_r - tau_nt2_r
839 f_l(6) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l - u_t2_l*tau_nt2_l
840 f_r(6) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r - u_t2_r*tau_nt2_r
845 u_l(1) = alpha_rho_l(1)
846 u_l(2) = alpha_rho_l(2)
852 u_l(8) = rho%L*tau_nn_l
853 u_l(9) = rho%L*tau_nt_l
854 u_l(10) = rho%L*tau_tt_l
855 u_l(11) = rho%L*tau_qq_l
857 u_r(1) = alpha_rho_r(1)
858 u_r(2) = alpha_rho_r(2)
864 u_r(8) = rho%R*tau_nn_r
865 u_r(9) = rho%R*tau_nt_r
866 u_r(10) = rho%R*tau_tt_r
867 u_r(11) = rho%R*tau_qq_r
869 f_l(1) = u_l(1)*u_n_l
870 f_l(2) = u_l(2)*u_n_l
871 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
872 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
873 f_l(5) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l
875 f_r(1) = u_r(1)*u_n_r
876 f_r(2) = u_r(2)*u_n_r
877 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
878 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
879 f_r(5) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r
884 a_l = rho%L*(s_l - u_n_l)
885 a_r = rho%R*(s_r - u_n_r)
888 s_m = ((ptot_r - ptot_l) + a_l*u_n_l - a_r*u_n_r)/(a_l - a_r + verysmall)
892 degenerate = (abs(denoma) < verysmall .or. .not. (s_l - verysmall <= s_m .and. s_m <= s_r + verysmall))
894 if (.not. degenerate)
then
897 ptot_star = ptot_l + a_l*(s_m - u_n_l)
902 rhol_star = rho%L*(s_l - u_n_l)/min(s_l - s_m, -verysmall)
903 rhor_star = rho%R*(s_r - u_n_r)/max(s_r - s_m, verysmall)
904 fac_l = (s_l - u_n_l)/min(s_l - s_m, -verysmall)
905 fac_r = (s_r - u_n_r)/max(s_r - s_m, verysmall)
907 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)
908 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)
910 if (riemann_hypo_adc)
then
914 dsigma = sigma_r - sigma_l
915 sigma_ref = max(max(abs(sigma_l), abs(sigma_r)), verysmall)
917 a_l_ref = sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L))
918 a_r_ref = sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R))
919 a_ref = max(max(a_l_ref, a_r_ref), verysmall)
922 dtau_nt = tau_nt_r - tau_nt_l
923 du_t2 = u_t2_r - u_t2_l
924 dtau_nt2 = tau_nt2_r - tau_nt2_l
926 sensor_ptot = (dsigma*dsigma)/((adc_kappa*sigma_ref)**2 + verysmall)
927 sensor_vt = (du_t*du_t + du_t2*du_t2)/((adc_kappa*a_ref)**2 + verysmall)
928 sensor_tnt = (dtau_nt*dtau_nt + dtau_nt2*dtau_nt2)/((adc_kappa*sigma_ref)**2 + verysmall)
930 sensor_combined = sensor_ptot + sensor_tnt + sensor_vt
932 phi = exp(-(sensor_combined**adc_power))
940# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
941#if defined(MFC_OpenACC)
942# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
944# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
945#elif defined(MFC_OpenMP)
946# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
948# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
951 do i = 1, eqn_idx%cont%end
952 alpha_rho_hat(i) = q_prim_vf(i)%sf(j + ipass - 1, k, l)
955 alpha_hat(i) = q_prim_vf(eqn_idx%E + i)%sf(j + ipass - 1, k, l)
958# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
959#if defined(MFC_OpenACC)
960# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
962# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
963#elif defined(MFC_OpenMP)
964# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
966# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
968 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
969 tau_e_hat(i) = q_prim_vf(eqn_idx%stress%beg - 1 + i)%sf(j + ipass - 1, k, l)
973 tau_nt2_hat = 0._wp; tau_t2t2_hat = 0._wp; tau_t1t2_hat = 0._wp
985 tau_qq_hat = tau_e_hat(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
992# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
993#if defined(MFC_OpenACC)
994# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
996# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
997#elif defined(MFC_OpenMP)
998# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1000# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1002 do i = 1, num_fluids
1003 rho_hat = rho_hat + alpha_rho_hat(i)
1008# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1009#if defined(MFC_OpenACC)
1010# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1012# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1013#elif defined(MFC_OpenMP)
1014# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1016# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1018 do i = 1, num_fluids
1019 g_hat = g_hat + alpha_hat(i)*
gs_rs(i)
1024 if (alt_soundspeed)
then
1025 pres_hat = q_prim_vf(eqn_idx%E)%sf(j + ipass - 1, k, l)
1026 blkmod1_hat = ((gammas(1) + 1._wp)*pres_hat + pi_infs(1))/gammas(1) + (4._wp/3._wp)*
gs_rs(1)
1027 blkmod2_hat = ((gammas(2) + 1._wp)*pres_hat + pi_infs(2))/gammas(2) + (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) = (s_r*f_l(i) - s_l*f_r(i) + s_l*s_r*(u_r(i) - u_l(i)))/(s_r - s_l + verysmall)
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# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1156# 677 "/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# 677 "/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# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1161# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1162# 677 "/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# 677 "/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# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1167# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1168# 677 "/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# 677 "/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# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1173# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1174# 677 "/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# 677 "/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# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1179# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1180# 677 "/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# 677 "/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# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1185# 680 "/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# 731 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 744 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1329# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 752 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1397# 731 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 738 "/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# 744 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1490# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 746 "/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# 752 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1558# 754 "/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 = (s_r*f_l(i) - s_l*f_r(i) + s_l*s_r*(u_r(i) - u_l(i)))/(s_r - s_l + verysmall)
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 = (s_r*f_l(3) - s_l*f_r(3) + s_l*s_r*(u_r(3) - u_l(3)))/(s_r - s_l + verysmall)
1610 p_face_hll = f_hll_c - rho_hll*u_n_hll_cons*u_n_hll_cons + tau_nn_hll
1615# 811 "/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# 811 "/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# 868 "/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# 898 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1760 if (ipass == 1)
then
1769# 898 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1770 if (ipass == 2)
then
1779# 908 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1788# 959 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1797# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1798#if defined(MFC_OpenACC)
1799# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1801# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1802#elif defined(MFC_OpenMP)
1803# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1805# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1807# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1810# 172 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1811# 173 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1812# 174 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1813 if (norm_dir == 2)
then
1816# 178 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1822# 184 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1823# 185 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1824# 186 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1825# 187 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1828# 190 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1832# 207 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1833# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1835# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1837# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1838#if defined(MFC_OpenACC)
1839# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1841# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1843# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1845# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1847# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1849# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1851# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1853# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1855# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1857# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1858#elif defined(MFC_OpenMP)
1859# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1861# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1863# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1865# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1867# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1869# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1871# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1873# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1875# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1877# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1879# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1881# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1883# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1885# 218 "/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# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1940#if defined(MFC_OpenACC)
1941# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1943# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1944#elif defined(MFC_OpenMP)
1945# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
1947# 270 "/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)
1999 rho%L = 0._wp; gamma%L = 0._wp; pi_inf%L = 0._wp; qv%L = 0._wp
2000 rho%R = 0._wp; gamma%R = 0._wp; pi_inf%R = 0._wp; qv%R = 0._wp
2002# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2003#if defined(MFC_OpenACC)
2004# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2006# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2007#elif defined(MFC_OpenMP)
2008# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2010# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2012 do i = 1, num_fluids
2013 rho%L = rho%L + alpha_rho_l(i)
2014 gamma%L = gamma%L + alpha_l(i)*gammas(i)
2015 pi_inf%L = pi_inf%L + alpha_l(i)*pi_infs(i)
2016 qv%L = qv%L + alpha_rho_l(i)*qvs(i)
2018 rho%R = rho%R + alpha_rho_r(i)
2019 gamma%R = gamma%R + alpha_r(i)*gammas(i)
2020 pi_inf%R = pi_inf%R + alpha_r(i)*pi_infs(i)
2021 qv%R = qv%R + alpha_rho_r(i)*qvs(i)
2024 g_l = 0._wp; g_r = 0._wp
2026# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2027#if defined(MFC_OpenACC)
2028# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2030# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2031#elif defined(MFC_OpenMP)
2032# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2034# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2036 do i = 1, num_fluids
2037 g_l = g_l + alpha_l(i)*
gs_rs(i)
2038 g_r = g_r + alpha_r(i)*
gs_rs(i)
2041 e%L = gamma%L*pres%L + pi_inf%L + 5e-1_wp*rho%L*vel_rms%L + qv%L
2042 e%R = gamma%R*pres%R + pi_inf%R + 5e-1_wp*rho%R*vel_rms%R + qv%R
2046 h%L = (e%L + pres%L)/rho%L
2047 h%R = (e%R + pres%R)/rho%R
2050# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2051#if defined(MFC_OpenACC)
2052# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2054# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2055#elif defined(MFC_OpenMP)
2056# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2058# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2060 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
2061 if (g_l > verysmall)
then
2062# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2063 e%L = e%L + (tau_e_l(i)*tau_e_l(i))/max(4._wp*g_l, verysmall)
2064# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2065 if ((n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5)))
then
2066# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2067 e%L = e%L + (tau_e_l(i)*tau_e_l(i))/max(4._wp*g_l, verysmall)
2068# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2070# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2072# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2073 if (g_r > verysmall)
then
2074# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2075 e%R = e%R + (tau_e_r(i)*tau_e_r(i))/max(4._wp*g_r, verysmall)
2076# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2077 if ((n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5)))
then
2078# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2079 e%R = e%R + (tau_e_r(i)*tau_e_r(i))/max(4._wp*g_r, verysmall)
2080# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2082# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2084# 356 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2094 s_l = min(u_n_l - sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)), &
2095 & u_n_r - sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)))
2096 s_r = max(u_n_r + sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)), &
2097 & u_n_l + sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)))
2099 if (p > 0 .and. .not. cyl_coord)
then
2101 u_l(1) = alpha_rho_l(1); u_r(1) = alpha_rho_r(1)
2102 u_l(2) = alpha_rho_l(2); u_r(2) = alpha_rho_r(2)
2103 u_l(3) = rho%L*u_n_l; u_r(3) = rho%R*u_n_r
2104 u_l(4) = rho%L*u_t_l; u_r(4) = rho%R*u_t_r
2105 u_l(5) = rho%L*u_t2_l; u_r(5) = rho%R*u_t2_r
2106 u_l(6) = e%L; u_r(6) = e%R
2107 u_l(7) = alpha_l(1); u_r(7) = alpha_r(1)
2108 u_l(8) = alpha_l(2); u_r(8) = alpha_r(2)
2109 u_l(9) = rho%L*tau_nn_l; u_r(9) = rho%R*tau_nn_r
2110 u_l(10) = rho%L*tau_nt_l; u_r(10) = rho%R*tau_nt_r
2111 u_l(11) = rho%L*tau_nt2_l; u_r(11) = rho%R*tau_nt2_r
2112 u_l(12) = rho%L*tau_tt_l; u_r(12) = rho%R*tau_tt_r
2113 u_l(13) = rho%L*tau_t2t2_l; u_r(13) = rho%R*tau_t2t2_r
2114 u_l(14) = rho%L*tau_t1t2_l; u_r(14) = rho%R*tau_t1t2_r
2116 f_l(1) = u_l(1)*u_n_l; f_r(1) = u_r(1)*u_n_r
2117 f_l(2) = u_l(2)*u_n_l; f_r(2) = u_r(2)*u_n_r
2118 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
2119 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
2120 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
2121 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
2122 f_l(5) = rho%L*u_n_l*u_t2_l - tau_nt2_l
2123 f_r(5) = rho%R*u_n_r*u_t2_r - tau_nt2_r
2124 f_l(6) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l - u_t2_l*tau_nt2_l
2125 f_r(6) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r - u_t2_r*tau_nt2_r
2130 u_l(1) = alpha_rho_l(1)
2131 u_l(2) = alpha_rho_l(2)
2132 u_l(3) = rho%L*u_n_l
2133 u_l(4) = rho%L*u_t_l
2137 u_l(8) = rho%L*tau_nn_l
2138 u_l(9) = rho%L*tau_nt_l
2139 u_l(10) = rho%L*tau_tt_l
2140 u_l(11) = rho%L*tau_qq_l
2142 u_r(1) = alpha_rho_r(1)
2143 u_r(2) = alpha_rho_r(2)
2144 u_r(3) = rho%R*u_n_r
2145 u_r(4) = rho%R*u_t_r
2149 u_r(8) = rho%R*tau_nn_r
2150 u_r(9) = rho%R*tau_nt_r
2151 u_r(10) = rho%R*tau_tt_r
2152 u_r(11) = rho%R*tau_qq_r
2154 f_l(1) = u_l(1)*u_n_l
2155 f_l(2) = u_l(2)*u_n_l
2156 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
2157 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
2158 f_l(5) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l
2160 f_r(1) = u_r(1)*u_n_r
2161 f_r(2) = u_r(2)*u_n_r
2162 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
2163 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
2164 f_r(5) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r
2169 a_l = rho%L*(s_l - u_n_l)
2170 a_r = rho%R*(s_r - u_n_r)
2171 denoma = (a_r - a_l)
2173 s_m = ((ptot_r - ptot_l) + a_l*u_n_l - a_r*u_n_r)/(a_l - a_r + verysmall)
2177 degenerate = (abs(denoma) < verysmall .or. .not. (s_l - verysmall <= s_m .and. s_m <= s_r + verysmall))
2179 if (.not. degenerate)
then
2182 ptot_star = ptot_l + a_l*(s_m - u_n_l)
2187 rhol_star = rho%L*(s_l - u_n_l)/min(s_l - s_m, -verysmall)
2188 rhor_star = rho%R*(s_r - u_n_r)/max(s_r - s_m, verysmall)
2189 fac_l = (s_l - u_n_l)/min(s_l - s_m, -verysmall)
2190 fac_r = (s_r - u_n_r)/max(s_r - s_m, verysmall)
2192 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)
2193 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)
2195 if (riemann_hypo_adc)
then
2199 dsigma = sigma_r - sigma_l
2200 sigma_ref = max(max(abs(sigma_l), abs(sigma_r)), verysmall)
2202 a_l_ref = sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L))
2203 a_r_ref = sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R))
2204 a_ref = max(max(a_l_ref, a_r_ref), verysmall)
2206 du_t = u_t_r - u_t_l
2207 dtau_nt = tau_nt_r - tau_nt_l
2208 du_t2 = u_t2_r - u_t2_l
2209 dtau_nt2 = tau_nt2_r - tau_nt2_l
2211 sensor_ptot = (dsigma*dsigma)/((adc_kappa*sigma_ref)**2 + verysmall)
2212 sensor_vt = (du_t*du_t + du_t2*du_t2)/((adc_kappa*a_ref)**2 + verysmall)
2213 sensor_tnt = (dtau_nt*dtau_nt + dtau_nt2*dtau_nt2)/((adc_kappa*sigma_ref)**2 + verysmall)
2215 sensor_combined = sensor_ptot + sensor_tnt + sensor_vt
2217 phi = exp(-(sensor_combined**adc_power))
2225# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2226#if defined(MFC_OpenACC)
2227# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2229# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2230#elif defined(MFC_OpenMP)
2231# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2233# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2236 do i = 1, eqn_idx%cont%end
2237 alpha_rho_hat(i) = q_prim_vf(i)%sf(j, k + ipass - 1, l)
2239 do i = 1, num_fluids
2240 alpha_hat(i) = q_prim_vf(eqn_idx%E + i)%sf(j, k + ipass - 1, l)
2243# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2244#if defined(MFC_OpenACC)
2245# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2247# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2248#elif defined(MFC_OpenMP)
2249# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2251# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2253 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
2254 tau_e_hat(i) = q_prim_vf(eqn_idx%stress%beg - 1 + i)%sf(j, k + ipass - 1, l)
2258 tau_nt2_hat = 0._wp; tau_t2t2_hat = 0._wp; tau_t1t2_hat = 0._wp
2270 tau_qq_hat = tau_e_hat(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
2277# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2278#if defined(MFC_OpenACC)
2279# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2281# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2282#elif defined(MFC_OpenMP)
2283# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2285# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2287 do i = 1, num_fluids
2288 rho_hat = rho_hat + alpha_rho_hat(i)
2293# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2294#if defined(MFC_OpenACC)
2295# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2297# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2298#elif defined(MFC_OpenMP)
2299# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2301# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2303 do i = 1, num_fluids
2304 g_hat = g_hat + alpha_hat(i)*
gs_rs(i)
2309 if (alt_soundspeed)
then
2310 pres_hat = q_prim_vf(eqn_idx%E)%sf(j, k + ipass - 1, l)
2311 blkmod1_hat = ((gammas(1) + 1._wp)*pres_hat + pi_infs(1))/gammas(1) + (4._wp/3._wp)*
gs_rs(1)
2312 blkmod2_hat = ((gammas(2) + 1._wp)*pres_hat + pi_infs(2))/gammas(2) + (4._wp/3._wp)*
gs_rs(2)
2313 k_hat = alpha_hat(1)*alpha_hat(2)*(blkmod2_hat - blkmod1_hat)/(alpha_hat(1)*blkmod2_hat &
2314 & + alpha_hat(2)*blkmod1_hat + verysmall)
2316 c_hat_1 = alpha_hat(1) + k_hat
2317 c_hat_2 = alpha_hat(2) - k_hat
2319 if (p > 0 .and. .not. cyl_coord)
then
2321 f_l(7) = u_l(7)*u_n_l - c_hat_1*u_n_l
2322 f_r(7) = u_r(7)*u_n_r - c_hat_1*u_n_r
2323 f_l(8) = u_l(8)*u_n_l - c_hat_2*u_n_l
2324 f_r(8) = u_r(8)*u_n_r - c_hat_2*u_n_r
2325 f_l(9) = u_l(9)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
2326 f_r(9) = u_r(9)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
2327 f_l(10) = u_l(10)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
2328 f_r(10) = u_r(10)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
2329 f_l(11) = u_l(11)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t2_l
2330 f_r(11) = u_r(11)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t2_r
2331 f_l(12) = u_l(12)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
2332 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
2333 f_r(12) = u_r(12)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
2334 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
2335 f_l(13) = u_l(13)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_l &
2336 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_l
2337 f_r(13) = u_r(13)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_r &
2338 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_r
2339 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 &
2340 & - rho_hat*tau_nt_hat*u_t2_l
2341 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 &
2342 & - rho_hat*tau_nt_hat*u_t2_r
2345 f_l(6) = u_l(6)*u_n_l - c_hat_1*u_n_l
2346 f_r(6) = u_r(6)*u_n_r - c_hat_1*u_n_r
2347 f_l(7) = u_l(7)*u_n_l - c_hat_2*u_n_l
2348 f_r(7) = u_r(7)*u_n_r - c_hat_2*u_n_r
2349 f_l(8) = u_l(8)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
2350 f_r(8) = u_r(8)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
2351 f_l(9) = u_l(9)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
2352 f_r(9) = u_r(9)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
2353 f_l(10) = u_l(10)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
2354 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
2355 f_r(10) = u_r(10)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
2356 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
2357 f_l(11) = u_l(11)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_l
2358 f_r(11) = u_r(11)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_r
2365 fan_fallback = degenerate
2366 shear_degenerate = .false.
2367 if (.not. fan_fallback)
then
2368 g_eff = g_hat + tau_nn_hat
2371 if (.not. (g_eff >= -g_eff_tol))
then
2373 fan_fallback = .true.
2374 else if (g_eff <= g_eff_tol)
then
2375 shear_degenerate = .true.
2379 c_nc = rho_hat*g_eff
2380 if (.not. (c_nc > 0._wp))
then
2383 shear_degenerate = (rho_hat > 0._wp)
2384 fan_fallback = .not. shear_degenerate
2388 sqrtc_nc = sqrt(c_nc)
2389 s_lstar = s_m - sqrtc_nc/rhol_star
2390 s_rstar = s_m + sqrtc_nc/rhor_star
2392 fan_fallback = (s_lstar < s_l .or. s_rstar > s_r)
2397 if (fan_fallback)
then
2399 if (s_l < 0._wp .and. s_r > 0._wp)
then
2401 f_hlld(i) = (s_r*f_l(i) - s_l*f_r(i) + s_l*s_r*(u_r(i) - u_l(i)))/(s_r - s_l + verysmall)
2403 else if (s_l >= 0._wp)
then
2404 f_hlld(1:ncomp) = f_l(1:ncomp)
2406 f_hlld(1:ncomp) = f_r(1:ncomp)
2411 ptot_star = 5e-1_wp*(ptot_l + ptot_r)
2414 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
2415 tau_nn_l_star = tau_nn_l
2416 tau_nn_r_star = tau_nn_r
2417 tau_qq_l_star = tau_qq_l
2418 tau_qq_r_star = tau_qq_r
2422 ptot_star = ptot_l + a_l*(s_m - u_n_l)
2424 if (shear_degenerate)
then
2426 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
2427 tau_nt_star = 5e-1_wp*(tau_nt_l + tau_nt_r)
2428 u_t2_star = 5e-1_wp*(u_t2_l + u_t2_r)
2429 tau_nt2_star = 5e-1_wp*(tau_nt2_l + tau_nt2_r)
2431 u_t_star = 5e-1_wp*((tau_nt_r - tau_nt_l)/sqrtc_nc + (u_t_r + u_t_l))
2432 tau_nt_star = 5e-1_wp*((u_t_r - u_t_l)*sqrtc_nc + (tau_nt_r + tau_nt_l))
2433 u_t2_star = 5e-1_wp*((tau_nt2_r - tau_nt2_l)/sqrtc_nc + (u_t2_r + u_t2_l))
2434 tau_nt2_star = 5e-1_wp*((u_t2_r - u_t2_l)*sqrtc_nc + (tau_nt2_r + tau_nt2_l))
2440# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2441# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2442 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))
2443# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2444 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))
2445# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2446# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2447# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2448 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))
2449# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2450 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))
2451# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2452# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2453# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2454 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))
2455# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2456 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))
2457# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2458# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2459# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2460 tau_t1t2_l_star = tau_t1t2_l + (rho_hat*tau_t1t2_hat*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
2461# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2462 tau_t1t2_r_star = tau_t1t2_r + (rho_hat*tau_t1t2_hat*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
2463# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2464# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2465# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2466 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))
2467# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2468 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))
2469# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2470# 680 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2472 if (shear_degenerate)
then
2474 tau_tt_l_starstar = tau_tt_l_star
2475 tau_tt_r_starstar = tau_tt_r_star
2476 tau_t2t2_l_starstar = tau_t2t2_l_star
2477 tau_t2t2_r_starstar = tau_t2t2_r_star
2478 tau_t1t2_l_starstar = tau_t1t2_l_star
2479 tau_t1t2_r_starstar = tau_t1t2_r_star
2480 e_l_starstar = e_l_star
2481 e_r_starstar = e_r_star
2483 tau_tt_l_starstar = tau_tt_l_star + 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_l)
2484 tau_tt_r_starstar = tau_tt_r_star - 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_r)
2485 tau_t2t2_l_starstar = tau_t2t2_l_star + 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
2487 tau_t2t2_r_starstar = tau_t2t2_r_star - 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
2489 tau_t1t2_l_starstar = tau_t1t2_l_star + rho_hat*(tau_nt2_hat*(u_t_star - u_t_l) &
2490 & + tau_nt_hat*(u_t2_star - u_t2_l))/sqrtc_nc
2491 tau_t1t2_r_starstar = tau_t1t2_r_star - rho_hat*(tau_nt2_hat*(u_t_star - u_t_r) &
2492 & + tau_nt_hat*(u_t2_star - u_t2_r))/sqrtc_nc
2493 e_l_starstar = e_l_star + (rhol_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_l*tau_nt_l) &
2494 & + (u_t2_star*tau_nt2_star - u_t2_l*tau_nt2_l))
2495 e_r_starstar = e_r_star - (rhor_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_r*tau_nt_r) &
2496 & + (u_t2_star*tau_nt2_star - u_t2_r*tau_nt2_r))
2499 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)
2500 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)
2502 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)
2503 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)
2518 f_hlld(1:ncomp) = f_l(1:ncomp)
2519 else if (zone == 5)
then
2520 f_hlld(1:ncomp) = f_r(1:ncomp)
2521# 731 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2522 else if (zone == 1 .or. zone == 2)
then
2527 if (p > 0 .and. .not. cyl_coord)
then
2528# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2531 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
2532 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
2534# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2537 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
2538 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
2540# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2541 us_c = rhol_star*s_m
2543 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
2544 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
2546# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2547 us_c = rhol_star*u_t_l
2548 uss_c = rhol_star*u_t_star
2549 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
2550 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
2552# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2553 us_c = rhol_star*u_t2_l
2554 uss_c = rhol_star*u_t2_star
2555 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
2556 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
2558# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2560 uss_c = e_l_starstar
2561 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
2562 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
2564# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2565 us_c = alpha1_l_star
2567 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
2568 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
2570# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2571 us_c = alpha2_l_star
2573 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
2574 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
2576# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2577 us_c = rhol_star*tau_nn_l_star
2579 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
2580 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
2582# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2583 us_c = rhol_star*tau_nt_l
2584 uss_c = rhol_star*tau_nt_star
2585 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
2586 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
2588# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2589 us_c = rhol_star*tau_nt2_l
2590 uss_c = rhol_star*tau_nt2_star
2591 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
2592 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
2594# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2595 us_c = rhol_star*tau_tt_l_star
2596 uss_c = rhol_star*tau_tt_l_starstar
2597 f_hlld(12) = f_l(12) + s_l*(us_c - u_l(12))
2598 if (zone == 2) f_hlld(12) = f_hlld(12) + s_lstar*(uss_c &
2600# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2601 us_c = rhol_star*tau_t2t2_l_star
2602 uss_c = rhol_star*tau_t2t2_l_starstar
2603 f_hlld(13) = f_l(13) + s_l*(us_c - u_l(13))
2604 if (zone == 2) f_hlld(13) = f_hlld(13) + s_lstar*(uss_c &
2606# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2607 us_c = rhol_star*tau_t1t2_l_star
2608 uss_c = rhol_star*tau_t1t2_l_starstar
2609 f_hlld(14) = f_l(14) + s_l*(us_c - u_l(14))
2610 if (zone == 2) f_hlld(14) = f_hlld(14) + s_lstar*(uss_c &
2612# 744 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2614# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2617 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
2618 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
2620# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2623 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
2624 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
2626# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2627 us_c = rhol_star*s_m
2629 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
2630 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
2632# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2633 us_c = rhol_star*u_t_l
2634 uss_c = rhol_star*u_t_star
2635 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
2636 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
2638# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2640 uss_c = e_l_starstar
2641 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
2642 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
2644# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2645 us_c = alpha1_l_star
2647 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
2648 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
2650# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2651 us_c = alpha2_l_star
2653 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
2654 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
2656# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2657 us_c = rhol_star*tau_nn_l_star
2659 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
2660 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
2662# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2663 us_c = rhol_star*tau_nt_l
2664 uss_c = rhol_star*tau_nt_star
2665 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
2666 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
2668# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2669 us_c = rhol_star*tau_tt_l_star
2670 uss_c = rhol_star*tau_tt_l_starstar
2671 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
2672 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
2674# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2675 us_c = rhol_star*tau_qq_l_star
2677 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
2678 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
2680# 752 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2682# 731 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2683 else if (zone == 3 .or. zone == 4)
then
2688 if (p > 0 .and. .not. cyl_coord)
then
2689# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2692 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
2693 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
2695# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2698 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
2699 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
2701# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2702 us_c = rhor_star*s_m
2704 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
2705 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
2707# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2708 us_c = rhor_star*u_t_r
2709 uss_c = rhor_star*u_t_star
2710 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
2711 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
2713# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2714 us_c = rhor_star*u_t2_r
2715 uss_c = rhor_star*u_t2_star
2716 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
2717 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
2719# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2721 uss_c = e_r_starstar
2722 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
2723 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
2725# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2726 us_c = alpha1_r_star
2728 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
2729 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
2731# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2732 us_c = alpha2_r_star
2734 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
2735 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
2737# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2738 us_c = rhor_star*tau_nn_r_star
2740 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
2741 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
2743# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2744 us_c = rhor_star*tau_nt_r
2745 uss_c = rhor_star*tau_nt_star
2746 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
2747 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
2749# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2750 us_c = rhor_star*tau_nt2_r
2751 uss_c = rhor_star*tau_nt2_star
2752 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
2753 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
2755# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2756 us_c = rhor_star*tau_tt_r_star
2757 uss_c = rhor_star*tau_tt_r_starstar
2758 f_hlld(12) = f_r(12) + s_r*(us_c - u_r(12))
2759 if (zone == 3) f_hlld(12) = f_hlld(12) + s_rstar*(uss_c &
2761# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2762 us_c = rhor_star*tau_t2t2_r_star
2763 uss_c = rhor_star*tau_t2t2_r_starstar
2764 f_hlld(13) = f_r(13) + s_r*(us_c - u_r(13))
2765 if (zone == 3) f_hlld(13) = f_hlld(13) + s_rstar*(uss_c &
2767# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2768 us_c = rhor_star*tau_t1t2_r_star
2769 uss_c = rhor_star*tau_t1t2_r_starstar
2770 f_hlld(14) = f_r(14) + s_r*(us_c - u_r(14))
2771 if (zone == 3) f_hlld(14) = f_hlld(14) + s_rstar*(uss_c &
2773# 744 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2775# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2778 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
2779 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
2781# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2784 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
2785 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
2787# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2788 us_c = rhor_star*s_m
2790 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
2791 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
2793# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2794 us_c = rhor_star*u_t_r
2795 uss_c = rhor_star*u_t_star
2796 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
2797 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
2799# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2801 uss_c = e_r_starstar
2802 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
2803 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
2805# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2806 us_c = alpha1_r_star
2808 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
2809 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
2811# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2812 us_c = alpha2_r_star
2814 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
2815 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
2817# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2818 us_c = rhor_star*tau_nn_r_star
2820 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
2821 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
2823# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2824 us_c = rhor_star*tau_nt_r
2825 uss_c = rhor_star*tau_nt_star
2826 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
2827 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
2829# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2830 us_c = rhor_star*tau_tt_r_star
2831 uss_c = rhor_star*tau_tt_r_starstar
2832 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
2833 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
2835# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2836 us_c = rhor_star*tau_qq_r_star
2838 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
2839 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
2841# 752 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2843# 754 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2848 if (riemann_hypo_adc)
then
2854 if (s_l < 0._wp .and. s_r > 0._wp)
then
2856 f_hll_c = (s_r*f_l(i) - s_l*f_r(i) + s_l*s_r*(u_r(i) - u_l(i)))/(s_r - s_l + verysmall)
2857 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
2862 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
2872 if (cyl_coord .and. (riemann_hypo_adc .or. fan_fallback))
then
2873 if (0._wp <= s_l)
then
2874 u_n_hll_trace = u_n_l; u_t_hll_trace = u_t_l
2875 p_face_hll = pres%L; tau_qq_face_hll = tau_qq_l
2876 else if (s_r <= 0._wp)
then
2877 u_n_hll_trace = u_n_r; u_t_hll_trace = u_t_r
2878 p_face_hll = pres%R; tau_qq_face_hll = tau_qq_r
2880 u_n_hll_trace = (s_r*u_n_l - s_l*u_n_r)/(s_r - s_l + verysmall)
2881 u_t_hll_trace = (s_r*u_t_l - s_l*u_t_r)/(s_r - s_l + verysmall)
2883 u_hll_c = (s_r*u_r(1) - s_l*u_l(1) - (f_r(1) - f_l(1)))/(s_r - s_l + verysmall)
2885 u_hll_c = (s_r*u_r(2) - s_l*u_l(2) - (f_r(2) - f_l(2)))/(s_r - s_l + verysmall)
2886 rho_hll = rho_hll + u_hll_c
2887 u_hll_c = (s_r*u_r(3) - s_l*u_l(3) - (f_r(3) - f_l(3)))/(s_r - s_l + verysmall)
2888 u_n_hll_cons = u_hll_c/(rho_hll + verysmall)
2889 u_hll_c = (s_r*u_r(8) - s_l*u_l(8) - (f_r(8) - f_l(8)))/(s_r - s_l + verysmall)
2890 tau_nn_hll = u_hll_c/(rho_hll + verysmall)
2891 u_hll_c = (s_r*u_r(11) - s_l*u_l(11) - (f_r(11) - f_l(11)))/(s_r - s_l + verysmall)
2892 tau_qq_face_hll = u_hll_c/(rho_hll + verysmall)
2894 f_hll_c = (s_r*f_l(3) - s_l*f_r(3) + s_l*s_r*(u_r(3) - u_l(3)))/(s_r - s_l + verysmall)
2895 p_face_hll = f_hll_c - rho_hll*u_n_hll_cons*u_n_hll_cons + tau_nn_hll
2900# 811 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2901 if (ipass == 1)
then
2902 if (p > 0 .and. .not. cyl_coord)
then
2914 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(9)
2915 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
2916 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(12)
2917 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(11)
2918 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(14)
2919 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
2920 else if (
dir_idx(1) == 2)
then
2921 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
2922 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
2923 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(9)
2924 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(14)
2925 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
2926 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
2928 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
2929 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(14)
2930 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(13)
2931 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(10)
2932 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
2933 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(9)
2944 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(9)
2946 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(8)
2947 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(10)
2949 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(10)
2950 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(8)
2953 flux_rsx_vf(j, k, l, eqn_idx%stress%end) = f_hlld(11)
2957# 811 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
2958 if (ipass == 2)
then
2959 if (p > 0 .and. .not. cyl_coord)
then
2977 else if (
dir_idx(1) == 2)
then
3014# 868 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3018 if (fan_fallback)
then
3021 u_n_face = u_n_hll_trace; u_t_face = u_t_hll_trace
3027 u_n_face = u_n_l; u_t_face = u_t_l
3028 else if (zone == 1)
then
3029 u_n_face = s_m; u_t_face = u_t_l
3030 else if (zone <= 3)
then
3031 u_n_face = s_m; u_t_face = u_t_star
3032 else if (zone == 4)
then
3033 u_n_face = s_m; u_t_face = u_t_r
3035 u_n_face = u_n_r; u_t_face = u_t_r
3039 if (riemann_hypo_adc)
then
3040 u_n_face = u_n_hll_trace + phi*(u_n_face - u_n_hll_trace)
3041 u_t_face = u_t_hll_trace + phi*(u_t_face - u_t_hll_trace)
3044# 898 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3045 if (ipass == 1)
then
3054# 898 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3055 if (ipass == 2)
then
3064# 908 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3073# 917 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3075# 920 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3076 if (ipass == 1)
then
3078# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3079#if defined(MFC_OpenACC)
3080# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3082# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3083#elif defined(MFC_OpenMP)
3084# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3086# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3091 if (fan_fallback)
then
3094 p_face = p_face_hll; tau_qq_face = tau_qq_face_hll
3098 p_face = pres%L; tau_qq_face = tau_qq_l
3099 else if (zone <= 2)
then
3100 p_face = ptot_star + tau_nn_l_star; tau_qq_face = tau_qq_l_star
3101 else if (zone <= 4)
then
3102 p_face = ptot_star + tau_nn_r_star; tau_qq_face = tau_qq_r_star
3104 p_face = pres%R; tau_qq_face = tau_qq_r
3108 if (riemann_hypo_adc)
then
3109 p_face = p_face_hll + phi*(p_face - p_face_hll)
3110 tau_qq_face = tau_qq_face_hll + phi*(tau_qq_face - tau_qq_face_hll)
3114 & eqn_idx%cont%end +
dir_idx(1)) - p_face + tau_qq_face
3118# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3119#if defined(MFC_OpenACC)
3120# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3122# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3123#elif defined(MFC_OpenMP)
3124# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3126# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3128 do i = eqn_idx%adv%beg, sys_size
3132# 920 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3133 if (ipass == 2)
then
3135# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3136#if defined(MFC_OpenACC)
3137# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3139# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3140#elif defined(MFC_OpenMP)
3141# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3143# 921 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3148 if (fan_fallback)
then
3151 p_face = p_face_hll; tau_qq_face = tau_qq_face_hll
3155 p_face = pres%L; tau_qq_face = tau_qq_l
3156 else if (zone <= 2)
then
3157 p_face = ptot_star + tau_nn_l_star; tau_qq_face = tau_qq_l_star
3158 else if (zone <= 4)
then
3159 p_face = ptot_star + tau_nn_r_star; tau_qq_face = tau_qq_r_star
3161 p_face = pres%R; tau_qq_face = tau_qq_r
3165 if (riemann_hypo_adc)
then
3166 p_face = p_face_hll + phi*(p_face - p_face_hll)
3167 tau_qq_face = tau_qq_face_hll + phi*(tau_qq_face - tau_qq_face_hll)
3171 & eqn_idx%cont%end +
dir_idx(1)) - p_face + tau_qq_face
3175# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3176#if defined(MFC_OpenACC)
3177# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3179# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3180#elif defined(MFC_OpenMP)
3181# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3183# 951 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3185 do i = eqn_idx%adv%beg, sys_size
3189# 957 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3191# 959 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3200# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3201#if defined(MFC_OpenACC)
3202# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3204# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3205#elif defined(MFC_OpenMP)
3206# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3208# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3210# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3213# 172 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3214# 173 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3215# 174 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3216 if (norm_dir == 3)
then
3219# 178 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3225# 184 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3226# 185 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3227# 186 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3228# 187 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3231# 190 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3235# 207 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3236# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3238# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3240# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3241#if defined(MFC_OpenACC)
3242# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3244# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3246# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3248# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3250# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3252# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3254# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3256# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3258# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3260# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3261#elif defined(MFC_OpenMP)
3262# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3264# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3266# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3268# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3270# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3272# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3274# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3276# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3278# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3280# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3282# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3284# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3286# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3288# 218 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3295 do i = 1, eqn_idx%cont%end
3296 alpha_rho_l(i) = ql_prim_rsx_vf(j, k, l, i)
3297 alpha_rho_r(i) = qr_prim_rsx_vf(j, k, l + 1, i)
3306 vel%L(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%cont%end + i)
3307 vel%R(i) = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%cont%end + i)
3310 vel_rms%L = vel%L(1)**2 + vel%L(2)**2 + vel%L(3)**2
3311 vel_rms%R = vel%R(1)**2 + vel%R(2)**2 + vel%R(3)**2
3313 do i = 1, num_fluids
3314 alpha_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%E + i)
3315 alpha_r(i) = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%E + i)
3322 do i = 1, num_fluids
3323 alpha_rho_l(i) = max(0._wp, alpha_rho_l(i))
3324 alpha_l(i) = min(max(0._wp, alpha_l(i)), 1._wp)
3325 alpha_l_sum = alpha_l_sum + alpha_l(i)
3327 alpha_l = alpha_l/max(alpha_l_sum, sgm_eps)
3329 do i = 1, num_fluids
3330 alpha_rho_r(i) = max(0._wp, alpha_rho_r(i))
3331 alpha_r(i) = min(max(0._wp, alpha_r(i)), 1._wp)
3332 alpha_r_sum = alpha_r_sum + alpha_r(i)
3334 alpha_r = alpha_r/max(alpha_r_sum, sgm_eps)
3337 pres%L = ql_prim_rsx_vf(j, k, l, eqn_idx%E)
3338 pres%R = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%E)
3342# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3343#if defined(MFC_OpenACC)
3344# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3346# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3347#elif defined(MFC_OpenMP)
3348# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3350# 270 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3352 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
3353 tau_e_l(i) = ql_prim_rsx_vf(j, k, l, eqn_idx%stress%beg - 1 + i)
3354 tau_e_r(i) = qr_prim_rsx_vf(j, k, l + 1, eqn_idx%stress%beg - 1 + i)
3357 u_t2_l = 0._wp; u_t2_r = 0._wp
3358 tau_nt2_l = 0._wp; tau_nt2_r = 0._wp
3359 tau_t2t2_l = 0._wp; tau_t2t2_r = 0._wp
3360 tau_t1t2_l = 0._wp; tau_t1t2_r = 0._wp
3367 else if (p == 0)
then
3383 tau_qq_l = tau_e_l(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
3384 tau_qq_r = tau_e_r(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
3386 tau_qq_l = 0._wp; tau_qq_r = 0._wp
3389 ptot_l = pres%L - tau_nn_l
3390 ptot_r = pres%R - tau_nn_r
3396 if (abs(ptot_r - ptot_l) < 1e-12_wp*max(abs(ptot_l), abs(ptot_r), 1._wp))
then
3397 ptot_l = 5e-1_wp*(ptot_l + ptot_r)
3402 rho%L = 0._wp; gamma%L = 0._wp; pi_inf%L = 0._wp; qv%L = 0._wp
3403 rho%R = 0._wp; gamma%R = 0._wp; pi_inf%R = 0._wp; qv%R = 0._wp
3405# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3406#if defined(MFC_OpenACC)
3407# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3409# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3410#elif defined(MFC_OpenMP)
3411# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3413# 323 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3415 do i = 1, num_fluids
3416 rho%L = rho%L + alpha_rho_l(i)
3417 gamma%L = gamma%L + alpha_l(i)*gammas(i)
3418 pi_inf%L = pi_inf%L + alpha_l(i)*pi_infs(i)
3419 qv%L = qv%L + alpha_rho_l(i)*qvs(i)
3421 rho%R = rho%R + alpha_rho_r(i)
3422 gamma%R = gamma%R + alpha_r(i)*gammas(i)
3423 pi_inf%R = pi_inf%R + alpha_r(i)*pi_infs(i)
3424 qv%R = qv%R + alpha_rho_r(i)*qvs(i)
3427 g_l = 0._wp; g_r = 0._wp
3429# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3430#if defined(MFC_OpenACC)
3431# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3433# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3434#elif defined(MFC_OpenMP)
3435# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3437# 337 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3439 do i = 1, num_fluids
3440 g_l = g_l + alpha_l(i)*
gs_rs(i)
3441 g_r = g_r + alpha_r(i)*
gs_rs(i)
3444 e%L = gamma%L*pres%L + pi_inf%L + 5e-1_wp*rho%L*vel_rms%L + qv%L
3445 e%R = gamma%R*pres%R + pi_inf%R + 5e-1_wp*rho%R*vel_rms%R + qv%R
3449 h%L = (e%L + pres%L)/rho%L
3450 h%R = (e%R + pres%R)/rho%R
3453# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3454#if defined(MFC_OpenACC)
3455# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3457# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3458#elif defined(MFC_OpenMP)
3459# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3461# 351 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3463 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
3464 if (g_l > verysmall)
then
3465# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3466 e%L = e%L + (tau_e_l(i)*tau_e_l(i))/max(4._wp*g_l, verysmall)
3467# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3468 if ((n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5)))
then
3469# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3470 e%L = e%L + (tau_e_l(i)*tau_e_l(i))/max(4._wp*g_l, verysmall)
3471# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3473# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3475# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3476 if (g_r > verysmall)
then
3477# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3478 e%R = e%R + (tau_e_r(i)*tau_e_r(i))/max(4._wp*g_r, verysmall)
3479# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3480 if ((n > 0 .and. p == 0 .and. i == 2) .or. (p > 0 .and. (i == 2 .or. i == 4 .or. i == 5)))
then
3481# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3482 e%R = e%R + (tau_e_r(i)*tau_e_r(i))/max(4._wp*g_r, verysmall)
3483# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3485# 353 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3487# 356 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3497 s_l = min(u_n_l - sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)), &
3498 & u_n_r - sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)))
3499 s_r = max(u_n_r + sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R)), &
3500 & u_n_l + sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L)))
3502 if (p > 0 .and. .not. cyl_coord)
then
3504 u_l(1) = alpha_rho_l(1); u_r(1) = alpha_rho_r(1)
3505 u_l(2) = alpha_rho_l(2); u_r(2) = alpha_rho_r(2)
3506 u_l(3) = rho%L*u_n_l; u_r(3) = rho%R*u_n_r
3507 u_l(4) = rho%L*u_t_l; u_r(4) = rho%R*u_t_r
3508 u_l(5) = rho%L*u_t2_l; u_r(5) = rho%R*u_t2_r
3509 u_l(6) = e%L; u_r(6) = e%R
3510 u_l(7) = alpha_l(1); u_r(7) = alpha_r(1)
3511 u_l(8) = alpha_l(2); u_r(8) = alpha_r(2)
3512 u_l(9) = rho%L*tau_nn_l; u_r(9) = rho%R*tau_nn_r
3513 u_l(10) = rho%L*tau_nt_l; u_r(10) = rho%R*tau_nt_r
3514 u_l(11) = rho%L*tau_nt2_l; u_r(11) = rho%R*tau_nt2_r
3515 u_l(12) = rho%L*tau_tt_l; u_r(12) = rho%R*tau_tt_r
3516 u_l(13) = rho%L*tau_t2t2_l; u_r(13) = rho%R*tau_t2t2_r
3517 u_l(14) = rho%L*tau_t1t2_l; u_r(14) = rho%R*tau_t1t2_r
3519 f_l(1) = u_l(1)*u_n_l; f_r(1) = u_r(1)*u_n_r
3520 f_l(2) = u_l(2)*u_n_l; f_r(2) = u_r(2)*u_n_r
3521 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
3522 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
3523 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
3524 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
3525 f_l(5) = rho%L*u_n_l*u_t2_l - tau_nt2_l
3526 f_r(5) = rho%R*u_n_r*u_t2_r - tau_nt2_r
3527 f_l(6) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l - u_t2_l*tau_nt2_l
3528 f_r(6) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r - u_t2_r*tau_nt2_r
3533 u_l(1) = alpha_rho_l(1)
3534 u_l(2) = alpha_rho_l(2)
3535 u_l(3) = rho%L*u_n_l
3536 u_l(4) = rho%L*u_t_l
3540 u_l(8) = rho%L*tau_nn_l
3541 u_l(9) = rho%L*tau_nt_l
3542 u_l(10) = rho%L*tau_tt_l
3543 u_l(11) = rho%L*tau_qq_l
3545 u_r(1) = alpha_rho_r(1)
3546 u_r(2) = alpha_rho_r(2)
3547 u_r(3) = rho%R*u_n_r
3548 u_r(4) = rho%R*u_t_r
3552 u_r(8) = rho%R*tau_nn_r
3553 u_r(9) = rho%R*tau_nt_r
3554 u_r(10) = rho%R*tau_tt_r
3555 u_r(11) = rho%R*tau_qq_r
3557 f_l(1) = u_l(1)*u_n_l
3558 f_l(2) = u_l(2)*u_n_l
3559 f_l(3) = rho%L*u_n_l*u_n_l + ptot_l
3560 f_l(4) = rho%L*u_n_l*u_t_l - tau_nt_l
3561 f_l(5) = (e%L + ptot_l)*u_n_l - u_t_l*tau_nt_l
3563 f_r(1) = u_r(1)*u_n_r
3564 f_r(2) = u_r(2)*u_n_r
3565 f_r(3) = rho%R*u_n_r*u_n_r + ptot_r
3566 f_r(4) = rho%R*u_n_r*u_t_r - tau_nt_r
3567 f_r(5) = (e%R + ptot_r)*u_n_r - u_t_r*tau_nt_r
3572 a_l = rho%L*(s_l - u_n_l)
3573 a_r = rho%R*(s_r - u_n_r)
3574 denoma = (a_r - a_l)
3576 s_m = ((ptot_r - ptot_l) + a_l*u_n_l - a_r*u_n_r)/(a_l - a_r + verysmall)
3580 degenerate = (abs(denoma) < verysmall .or. .not. (s_l - verysmall <= s_m .and. s_m <= s_r + verysmall))
3582 if (.not. degenerate)
then
3585 ptot_star = ptot_l + a_l*(s_m - u_n_l)
3590 rhol_star = rho%L*(s_l - u_n_l)/min(s_l - s_m, -verysmall)
3591 rhor_star = rho%R*(s_r - u_n_r)/max(s_r - s_m, verysmall)
3592 fac_l = (s_l - u_n_l)/min(s_l - s_m, -verysmall)
3593 fac_r = (s_r - u_n_r)/max(s_r - s_m, verysmall)
3595 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)
3596 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)
3598 if (riemann_hypo_adc)
then
3602 dsigma = sigma_r - sigma_l
3603 sigma_ref = max(max(abs(sigma_l), abs(sigma_r)), verysmall)
3605 a_l_ref = sqrt(max(verysmall, c%L*c%L + ((4._wp/3._wp)*g_l + tau_nn_l)/rho%L))
3606 a_r_ref = sqrt(max(verysmall, c%R*c%R + ((4._wp/3._wp)*g_r + tau_nn_r)/rho%R))
3607 a_ref = max(max(a_l_ref, a_r_ref), verysmall)
3609 du_t = u_t_r - u_t_l
3610 dtau_nt = tau_nt_r - tau_nt_l
3611 du_t2 = u_t2_r - u_t2_l
3612 dtau_nt2 = tau_nt2_r - tau_nt2_l
3614 sensor_ptot = (dsigma*dsigma)/((adc_kappa*sigma_ref)**2 + verysmall)
3615 sensor_vt = (du_t*du_t + du_t2*du_t2)/((adc_kappa*a_ref)**2 + verysmall)
3616 sensor_tnt = (dtau_nt*dtau_nt + dtau_nt2*dtau_nt2)/((adc_kappa*sigma_ref)**2 + verysmall)
3618 sensor_combined = sensor_ptot + sensor_tnt + sensor_vt
3620 phi = exp(-(sensor_combined**adc_power))
3628# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3629#if defined(MFC_OpenACC)
3630# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3632# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3633#elif defined(MFC_OpenMP)
3634# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3636# 495 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3639 do i = 1, eqn_idx%cont%end
3640 alpha_rho_hat(i) = q_prim_vf(i)%sf(j, k, l + ipass - 1)
3642 do i = 1, num_fluids
3643 alpha_hat(i) = q_prim_vf(eqn_idx%E + i)%sf(j, k, l + ipass - 1)
3646# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3647#if defined(MFC_OpenACC)
3648# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3650# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3651#elif defined(MFC_OpenMP)
3652# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3654# 503 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3656 do i = 1, eqn_idx%stress%end - eqn_idx%stress%beg + 1
3657 tau_e_hat(i) = q_prim_vf(eqn_idx%stress%beg - 1 + i)%sf(j, k, l + ipass - 1)
3661 tau_nt2_hat = 0._wp; tau_t2t2_hat = 0._wp; tau_t1t2_hat = 0._wp
3673 tau_qq_hat = tau_e_hat(eqn_idx%stress%end - eqn_idx%stress%beg + 1)
3680# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3681#if defined(MFC_OpenACC)
3682# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3684# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3685#elif defined(MFC_OpenMP)
3686# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3688# 527 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3690 do i = 1, num_fluids
3691 rho_hat = rho_hat + alpha_rho_hat(i)
3696# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3697#if defined(MFC_OpenACC)
3698# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3700# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3701#elif defined(MFC_OpenMP)
3702# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3704# 533 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3706 do i = 1, num_fluids
3707 g_hat = g_hat + alpha_hat(i)*
gs_rs(i)
3712 if (alt_soundspeed)
then
3713 pres_hat = q_prim_vf(eqn_idx%E)%sf(j, k, l + ipass - 1)
3714 blkmod1_hat = ((gammas(1) + 1._wp)*pres_hat + pi_infs(1))/gammas(1) + (4._wp/3._wp)*
gs_rs(1)
3715 blkmod2_hat = ((gammas(2) + 1._wp)*pres_hat + pi_infs(2))/gammas(2) + (4._wp/3._wp)*
gs_rs(2)
3716 k_hat = alpha_hat(1)*alpha_hat(2)*(blkmod2_hat - blkmod1_hat)/(alpha_hat(1)*blkmod2_hat &
3717 & + alpha_hat(2)*blkmod1_hat + verysmall)
3719 c_hat_1 = alpha_hat(1) + k_hat
3720 c_hat_2 = alpha_hat(2) - k_hat
3722 if (p > 0 .and. .not. cyl_coord)
then
3724 f_l(7) = u_l(7)*u_n_l - c_hat_1*u_n_l
3725 f_r(7) = u_r(7)*u_n_r - c_hat_1*u_n_r
3726 f_l(8) = u_l(8)*u_n_l - c_hat_2*u_n_l
3727 f_r(8) = u_r(8)*u_n_r - c_hat_2*u_n_r
3728 f_l(9) = u_l(9)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
3729 f_r(9) = u_r(9)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
3730 f_l(10) = u_l(10)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
3731 f_r(10) = u_r(10)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
3732 f_l(11) = u_l(11)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t2_l
3733 f_r(11) = u_r(11)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t2_r
3734 f_l(12) = u_l(12)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
3735 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
3736 f_r(12) = u_r(12)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
3737 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
3738 f_l(13) = u_l(13)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_l &
3739 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_l
3740 f_r(13) = u_r(13)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_t2t2_hat)*u_n_r &
3741 & - 2._wp*rho_hat*tau_nt2_hat*u_t2_r
3742 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 &
3743 & - rho_hat*tau_nt_hat*u_t2_l
3744 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 &
3745 & - rho_hat*tau_nt_hat*u_t2_r
3748 f_l(6) = u_l(6)*u_n_l - c_hat_1*u_n_l
3749 f_r(6) = u_r(6)*u_n_r - c_hat_1*u_n_r
3750 f_l(7) = u_l(7)*u_n_l - c_hat_2*u_n_l
3751 f_r(7) = u_r(7)*u_n_r - c_hat_2*u_n_r
3752 f_l(8) = u_l(8)*u_n_l - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_l
3753 f_r(8) = u_r(8)*u_n_r - rho_hat*(4._wp/3._wp*g_hat + tau_nn_hat)*u_n_r
3754 f_l(9) = u_l(9)*u_n_l - rho_hat*(g_hat + tau_nn_hat)*u_t_l
3755 f_r(9) = u_r(9)*u_n_r - rho_hat*(g_hat + tau_nn_hat)*u_t_r
3756 f_l(10) = u_l(10)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_l &
3757 & - 2._wp*rho_hat*tau_nt_hat*u_t_l
3758 f_r(10) = u_r(10)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_tt_hat)*u_n_r &
3759 & - 2._wp*rho_hat*tau_nt_hat*u_t_r
3760 f_l(11) = u_l(11)*u_n_l + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_l
3761 f_r(11) = u_r(11)*u_n_r + rho_hat*(2._wp/3._wp*g_hat + tau_qq_hat)*u_n_r
3768 fan_fallback = degenerate
3769 shear_degenerate = .false.
3770 if (.not. fan_fallback)
then
3771 g_eff = g_hat + tau_nn_hat
3774 if (.not. (g_eff >= -g_eff_tol))
then
3776 fan_fallback = .true.
3777 else if (g_eff <= g_eff_tol)
then
3778 shear_degenerate = .true.
3782 c_nc = rho_hat*g_eff
3783 if (.not. (c_nc > 0._wp))
then
3786 shear_degenerate = (rho_hat > 0._wp)
3787 fan_fallback = .not. shear_degenerate
3791 sqrtc_nc = sqrt(c_nc)
3792 s_lstar = s_m - sqrtc_nc/rhol_star
3793 s_rstar = s_m + sqrtc_nc/rhor_star
3795 fan_fallback = (s_lstar < s_l .or. s_rstar > s_r)
3800 if (fan_fallback)
then
3802 if (s_l < 0._wp .and. s_r > 0._wp)
then
3804 f_hlld(i) = (s_r*f_l(i) - s_l*f_r(i) + s_l*s_r*(u_r(i) - u_l(i)))/(s_r - s_l + verysmall)
3806 else if (s_l >= 0._wp)
then
3807 f_hlld(1:ncomp) = f_l(1:ncomp)
3809 f_hlld(1:ncomp) = f_r(1:ncomp)
3814 ptot_star = 5e-1_wp*(ptot_l + ptot_r)
3817 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
3818 tau_nn_l_star = tau_nn_l
3819 tau_nn_r_star = tau_nn_r
3820 tau_qq_l_star = tau_qq_l
3821 tau_qq_r_star = tau_qq_r
3825 ptot_star = ptot_l + a_l*(s_m - u_n_l)
3827 if (shear_degenerate)
then
3829 u_t_star = 5e-1_wp*(u_t_l + u_t_r)
3830 tau_nt_star = 5e-1_wp*(tau_nt_l + tau_nt_r)
3831 u_t2_star = 5e-1_wp*(u_t2_l + u_t2_r)
3832 tau_nt2_star = 5e-1_wp*(tau_nt2_l + tau_nt2_r)
3834 u_t_star = 5e-1_wp*((tau_nt_r - tau_nt_l)/sqrtc_nc + (u_t_r + u_t_l))
3835 tau_nt_star = 5e-1_wp*((u_t_r - u_t_l)*sqrtc_nc + (tau_nt_r + tau_nt_l))
3836 u_t2_star = 5e-1_wp*((tau_nt2_r - tau_nt2_l)/sqrtc_nc + (u_t2_r + u_t2_l))
3837 tau_nt2_star = 5e-1_wp*((u_t2_r - u_t2_l)*sqrtc_nc + (tau_nt2_r + tau_nt2_l))
3843# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3844# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3845 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))
3846# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3847 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))
3848# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3849# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3850# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3851 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))
3852# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3853 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))
3854# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3855# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3856# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3857 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))
3858# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3859 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))
3860# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3861# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3862# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3863 tau_t1t2_l_star = tau_t1t2_l + (rho_hat*tau_t1t2_hat*(u_n_l - s_m))/(rho%L*(u_n_l - s_l))
3864# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3865 tau_t1t2_r_star = tau_t1t2_r + (rho_hat*tau_t1t2_hat*(u_n_r - s_m))/(rho%R*(u_n_r - s_r))
3866# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3867# 676 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3868# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3869 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))
3870# 677 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3871 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))
3872# 679 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3873# 680 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3875 if (shear_degenerate)
then
3877 tau_tt_l_starstar = tau_tt_l_star
3878 tau_tt_r_starstar = tau_tt_r_star
3879 tau_t2t2_l_starstar = tau_t2t2_l_star
3880 tau_t2t2_r_starstar = tau_t2t2_r_star
3881 tau_t1t2_l_starstar = tau_t1t2_l_star
3882 tau_t1t2_r_starstar = tau_t1t2_r_star
3883 e_l_starstar = e_l_star
3884 e_r_starstar = e_r_star
3886 tau_tt_l_starstar = tau_tt_l_star + 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_l)
3887 tau_tt_r_starstar = tau_tt_r_star - 2._wp*rho_hat*tau_nt_hat/sqrtc_nc*(u_t_star - u_t_r)
3888 tau_t2t2_l_starstar = tau_t2t2_l_star + 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
3890 tau_t2t2_r_starstar = tau_t2t2_r_star - 2._wp*rho_hat*tau_nt2_hat/sqrtc_nc*(u_t2_star &
3892 tau_t1t2_l_starstar = tau_t1t2_l_star + rho_hat*(tau_nt2_hat*(u_t_star - u_t_l) &
3893 & + tau_nt_hat*(u_t2_star - u_t2_l))/sqrtc_nc
3894 tau_t1t2_r_starstar = tau_t1t2_r_star - rho_hat*(tau_nt2_hat*(u_t_star - u_t_r) &
3895 & + tau_nt_hat*(u_t2_star - u_t2_r))/sqrtc_nc
3896 e_l_starstar = e_l_star + (rhol_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_l*tau_nt_l) &
3897 & + (u_t2_star*tau_nt2_star - u_t2_l*tau_nt2_l))
3898 e_r_starstar = e_r_star - (rhor_star/sqrtc_nc)*((u_t_star*tau_nt_star - u_t_r*tau_nt_r) &
3899 & + (u_t2_star*tau_nt2_star - u_t2_r*tau_nt2_r))
3902 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)
3903 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)
3905 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)
3906 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)
3921 f_hlld(1:ncomp) = f_l(1:ncomp)
3922 else if (zone == 5)
then
3923 f_hlld(1:ncomp) = f_r(1:ncomp)
3924# 731 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3925 else if (zone == 1 .or. zone == 2)
then
3930 if (p > 0 .and. .not. cyl_coord)
then
3931# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3934 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
3935 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
3937# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3940 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
3941 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
3943# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3944 us_c = rhol_star*s_m
3946 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
3947 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
3949# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3950 us_c = rhol_star*u_t_l
3951 uss_c = rhol_star*u_t_star
3952 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
3953 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
3955# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3956 us_c = rhol_star*u_t2_l
3957 uss_c = rhol_star*u_t2_star
3958 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
3959 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
3961# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3963 uss_c = e_l_starstar
3964 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
3965 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
3967# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3968 us_c = alpha1_l_star
3970 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
3971 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
3973# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3974 us_c = alpha2_l_star
3976 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
3977 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
3979# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3980 us_c = rhol_star*tau_nn_l_star
3982 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
3983 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
3985# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3986 us_c = rhol_star*tau_nt_l
3987 uss_c = rhol_star*tau_nt_star
3988 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
3989 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
3991# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3992 us_c = rhol_star*tau_nt2_l
3993 uss_c = rhol_star*tau_nt2_star
3994 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
3995 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
3997# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
3998 us_c = rhol_star*tau_tt_l_star
3999 uss_c = rhol_star*tau_tt_l_starstar
4000 f_hlld(12) = f_l(12) + s_l*(us_c - u_l(12))
4001 if (zone == 2) f_hlld(12) = f_hlld(12) + s_lstar*(uss_c &
4003# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4004 us_c = rhol_star*tau_t2t2_l_star
4005 uss_c = rhol_star*tau_t2t2_l_starstar
4006 f_hlld(13) = f_l(13) + s_l*(us_c - u_l(13))
4007 if (zone == 2) f_hlld(13) = f_hlld(13) + s_lstar*(uss_c &
4009# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4010 us_c = rhol_star*tau_t1t2_l_star
4011 uss_c = rhol_star*tau_t1t2_l_starstar
4012 f_hlld(14) = f_l(14) + s_l*(us_c - u_l(14))
4013 if (zone == 2) f_hlld(14) = f_hlld(14) + s_lstar*(uss_c &
4015# 744 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4017# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4020 f_hlld(1) = f_l(1) + s_l*(us_c - u_l(1))
4021 if (zone == 2) f_hlld(1) = f_hlld(1) + s_lstar*(uss_c &
4023# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4026 f_hlld(2) = f_l(2) + s_l*(us_c - u_l(2))
4027 if (zone == 2) f_hlld(2) = f_hlld(2) + s_lstar*(uss_c &
4029# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4030 us_c = rhol_star*s_m
4032 f_hlld(3) = f_l(3) + s_l*(us_c - u_l(3))
4033 if (zone == 2) f_hlld(3) = f_hlld(3) + s_lstar*(uss_c &
4035# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4036 us_c = rhol_star*u_t_l
4037 uss_c = rhol_star*u_t_star
4038 f_hlld(4) = f_l(4) + s_l*(us_c - u_l(4))
4039 if (zone == 2) f_hlld(4) = f_hlld(4) + s_lstar*(uss_c &
4041# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4043 uss_c = e_l_starstar
4044 f_hlld(5) = f_l(5) + s_l*(us_c - u_l(5))
4045 if (zone == 2) f_hlld(5) = f_hlld(5) + s_lstar*(uss_c &
4047# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4048 us_c = alpha1_l_star
4050 f_hlld(6) = f_l(6) + s_l*(us_c - u_l(6))
4051 if (zone == 2) f_hlld(6) = f_hlld(6) + s_lstar*(uss_c &
4053# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4054 us_c = alpha2_l_star
4056 f_hlld(7) = f_l(7) + s_l*(us_c - u_l(7))
4057 if (zone == 2) f_hlld(7) = f_hlld(7) + s_lstar*(uss_c &
4059# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4060 us_c = rhol_star*tau_nn_l_star
4062 f_hlld(8) = f_l(8) + s_l*(us_c - u_l(8))
4063 if (zone == 2) f_hlld(8) = f_hlld(8) + s_lstar*(uss_c &
4065# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4066 us_c = rhol_star*tau_nt_l
4067 uss_c = rhol_star*tau_nt_star
4068 f_hlld(9) = f_l(9) + s_l*(us_c - u_l(9))
4069 if (zone == 2) f_hlld(9) = f_hlld(9) + s_lstar*(uss_c &
4071# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4072 us_c = rhol_star*tau_tt_l_star
4073 uss_c = rhol_star*tau_tt_l_starstar
4074 f_hlld(10) = f_l(10) + s_l*(us_c - u_l(10))
4075 if (zone == 2) f_hlld(10) = f_hlld(10) + s_lstar*(uss_c &
4077# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4078 us_c = rhol_star*tau_qq_l_star
4080 f_hlld(11) = f_l(11) + s_l*(us_c - u_l(11))
4081 if (zone == 2) f_hlld(11) = f_hlld(11) + s_lstar*(uss_c &
4083# 752 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4085# 731 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4086 else if (zone == 3 .or. zone == 4)
then
4091 if (p > 0 .and. .not. cyl_coord)
then
4092# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4095 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
4096 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
4098# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4101 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
4102 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
4104# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4105 us_c = rhor_star*s_m
4107 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
4108 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
4110# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4111 us_c = rhor_star*u_t_r
4112 uss_c = rhor_star*u_t_star
4113 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
4114 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
4116# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4117 us_c = rhor_star*u_t2_r
4118 uss_c = rhor_star*u_t2_star
4119 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
4120 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
4122# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4124 uss_c = e_r_starstar
4125 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
4126 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
4128# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4129 us_c = alpha1_r_star
4131 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
4132 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
4134# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4135 us_c = alpha2_r_star
4137 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
4138 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
4140# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4141 us_c = rhor_star*tau_nn_r_star
4143 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
4144 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
4146# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4147 us_c = rhor_star*tau_nt_r
4148 uss_c = rhor_star*tau_nt_star
4149 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
4150 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
4152# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4153 us_c = rhor_star*tau_nt2_r
4154 uss_c = rhor_star*tau_nt2_star
4155 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
4156 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
4158# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4159 us_c = rhor_star*tau_tt_r_star
4160 uss_c = rhor_star*tau_tt_r_starstar
4161 f_hlld(12) = f_r(12) + s_r*(us_c - u_r(12))
4162 if (zone == 3) f_hlld(12) = f_hlld(12) + s_rstar*(uss_c &
4164# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4165 us_c = rhor_star*tau_t2t2_r_star
4166 uss_c = rhor_star*tau_t2t2_r_starstar
4167 f_hlld(13) = f_r(13) + s_r*(us_c - u_r(13))
4168 if (zone == 3) f_hlld(13) = f_hlld(13) + s_rstar*(uss_c &
4170# 738 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4171 us_c = rhor_star*tau_t1t2_r_star
4172 uss_c = rhor_star*tau_t1t2_r_starstar
4173 f_hlld(14) = f_r(14) + s_r*(us_c - u_r(14))
4174 if (zone == 3) f_hlld(14) = f_hlld(14) + s_rstar*(uss_c &
4176# 744 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4178# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4181 f_hlld(1) = f_r(1) + s_r*(us_c - u_r(1))
4182 if (zone == 3) f_hlld(1) = f_hlld(1) + s_rstar*(uss_c &
4184# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4187 f_hlld(2) = f_r(2) + s_r*(us_c - u_r(2))
4188 if (zone == 3) f_hlld(2) = f_hlld(2) + s_rstar*(uss_c &
4190# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4191 us_c = rhor_star*s_m
4193 f_hlld(3) = f_r(3) + s_r*(us_c - u_r(3))
4194 if (zone == 3) f_hlld(3) = f_hlld(3) + s_rstar*(uss_c &
4196# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4197 us_c = rhor_star*u_t_r
4198 uss_c = rhor_star*u_t_star
4199 f_hlld(4) = f_r(4) + s_r*(us_c - u_r(4))
4200 if (zone == 3) f_hlld(4) = f_hlld(4) + s_rstar*(uss_c &
4202# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4204 uss_c = e_r_starstar
4205 f_hlld(5) = f_r(5) + s_r*(us_c - u_r(5))
4206 if (zone == 3) f_hlld(5) = f_hlld(5) + s_rstar*(uss_c &
4208# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4209 us_c = alpha1_r_star
4211 f_hlld(6) = f_r(6) + s_r*(us_c - u_r(6))
4212 if (zone == 3) f_hlld(6) = f_hlld(6) + s_rstar*(uss_c &
4214# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4215 us_c = alpha2_r_star
4217 f_hlld(7) = f_r(7) + s_r*(us_c - u_r(7))
4218 if (zone == 3) f_hlld(7) = f_hlld(7) + s_rstar*(uss_c &
4220# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4221 us_c = rhor_star*tau_nn_r_star
4223 f_hlld(8) = f_r(8) + s_r*(us_c - u_r(8))
4224 if (zone == 3) f_hlld(8) = f_hlld(8) + s_rstar*(uss_c &
4226# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4227 us_c = rhor_star*tau_nt_r
4228 uss_c = rhor_star*tau_nt_star
4229 f_hlld(9) = f_r(9) + s_r*(us_c - u_r(9))
4230 if (zone == 3) f_hlld(9) = f_hlld(9) + s_rstar*(uss_c &
4232# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4233 us_c = rhor_star*tau_tt_r_star
4234 uss_c = rhor_star*tau_tt_r_starstar
4235 f_hlld(10) = f_r(10) + s_r*(us_c - u_r(10))
4236 if (zone == 3) f_hlld(10) = f_hlld(10) + s_rstar*(uss_c &
4238# 746 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4239 us_c = rhor_star*tau_qq_r_star
4241 f_hlld(11) = f_r(11) + s_r*(us_c - u_r(11))
4242 if (zone == 3) f_hlld(11) = f_hlld(11) + s_rstar*(uss_c &
4244# 752 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4246# 754 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4251 if (riemann_hypo_adc)
then
4257 if (s_l < 0._wp .and. s_r > 0._wp)
then
4259 f_hll_c = (s_r*f_l(i) - s_l*f_r(i) + s_l*s_r*(u_r(i) - u_l(i)))/(s_r - s_l + verysmall)
4260 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
4265 f_hlld(i) = f_hll_c + phi*(f_hlld(i) - f_hll_c)
4275 if (cyl_coord .and. (riemann_hypo_adc .or. fan_fallback))
then
4276 if (0._wp <= s_l)
then
4277 u_n_hll_trace = u_n_l; u_t_hll_trace = u_t_l
4278 p_face_hll = pres%L; tau_qq_face_hll = tau_qq_l
4279 else if (s_r <= 0._wp)
then
4280 u_n_hll_trace = u_n_r; u_t_hll_trace = u_t_r
4281 p_face_hll = pres%R; tau_qq_face_hll = tau_qq_r
4283 u_n_hll_trace = (s_r*u_n_l - s_l*u_n_r)/(s_r - s_l + verysmall)
4284 u_t_hll_trace = (s_r*u_t_l - s_l*u_t_r)/(s_r - s_l + verysmall)
4286 u_hll_c = (s_r*u_r(1) - s_l*u_l(1) - (f_r(1) - f_l(1)))/(s_r - s_l + verysmall)
4288 u_hll_c = (s_r*u_r(2) - s_l*u_l(2) - (f_r(2) - f_l(2)))/(s_r - s_l + verysmall)
4289 rho_hll = rho_hll + u_hll_c
4290 u_hll_c = (s_r*u_r(3) - s_l*u_l(3) - (f_r(3) - f_l(3)))/(s_r - s_l + verysmall)
4291 u_n_hll_cons = u_hll_c/(rho_hll + verysmall)
4292 u_hll_c = (s_r*u_r(8) - s_l*u_l(8) - (f_r(8) - f_l(8)))/(s_r - s_l + verysmall)
4293 tau_nn_hll = u_hll_c/(rho_hll + verysmall)
4294 u_hll_c = (s_r*u_r(11) - s_l*u_l(11) - (f_r(11) - f_l(11)))/(s_r - s_l + verysmall)
4295 tau_qq_face_hll = u_hll_c/(rho_hll + verysmall)
4297 f_hll_c = (s_r*f_l(3) - s_l*f_r(3) + s_l*s_r*(u_r(3) - u_l(3)))/(s_r - s_l + verysmall)
4298 p_face_hll = f_hll_c - rho_hll*u_n_hll_cons*u_n_hll_cons + tau_nn_hll
4303# 811 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4304 if (ipass == 1)
then
4305 if (p > 0 .and. .not. cyl_coord)
then
4317 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(9)
4318 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
4319 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(12)
4320 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(11)
4321 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(14)
4322 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
4323 else if (
dir_idx(1) == 2)
then
4324 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
4325 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(10)
4326 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(9)
4327 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(14)
4328 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
4329 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(13)
4331 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(12)
4332 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(14)
4333 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(13)
4334 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 3) = f_hlld(10)
4335 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 4) = f_hlld(11)
4336 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 5) = f_hlld(9)
4347 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 1) = f_hlld(9)
4349 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(8)
4350 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(10)
4352 flux_rsx_vf(j, k, l, eqn_idx%stress%beg) = f_hlld(10)
4353 flux_rsx_vf(j, k, l, eqn_idx%stress%beg + 2) = f_hlld(8)
4356 flux_rsx_vf(j, k, l, eqn_idx%stress%end) = f_hlld(11)
4360# 811 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4361 if (ipass == 2)
then
4362 if (p > 0 .and. .not. cyl_coord)
then
4380 else if (
dir_idx(1) == 2)
then
4417# 868 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4421 if (fan_fallback)
then
4424 u_n_face = u_n_hll_trace; u_t_face = u_t_hll_trace
4430 u_n_face = u_n_l; u_t_face = u_t_l
4431 else if (zone == 1)
then
4432 u_n_face = s_m; u_t_face = u_t_l
4433 else if (zone <= 3)
then
4434 u_n_face = s_m; u_t_face = u_t_star
4435 else if (zone == 4)
then
4436 u_n_face = s_m; u_t_face = u_t_r
4438 u_n_face = u_n_r; u_t_face = u_t_r
4442 if (riemann_hypo_adc)
then
4443 u_n_face = u_n_hll_trace + phi*(u_n_face - u_n_hll_trace)
4444 u_t_face = u_t_hll_trace + phi*(u_t_face - u_t_hll_trace)
4447# 898 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4448 if (ipass == 1)
then
4457# 898 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4458 if (ipass == 2)
then
4467# 908 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4476# 959 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4485# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4486#if defined(MFC_OpenACC)
4487# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4489# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4490#elif defined(MFC_OpenMP)
4491# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4493# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4495# 966 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"
4498# 969 "/home/runner/work/MFC/MFC/src/simulation/m_riemann_solver_hypo_hlld.fpp"