define real flux
define character testo*28
!
poly 
mean
if (poly$sum.eq.0) then
 say "Geminga... no valid map found."
 return
endif

say " "
for beam_size 9.75 to 11.5 by 0.25 
 let flux poly$sum|(beam_size*sec)^2|pi*4*log(2)
 if (aaa%unit.eq." mJy") then
   if (flux.lt.1000) then
     say "Total Flux: " flux " mJy  (beam: " beam_size " x " beam_size ")" /format a12 f7.2 a13 f4.1 a3 f4.1 a1
    else
     say "Total Flux: " flux|1000 " Jy   (beam: " beam_size " x " beam_size ")" /format a12 f7.2 a13 f4.1 a3 f4.1 a1
    endif
 endif 
 if (aaa%unit.eq." Jy") then
       say "Total Flux: " flux " Jy   (beam: " beam_size " x " beam_size ")" /format a12 f7.2 a13 f4.1 a3 f4.1 a1
 endif
 if (aaa%unit.eq." Counts") then
   if (flux.gt.1000000) then
     say "Total Flux: " flux|1e6 " MegaCounts (beam: " beam_size " x " beam_size ")" /format a12 f7.2 a19 f4.1 a3 f4.1 a1
   else
     say "Total Flux: " flux " Counts (beam: " beam_size " x " beam_size ")" /format a12 i7 a15 f4.1 a3 f4.1 a1
   endif
 endif
next

if (beam_az.ne.0.and.beam_el.ne.0) then
  say " "
  let testo " as from last pointing scan)" 
  let flux poly$sum|(beam_az*beam_el*sec^2)|pi*4*log(2)
  if (aaa%unit.eq." mJy") then
    if (flux.lt.1000) then
      say "Total Flux: " flux " mJy  (beam: " beam_az " x " beam_el 'testo' /format a12 f7.2 a13 f4.1 a3 f4.1 a28
     else
      say "Total Flux: " flux|1000 " Jy   (beam: " beam_az " x " beam_el 'testo' /format a12 f7.2 a13 f4.1 a3 f4.1 a28
     endif
  endif 
  if (aaa%unit.eq." Jy") then
        say "Total Flux: " flux " Jy   (beam: " beam_az " x " beam_el 'testo' /format a12 f7.2 a13 f4.1 a3 f4.1 a28
  endif
  if (aaa%unit.eq." Counts") then
    if (flux.gt.1000000) then
      say "Total Flux: " flux|1e6 " MegaCounts (beam: " beam_az " x " beam_el 'testo' /format a12 f7.2 a19 f4.1 a3 f4.1 a28
    else
      say "Total Flux: " flux " Counts (beam: " beam_az " x " beam_el 'testo' /format a12 i7 a15 f4.1 a3 f4.1 a28
    endif
  endif
endif