GOTO, weiter

mkmov, mov, 'AR*'

  ;;; 57 is bad - replace
mov(*, *, 57) = total(mov(*, *, [56, 58]), 3)/2

enhance, mov

rubb, mov, av=5, cell=24, step=16, range=3, plot=3, /same

mov(*, *, 57) = total(mov(*, *, [56, 58]), 3)/2

save,mov

f = findfile('AR*')
tim = 3600l*strmid(f, 15, 2)+60l*strmid(f, 17, 2)+strmid(f, 19, 2)
tim0 = tim(0)
tim = tim-tim0

sh = fltarr(2, 65)
FOR i=0, 64 DO sh(*, i) = shc(ref, mov(*, *, i), /n2, /int)
FOR i=0, 64 DO mov(*, *, i) = f_shift(mov(*, *, i), sh(0, i), sh(1, i))


mov1 = intarr(333, 235, 64)
FOR i=0, 63 DO mov1(*, *, i) = binsum(mov(*, *, i), 2)
FOR i=0, 63 DO mov1(*, *, i) = 1000./avg(mov1(120:254, *, i))*mov1(*, *, i)


pmov = fltarr(333, 235, 33)
FOR i=0, 332 DO $
  FOR j=0, 234 DO $
  pmov(i, j, *) = power(reform(mov1(i, j, *)), win='cosine')
;; pmov(i, j, *) = power(reform(mov1(i, j, *)), win='Hanning')


 ;;;
 ;;; How reliable are the power values?

prel = intarr(333, 235, 33)
FOR i=0, 332 DO BEGIN
    FOR j=0, 234 DO BEGIN
        FOR k=0, 99 DO BEGIN
            tmp = power(reform(mov1(i, j, sort(randomu(seed, 100)))), $
                        win='cosine')
            prel(i, j, *) = prel(i, j, *) + (pmov(i, j, *) GT tmp)
        ENDFOR
    ENDFOR
ENDFOR


stop

f = findgen(33)/32*(0.5/25)
pic=total(mov1,3)
pic=pic/avg(pic)
tvimg,histo_opt(total(pmov(*,*,8:10),3)),/asp,subtit=""
contour,pic,/over,lev=[.4,.6,1.2],c_line=[0,0,4]
p=tvrd()
write_jpeg,'5mhz.jpg',p
tvimg,histo_opt(total(pmov(*,*,5:7),3)),/asp,subtit=""
contour,pic,/over,lev=[.4,.6,1.2],c_line=[0,0,4]
p=tvrd()
write_jpeg,'3mhz.jpg',p
tvimg,pic,/asp,subtit=''
p=tvrd()
write_jpeg,'wl.jpg',p

show, pic, /cent
winplot
erase
contour,pic,/over,lev=[.4,.6],c_line=[0,0]
m1 = tvrd()
erase
contour,pic,/over,lev=1.2,c_line=4
m2 = tvrd()






s = size(pic)
window, 1, xs=3*s(1), ys=4*s(2)

 ;;; following does not work: power subtracts mean.  Use pic = average
 ; norm = sqrt(pmov(*, *, 0))
norm = pic
Weiter:


off = 0
FOR i=off, off+11 DO BEGIN
    x = s(1)*((i-off) MOD 3)
    y = s(2)*((i-off)/3)
    tvscl, histo_opt(pmov(*, *, i)/norm ), x, y
    xyouts, x+10, y+10, /dev, $
      string(1000*f(i), 1./f(i)/60, form="(f5.2,'mHz, ',f4.1,'min')")
ENDFOR

showmask, m1, m2


END