9574 & flux_src_n_vf_arg, Kterm_arg)
9575 integer,
intent(in) :: current_idir
9576 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf_arg
9577 type(vector_field),
intent(in) :: q_cons_vf_arg
9578 type(vector_field),
intent(in) :: q_prim_vf_arg
9579 type(vector_field),
intent(in) :: flux_src_n_vf_arg
9580 real(wp),
allocatable,
dimension(:,:,:),
intent(in) :: Kterm_arg
9581 integer :: j_adv, k_idx, l_idx, q_idx
9582 real(wp) :: local_inv_ds, local_term_coeff, local_flux1, local_flux2
9583 real(wp) :: local_k_term_val
9593 select case (current_idir)
9595 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9598# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9600# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9601#if defined(MFC_OpenACC)
9602# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9604# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9605#elif defined(MFC_OpenMP)
9606# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9608# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9610# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9612# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9614# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9616# 1440 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9618# 1442 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9619 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9623 local_inv_ds = 1._wp/dx(k_idx)
9624 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(k_idx, l_idx, q_idx)
9625 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
9626 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9627 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9628 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9634# 1456 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9635#if defined(MFC_OpenACC)
9636# 1456 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9638# 1456 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9639#elif defined(MFC_OpenMP)
9640# 1456 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9642# 1456 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9644# 1456 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9646 if (alt_soundspeed)
then
9648# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9650# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9651#if defined(MFC_OpenACC)
9652# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9654# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9655#elif defined(MFC_OpenMP)
9656# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9658# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9660# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9662# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9664# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9666# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9668# 1460 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9669 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9670 local_inv_ds = 1._wp/dx(k_idx)
9671 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9673 local_flux2 =
nc_iface_vel_n(1)%vf(1)%sf(k_idx - 1, l_idx, q_idx)
9674 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, &
9675 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9676 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, &
9677 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9678 end do;
end do;
end do
9680# 1470 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9681#if defined(MFC_OpenACC)
9682# 1470 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9684# 1470 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9685#elif defined(MFC_OpenMP)
9686# 1470 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9688# 1470 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9690# 1470 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9693 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9696# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9698# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9699#if defined(MFC_OpenACC)
9700# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9702# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9703#elif defined(MFC_OpenMP)
9704# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9706# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9708# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9710# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9712# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9714# 1474 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9716# 1476 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9717 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9718 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9719 local_inv_ds = 1._wp/dx(k_idx)
9720 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9721 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
9722 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
9723 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9724 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9725 end do;
end do;
end do
9728# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9729#if defined(MFC_OpenACC)
9730# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9732# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9733#elif defined(MFC_OpenMP)
9734# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9736# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9738# 1486 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9740 if (alt_soundspeed)
then
9742# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9744# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9745#if defined(MFC_OpenACC)
9746# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9748# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9749#elif defined(MFC_OpenMP)
9750# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9752# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9754# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9756# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9758# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9760# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9762# 1490 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9763 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9764 local_inv_ds = 1._wp/dx(k_idx)
9765 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9766 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
9767 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
9768 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, &
9769 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9770 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, &
9771 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9772 end do;
end do;
end do
9774# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9775#if defined(MFC_OpenACC)
9776# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9778# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9779#elif defined(MFC_OpenMP)
9780# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9782# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9784# 1500 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9789 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9792# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9794# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9795#if defined(MFC_OpenACC)
9796# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9798# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9799#elif defined(MFC_OpenMP)
9800# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9802# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9804# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9806# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9808# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9810# 1506 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9812# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9813 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9817 local_inv_ds = 1._wp/dy(k_idx)
9818 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(q_idx, k_idx, l_idx)
9819 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
9820 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
9821 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
9822 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9828# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9829#if defined(MFC_OpenACC)
9830# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9832# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9833#elif defined(MFC_OpenMP)
9834# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9836# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9838# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9840 if (alt_soundspeed)
then
9842# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9844# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9845#if defined(MFC_OpenACC)
9846# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9848# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9849#elif defined(MFC_OpenMP)
9850# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9852# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9854# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9856# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9858# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9860# 1524 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9862# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9863 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9864 local_inv_ds = 1._wp/dy(k_idx)
9865 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
9867 local_flux2 =
nc_iface_vel_n(2)%vf(2)%sf(q_idx, k_idx - 1, l_idx)
9868 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, &
9869 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9870 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, &
9871 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9872 end do;
end do;
end do
9874# 1536 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9875#if defined(MFC_OpenACC)
9876# 1536 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9878# 1536 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9879#elif defined(MFC_OpenMP)
9880# 1536 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9882# 1536 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9884# 1536 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9888# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9890# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9891#if defined(MFC_OpenACC)
9892# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9894# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9895#elif defined(MFC_OpenMP)
9896# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9898# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9900# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9902# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9904# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9906# 1538 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9908 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9909 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
9910 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
9911 & k_idx, l_idx) + local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
9912 & k_idx, l_idx)/y_cc(k_idx)
9913 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
9914 & k_idx, l_idx) - local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
9915 & k_idx, l_idx)/y_cc(k_idx)
9916 end do;
end do;
end do
9918# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9919#if defined(MFC_OpenACC)
9920# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9922# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9923#elif defined(MFC_OpenMP)
9924# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9926# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9928# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9932 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9935# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9937# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9938#if defined(MFC_OpenACC)
9939# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9941# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9942#elif defined(MFC_OpenMP)
9943# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9945# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9947# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9949# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9951# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9953# 1553 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9955# 1555 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9956 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9957 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9958 local_inv_ds = 1._wp/dy(k_idx)
9959 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
9960 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
9961 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
9962 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
9963 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9964 end do;
end do;
end do
9967# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9968#if defined(MFC_OpenACC)
9969# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9971# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9972#elif defined(MFC_OpenMP)
9973# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9975# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9977# 1565 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9979 if (alt_soundspeed)
then
9981# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9983# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9984#if defined(MFC_OpenACC)
9985# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9987# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9988#elif defined(MFC_OpenMP)
9989# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9991# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9993# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9995# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9997# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9999# 1567 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10001# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10002 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10003 local_inv_ds = 1._wp/dy(k_idx)
10004 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10005 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10006 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10007 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, &
10008 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10009 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, &
10010 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10011 end do;
end do;
end do
10013# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10014#if defined(MFC_OpenACC)
10015# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10017# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10018#elif defined(MFC_OpenMP)
10019# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10021# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10023# 1579 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10025 if (cyl_coord)
then
10027# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10029# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10030#if defined(MFC_OpenACC)
10031# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10033# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10034#elif defined(MFC_OpenMP)
10035# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10037# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10039# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10041# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10043# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10045# 1581 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10047# 1583 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10048 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10049 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10050 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10051 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10052 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
10053 & k_idx, l_idx) + (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10054 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
10055 & k_idx, l_idx) - (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10056 end do;
end do;
end do
10058# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10059#if defined(MFC_OpenACC)
10060# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10062# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10063#elif defined(MFC_OpenMP)
10064# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10066# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10068# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10074 if (grid_geometry == 3)
then
10075 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10077# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10079# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10080#if defined(MFC_OpenACC)
10081# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10083# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10084#elif defined(MFC_OpenMP)
10085# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10087# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10089# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10091# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10093# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10095# 1599 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10097# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10098 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10102 local_inv_ds = 1._wp/dz(k_idx)
10103 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10104 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10105 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10106 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10107 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10113# 1615 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10114#if defined(MFC_OpenACC)
10115# 1615 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10117# 1615 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10118#elif defined(MFC_OpenMP)
10119# 1615 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10121# 1615 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10123# 1615 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10125 if (alt_soundspeed)
then
10127# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10129# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10130#if defined(MFC_OpenACC)
10131# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10133# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10134#elif defined(MFC_OpenMP)
10135# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10137# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10139# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10141# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10143# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10145# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10147# 1619 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10148 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10149 local_inv_ds = 1._wp/dz(k_idx)
10150 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10152 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10153 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10154 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10155 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10156 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10157 end do;
end do;
end do
10159# 1629 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10160#if defined(MFC_OpenACC)
10161# 1629 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10163# 1629 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10164#elif defined(MFC_OpenMP)
10165# 1629 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10167# 1629 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10169# 1629 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10172 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10175# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10177# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10178#if defined(MFC_OpenACC)
10179# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10181# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10182#elif defined(MFC_OpenMP)
10183# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10185# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10187# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10189# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10191# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10193# 1633 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10195# 1635 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10196 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10197 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10198 local_inv_ds = 1._wp/dz(k_idx)
10199 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10200 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10201 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10202 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10203 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10204 end do;
end do;
end do
10207# 1645 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10208#if defined(MFC_OpenACC)
10209# 1645 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10211# 1645 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10212#elif defined(MFC_OpenMP)
10213# 1645 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10215# 1645 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10217# 1645 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10219 if (alt_soundspeed)
then
10221# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10223# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10224#if defined(MFC_OpenACC)
10225# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10227# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10228#elif defined(MFC_OpenMP)
10229# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10231# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10233# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10235# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10237# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10239# 1647 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10241# 1649 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10242 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10243 local_inv_ds = 1._wp/dz(k_idx)
10244 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10245 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10246 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10247 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10248 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10249 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10250 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10251 end do;
end do;
end do
10253# 1659 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10254#if defined(MFC_OpenACC)
10255# 1659 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10257# 1659 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10258#elif defined(MFC_OpenMP)
10259# 1659 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10261# 1659 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10263# 1659 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10268 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10271# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10273# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10274#if defined(MFC_OpenACC)
10275# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10277# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10278#elif defined(MFC_OpenMP)
10279# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10281# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10283# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10285# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10287# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10289# 1665 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10291# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10292 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10296 local_inv_ds = 1._wp/dz(k_idx)
10297 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10298 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10299 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10300 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10301 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10307# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10308#if defined(MFC_OpenACC)
10309# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10311# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10312#elif defined(MFC_OpenMP)
10313# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10315# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10317# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10319 if (alt_soundspeed)
then
10321# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10323# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10324#if defined(MFC_OpenACC)
10325# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10327# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10328#elif defined(MFC_OpenMP)
10329# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10331# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10333# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10335# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10337# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10339# 1683 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10341# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10342 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10343 local_inv_ds = 1._wp/dz(k_idx)
10344 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10346 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10347 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10348 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10349 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10350 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10351 end do;
end do;
end do
10353# 1695 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10354#if defined(MFC_OpenACC)
10355# 1695 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10357# 1695 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10358#elif defined(MFC_OpenMP)
10359# 1695 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10361# 1695 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10363# 1695 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10366 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10369# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10371# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10372#if defined(MFC_OpenACC)
10373# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10375# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10376#elif defined(MFC_OpenMP)
10377# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10379# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10381# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10383# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10385# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10387# 1699 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10389# 1701 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10390 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10391 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10392 local_inv_ds = 1._wp/dz(k_idx)
10393 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10394 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10395 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10396 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10397 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10398 end do;
end do;
end do
10401# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10402#if defined(MFC_OpenACC)
10403# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10405# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10406#elif defined(MFC_OpenMP)
10407# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10409# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10411# 1711 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10413 if (alt_soundspeed)
then
10415# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10417# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10418#if defined(MFC_OpenACC)
10419# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10421# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10422#elif defined(MFC_OpenMP)
10423# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10425# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10427# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10429# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10431# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10433# 1713 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10435# 1715 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10436 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10437 local_inv_ds = 1._wp/dz(k_idx)
10438 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10439 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10440 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10441 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10442 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10443 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10444 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10445 end do;
end do;
end do
10447# 1725 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10448#if defined(MFC_OpenACC)
10449# 1725 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10451# 1725 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10452#elif defined(MFC_OpenMP)
10453# 1725 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10455# 1725 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10457# 1725 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"