9600 & flux_src_n_vf_arg, Kterm_arg)
9601 integer,
intent(in) :: current_idir
9602 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf_arg
9603 type(vector_field),
intent(in) :: q_cons_vf_arg
9604 type(vector_field),
intent(in) :: q_prim_vf_arg
9605 type(vector_field),
intent(in) :: flux_src_n_vf_arg
9606 real(wp),
allocatable,
dimension(:,:,:),
intent(in) :: Kterm_arg
9607 integer :: j_adv, k_idx, l_idx, q_idx
9608 real(wp) :: local_inv_ds, local_term_coeff, local_flux1, local_flux2
9609 real(wp) :: local_k_term_val
9619 select case (current_idir)
9621 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9624# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9626# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9627#if defined(MFC_OpenACC)
9628# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9630# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9631#elif defined(MFC_OpenMP)
9632# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9634# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9636# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9638# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9640# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9642# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9644# 1444 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9645 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9649 local_inv_ds = 1._wp/dx(k_idx)
9650 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(k_idx, l_idx, q_idx)
9651 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
9652 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9653 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9654 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9660# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9661#if defined(MFC_OpenACC)
9662# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9664# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9665#elif defined(MFC_OpenMP)
9666# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9668# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9670# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9672 if (alt_soundspeed)
then
9674# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9676# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9677#if defined(MFC_OpenACC)
9678# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9680# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9681#elif defined(MFC_OpenMP)
9682# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9684# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9686# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9688# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9690# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9692# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9694# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9695 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9696 local_inv_ds = 1._wp/dx(k_idx)
9697 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9699 local_flux2 =
nc_iface_vel_n(1)%vf(1)%sf(k_idx - 1, l_idx, q_idx)
9700 rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, l_idx, &
9701 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9702 rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, l_idx, &
9703 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9704 end do;
end do;
end do
9706# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9707#if defined(MFC_OpenACC)
9708# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9710# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9711#elif defined(MFC_OpenMP)
9712# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9714# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9716# 1472 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9719 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9722# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9724# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9725#if defined(MFC_OpenACC)
9726# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9728# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9729#elif defined(MFC_OpenMP)
9730# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9732# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9734# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9736# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9738# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9740# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9742# 1478 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9743 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9744 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9745 local_inv_ds = 1._wp/dx(k_idx)
9746 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9747 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
9748 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
9749 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9750 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9751 end do;
end do;
end do
9754# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9755#if defined(MFC_OpenACC)
9756# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9758# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9759#elif defined(MFC_OpenMP)
9760# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9762# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9764# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9766 if (alt_soundspeed)
then
9768# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9770# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9771#if defined(MFC_OpenACC)
9772# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9774# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9775#elif defined(MFC_OpenMP)
9776# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9778# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9780# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9782# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9784# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9786# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9788# 1492 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9789 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9790 local_inv_ds = 1._wp/dx(k_idx)
9791 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9792 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
9793 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
9794 rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(k_idx, l_idx, &
9795 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9796 rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(k_idx, l_idx, &
9797 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9798 end do;
end do;
end do
9800# 1502 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9801#if defined(MFC_OpenACC)
9802# 1502 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9804# 1502 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9805#elif defined(MFC_OpenMP)
9806# 1502 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9808# 1502 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9810# 1502 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9815 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9818# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9820# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9821#if defined(MFC_OpenACC)
9822# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9824# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9825#elif defined(MFC_OpenMP)
9826# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9828# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9830# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9832# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9834# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9836# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9838# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9839 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9843 local_inv_ds = 1._wp/dy(k_idx)
9844 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(q_idx, k_idx, l_idx)
9845 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
9846 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
9847 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
9848 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9854# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9855#if defined(MFC_OpenACC)
9856# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9858# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9859#elif defined(MFC_OpenMP)
9860# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9862# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9864# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9866 if (alt_soundspeed)
then
9868# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9870# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9871#if defined(MFC_OpenACC)
9872# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9874# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9875#elif defined(MFC_OpenMP)
9876# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9878# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9880# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9882# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9884# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9886# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9888# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9889 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9890 local_inv_ds = 1._wp/dy(k_idx)
9891 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
9893 local_flux2 =
nc_iface_vel_n(2)%vf(2)%sf(q_idx, k_idx - 1, l_idx)
9894 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, &
9895 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9896 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, &
9897 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9898 end do;
end do;
end do
9900# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9901#if defined(MFC_OpenACC)
9902# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9904# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9905#elif defined(MFC_OpenMP)
9906# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9908# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9910# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9914# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9916# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9917#if defined(MFC_OpenACC)
9918# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9920# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9921#elif defined(MFC_OpenMP)
9922# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9924# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9926# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9928# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9930# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9932# 1540 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9934 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9935 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
9936 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
9937 & k_idx, l_idx) + local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
9938 & k_idx, l_idx)/y_cc(k_idx)
9939 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
9940 & k_idx, l_idx) - local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
9941 & k_idx, l_idx)/y_cc(k_idx)
9942 end do;
end do;
end do
9944# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9945#if defined(MFC_OpenACC)
9946# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9948# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9949#elif defined(MFC_OpenMP)
9950# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9952# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9954# 1550 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9958 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9961# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9963# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9964#if defined(MFC_OpenACC)
9965# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9967# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9968#elif defined(MFC_OpenMP)
9969# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9971# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9973# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9975# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9977# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9979# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9981# 1557 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9982 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9983 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9984 local_inv_ds = 1._wp/dy(k_idx)
9985 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
9986 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
9987 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
9988 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
9989 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9990 end do;
end do;
end do
9993# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9994#if defined(MFC_OpenACC)
9995# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9997# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9998#elif defined(MFC_OpenMP)
9999# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10001# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10003# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10005 if (alt_soundspeed)
then
10007# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10009# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10010#if defined(MFC_OpenACC)
10011# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10013# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10014#elif defined(MFC_OpenMP)
10015# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10017# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10019# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10021# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10023# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10025# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10027# 1571 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10028 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10029 local_inv_ds = 1._wp/dy(k_idx)
10030 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10031 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10032 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10033 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, &
10034 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10035 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, &
10036 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10037 end do;
end do;
end do
10039# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10040#if defined(MFC_OpenACC)
10041# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10043# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10044#elif defined(MFC_OpenMP)
10045# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10047# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10049# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10051 if (cyl_coord)
then
10053# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10055# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10056#if defined(MFC_OpenACC)
10057# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10059# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10060#elif defined(MFC_OpenMP)
10061# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10063# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10065# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10067# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10069# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10071# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10073# 1585 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10074 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10075 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10076 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10077 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10078 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
10079 & k_idx, l_idx) + (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10080 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
10081 & k_idx, l_idx) - (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10082 end do;
end do;
end do
10084# 1594 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10085#if defined(MFC_OpenACC)
10086# 1594 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10088# 1594 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10089#elif defined(MFC_OpenMP)
10090# 1594 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10092# 1594 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10094# 1594 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10100 if (grid_geometry == 3)
then
10101 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10103# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10105# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10106#if defined(MFC_OpenACC)
10107# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10109# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10110#elif defined(MFC_OpenMP)
10111# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10113# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10115# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10117# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10119# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10121# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10123# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10124 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10128 local_inv_ds = 1._wp/dz(k_idx)
10129 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10130 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10131 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10132 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10133 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10139# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10140#if defined(MFC_OpenACC)
10141# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10143# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10144#elif defined(MFC_OpenMP)
10145# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10147# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10149# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10151 if (alt_soundspeed)
then
10153# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10155# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10156#if defined(MFC_OpenACC)
10157# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10159# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10160#elif defined(MFC_OpenMP)
10161# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10163# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10165# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10167# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10169# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10171# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10173# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10174 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10175 local_inv_ds = 1._wp/dz(k_idx)
10176 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10178 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10179 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10180 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10181 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10182 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10183 end do;
end do;
end do
10185# 1631 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10186#if defined(MFC_OpenACC)
10187# 1631 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10189# 1631 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10190#elif defined(MFC_OpenMP)
10191# 1631 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10193# 1631 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10195# 1631 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10198 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10201# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10203# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10204#if defined(MFC_OpenACC)
10205# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10207# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10208#elif defined(MFC_OpenMP)
10209# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10211# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10213# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10215# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10217# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10219# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10221# 1637 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10222 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10223 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10224 local_inv_ds = 1._wp/dz(k_idx)
10225 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10226 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10227 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10228 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10229 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10230 end do;
end do;
end do
10233# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10234#if defined(MFC_OpenACC)
10235# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10237# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10238#elif defined(MFC_OpenMP)
10239# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10241# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10243# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10245 if (alt_soundspeed)
then
10247# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10249# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10250#if defined(MFC_OpenACC)
10251# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10253# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10254#elif defined(MFC_OpenMP)
10255# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10257# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10259# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10261# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10263# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10265# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10267# 1651 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10268 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10269 local_inv_ds = 1._wp/dz(k_idx)
10270 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10271 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10272 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10273 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10274 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10275 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10276 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10277 end do;
end do;
end do
10279# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10280#if defined(MFC_OpenACC)
10281# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10283# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10284#elif defined(MFC_OpenMP)
10285# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10287# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10289# 1661 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10294 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10297# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10299# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10300#if defined(MFC_OpenACC)
10301# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10303# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10304#elif defined(MFC_OpenMP)
10305# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10307# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10309# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10311# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10313# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10315# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10317# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10318 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10322 local_inv_ds = 1._wp/dz(k_idx)
10323 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10324 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10325 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10326 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10327 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10333# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10334#if defined(MFC_OpenACC)
10335# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10337# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10338#elif defined(MFC_OpenMP)
10339# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10341# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10343# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10345 if (alt_soundspeed)
then
10347# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10349# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10350#if defined(MFC_OpenACC)
10351# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10353# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10354#elif defined(MFC_OpenMP)
10355# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10357# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10359# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10361# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10363# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10365# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10367# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10368 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10369 local_inv_ds = 1._wp/dz(k_idx)
10370 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10372 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10373 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10374 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10375 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10376 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10377 end do;
end do;
end do
10379# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10380#if defined(MFC_OpenACC)
10381# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10383# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10384#elif defined(MFC_OpenMP)
10385# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10387# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10389# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10392 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10395# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10397# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10398#if defined(MFC_OpenACC)
10399# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10401# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10402#elif defined(MFC_OpenMP)
10403# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10405# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10407# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10409# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10411# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10413# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10415# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10416 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10417 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10418 local_inv_ds = 1._wp/dz(k_idx)
10419 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10420 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10421 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10422 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10423 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10424 end do;
end do;
end do
10427# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10428#if defined(MFC_OpenACC)
10429# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10431# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10432#elif defined(MFC_OpenMP)
10433# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10435# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10437# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10439 if (alt_soundspeed)
then
10441# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10443# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10444#if defined(MFC_OpenACC)
10445# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10447# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10448#elif defined(MFC_OpenMP)
10449# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10451# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10453# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10455# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10457# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10459# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10461# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10462 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10463 local_inv_ds = 1._wp/dz(k_idx)
10464 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10465 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10466 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10467 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10468 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10469 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10470 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10471 end do;
end do;
end do
10473# 1727 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10474#if defined(MFC_OpenACC)
10475# 1727 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10477# 1727 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10478#elif defined(MFC_OpenMP)
10479# 1727 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10481# 1727 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10483# 1727 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"