7154 & flux_src_n_vf_arg, Kterm_arg)
7155 integer,
intent(in) :: current_idir
7156 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf_arg
7157 type(vector_field),
intent(in) :: q_cons_vf_arg
7158 type(vector_field),
intent(in) :: q_prim_vf_arg
7159 type(vector_field),
intent(in) :: flux_src_n_vf_arg
7161 real(wp),
allocatable,
dimension(:,:,:),
intent(in) :: Kterm_arg
7162 integer :: j_adv, k_idx, l_idx, q_idx
7163 real(wp) :: local_inv_ds, local_term_coeff, local_flux1, local_flux2
7164 real(wp) :: local_q_cons_val, local_k_term_val
7165 logical :: use_standard_riemann
7167 select case (current_idir)
7169 use_standard_riemann = (riemann_solver == 1 .or. riemann_solver == 4)
7170 if (use_standard_riemann)
then
7172# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7174# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7175#if defined(MFC_OpenACC)
7176# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7178# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7179#elif defined(MFC_OpenMP)
7180# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7182# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7184# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7186# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7188# 1118 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7190# 1120 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7191 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7195 local_inv_ds = 1._wp/dx(k_idx)
7196 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(k_idx, l_idx, q_idx)
7197 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
7198 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
7199 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
7200 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7206# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7207#if defined(MFC_OpenACC)
7208# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7210# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7211#elif defined(MFC_OpenMP)
7212# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7214# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7216# 1134 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7219 if (alt_soundspeed)
then
7220 if (bubbles_euler .neqv. .true.)
then
7222# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7224# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7225#if defined(MFC_OpenACC)
7226# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7228# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7229#elif defined(MFC_OpenMP)
7230# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7232# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7234# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7236# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7238# 1138 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7240# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7241 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
7242 local_inv_ds = 1._wp/dx(k_idx)
7243 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx)
7244 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
7245 local_term_coeff = local_q_cons_val - local_k_term_val
7246 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx)
7247 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(k_idx - 1, l_idx, q_idx)
7248 rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, &
7249 & l_idx, q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7250 end do;
end do;
end do
7252# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7253#if defined(MFC_OpenACC)
7254# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7256# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7257#elif defined(MFC_OpenMP)
7258# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7260# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7262# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7266# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7268# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7269#if defined(MFC_OpenACC)
7270# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7272# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7273#elif defined(MFC_OpenMP)
7274# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7276# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7278# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7280# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7282# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7284# 1154 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7285 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
7286 local_inv_ds = 1._wp/dx(k_idx)
7287 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
7288 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
7289 local_term_coeff = local_q_cons_val + local_k_term_val
7290 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
7291 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
7292 rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, &
7293 & l_idx, q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7294 end do;
end do;
end do
7296# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7297#if defined(MFC_OpenACC)
7298# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7300# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7301#elif defined(MFC_OpenMP)
7302# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7304# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7306# 1164 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7311# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7313# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7314#if defined(MFC_OpenACC)
7315# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7317# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7318#elif defined(MFC_OpenMP)
7319# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7321# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7323# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7325# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7327# 1167 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7329# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7330 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7331 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
7332 local_inv_ds = 1._wp/dx(k_idx)
7333 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
7334 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
7335 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
7336 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
7337 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7338 end do;
end do;
end do
7341# 1179 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7342#if defined(MFC_OpenACC)
7343# 1179 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7345# 1179 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7346#elif defined(MFC_OpenMP)
7347# 1179 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7349# 1179 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7351# 1179 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7357 use_standard_riemann = (riemann_solver == 1 .or. riemann_solver == 4)
7358 if (use_standard_riemann)
then
7360# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7362# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7363#if defined(MFC_OpenACC)
7364# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7366# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7367#elif defined(MFC_OpenMP)
7368# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7370# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7372# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7374# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7376# 1186 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7378# 1188 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7379 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7383 local_inv_ds = 1._wp/dy(k_idx)
7384 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(q_idx, k_idx, l_idx)
7385 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
7386 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
7387 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
7388 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7394# 1202 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7395#if defined(MFC_OpenACC)
7396# 1202 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7398# 1202 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7399#elif defined(MFC_OpenMP)
7400# 1202 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7402# 1202 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7404# 1202 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7407 if (alt_soundspeed)
then
7408 if (bubbles_euler .neqv. .true.)
then
7410# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7412# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7413#if defined(MFC_OpenACC)
7414# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7416# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7417#elif defined(MFC_OpenMP)
7418# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7420# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7422# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7424# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7426# 1206 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7428# 1208 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7429 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
7430 local_inv_ds = 1._wp/dy(k_idx)
7431 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx)
7432 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
7433 local_term_coeff = local_q_cons_val - local_k_term_val
7434 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx)
7435 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(q_idx, k_idx - 1, l_idx)
7436 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
7437 & k_idx, l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7439 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
7440 & k_idx, l_idx) - (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
7442 end do;
end do;
end do
7444# 1222 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7445#if defined(MFC_OpenACC)
7446# 1222 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7448# 1222 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7449#elif defined(MFC_OpenMP)
7450# 1222 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7452# 1222 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7454# 1222 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7458# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7460# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7461#if defined(MFC_OpenACC)
7462# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7464# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7465#elif defined(MFC_OpenMP)
7466# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7468# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7470# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7472# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7474# 1224 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7476# 1226 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7477 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
7478 local_inv_ds = 1._wp/dy(k_idx)
7479 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
7480 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
7481 local_term_coeff = local_q_cons_val + local_k_term_val
7482 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
7483 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
7484 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
7485 & k_idx, l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7487 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
7488 & k_idx, l_idx) + (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
7490 end do;
end do;
end do
7492# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7493#if defined(MFC_OpenACC)
7494# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7496# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7497#elif defined(MFC_OpenMP)
7498# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7500# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7502# 1240 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7507# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7509# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7510#if defined(MFC_OpenACC)
7511# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7513# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7514#elif defined(MFC_OpenMP)
7515# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7517# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7519# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7521# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7523# 1243 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7525# 1245 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7526 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7527 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
7528 local_inv_ds = 1._wp/dy(k_idx)
7529 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
7530 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
7531 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
7532 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
7533 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7534 end do;
end do;
end do
7537# 1255 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7538#if defined(MFC_OpenACC)
7539# 1255 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7541# 1255 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7542#elif defined(MFC_OpenMP)
7543# 1255 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7545# 1255 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7547# 1255 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7553 if (grid_geometry == 3)
then
7554 use_standard_riemann = (riemann_solver == 1)
7556 use_standard_riemann = (riemann_solver == 1 .or. riemann_solver == 4)
7559 if (use_standard_riemann)
then
7561# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7563# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7564#if defined(MFC_OpenACC)
7565# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7567# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7568#elif defined(MFC_OpenMP)
7569# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7571# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7573# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7575# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7577# 1267 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7579# 1269 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7580 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7584 local_inv_ds = 1._wp/dz(k_idx)
7585 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
7586 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
7587 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
7588 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
7589 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7595# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7596#if defined(MFC_OpenACC)
7597# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7599# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7600#elif defined(MFC_OpenMP)
7601# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7603# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7605# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7608 if (alt_soundspeed)
then
7609 if (bubbles_euler .neqv. .true.)
then
7611# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7613# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7614#if defined(MFC_OpenACC)
7615# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7617# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7618#elif defined(MFC_OpenMP)
7619# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7621# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7623# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7625# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7627# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7629# 1289 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7630 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
7631 local_inv_ds = 1._wp/dz(k_idx)
7632 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx)
7633 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
7634 local_term_coeff = local_q_cons_val - local_k_term_val
7635 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx)
7636 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx - 1)
7637 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
7638 & q_idx, k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7639 end do;
end do;
end do
7641# 1299 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7642#if defined(MFC_OpenACC)
7643# 1299 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7645# 1299 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7646#elif defined(MFC_OpenMP)
7647# 1299 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7649# 1299 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7651# 1299 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7655# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7657# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7658#if defined(MFC_OpenACC)
7659# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7661# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7662#elif defined(MFC_OpenMP)
7663# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7665# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7667# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7669# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7671# 1301 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7673# 1303 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7674 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
7675 local_inv_ds = 1._wp/dz(k_idx)
7676 local_q_cons_val = q_cons_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
7677 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
7678 local_term_coeff = local_q_cons_val + local_k_term_val
7679 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
7680 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
7681 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
7682 & q_idx, k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7683 end do;
end do;
end do
7685# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7686#if defined(MFC_OpenACC)
7687# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7689# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7690#elif defined(MFC_OpenMP)
7691# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7693# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7695# 1313 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7700# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7702# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7703#if defined(MFC_OpenACC)
7704# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7706# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7707#elif defined(MFC_OpenMP)
7708# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7710# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7712# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7714# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7716# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7718# 1318 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7719 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
7720 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
7721 local_inv_ds = 1._wp/dz(k_idx)
7722 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
7723 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
7724 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
7725 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
7726 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
7727 end do;
end do;
end do
7730# 1328 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7731#if defined(MFC_OpenACC)
7732# 1328 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7734# 1328 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7735#elif defined(MFC_OpenMP)
7736# 1328 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7738# 1328 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
7740# 1328 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"