!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
define integer level
if pro%narg.eq.0 then
  let level 1
else if "&1".eq."BRIEF" then
  let level -1
else if ("&1".eq."?") then
  let level 1
else 
  let level 2
endif
!
if (level.ge.0) then
  say "CLEAN deconvolves the DIRTY image map using the dirty BEAM"
  say " "
  say "* The outputs are the CLEAN and RESIDUAL images, and the "
  say "  Clean Component Table, CCT"
  say " "
endif
!
define character r*5 n*5 b*5
r = 'ansi%red'
n = 'ansi%none'
b = 'ansi%blue'
!
define double temp
let umet 'method' /upper /new chara*6 
say "  Cleaning method                "'b'"CLEAN_METHOD    "'n'"[" 'r''method''n' "]     METHOD "
let temp 'nint(gain*1000)/1000'
say "  Loop gain                      "'b'"CLEAN_GAIN      "'n'"[" 'r''temp''n' "]        GAIN   "
say "  Stopping level                 "'b'"CLEAN_STOP      "'n'"[" 'r''clean_stop[1]'" "'clean_stop[2]''n' "]  "
let temp 'nint(fres*10000)/10000'
say "    Fractional residual            "'b'"CLEAN_FRES    "'n'"[" 'r''temp''n' "]          FRES"
say "    Absolute residual (Jy)         "'b'"CLEAN_ARES    "'n'"[" 'r''ares''n' "]          ARES" 
say "    Number of Clean components     "'b'"CLEAN_NITER   "'n'"[" 'r''niter''n' "]          NITER "
if (umet.ne."SDI") then
  say "    Minimum number of components   "'b'"CLEAN_NKEEP   "'n'"[" 'r''clean_nkeep''n' "]   "
endif
if (umet.eq."CLARK") then
  say "  Number of major cycles         "'b'"CLEAN_NCYCLE    "'n'"[" 'r''clean_ncycle''n' "]    "
  if (level.gt.1) then
    say "  Number of components per cycle "'b'"CLEAN_NGOAL     "'n'"[" 'r''clean_ngoal''n' "]     "
    say "  Caution factor                 "'b'"CLEAN_WORRY     "'n'"[" 'r''clean_worry''n' "]     "
    say "  Speed-up factor                "'b'"CLEAN_SPEEDY    "'n'"[" 'r''clean_speedy''n' "]    "
  endif
else if (umet.eq."MULTI") then
  say "  MultiScale Smoothing factor    "'b'"CLEAN_SMOOTH    "'n'"[" 'r''clean_smooth''n'"]     "
  say "  MultiScale Inflation factor    "'b'"CLEAN_INFLATE   "'n'"[" 'r''clean_inflate''n'"]    "
else if (umet.eq."MRC") then
  say "  MRC Smoothing ratio            "'b'"CLEAN_RATIO     "'n'"[" 'r''clean_ratio''n' "]     "
endif
say "  Number of positive components  "'b'"CLEAN_POSITIVE  "'n'"[" 'r''clean_positive''n' "]    "
say " "
say "  Clean beam size  (User)        "'b'"BEAM_SIZE       [ "'r''beam_size[1]'""" by "'beam_size[2]'""" at PA "'beam_size[3]''n'"° ]"
if exist(beam_fitted) then
  define double bf /like beam_fitted
  let bf 0.001*nint(beam_fitted*1000)
  let bf[3] 0.1*nint(10*bf[3]) 
  let temp 0.01*nint(beam_jvm*100)
  say "  Clean beam size  (Fitted)      "'b'"BEAM_FITTED     [ "'r''bf[1]'""" by "'bf[2]'""" at PA "'bf[3]''n'"° ], JvM "'temp'
endif
! 
if (mosaic) then
  say "* Mosaic mode is ON "
  let temp 'int(clean_search*1000)/1000'
  say "  Threshold to search components "'b'"CLEAN_SEARCH    "'n'"[ "'r''temp''n'" ]"
  let temp 'int(clean_truncate*1000)/1000'
  say "  Threshold for restoration      "'b'"CLEAN_TRUNCATE  "'n'"[ "'r''temp''n'" ]"
else
  if (level.gt.1) say "* Mosaic mode is OFF "
endif
!
if (level.lt.2) return
say "  BLC [" 'r''blc[1]' 'blc[2]''n' "]             TRC [" 'r''trc[1]' 'trc[2]''n' "]            BEAM_PATCH [" 'r''beam_patch[1]' 'beam_patch[2]''n' "]"
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
