histogram_get_range Subroutine

subroutine histogram_get_range(h,index,lower,upper)
  type(histogram),intent(in) :: h
  integer,intent(in)         :: index
  real(8),intent(out)        :: lower,upper
  if(index>=h%n)then
     print*,"error: *i lies outside valid range=0...n-1"
     stop
  endif
  lower=h%range(index)
  upper=h%range(index+1)
end subroutine histogram_get_range