if .not.exist(FLUX%VALUES) then
  message I "SHOW FLUX" "No Flux values, use FLUX command before"
  return error
endif
!
if (flux%nc.le.1) then
  message I "SHOW FLUX" "Only 1 channel"
  for i 1 to flux_nf
    SAY "Field "'i'",  Flux "'flux%values[1,i]'" Jy"
  next
  return
endif
!
gtvl\destroy all
define integer nr rf rl
rf = 1
rl = 1
if (pro%narg.ge.1) then
  let rf min('pro%arg[1]',flux%nf)
  if pro%narg.ge.2 then
    let rl min('pro%arg[2]',flux%nf)
    if (rl.lt.rf) then
      let nr rl
      let rl rf
      let rf nr
    endif
  else
    let rl rf
  endif
  !
endif
let nr rl-rf+1
!
@ window_tools
@ window_init
let aspect 2
@ window_xy 1 'nr' 1 1
greg1\limits /var flux%velocities flux%values[1]
for i rf to rl
  lim = = < > /var flux%velocities flux%values[i]
next
greg1\axis xl
!
for i rf to rl
  @ window_xy 1 'nr' 1 'i-rf+1'
  greg1\histogram flux%velocities flux%values[i]
  greg1\axis yl
  greg1\axis yr
  greg1\axis xu
  greg1\axis xl /label n
  if (nr.gt.1)  draw text 1 -1 'i' /char 7
next
@ window_xy 1 'nr' 1 1
greg1\label "Velocity (km/s)" /X
greg1\label "Flux Density (Jy) " /Y
!
