9861 & flux_src_n_vf_arg, Kterm_arg)
9862 integer,
intent(in) :: current_idir
9863 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf_arg
9864 type(vector_field),
intent(in) :: q_cons_vf_arg
9865 type(vector_field),
intent(in) :: q_prim_vf_arg
9866 type(vector_field),
intent(in) :: flux_src_n_vf_arg
9867 real(wp),
allocatable,
dimension(:,:,:),
intent(in) :: Kterm_arg
9868 integer :: j_adv, k_idx, l_idx, q_idx
9869 real(wp) :: local_inv_ds, local_term_coeff, local_flux1, local_flux2
9870 real(wp) :: local_k_term_val
9880 select case (current_idir)
9882 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9885# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9887# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9888#if defined(MFC_OpenACC)
9889# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9891# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9892#elif defined(MFC_OpenMP)
9893# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9895# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9897# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9899# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9901# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9903# 1465 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9905# 1467 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9906 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9910 local_inv_ds = 1._wp/dx(k_idx)
9911 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(k_idx, l_idx, q_idx)
9912 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
9913 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9914 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9915 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9921# 1481 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9922#if defined(MFC_OpenACC)
9923# 1481 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9925# 1481 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9926#elif defined(MFC_OpenMP)
9927# 1481 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9929# 1481 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9931# 1481 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9933# 1482 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9934 if (alt_soundspeed)
then
9936# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9938# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9939#if defined(MFC_OpenACC)
9940# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9942# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9943#elif defined(MFC_OpenMP)
9944# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9946# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9948# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9950# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9952# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9954# 1483 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9956# 1485 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9957 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9958 local_inv_ds = 1._wp/dx(k_idx)
9959 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9961 local_flux2 =
nc_iface_vel_n(1)%vf(1)%sf(k_idx - 1, l_idx, q_idx)
9962 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, &
9963 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9964 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, &
9965 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9966 end do;
end do;
end do
9968# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9969#if defined(MFC_OpenACC)
9970# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9972# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9973#elif defined(MFC_OpenMP)
9974# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9976# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9978# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9980# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9982 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9985# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9987# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9988#if defined(MFC_OpenACC)
9989# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9991# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9992#elif defined(MFC_OpenMP)
9993# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9995# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9997# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9999# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10001# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10003# 1499 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10005# 1501 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10006 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10007 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
10008 local_inv_ds = 1._wp/dx(k_idx)
10009 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
10010 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
10011 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
10012 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
10013 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10014 end do;
end do;
end do
10017# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10018#if defined(MFC_OpenACC)
10019# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10021# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10022#elif defined(MFC_OpenMP)
10023# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10025# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10027# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10029# 1512 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10030 if (alt_soundspeed)
then
10032# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10034# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10035#if defined(MFC_OpenACC)
10036# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10038# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10039#elif defined(MFC_OpenMP)
10040# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10042# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10044# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10046# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10048# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10050# 1513 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10052# 1515 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10053 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
10054 local_inv_ds = 1._wp/dx(k_idx)
10055 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
10056 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
10057 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
10058 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, &
10059 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10060 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, &
10061 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10062 end do;
end do;
end do
10064# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10065#if defined(MFC_OpenACC)
10066# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10068# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10069#elif defined(MFC_OpenMP)
10070# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10072# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10074# 1525 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10076# 1526 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10080 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10083# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10085# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10086#if defined(MFC_OpenACC)
10087# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10089# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10090#elif defined(MFC_OpenMP)
10091# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10093# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10095# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10097# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10099# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10101# 1531 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10103# 1533 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10104 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10108 local_inv_ds = 1._wp/dy(k_idx)
10109 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(q_idx, k_idx, l_idx)
10110 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
10111 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
10112 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
10113 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10119# 1547 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10120#if defined(MFC_OpenACC)
10121# 1547 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10123# 1547 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10124#elif defined(MFC_OpenMP)
10125# 1547 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10127# 1547 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10129# 1547 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10131# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10132 if (alt_soundspeed)
then
10134# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10136# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10137#if defined(MFC_OpenACC)
10138# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10140# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10141#elif defined(MFC_OpenMP)
10142# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10144# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10146# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10148# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10150# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10152# 1549 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10154# 1551 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10155 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10156 local_inv_ds = 1._wp/dy(k_idx)
10157 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10159 local_flux2 =
nc_iface_vel_n(2)%vf(2)%sf(q_idx, k_idx - 1, l_idx)
10160 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, &
10161 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10162 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, &
10163 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10164 end do;
end do;
end do
10166# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10167#if defined(MFC_OpenACC)
10168# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10170# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10171#elif defined(MFC_OpenMP)
10172# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10174# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10176# 1561 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10178# 1562 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10179 if (cyl_coord)
then
10181# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10183# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10184#if defined(MFC_OpenACC)
10185# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10187# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10188#elif defined(MFC_OpenMP)
10189# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10191# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10193# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10195# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10197# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10199# 1563 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10201# 1564 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10202 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10203 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10204 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
10205 & k_idx, l_idx) + local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
10206 & k_idx, l_idx)/y_cc(k_idx)
10207 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
10208 & k_idx, l_idx) - local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
10209 & k_idx, l_idx)/y_cc(k_idx)
10210 end do;
end do;
end do
10212# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10213#if defined(MFC_OpenACC)
10214# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10216# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10217#elif defined(MFC_OpenMP)
10218# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10220# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10222# 1573 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10224# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10227 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10230# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10232# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10233#if defined(MFC_OpenACC)
10234# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10236# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10237#elif defined(MFC_OpenMP)
10238# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10240# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10242# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10244# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10246# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10248# 1578 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10250# 1580 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10251 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10252 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10253 local_inv_ds = 1._wp/dy(k_idx)
10254 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
10255 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10256 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10257 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
10258 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10259 end do;
end do;
end do
10262# 1590 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10263#if defined(MFC_OpenACC)
10264# 1590 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10266# 1590 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10267#elif defined(MFC_OpenMP)
10268# 1590 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10270# 1590 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10272# 1590 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10274# 1591 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10275 if (alt_soundspeed)
then
10277# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10279# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10280#if defined(MFC_OpenACC)
10281# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10283# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10284#elif defined(MFC_OpenMP)
10285# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10287# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10289# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10291# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10293# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10295# 1592 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10297# 1594 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10298 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10299 local_inv_ds = 1._wp/dy(k_idx)
10300 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10301 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10302 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10303 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, &
10304 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10305 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, &
10306 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10307 end do;
end do;
end do
10309# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10310#if defined(MFC_OpenACC)
10311# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10313# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10314#elif defined(MFC_OpenMP)
10315# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10317# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10319# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10321# 1605 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10322 if (cyl_coord)
then
10324# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10326# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10327#if defined(MFC_OpenACC)
10328# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10330# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10331#elif defined(MFC_OpenMP)
10332# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10334# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10336# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10338# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10340# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10342# 1606 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10344# 1608 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10345 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10346 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10347 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10348 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10349 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
10350 & k_idx, l_idx) + (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10351 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
10352 & k_idx, l_idx) - (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10353 end do;
end do;
end do
10355# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10356#if defined(MFC_OpenACC)
10357# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10359# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10360#elif defined(MFC_OpenMP)
10361# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10363# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10365# 1617 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10367# 1618 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10372 if (grid_geometry == 3)
then
10373 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10375# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10377# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10378#if defined(MFC_OpenACC)
10379# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10381# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10382#elif defined(MFC_OpenMP)
10383# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10385# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10387# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10389# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10391# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10393# 1624 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10395# 1626 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10396 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10400 local_inv_ds = 1._wp/(dz(k_idx)*y_cc(q_idx))
10401 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10402 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10403 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10404 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10405 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10411# 1640 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10412#if defined(MFC_OpenACC)
10413# 1640 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10415# 1640 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10416#elif defined(MFC_OpenMP)
10417# 1640 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10419# 1640 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10421# 1640 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10423# 1641 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10424 if (alt_soundspeed)
then
10426# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10428# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10429#if defined(MFC_OpenACC)
10430# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10432# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10433#elif defined(MFC_OpenMP)
10434# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10436# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10438# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10440# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10442# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10444# 1642 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10446# 1644 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10447 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10448 local_inv_ds = 1._wp/(dz(k_idx)*y_cc(q_idx))
10449 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10451 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10452 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10453 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10454 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10455 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10456 end do;
end do;
end do
10458# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10459#if defined(MFC_OpenACC)
10460# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10462# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10463#elif defined(MFC_OpenMP)
10464# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10466# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10468# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10470# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10472 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10475# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10477# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10478#if defined(MFC_OpenACC)
10479# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10481# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10482#elif defined(MFC_OpenMP)
10483# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10485# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10487# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10489# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10491# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10493# 1658 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10495# 1660 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10496 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10497 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10498 local_inv_ds = 1._wp/(dz(k_idx)*y_cc(q_idx))
10499 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10500 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10501 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10502 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10503 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10504 end do;
end do;
end do
10507# 1670 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10508#if defined(MFC_OpenACC)
10509# 1670 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10511# 1670 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10512#elif defined(MFC_OpenMP)
10513# 1670 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10515# 1670 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10517# 1670 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10519# 1671 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10520 if (alt_soundspeed)
then
10522# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10524# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10525#if defined(MFC_OpenACC)
10526# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10528# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10529#elif defined(MFC_OpenMP)
10530# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10532# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10534# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10536# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10538# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10540# 1672 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10542# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10543 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10544 local_inv_ds = 1._wp/(dz(k_idx)*y_cc(q_idx))
10545 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10546 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10547 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10548 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10549 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10550 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10551 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10552 end do;
end do;
end do
10554# 1684 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10555#if defined(MFC_OpenACC)
10556# 1684 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10558# 1684 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10559#elif defined(MFC_OpenMP)
10560# 1684 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10562# 1684 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10564# 1684 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10566# 1685 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10570 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10573# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10575# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10576#if defined(MFC_OpenACC)
10577# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10579# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10580#elif defined(MFC_OpenMP)
10581# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10583# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10585# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10587# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10589# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10591# 1690 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10593# 1692 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10594 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10598 local_inv_ds = 1._wp/dz(k_idx)
10599 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10600 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10601 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10602 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10603 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10609# 1706 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10610#if defined(MFC_OpenACC)
10611# 1706 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10613# 1706 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10614#elif defined(MFC_OpenMP)
10615# 1706 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10617# 1706 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10619# 1706 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10621# 1707 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10622 if (alt_soundspeed)
then
10624# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10626# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10627#if defined(MFC_OpenACC)
10628# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10630# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10631#elif defined(MFC_OpenMP)
10632# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10634# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10636# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10638# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10640# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10642# 1708 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10644# 1710 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10645 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10646 local_inv_ds = 1._wp/dz(k_idx)
10647 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10649 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10650 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10651 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10652 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10653 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10654 end do;
end do;
end do
10656# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10657#if defined(MFC_OpenACC)
10658# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10660# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10661#elif defined(MFC_OpenMP)
10662# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10664# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10666# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10668# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10670 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10673# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10675# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10676#if defined(MFC_OpenACC)
10677# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10679# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10680#elif defined(MFC_OpenMP)
10681# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10683# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10685# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10687# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10689# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10691# 1724 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10693# 1726 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10694 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10695 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10696 local_inv_ds = 1._wp/dz(k_idx)
10697 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10698 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10699 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10700 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10701 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10702 end do;
end do;
end do
10705# 1736 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10706#if defined(MFC_OpenACC)
10707# 1736 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10709# 1736 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10710#elif defined(MFC_OpenMP)
10711# 1736 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10713# 1736 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10715# 1736 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10717# 1737 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10718 if (alt_soundspeed)
then
10720# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10722# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10723#if defined(MFC_OpenACC)
10724# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10726# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10727#elif defined(MFC_OpenMP)
10728# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10730# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10732# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10734# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10736# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10738# 1738 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10740# 1740 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10741 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10742 local_inv_ds = 1._wp/dz(k_idx)
10743 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10744 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10745 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10746 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10747 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10748 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10749 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10750 end do;
end do;
end do
10752# 1750 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10753#if defined(MFC_OpenACC)
10754# 1750 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10756# 1750 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10757#elif defined(MFC_OpenMP)
10758# 1750 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10760# 1750 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10762# 1750 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10764# 1751 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"