7181 & flux_src_n_vf_arg, Kterm_arg)
7182 integer,
intent(in) :: current_idir
7183 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf_arg
7184 type(vector_field),
intent(in) :: q_cons_vf_arg
7185 type(vector_field),
intent(in) :: q_prim_vf_arg
7186 type(vector_field),
intent(in) :: flux_src_n_vf_arg
7188 real(wp),
allocatable,
dimension(:,:,:),
intent(in) :: Kterm_arg
7189 integer :: j_adv, k_idx, l_idx, q_idx
7190 real(wp) :: local_inv_ds, local_term_coeff, local_flux1, local_flux2
7191 real(wp) :: local_q_cons_val, local_k_term_val
7192 logical :: use_standard_riemann
7194 select case (current_idir)
7196 use_standard_riemann = (riemann_solver == riemann_solver_hll .or. riemann_solver == riemann_solver_hlld)
7197 if (use_standard_riemann)
then
7199# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7201# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7202#if defined(MFC_OpenACC)
7203# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7205# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7206#elif defined(MFC_OpenMP)
7207# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7209# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7211# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7213# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7215# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7217# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7219# 1121 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7220 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7224 local_inv_ds = 1._wp/dx(k_idx)
7225 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(k_idx, l_idx, q_idx)
7226 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
7227 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
7228 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
7229 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7235# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7236#if defined(MFC_OpenACC)
7237# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7239# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7240#elif defined(MFC_OpenMP)
7241# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7243# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7245# 1135 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7248 if (alt_soundspeed)
then
7249 if (bubbles_euler .neqv. .true.)
then
7251# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7253# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7254#if defined(MFC_OpenACC)
7255# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7257# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7258#elif defined(MFC_OpenMP)
7259# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7261# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7263# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7265# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7267# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7269# 1139 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7271# 1141 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7272 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
7273 local_inv_ds = 1._wp/dx(k_idx)
7274 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx)
7275 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
7276 local_term_coeff = local_q_cons_val - local_k_term_val
7277 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx)
7278 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(k_idx - 1, l_idx, q_idx)
7279 rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, &
7280 & l_idx, q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7281 end do;
end do;
end do
7283# 1151 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7284#if defined(MFC_OpenACC)
7285# 1151 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7287# 1151 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7288#elif defined(MFC_OpenMP)
7289# 1151 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7291# 1151 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7293# 1151 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7297# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7299# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7300#if defined(MFC_OpenACC)
7301# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7303# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7304#elif defined(MFC_OpenMP)
7305# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7307# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7309# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7311# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7313# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7315# 1153 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7317# 1155 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7318 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
7319 local_inv_ds = 1._wp/dx(k_idx)
7320 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
7321 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
7322 local_term_coeff = local_q_cons_val + local_k_term_val
7323 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
7324 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
7325 rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, &
7326 & l_idx, q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7327 end do;
end do;
end do
7329# 1165 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7330#if defined(MFC_OpenACC)
7331# 1165 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7333# 1165 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7334#elif defined(MFC_OpenMP)
7335# 1165 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7337# 1165 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7339# 1165 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7344# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7346# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7347#if defined(MFC_OpenACC)
7348# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7350# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7351#elif defined(MFC_OpenMP)
7352# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7354# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7356# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7358# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7360# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7362# 1168 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7364# 1170 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7365 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7366 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
7367 local_inv_ds = 1._wp/dx(k_idx)
7368 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
7369 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
7370 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
7371 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
7372 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7373 end do;
end do;
end do
7376# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7377#if defined(MFC_OpenACC)
7378# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7380# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7381#elif defined(MFC_OpenMP)
7382# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7384# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7386# 1180 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7392 use_standard_riemann = (riemann_solver == riemann_solver_hll .or. riemann_solver == riemann_solver_hlld)
7393 if (use_standard_riemann)
then
7395# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7397# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7398#if defined(MFC_OpenACC)
7399# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7401# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7402#elif defined(MFC_OpenMP)
7403# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7405# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7407# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7409# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7411# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7413# 1187 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7415# 1189 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7416 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7420 local_inv_ds = 1._wp/dy(k_idx)
7421 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(q_idx, k_idx, l_idx)
7422 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
7423 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
7424 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
7425 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7431# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7432#if defined(MFC_OpenACC)
7433# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7435# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7436#elif defined(MFC_OpenMP)
7437# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7439# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7441# 1203 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7444 if (alt_soundspeed)
then
7445 if (bubbles_euler .neqv. .true.)
then
7447# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7449# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7450#if defined(MFC_OpenACC)
7451# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7453# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7454#elif defined(MFC_OpenMP)
7455# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7457# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7459# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7461# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7463# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7465# 1207 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7467# 1209 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7468 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
7469 local_inv_ds = 1._wp/dy(k_idx)
7470 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx)
7471 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
7472 local_term_coeff = local_q_cons_val - local_k_term_val
7473 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx)
7474 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(q_idx, k_idx - 1, l_idx)
7475 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
7476 & k_idx, l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7478 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
7479 & k_idx, l_idx) - (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
7481 end do;
end do;
end do
7483# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7484#if defined(MFC_OpenACC)
7485# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7487# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7488#elif defined(MFC_OpenMP)
7489# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7491# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7493# 1223 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7497# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7499# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7500#if defined(MFC_OpenACC)
7501# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7503# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7504#elif defined(MFC_OpenMP)
7505# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7507# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7509# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7511# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7513# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7515# 1225 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7517# 1227 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7518 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
7519 local_inv_ds = 1._wp/dy(k_idx)
7520 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
7521 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
7522 local_term_coeff = local_q_cons_val + local_k_term_val
7523 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
7524 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
7525 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
7526 & k_idx, l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7528 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
7529 & k_idx, l_idx) + (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
7531 end do;
end do;
end do
7533# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7534#if defined(MFC_OpenACC)
7535# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7537# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7538#elif defined(MFC_OpenMP)
7539# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7541# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7543# 1241 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7548# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7550# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7551#if defined(MFC_OpenACC)
7552# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7554# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7555#elif defined(MFC_OpenMP)
7556# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7558# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7560# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7562# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7564# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7566# 1244 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7568# 1246 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7569 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7570 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
7571 local_inv_ds = 1._wp/dy(k_idx)
7572 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
7573 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
7574 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
7575 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
7576 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7577 end do;
end do;
end do
7580# 1256 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7581#if defined(MFC_OpenACC)
7582# 1256 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7584# 1256 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7585#elif defined(MFC_OpenMP)
7586# 1256 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7588# 1256 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7590# 1256 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7596 if (grid_geometry == 3)
then
7597 use_standard_riemann = (riemann_solver == riemann_solver_hll)
7599 use_standard_riemann = (riemann_solver == riemann_solver_hll .or. riemann_solver == riemann_solver_hlld)
7602 if (use_standard_riemann)
then
7604# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7606# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7607#if defined(MFC_OpenACC)
7608# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7610# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7611#elif defined(MFC_OpenMP)
7612# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7614# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7616# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7618# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7620# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7622# 1268 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7624# 1270 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7625 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7629 local_inv_ds = 1._wp/dz(k_idx)
7630 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
7631 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
7632 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
7633 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
7634 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7640# 1284 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7641#if defined(MFC_OpenACC)
7642# 1284 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7644# 1284 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7645#elif defined(MFC_OpenMP)
7646# 1284 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7648# 1284 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7650# 1284 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7653 if (alt_soundspeed)
then
7654 if (bubbles_euler .neqv. .true.)
then
7656# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7658# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7659#if defined(MFC_OpenACC)
7660# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7662# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7663#elif defined(MFC_OpenMP)
7664# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7666# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7668# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7670# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7672# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7674# 1288 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7676# 1290 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7677 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
7678 local_inv_ds = 1._wp/dz(k_idx)
7679 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx)
7680 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
7681 local_term_coeff = local_q_cons_val - local_k_term_val
7682 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx)
7683 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx - 1)
7684 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
7685 & q_idx, k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7686 end do;
end do;
end do
7688# 1300 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7689#if defined(MFC_OpenACC)
7690# 1300 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7692# 1300 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7693#elif defined(MFC_OpenMP)
7694# 1300 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7696# 1300 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7698# 1300 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7702# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7704# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7705#if defined(MFC_OpenACC)
7706# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7708# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7709#elif defined(MFC_OpenMP)
7710# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7712# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7714# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7716# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7718# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7720# 1302 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7722# 1304 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7723 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
7724 local_inv_ds = 1._wp/dz(k_idx)
7725 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
7726 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
7727 local_term_coeff = local_q_cons_val + local_k_term_val
7728 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
7729 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
7730 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
7731 & q_idx, k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7732 end do;
end do;
end do
7734# 1314 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7735#if defined(MFC_OpenACC)
7736# 1314 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7738# 1314 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7739#elif defined(MFC_OpenMP)
7740# 1314 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7742# 1314 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7744# 1314 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7749# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7751# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7752#if defined(MFC_OpenACC)
7753# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7755# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7756#elif defined(MFC_OpenMP)
7757# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7759# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7761# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7763# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7765# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7767# 1317 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7769# 1319 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7770 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7771 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
7772 local_inv_ds = 1._wp/dz(k_idx)
7773 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
7774 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
7775 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
7776 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
7777 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7778 end do;
end do;
end do
7781# 1329 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7782#if defined(MFC_OpenACC)
7783# 1329 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7785# 1329 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7786#elif defined(MFC_OpenMP)
7787# 1329 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7789# 1329 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7791# 1329 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"