on error return 
if (&1major.eq.0) then 
  say "No clean beam, flux density estimate impossible" 
  return 
endif 
polygon &2 
define real flux 
mean 
let flux poly$sum|pi|&1major|&1minor*4*log(2) 
if (flux.lt.0.01) then 
  say "Total Flux "'0.1*nint(flux*10000.)'" mJy" 
else if (flux.lt.1.0) then 
  say "Total Flux "'nint(flux*1000.)'" mJy" 
else 
  say "Total Flux "'0.01*nint(flux*100.)'" Jy" 
endif