define real old_rms noise_gain /like chan_rms
define real effective_gain effective_chan
define inte n1 n2
define char planet*40[7] dummy*4

let planet[1] "NEPT"
let planet[2] "URAN"
let planet[3] "SATU"
let planet[4] "JUPI"
let planet[5] "MARS"
let planet[6] "VENU"
let planet[7] "MERC"

let dummy 'source'

for i 1 to 7
  if (dummy.eq.planet[i]) then
    say "Skynoise removal is (in general) not feasible on a planet"
    say "Do it by hand ..."
    return
  endif
next

if (.not.exist(chan_flag)) then
  var flag
endif
let old_rms chan_rms
compute effective_gain mean chan_rms
if (effective_gain.eq.0) then
  say "I-BASE,["'scan'"] Remove baseline first"
  return
endif
base /noise
! base /corr
let noise_gain old_rms|chan_rms /where chan_rms.ne.0.and.old_rms.ne.0.and.chan_flag.eq..false.
compute effective_gain mean noise_gain
let noise_gain 1 /where noise_gain.ne.0
compute effective_chan mean noise_gain
SAY " "
SAY "Command was: BASELINE /NOISE"
let n1 nint(nchan*effective_chan)
let n2 nint(1000*(effective_gain|effective_chan-1))|10
SAY "Mean gain over "'n1'" channels is: "'n2'"%"
delete /var effective_gain effective_chan old_rms