; file: tinker.pro for ~rutten/data/DOT/2007-04-12/QuietSun/Cubes
; last: Oct 15 2010 

ximovie,'QuietSun-ha-0500-20070412-cube-1034x776x180-dt30',1034,776,/int,mag=0.8,afile2='QuietSun-gb-0000-20070412-cube-1034x776x180-dt30'

ximovie,'QuietSun-rc-0000-20070412-cube-1034x776x180-dt30',1034,776,/int,mag=0.8,afile2='QuietSun-gb-0000-20070412-cube-1034x776x180-dt30'
;# perfect alignment

ximovie,'QuietSun-ha-0500-20070412-cube-1034x776x180-dt30',1034,776,/int,mag=0.8,afile2='QuietSun-ha+0500-20070412-cube-1034x776x180-dt30'
;# enormous difference between the two  @ make sum-dopp cubes
;# alignment seems good (should be since KvdL)

ximovie,'QuietSun-ca-flip-20070412-cube-1034x776x90-dt60-even',1034,776,/int,mag=0.8,afile2='QuietSun-ca-flip-20070412-cube-1034x776x90-dt60-odd'

; load for processing 
caw=readdotcube('QuietSun-ca-flip-20070412-cube-1034x776x90-dt60-odd')
cac=readdotcube('QuietSun-ca-flip-20070412-cube-1034x776x90-dt60-even')
gb=readdotcube('QuietSun-gb-0000-20070412-cube-1034x776x180-dt30')
hawing=readdotcube('QuietSun-ha-wing-20070412-cube-1034x776x180-dt30')
hadopp=readdotcube('QuietSun-ha-dopp-20070412-cube-1034x776x180-dt30')

; make small ones for comparison
nt=180
gb2=caw
for it=1,nt-1,2 do gb2(*,*,it/2)=gb(*,*,it)
;; odd = synchronous starting at it=1 (Pit custum) 

nt=180
had2=cac
for it=0,nt-1,2 do had2(*,*,it/2)=hadopp(*,*,it)
;; even = synchronous starting at it=0 

; inspect
cubeximovie,gb2,cube2=caw,mag=0.7
  ;# perfect alignment
cubeximovie,caw,cube2=had2,mag=0.7

ximovie,'QuietSun-ha-wing-20070412-cube-1034x776x180-dt30',1034,776,/int,mag=0.8,afile2='QuietSun-ha-dopp-20070412-cube-1034x776x180-dt30'

xslice,hadopp,root=root
  ;# oscillations everywhere in QS, extended

; make confusograms
plotconfusogram,hawing,hadopp,30,'hawing-hadopp-confus.ps'
plotconfusogram,gb2,caw,60,'gb-caw-confus.ps'
plotconfusogram,gb2,cac,60,'gb-cac-confus.ps'
  ;# wrong timing
plotconfusogram,gb,hadopp,30,'gb-hadopp-confus.ps'
  ;# no coherence in 3-min peak but smooth phase diff decay
plotconfusogram,caw,had2,60,'caw-hadopp-confus.ps'
plotconfusogram,cac,had2,60,'cac-hadopp-confus.ps'  ; most interesting one
plotconfusogram,caw,cac,60,'caw-cac-confus.ps'   
  ; todo = undo dT=30s sampling shift

plotkophasediff,gb2,caw,0.071,60,'ko-phasediff-gb-caw.ps',contours=1,$
    kmax=0.3,minphasediff=-180,maxphasediff=90
plotkophasediff,cac,had2,0.071,60,'ko-phasediff-cac-hadopp.ps',contours=0,$
    kmax=0.2,minphasediff=-180,maxphasediff=30,lamb=7,/fundamental
plotkopower,hadopp,0.071,30,'ko-power-hadopp.ps',maxpower=3,fmax=0.5,kmax=0.2,/fundamental,/contours

spawn,'date'