get_residue Subroutine

  1. ausm.f90
  2. ausm
  3. get_residue

public subroutine get_residue()

Compute the residue for the AUSM scheme

Arguments

None

Calls

proc~~get_residue~5~~CallsGraph proc~get_residue~5 get_residue proc~dmsg dmsg proc~get_residue~5->proc~dmsg

Contents

Source Code

get_residue

Source Code

        subroutine get_residue()
            !< Compute the residue for the AUSM scheme
            !-----------------------------------------------------------

            implicit none
            
            integer :: i, j, k, l

            call dmsg(1, 'ausm', 'compute_residue')

            do l = 1, n_var
             do k = 1, kmx - 1
              do j = 1, jmx - 1
               do i = 1, imx - 1
               residue(i, j, k, l) = (F(i+1, j, k, l) - F(i, j, k, l)) &
                                   + (G(i, j+1, k, l) - G(i, j, k, l)) &
                                   + (H(i, j, k+1, l) - H(i, j, k, l))
               !print*, i, j, l, F(i+1, j,k,l), F(i,j,k,l), G(i,j+1,k,l)&
               !       , G(i,j,k,l), H(i,j,k+1,l), H(i,j,k,l)
               end do
              end do
             end do
            end do
        
        end subroutine get_residue