FIXME : Add documentation
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | n | |||
real(kind=pgkind), | intent(in) | :: | xpts(*) | |||
real(kind=pgkind), | intent(in) | :: | ypts(*) | |||
real(kind=pgkind), | intent(in) | :: | blnk |
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