9814 & flux_src_n_vf_arg, Kterm_arg)
9815 integer,
intent(in) :: current_idir
9816 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf_arg
9817 type(vector_field),
intent(in) :: q_cons_vf_arg
9818 type(vector_field),
intent(in) :: q_prim_vf_arg
9819 type(vector_field),
intent(in) :: flux_src_n_vf_arg
9820 real(wp),
allocatable,
dimension(:,:,:),
intent(in) :: Kterm_arg
9821 integer :: j_adv, k_idx, l_idx, q_idx
9822 real(wp) :: local_inv_ds, local_term_coeff, local_flux1, local_flux2
9823 real(wp) :: local_k_term_val
9833 select case (current_idir)
9835 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9838# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9840# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9841#if defined(MFC_OpenACC)
9842# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9844# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9845#elif defined(MFC_OpenMP)
9846# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9848# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9850# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9852# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9854# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9856# 1462 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9858# 1464 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9859 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9863 local_inv_ds = 1._wp/dx(k_idx)
9864 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(k_idx, l_idx, q_idx)
9865 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
9866 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9867 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9868 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9874# 1478 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9875#if defined(MFC_OpenACC)
9876# 1478 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9878# 1478 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9879#elif defined(MFC_OpenMP)
9880# 1478 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9882# 1478 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9884# 1478 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9886# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9887 if (alt_soundspeed)
then
9889# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9891# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9892#if defined(MFC_OpenACC)
9893# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9895# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9896#elif defined(MFC_OpenMP)
9897# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9899# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9901# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9903# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9905# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9907# 1480 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9909# 1482 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9910 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9911 local_inv_ds = 1._wp/dx(k_idx)
9912 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9914 local_flux2 =
nc_iface_vel_n(1)%vf(1)%sf(k_idx - 1, l_idx, q_idx)
9915 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, &
9916 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9917 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, &
9918 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9919 end do;
end do;
end do
9921# 1492 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9922#if defined(MFC_OpenACC)
9923# 1492 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9925# 1492 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9926#elif defined(MFC_OpenMP)
9927# 1492 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9929# 1492 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9931# 1492 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9933# 1493 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9935 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9938# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9940# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9941#if defined(MFC_OpenACC)
9942# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9944# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9945#elif defined(MFC_OpenMP)
9946# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9948# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9950# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9952# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9954# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9956# 1496 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9958# 1498 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9959 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9960 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9961 local_inv_ds = 1._wp/dx(k_idx)
9962 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9963 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
9964 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
9965 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9966 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9967 end do;
end do;
end do
9970# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9971#if defined(MFC_OpenACC)
9972# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9974# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9975#elif defined(MFC_OpenMP)
9976# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9978# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9980# 1508 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9982# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9983 if (alt_soundspeed)
then
9985# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9987# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9988#if defined(MFC_OpenACC)
9989# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9991# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9992#elif defined(MFC_OpenMP)
9993# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9995# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9997# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9999# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10001# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10003# 1510 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10005# 1512 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10006 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
10007 local_inv_ds = 1._wp/dx(k_idx)
10008 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
10009 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
10010 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
10011 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, &
10012 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10013 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, &
10014 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10015 end do;
end do;
end do
10017# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10018#if defined(MFC_OpenACC)
10019# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10021# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10022#elif defined(MFC_OpenMP)
10023# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10025# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10027# 1522 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10029# 1523 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10033 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10036# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10038# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10039#if defined(MFC_OpenACC)
10040# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10042# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10043#elif defined(MFC_OpenMP)
10044# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10046# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10048# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10050# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10052# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10054# 1528 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10056# 1530 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10057 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10061 local_inv_ds = 1._wp/dy(k_idx)
10062 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(q_idx, k_idx, l_idx)
10063 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
10064 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
10065 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
10066 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10072# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10073#if defined(MFC_OpenACC)
10074# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10076# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10077#elif defined(MFC_OpenMP)
10078# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10080# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10082# 1544 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10084# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10085 if (alt_soundspeed)
then
10087# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10089# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10090#if defined(MFC_OpenACC)
10091# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10093# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10094#elif defined(MFC_OpenMP)
10095# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10097# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10099# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10101# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10103# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10105# 1546 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10107# 1548 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10108 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10109 local_inv_ds = 1._wp/dy(k_idx)
10110 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10112 local_flux2 =
nc_iface_vel_n(2)%vf(2)%sf(q_idx, k_idx - 1, l_idx)
10113 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, &
10114 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10115 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, &
10116 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10117 end do;
end do;
end do
10119# 1558 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10120#if defined(MFC_OpenACC)
10121# 1558 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10123# 1558 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10124#elif defined(MFC_OpenMP)
10125# 1558 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10127# 1558 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10129# 1558 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10131# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10132 if (cyl_coord)
then
10134# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10136# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10137#if defined(MFC_OpenACC)
10138# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10140# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10141#elif defined(MFC_OpenMP)
10142# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10144# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10146# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10148# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10150# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10152# 1560 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10154# 1561 "/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_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10157 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
10158 & k_idx, l_idx) + local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
10159 & k_idx, l_idx)/y_cc(k_idx)
10160 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
10161 & k_idx, l_idx) - local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
10162 & k_idx, l_idx)/y_cc(k_idx)
10163 end do;
end do;
end do
10165# 1570 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10166#if defined(MFC_OpenACC)
10167# 1570 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10169# 1570 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10170#elif defined(MFC_OpenMP)
10171# 1570 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10173# 1570 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10175# 1570 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10177# 1571 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10180 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10183# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10185# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10186#if defined(MFC_OpenACC)
10187# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10189# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10190#elif defined(MFC_OpenMP)
10191# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10193# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10195# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10197# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10199# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10201# 1575 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10203# 1577 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10204 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10205 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10206 local_inv_ds = 1._wp/dy(k_idx)
10207 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
10208 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10209 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10210 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
10211 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10212 end do;
end do;
end do
10215# 1587 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10216#if defined(MFC_OpenACC)
10217# 1587 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10219# 1587 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10220#elif defined(MFC_OpenMP)
10221# 1587 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10223# 1587 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10225# 1587 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10227# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10228 if (alt_soundspeed)
then
10230# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10232# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10233#if defined(MFC_OpenACC)
10234# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10236# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10237#elif defined(MFC_OpenMP)
10238# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10240# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10242# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10244# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10246# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10248# 1589 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10250# 1591 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10251 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10252 local_inv_ds = 1._wp/dy(k_idx)
10253 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10254 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10255 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10256 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, &
10257 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10258 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, &
10259 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10260 end do;
end do;
end do
10262# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10263#if defined(MFC_OpenACC)
10264# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10266# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10267#elif defined(MFC_OpenMP)
10268# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10270# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10272# 1601 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10274# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10275 if (cyl_coord)
then
10277# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10279# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10280#if defined(MFC_OpenACC)
10281# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10283# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10284#elif defined(MFC_OpenMP)
10285# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10287# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10289# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10291# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10293# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10295# 1603 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10297# 1605 "/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_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10300 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10301 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10302 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
10303 & k_idx, l_idx) + (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10304 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
10305 & k_idx, l_idx) - (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10306 end do;
end do;
end do
10308# 1614 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10309#if defined(MFC_OpenACC)
10310# 1614 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10312# 1614 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10313#elif defined(MFC_OpenMP)
10314# 1614 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10316# 1614 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10318# 1614 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10320# 1615 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10325 if (grid_geometry == 3)
then
10326 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10328# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10330# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10331#if defined(MFC_OpenACC)
10332# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10334# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10335#elif defined(MFC_OpenMP)
10336# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10338# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10340# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10342# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10344# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10346# 1621 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10348# 1623 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10349 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10353 local_inv_ds = 1._wp/dz(k_idx)
10354 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10355 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10356 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10357 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10358 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10364# 1637 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10365#if defined(MFC_OpenACC)
10366# 1637 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10368# 1637 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10369#elif defined(MFC_OpenMP)
10370# 1637 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10372# 1637 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10374# 1637 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10376# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10377 if (alt_soundspeed)
then
10379# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10381# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10382#if defined(MFC_OpenACC)
10383# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10385# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10386#elif defined(MFC_OpenMP)
10387# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10389# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10391# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10393# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10395# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10397# 1639 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10399# 1641 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10400 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10401 local_inv_ds = 1._wp/dz(k_idx)
10402 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10404 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10405 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10406 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10407 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10408 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10409 end do;
end do;
end do
10411# 1651 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10412#if defined(MFC_OpenACC)
10413# 1651 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10415# 1651 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10416#elif defined(MFC_OpenMP)
10417# 1651 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10419# 1651 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10421# 1651 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10423# 1652 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10425 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10428# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10430# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10431#if defined(MFC_OpenACC)
10432# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10434# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10435#elif defined(MFC_OpenMP)
10436# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10438# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10440# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10442# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10444# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10446# 1655 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10448# 1657 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10449 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10450 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10451 local_inv_ds = 1._wp/dz(k_idx)
10452 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10453 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10454 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10455 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10456 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10457 end do;
end do;
end do
10460# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10461#if defined(MFC_OpenACC)
10462# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10464# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10465#elif defined(MFC_OpenMP)
10466# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10468# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10470# 1667 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10472# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10473 if (alt_soundspeed)
then
10475# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10477# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10478#if defined(MFC_OpenACC)
10479# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10481# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10482#elif defined(MFC_OpenMP)
10483# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10485# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10487# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10489# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10491# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10493# 1669 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10495# 1671 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10496 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10497 local_inv_ds = 1._wp/dz(k_idx)
10498 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10499 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10500 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10501 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10502 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10503 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10504 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10505 end do;
end do;
end do
10507# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10508#if defined(MFC_OpenACC)
10509# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10511# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10512#elif defined(MFC_OpenMP)
10513# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10515# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10517# 1681 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10519# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10523 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10526# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10528# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10529#if defined(MFC_OpenACC)
10530# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10532# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10533#elif defined(MFC_OpenMP)
10534# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10536# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10538# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10540# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10542# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10544# 1687 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10546# 1689 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10547 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10551 local_inv_ds = 1._wp/dz(k_idx)
10552 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10553 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10554 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10555 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10556 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10562# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10563#if defined(MFC_OpenACC)
10564# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10566# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10567#elif defined(MFC_OpenMP)
10568# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10570# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10572# 1703 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10574# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10575 if (alt_soundspeed)
then
10577# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10579# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10580#if defined(MFC_OpenACC)
10581# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10583# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10584#elif defined(MFC_OpenMP)
10585# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10587# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10589# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10591# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10593# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10595# 1705 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10597# 1707 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10598 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10599 local_inv_ds = 1._wp/dz(k_idx)
10600 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10602 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10603 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10604 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10605 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10606 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10607 end do;
end do;
end do
10609# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10610#if defined(MFC_OpenACC)
10611# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10613# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10614#elif defined(MFC_OpenMP)
10615# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10617# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10619# 1717 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10621# 1718 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10623 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10626# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10628# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10629#if defined(MFC_OpenACC)
10630# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10632# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10633#elif defined(MFC_OpenMP)
10634# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10636# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10638# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10640# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10642# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10644# 1721 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10646# 1723 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10647 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10648 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10649 local_inv_ds = 1._wp/dz(k_idx)
10650 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10651 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10652 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10653 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10654 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10655 end do;
end do;
end do
10658# 1733 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10659#if defined(MFC_OpenACC)
10660# 1733 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10662# 1733 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10663#elif defined(MFC_OpenMP)
10664# 1733 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10666# 1733 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10668# 1733 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10670# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10671 if (alt_soundspeed)
then
10673# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10675# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10676#if defined(MFC_OpenACC)
10677# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10679# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10680#elif defined(MFC_OpenMP)
10681# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10683# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10685# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10687# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10689# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10691# 1735 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10693# 1737 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10694 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10695 local_inv_ds = 1._wp/dz(k_idx)
10696 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10697 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10698 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10699 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10700 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10701 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10702 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10703 end do;
end do;
end do
10705# 1747 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10706#if defined(MFC_OpenACC)
10707# 1747 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10709# 1747 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10710#elif defined(MFC_OpenMP)
10711# 1747 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10713# 1747 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10715# 1747 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10717# 1748 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"