@@ -818,30 +818,6 @@ contains
818818 ! ===============================================================
819819 call nvtxStartRange(" RHS-Riemann" )
820820
821- ! Computing Riemann Solver Flux and Source Flux =================
822- ! do l = 0 , p
823- ! do k = 0 , n
824- ! do j = 0 , m
825- ! do q = 1 , sys_size
826- ! if (isnan(q_prim_qp%vf(q)%sf(j,k,l))) then
827- ! write (* ,* ) t_step,id,q,j,k,l,' q_prim_qp'
828- ! end if
829- ! end do
830- ! do q = momxb,momxe
831- ! if (isnan(dqR_prim_dx_n(id)%vf(q)%sf(j,k,l))) then
832- ! write (* ,* ) t_step,id,q,j,k,l,' dqR_prim_dx_n'
833- ! end if
834- ! if (isnan(dqR_prim_dy_n(id)%vf(q)%sf(j,k,l))) then
835- ! write (* ,* ) t_step,id,q,j,k,l,' dqR_prim_dy_n'
836- ! end if
837- ! if (isnan(dqR_prim_dz_n(id)%vf(q)%sf(j,k,l))) then
838- ! write (* ,* ) t_step,id,q,j,k,l,' dqR_prim_dz_n'
839- ! end if
840- ! end do
841- ! end do
842- ! end do
843- ! end do
844-
845821 call s_riemann_solver(qR_rsx_vf, qR_rsy_vf, qR_rsz_vf, &
846822 dqR_prim_dx_n(id)%vf, &
847823 dqR_prim_dy_n(id)%vf, &
@@ -858,22 +834,7 @@ contains
858834 flux_gsrc_n(id)%vf, &
859835 id, ix, iy, iz)
860836 call nvtxEndRange
861- ! do l = 0 , p
862- ! do k = 0 , n
863- ! do j = 0 , m
864- ! do q = 1 , sys_size
865- ! if (isnan(flux_n(id)%vf(q)%sf(j,k,l))) then
866- ! write (* ,* ) t_step,id,q,j,k,l,' flux_n'
867- ! end if
868- ! end do
869- ! do q = momxb,advxe
870- ! if (isnan(flux_src_n(id)%vf(q)%sf(j,k,l))) then
871- ! write (* ,* ) t_step,id,q,j,k,l,' flux_src_n'
872- ! end if
873- ! end do
874- ! end do
875- ! end do
876- ! end do
837+
877838 ! ===============================================================
878839
879840
0 commit comments