dedge Function

private function dedge(a, r, n, iside)

FIXME : Add documentation

Arguments

Type IntentOptional Attributes Name
real, intent(in), dimension(n) :: a
real, intent(in), dimension(n) :: r
integer, intent(in) :: n
integer, intent(in) :: iside

Return Value real


Contents

Source Code

dedge

Source Code

  real function dedge(a,r,n,iside)
    implicit none
    integer, intent(in) :: n, iside
    real, dimension(n), intent(in) :: a, r
!
! Not quite right for non-uniform r mesh
!
    if(iside.eq.1) then
       dedge=-(3.*a(1)-4.*a(2)+a(3))/(r(3)-r(1))
    else
       dedge=(3.*a(iside)-4.*a(iside-1)+a(iside-2))/(r(iside)-r(iside-2))
    endif

    return
  end function dedge