get_residue Subroutine

  1. slau.f90
  2. slau
  3. get_residue

public subroutine get_residue()

Compute the residue for the slau scheme

Arguments

None

Calls

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

Contents

Source Code

get_residue

Source Code

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

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

            call dmsg(1, 'slau', '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