; IDL Version 8.1, Mac OS X (darwin x86_64 m64)
; Journal File for Bob@steinr.pa.msu.edu
; Working directory: /Volumes/d1/96-20/Slices
; Date: Sat Oct  8 16:20:40 2011
 
journ
restore,/ver,'vv_hd96_1355.save'
help,vy01
t=82
tx=81
vy01h=vy01
 for t=0,69 do begin tmp=seeing(reform(vy01[*,*,t]),mtf='airy',ap=50,size=48,ora=0.344) & vy01h[*,*,t]=tmp
;     0.113530    0.0513644
retall
 for t=0,81 do begin tmp=seeing(reform(vy01[*,*,t]),mtf='airy',ap=50,size=48,ora=0.344) & vy01h[*,*,t]=tmp
;     0.113530    0.0513644
;     0.112851    0.0503663
;     0.117690    0.0531598
;     0.120083    0.0535228
;     0.121071    0.0533047
;     0.121913    0.0522714
;     0.124006    0.0528568
;     0.125836    0.0537168
;     0.127027    0.0541416
;     0.126836    0.0535822
;     0.126648    0.0533263
;     0.126435    0.0533657
;     0.125664    0.0528824
;     0.125909    0.0531536
;     0.125418    0.0532273
;     0.124502    0.0527260
;     0.123905    0.0525969
;     0.123280    0.0523749
;     0.122787    0.0523009
;     0.122493    0.0524607
;     0.122062    0.0523680
;     0.121527    0.0521574
;     0.121267    0.0522468
;     0.120774    0.0518857
;     0.121102    0.0519195
;     0.120900    0.0516177
;     0.121120    0.0520090
;     0.121113    0.0519912
;     0.121229    0.0520933
;     0.121439    0.0523035
;     0.120960    0.0517836
;     0.121370    0.0519944
;     0.121224    0.0519527
;     0.121089    0.0519958
;     0.120839    0.0518696
;     0.120824    0.0516652
;     0.120622    0.0513579
;     0.120737    0.0518804
;     0.120660    0.0515253
;     0.120641    0.0515568
;     0.120761    0.0517375
;     0.120476    0.0513113
;     0.120712    0.0515012
;     0.120755    0.0514713
;     0.120452    0.0513272
;     0.120248    0.0515722
;     0.120783    0.0517491
;     0.120939    0.0518373
;     0.121035    0.0515775
;     0.121188    0.0517136
;     0.121248    0.0517767
;     0.121009    0.0512681
;     0.121328    0.0516657
;     0.121298    0.0515604
;     0.121149    0.0515334
;     0.121320    0.0517709
;     0.121613    0.0519907
;     0.121514    0.0517366
;     0.121215    0.0515029
;     0.121441    0.0521185
;     0.120797    0.0515633
;     0.121037    0.0517175
;     0.121116    0.0518389
;     0.120629    0.0511762
;     0.120440    0.0513087
;     0.120717    0.0517408
;     0.120603    0.0512440
;     0.120487    0.0513630
;     0.120349    0.0513161
;     0.120439    0.0509932
;     0.120700    0.0515617
;     0.120430    0.0513143
;     0.120533    0.0513737
;     0.120762    0.0514511
;     0.120654    0.0512223
;     0.121112    0.0515695
;     0.120837    0.0509780
;     0.120944    0.0512301
;     0.121538    0.0519194
;     0.121244    0.0515070
;     0.120491    0.0513621
;     0.120650    0.0513797
help;,vy01
help,vy01
vy01p=fltarr(2016,82)
for t=0,81 do begin sizesp,reform(vy01h[*,*,t]),size=48,wavenumbers=k,spectrum=tmp & vy01p[*,t]=tmp
; % Array subscript for VY01P must have same size as source expression.
help,tmp
help,vy01h
vy01p=fltarr(1007,82)
for t=0,81 do begin sizesp,reform(vy01h[*,*,t]),size=48,wavenumbers=k,spectrum=tmp & vy01p[*,t]=tmp
stat,v01p
; % MIN: Variable is undefined: A.
stat,vy01p
; % MIN: Variable is undefined: A.
retall
stat,vy01p
;aver =    2.052e-05
;rms  =    6.104e-05     rms/aver =        2.974
;min  =    2.977e-19     min/aver =    1.450e-14        389   1
;max  =    0.0004110     max/aver =        20.03         35   4
 vy01pm=haver(vy01p)
help,vy01pm
plot,k,vy01pm,/xlog,/ylog
print,(k[2]-k[1])
;     0.130900
plot,k,.5*vy01pm*2.7e5/(k[2]-k[1]),/xlog,/ylog,yr=[10,1000],xr=[0.1,100],ytit='EK',xtit='k (1/Mm)',chars=2
oplot,k,k^(0.3)*1.9e2,lin=1
oplot,k,k^(.5)*1.5e2,lin=1
oplot,k,k^(.5)*1.3e2,lin=1
oplot,k,k^(.5)*1.4e2,lin=1
plot,k,.5*vy01pm*2.7e5/(k[2]-k[1]),/xlog,/ylog,yr=[10,1000],xr=[0.1,100],ytit='EK',xtit='k (1/Mm)',chars=2
oplot,k,k^(0.3)*1.9e2,lin=1
oplot,k,k^(.5)*1.4e2,lin=1
oplot,k,k^(0.3)*1.8e2,lin=1
plot,k,.5*vy01pm*2.7e5/(k[2]-k[1]),/xlog,/ylog,yr=[10,1000],xr=[0.1,100],ytit='EK',xtit='k (1/Mm)',chars=2
oplot,k,k^(0.3)*1.8e2,lin=1
oplot,k,k^(.5)*1.4e2,lin=1
xyouts,.2,90,'0.5'
xyouts,.2,40,'0.5'
xyouts,.2,45,'0.5'
xyouts,.2,90,'0.3'
xyouts,.2,120,'0.3'
xyouts,.2,130,'0.3'
ps
plot,k,.5*vy01pm*2.7e5/(k[2]-k[1]),/xlog,/ylog,yr=[10,1000],xr=[0.1,100],ytit='EK',xtit='k (1/Mm)',chars=2
oplot,k,k^(0.3)*1.8e2,lin=1
oplot,k,k^(.5)*1.4e2,lin=1
xyouts,.2,45,'0.5'
xyouts,.2,130,'0.3'
psend
x
$mv idl.ps Ek_spec_96_t1355.ps