9607 & flux_src_n_vf_arg, Kterm_arg)
9608 integer,
intent(in) :: current_idir
9609 type(scalar_field),
dimension(sys_size),
intent(inout) :: rhs_vf_arg
9610 type(vector_field),
intent(in) :: q_cons_vf_arg
9611 type(vector_field),
intent(in) :: q_prim_vf_arg
9612 type(vector_field),
intent(in) :: flux_src_n_vf_arg
9613 real(wp),
allocatable,
dimension(:,:,:),
intent(in) :: Kterm_arg
9614 integer :: j_adv, k_idx, l_idx, q_idx
9615 real(wp) :: local_inv_ds, local_term_coeff, local_flux1, local_flux2
9616 real(wp) :: local_k_term_val
9626 select case (current_idir)
9628 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9631# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9633# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9634#if defined(MFC_OpenACC)
9635# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9637# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9638#elif defined(MFC_OpenMP)
9639# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9641# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9643# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9645# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9647# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9649# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9651# 1463 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9652 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9656 local_inv_ds = 1._wp/dx(k_idx)
9657 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(k_idx, l_idx, q_idx)
9658 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx - 1, l_idx, q_idx)
9659 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9660 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9661 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9667# 1477 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9668#if defined(MFC_OpenACC)
9669# 1477 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9671# 1477 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9672#elif defined(MFC_OpenMP)
9673# 1477 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9675# 1477 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9677# 1477 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9679 if (alt_soundspeed)
then
9681# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9683# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9684#if defined(MFC_OpenACC)
9685# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9687# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9688#elif defined(MFC_OpenMP)
9689# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9691# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9693# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9695# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9697# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9699# 1479 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9701# 1481 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9702 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9703 local_inv_ds = 1._wp/dx(k_idx)
9704 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9706 local_flux2 =
nc_iface_vel_n(1)%vf(1)%sf(k_idx - 1, l_idx, q_idx)
9707 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, &
9708 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9709 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, &
9710 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9711 end do;
end do;
end do
9713# 1491 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9714#if defined(MFC_OpenACC)
9715# 1491 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9717# 1491 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9718#elif defined(MFC_OpenMP)
9719# 1491 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9721# 1491 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9723# 1491 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9726 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9729# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9731# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9732#if defined(MFC_OpenACC)
9733# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9735# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9736#elif defined(MFC_OpenMP)
9737# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9739# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9741# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9743# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9745# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9747# 1495 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9749# 1497 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9750 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9751 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9752 local_inv_ds = 1._wp/dx(k_idx)
9753 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(k_idx, l_idx, q_idx)
9754 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
9755 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
9756 rhs_vf_arg(j_adv)%sf(k_idx, l_idx, q_idx) = rhs_vf_arg(j_adv)%sf(k_idx, l_idx, &
9757 & q_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9758 end do;
end do;
end do
9761# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9762#if defined(MFC_OpenACC)
9763# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9765# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9766#elif defined(MFC_OpenMP)
9767# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9769# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9771# 1507 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9773 if (alt_soundspeed)
then
9775# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9777# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9778#if defined(MFC_OpenACC)
9779# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9781# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9782#elif defined(MFC_OpenMP)
9783# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9785# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9787# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9789# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9791# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9793# 1509 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9795# 1511 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9796 do q_idx = 0, p;
do l_idx = 0, n;
do k_idx = 0, m
9797 local_inv_ds = 1._wp/dx(k_idx)
9798 local_k_term_val = kterm_arg(k_idx, l_idx, q_idx)
9799 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx, l_idx, q_idx)
9800 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(k_idx - 1, l_idx, q_idx)
9801 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, &
9802 & q_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9803 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, &
9804 & q_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9805 end do;
end do;
end do
9807# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9808#if defined(MFC_OpenACC)
9809# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9811# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9812#elif defined(MFC_OpenMP)
9813# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9815# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9817# 1521 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9822 if (adv_src_mode == adv_src_mode_alpha_iface)
then
9825# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9827# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9828#if defined(MFC_OpenACC)
9829# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9831# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9832#elif defined(MFC_OpenMP)
9833# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9835# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9837# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9839# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9841# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9843# 1527 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9845# 1529 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9846 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9850 local_inv_ds = 1._wp/dy(k_idx)
9851 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(q_idx, k_idx, l_idx)
9852 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx - 1, l_idx)
9853 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
9854 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
9855 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9861# 1543 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9862#if defined(MFC_OpenACC)
9863# 1543 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9865# 1543 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9866#elif defined(MFC_OpenMP)
9867# 1543 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9869# 1543 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9871# 1543 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9873 if (alt_soundspeed)
then
9875# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9877# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9878#if defined(MFC_OpenACC)
9879# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9881# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9882#elif defined(MFC_OpenMP)
9883# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9885# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9887# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9889# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9891# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9893# 1545 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9895# 1547 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9896 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9897 local_inv_ds = 1._wp/dy(k_idx)
9898 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
9900 local_flux2 =
nc_iface_vel_n(2)%vf(2)%sf(q_idx, k_idx - 1, l_idx)
9901 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, &
9902 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9903 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, &
9904 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
9905 end do;
end do;
end do
9907# 1557 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9908#if defined(MFC_OpenACC)
9909# 1557 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9911# 1557 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9912#elif defined(MFC_OpenMP)
9913# 1557 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9915# 1557 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9917# 1557 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9921# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9923# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9924#if defined(MFC_OpenACC)
9925# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9927# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9928#elif defined(MFC_OpenMP)
9929# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9931# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9933# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9935# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9937# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9939# 1559 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9941 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9942 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
9943 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
9944 & k_idx, l_idx) + local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
9945 & k_idx, l_idx)/y_cc(k_idx)
9946 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
9947 & k_idx, l_idx) - local_k_term_val*q_prim_vf_arg%vf(eqn_idx%mom%beg + 1)%sf(q_idx, &
9948 & k_idx, l_idx)/y_cc(k_idx)
9949 end do;
end do;
end do
9951# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9952#if defined(MFC_OpenACC)
9953# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9955# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9956#elif defined(MFC_OpenMP)
9957# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9959# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9961# 1569 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9965 else if (adv_src_mode == adv_src_mode_vel_iface)
then
9968# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9970# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9971#if defined(MFC_OpenACC)
9972# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9974# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9975#elif defined(MFC_OpenMP)
9976# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9978# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9980# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9982# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9984# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9986# 1574 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9988# 1576 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
9989 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
9990 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
9991 local_inv_ds = 1._wp/dy(k_idx)
9992 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(q_idx, k_idx, l_idx)
9993 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
9994 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
9995 rhs_vf_arg(j_adv)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(j_adv)%sf(q_idx, k_idx, &
9996 & l_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
9997 end do;
end do;
end do
10000# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10001#if defined(MFC_OpenACC)
10002# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10004# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10005#elif defined(MFC_OpenMP)
10006# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10008# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10010# 1586 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10012 if (alt_soundspeed)
then
10014# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10016# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10017#if defined(MFC_OpenACC)
10018# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10020# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10021#elif defined(MFC_OpenMP)
10022# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10024# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10026# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10028# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10030# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10032# 1588 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10034# 1590 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10035 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10036 local_inv_ds = 1._wp/dy(k_idx)
10037 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10038 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10039 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10040 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, &
10041 & l_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10042 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, &
10043 & l_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10044 end do;
end do;
end do
10046# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10047#if defined(MFC_OpenACC)
10048# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10050# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10051#elif defined(MFC_OpenMP)
10052# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10054# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10056# 1600 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10058 if (cyl_coord)
then
10060# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10062# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10063#if defined(MFC_OpenACC)
10064# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10066# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10067#elif defined(MFC_OpenMP)
10068# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10070# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10072# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10074# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10076# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10078# 1602 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10080# 1604 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10081 do l_idx = 0, p;
do k_idx = 0, n;
do q_idx = 0, m
10082 local_k_term_val = kterm_arg(q_idx, k_idx, l_idx)
10083 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx)
10084 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(q_idx, k_idx - 1, l_idx)
10085 rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(q_idx, &
10086 & k_idx, l_idx) + (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10087 rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, k_idx, l_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(q_idx, &
10088 & k_idx, l_idx) - (local_k_term_val/(2._wp*y_cc(k_idx)))*(local_flux1 + local_flux2)
10089 end do;
end do;
end do
10091# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10092#if defined(MFC_OpenACC)
10093# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10095# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10096#elif defined(MFC_OpenMP)
10097# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10099# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10101# 1613 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10107 if (grid_geometry == 3)
then
10108 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10110# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10112# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10113#if defined(MFC_OpenACC)
10114# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10116# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10117#elif defined(MFC_OpenMP)
10118# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10120# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10122# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10124# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10126# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10128# 1620 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10130# 1622 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10131 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10135 local_inv_ds = 1._wp/dz(k_idx)
10136 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10137 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10138 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10139 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10140 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10146# 1636 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10147#if defined(MFC_OpenACC)
10148# 1636 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10150# 1636 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10151#elif defined(MFC_OpenMP)
10152# 1636 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10154# 1636 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10156# 1636 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10158 if (alt_soundspeed)
then
10160# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10162# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10163#if defined(MFC_OpenACC)
10164# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10166# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10167#elif defined(MFC_OpenMP)
10168# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10170# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10172# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10174# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10176# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10178# 1638 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10180# 1640 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10181 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10182 local_inv_ds = 1._wp/dz(k_idx)
10183 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10185 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10186 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10187 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10188 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10189 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10190 end do;
end do;
end do
10192# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10193#if defined(MFC_OpenACC)
10194# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10196# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10197#elif defined(MFC_OpenMP)
10198# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10200# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10202# 1650 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10205 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10208# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10210# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10211#if defined(MFC_OpenACC)
10212# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10214# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10215#elif defined(MFC_OpenMP)
10216# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10218# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10220# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10222# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10224# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10226# 1654 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10228# 1656 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10229 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10230 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10231 local_inv_ds = 1._wp/dz(k_idx)
10232 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10233 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10234 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10235 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10236 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10237 end do;
end do;
end do
10240# 1666 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10241#if defined(MFC_OpenACC)
10242# 1666 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10244# 1666 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10245#elif defined(MFC_OpenMP)
10246# 1666 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10248# 1666 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10250# 1666 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10252 if (alt_soundspeed)
then
10254# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10256# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10257#if defined(MFC_OpenACC)
10258# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10260# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10261#elif defined(MFC_OpenMP)
10262# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10264# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10266# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10268# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10270# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10272# 1668 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10274# 1670 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10275 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10276 local_inv_ds = 1._wp/dz(k_idx)
10277 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10278 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10279 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10280 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10281 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10282 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10283 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10284 end do;
end do;
end do
10286# 1680 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10287#if defined(MFC_OpenACC)
10288# 1680 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10290# 1680 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10291#elif defined(MFC_OpenMP)
10292# 1680 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10294# 1680 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10296# 1680 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10301 if (adv_src_mode == adv_src_mode_alpha_iface)
then
10304# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10306# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10307#if defined(MFC_OpenACC)
10308# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10310# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10311#elif defined(MFC_OpenMP)
10312# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10314# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10316# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10318# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10320# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10322# 1686 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10324# 1688 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10325 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10329 local_inv_ds = 1._wp/dz(k_idx)
10330 local_term_coeff = q_prim_vf_arg%vf(eqn_idx%cont%end + current_idir)%sf(l_idx, q_idx, k_idx)
10331 local_flux1 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx - 1)
10332 local_flux2 = flux_src_n_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10333 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10334 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10340# 1702 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10341#if defined(MFC_OpenACC)
10342# 1702 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10344# 1702 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10345#elif defined(MFC_OpenMP)
10346# 1702 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10348# 1702 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10350# 1702 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10352 if (alt_soundspeed)
then
10354# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10356# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10357#if defined(MFC_OpenACC)
10358# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10360# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10361#elif defined(MFC_OpenMP)
10362# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10364# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10366# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10368# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10370# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10372# 1704 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10374# 1706 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10375 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10376 local_inv_ds = 1._wp/dz(k_idx)
10377 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10379 local_flux2 =
nc_iface_vel_n(3)%vf(3)%sf(l_idx, q_idx, k_idx - 1)
10380 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10381 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10382 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10383 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10384 end do;
end do;
end do
10386# 1716 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10387#if defined(MFC_OpenACC)
10388# 1716 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10390# 1716 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10391#elif defined(MFC_OpenMP)
10392# 1716 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10394# 1716 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10396# 1716 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10399 else if (adv_src_mode == adv_src_mode_vel_iface)
then
10402# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10404# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10405#if defined(MFC_OpenACC)
10406# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10408# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10409#elif defined(MFC_OpenMP)
10410# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10412# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10414# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10416# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10418# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10420# 1720 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10422# 1722 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10423 do j_adv = eqn_idx%adv%beg, eqn_idx%adv%end
10424 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10425 local_inv_ds = 1._wp/dz(k_idx)
10426 local_term_coeff = q_cons_vf_arg%vf(j_adv)%sf(l_idx, q_idx, k_idx)
10427 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10428 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10429 rhs_vf_arg(j_adv)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(j_adv)%sf(l_idx, q_idx, &
10430 & k_idx) + local_inv_ds*local_term_coeff*(local_flux1 - local_flux2)
10431 end do;
end do;
end do
10434# 1732 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10435#if defined(MFC_OpenACC)
10436# 1732 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10438# 1732 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10439#elif defined(MFC_OpenMP)
10440# 1732 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10442# 1732 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10444# 1732 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10446 if (alt_soundspeed)
then
10448# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10450# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10451#if defined(MFC_OpenACC)
10452# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10454# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10455#elif defined(MFC_OpenMP)
10456# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10458# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10460# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10462# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10464# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10466# 1734 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10468# 1736 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10469 do k_idx = 0, p;
do q_idx = 0, n;
do l_idx = 0, m
10470 local_inv_ds = 1._wp/dz(k_idx)
10471 local_k_term_val = kterm_arg(l_idx, q_idx, k_idx)
10472 local_flux1 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx)
10473 local_flux2 = flux_src_n_vf_arg%vf(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx - 1)
10474 rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%beg)%sf(l_idx, &
10475 & q_idx, k_idx) + local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10476 rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, q_idx, k_idx) = rhs_vf_arg(eqn_idx%adv%end)%sf(l_idx, &
10477 & q_idx, k_idx) - local_k_term_val*local_inv_ds*(local_flux1 - local_flux2)
10478 end do;
end do;
end do
10480# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10481#if defined(MFC_OpenACC)
10482# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10484# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10485#elif defined(MFC_OpenMP)
10486# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10488# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"
10490# 1746 "/home/runner/work/MFC/MFC/src/simulation/m_rhs.fpp"