get_residue Subroutine

  1. van_leer.f90
  2. van_leer
  3. get_residue

public subroutine get_residue()

Compute the residue using the Van-Leer scheme

Arguments

None

Calls

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

Contents

Source Code

get_residue

Source Code

        subroutine get_residue()
            !< Compute the residue using the Van-Leer scheme
            !-----------------------------------------------------------

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

            call dmsg(1, 'van_leer', '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)
               end do
              end do
             end do
            end do
        
        end subroutine get_residue