get_residue Subroutine

  1. ldfss0.f90
  2. ldfss0
  3. get_residue

public subroutine get_residue()

Compute the residue for the ldfss0 scheme

Arguments

None

Calls

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

Contents

Source Code

get_residue

Source Code

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

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

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