pglineb Subroutine

private subroutine pglineb(n, xpts, ypts, blnk)

FIXME : Add documentation

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: n
real(kind=pgkind), intent(in) :: xpts(*)
real(kind=pgkind), intent(in) :: ypts(*)
real(kind=pgkind), intent(in) :: blnk

Contents

Source Code


Source Code

  subroutine pglineb(n,xpts,ypts,blnk)
    implicit none

    integer :: i
    integer, intent(in) :: n
    real (kind=pgkind), intent(in) :: xpts(*), ypts(*)
    integer :: n2
    real (kind=pgkind) :: xpts2(n), ypts2(n)
    real (kind=pgkind), intent(in) :: blnk
    real, parameter :: epsilon=1.e-20


    n2=0
    do i = 1, n
       if(abs(ypts(i)-blnk).lt.epsilon) cycle
       n2=n2+1
       xpts2(n2)=xpts(i)
       ypts2(n2)=ypts(i)
    end do

    call pgline(n2,xpts2,ypts2)

    return

  end subroutine pglineb