; IDL Version 8.1, Mac OS X (darwin x86_64 m64)
; Journal File for Bob@steinr.pa.msu.edu
; Working directory: /Volumes/d2/MHDV/mhd12v/Slices
; Date: Thu Jun  9 13:51:25 2011
 
journ
window,xs=1200,ys=1200
restore,/ver,'bb_mhd12_6km_232.save'
print,timet
;      232.016
help,by01
tvscl,by01
ploghist,by01,/prob
; % PLOT_IO: Warning: Infinite plot range.
ploghist,by01,/prob,chars=2,yr=[1.e-6,1]
ploghist,by01,/prob,chars=2,yr=[1.e-6,1],n=500
!x.style=
; % Syntax error.
!x.type=1
ploghist,by01,/prob,chars=2,yr=[1.e-6,1],n=500
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,1],n=500
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,1],n=1000
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,1],n=500
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,.1],n=500
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,.1],n=500,xr=[1.e-4,1]
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,.1],n=500
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,.1],n=500,xr=[1.e-4,10]
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,.1],n=500,xr=[1.e-4,1
; % Syntax error.
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,.1],n=500,xr=[-4,1]
ploghist,alog10(abs(by01)),/prob,chars=2,yr=[1.e-6,.1],n=500,xr=[-5,1]
ploghist,abs(by01),/prob,chars=2,n=500
; % PLOT_IO: Warning: Infinite plot range.
ploghist,abs(by01),/prob,chars=2,n=500,yr=[1.e-5,1]
ploghist,abs(by01),/prob,chars=2,n=500,yr=[1.e-6,1]
ploghist,abs(by01),/prob,chars=2,n=500,yr=[1.e-7,1]
plhist,abs(by01),/prob,chars=2,n=500,yr=[1.e-7,1]
ploghist,abs(by01),/prob,chars=2,n=500,yr=[1.e-7,1]
.r ploghist
ploghist,abs(by01),/prob,chars=2,n=500,yr=[1.e-7,1]
ploghist,abs(by01),/prob,chars=2,n=500,yr=[1.e-5,1],xr=[1.e-4,10]
ploghist,abs(by01),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-4,10]
bh01=sqrt(bx01^2+by01^2)
ploghist,bh01,/prob,chars=2,n=1000,oplt=2
help,bh01
bhnst=seeing(bh01,size=12,ap=160)
;     0.249462     0.215570
ploghist,bhnst,/prob,chars=2,n=1000,oplt=3
bynst=seeing(by01,size=12,ap=160)
;     0.246014     0.211855
ploghist,abs(bynst),/prob,chars=2,n=1000,oplt=4
x=findgen(2000)/700.
x=x[1:*]
oplot,x,x^(-1.5)*1.e-4,lin=1
oplot,x,x^(-1.5)*3.e-5,lin=1
oplot,x,x^(-1.5)*2.e-5,lin=1
oplot,x,x^(-1.5)*2.5e-5,lin=1
ploghist,abs(by01),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-4,10]
oplot,x,x^(-1.5)*2.5e-5,lin=1
bysst=seeing(by01,size=12,ap=100)
;     0.246014     0.197700
ploghist,abs(bysst),/prob,chars=2,n=1000,oplt=4
bhsst=seeing(bh01,size=12,ap=100)
;     0.249462     0.201475
ploghist,abs(bhsst),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bh01),/prob,chars=2,n=1000,oplt=2
ploghist,abs(by01),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-4,3]
!x.style=1
ploghist,abs(by01),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-4,3]
ploghist,abs(bysst),/prob,chars=2,n=1000,oplt=4
ploghist,abs(by01),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]
ploghist,abs(bysst),/prob,chars=2,n=1000,oplt=4
oplot,x,x^(-1.5)*2.5e-5,lin=1
ploghist,abs(bhsst),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bh01),/prob,chars=2,n=1000,oplt=2
!p.thick=3
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3],xtit='!5B (G)',ytit='Probability'
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]*1.e3,xtit='!5B (G)',ytit='Probability'
ploghist,abs(bh01*1.e3),/prob,chars=2,n=1000,oplt=2
ploghist,abs(bysst),/prob,chars=2,n=1000,oplt=3
ploghist,abs(byss1.e3),/prob,chars=2,n=1000,oplt=3
; % Expression must be a structure in this context: BYSS1.
ploghist,abs(bysst*1.e3),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bhsst*1.e3),/prob,chars=2,n=1000,oplt=3
label,x=.3,y=.4,lin=0,'!5Vertical'
label,lin=2,'!5Horizontal'
label,lin=3,'!5SST Vertial'
label,lin=4,'!5SST Horizontal
ps
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]*1.e3,xtit='!5B (G)',ytit='Probability'
ploghist,abs(bh01*1.e3),/prob,chars=2,n=1000,oplt=2
ploghist,abs(bysst*1.e3),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bhsst*1.e3),/prob,chars=2,n=1000,oplt=4
label,x=.3,y=.4,lin=0,'!5Vertical'
label,lin=2,'!5Horizontal'
label,lin=3,'!5SST Vertical'
label,lin=4,'!5SST Horizontal
psend
x
ps
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]*1.e3,xtit='!5B (G)',ytit='Probability'
ploghist,abs(bh01*1.e3),/prob,chars=2,n=1000,oplt=2
ploghist,abs(bysst*1.e3),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bhsst*1.e3),/prob,chars=2,n=1000,oplt=4
label,x=.4,y=.5,lin=0,'!5Vertical'
label,lin=2,'!5Horizontal'
label,lin=3,'!5SST Vertical'
label,lin=4,'!5SST Horizontal
psend
x
$mv idl.ps Bvh_hist_12v_6km.ps
help,bh01
restore,/ver,'int_mhd12_6km_232.save'
help,int
tvscl,congrid(int,800,800)
intsst=seeing(int,size=12,ap=100)
;    0.0398176    0.0315356
intnst=seeing(int,size=12,ap=160)
;    0.0398176    0.0342180
inthin=seeing(int,size=12,ap=50)
;    0.0398176    0.0257494
intatst=seeing(int,size=12,ap=400)
;    0.0398176    0.0373459
ploghist,int
; % PLOT_OO: Warning: Infinite plot range.
.r ploghist
ploghist,int,/prob
; % PLOT_IO: Warning: Infinite plot range.
intm=aver(int)
int1=int/intm
ploghist,int1,/prob,n=1000,yr=[1.e-6,1]
ploghist,int1,/prob,n=1000,yr=[1.e-8,.01]
ploghist,int1,/prob,n=1000,yr=[1.e-6,.01]
ploghist,int1,/prob,n=1000,yr=[1.e-6,.01],xr=[.5,2.1]
ploghist,int1,/prob,n=1000,yr=[1.e-6,.01],xr=[.5,2.1],chars=2
ploghist,int1,/prob,n=1000,yr=[1.e-6,.01],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
.r ploghist.pro
.r ploghist.pro
; % No commands containing 'By01' exist in the command buffer.
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]*1.e3,xtit='!5B (G)',ytit='Probability'
ploghist,abs(bh01*1.e3),/prob,chars=2,n=1000,oplt=2
ploghist,abs(bysst*1.e3),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bhsst*1.e3),/prob,chars=2,n=1000,oplt=4
label,x=.4,y=.5,lin=0,'!5Vertical'
label,lin=2,'!5Horizontal'
label,lin=3,'!5SST Vertical'
label,lin=4,'!5SST Horizontal
ps
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]*1.e3,xtit='!5B (G)',ytit='Probability'
ploghist,abs(bh01*1.e3),/prob,chars=2,n=1000,oplt=2
ploghist,abs(bysst*1.e3),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bhsst*1.e3),/prob,chars=2,n=1000,oplt=4
label,x=.4,y=.5,lin=0,'!5Vertical'
label,lin=2,'!5Horizontal'
label,lin=3,'!5SST Vertical'
label,lin=4,'!5SST Horizontal
ps
x
$mv idl.ps Bvh_hist_12v_6km.ps
ps
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]*1.e3,xtit='!5B (G)',ytit='Probability'
ploghist,abs(bh01*1.e3),/prob,chars=2,n=1000,oplt=2
ploghist,abs(bysst*1.e3),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bhsst*1.e3),/prob,chars=2,n=1000,oplt=4
label,x=.4,y=.5,lin=0,'!5Vertical'
label,lin=2,'!5Horizontal'
label,lin=3,'!5SST Vertical'
label,lin=4,'!5SST Horizontal'
psend
x
$mv idl.ps Bvh_hist_12v_6km.ps
ps
ploghist,abs(by01*1.e3),/prob,chars=2,n=1000,yr=[1.e-5,1],xr=[1.e-3,3]*1.e3,xtit='!5B (G)',ytit='Probability'
ploghist,abs(bh01*1.e3),/prob,chars=2,n=1000,oplt=2
ploghist,abs(bysst*1.e3),/prob,chars=2,n=1000,oplt=3
ploghist,abs(bhsst*1.e3),/prob,chars=2,n=1000,oplt=4
label,x=.4,y=.5,lin=0,'!5Vertical'
label,lin=2,'!5Horizontal'
label,lin=3,'!5SST Vertical'
label,lin=4,'!5SST Horizontal'
psend
x
$mv idl.ps Bvh_hist_12v_6km.ps
.r ploghist
ploghist,int1,/prob,n=1000,yr=[1.e-6,.01],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,int1,/prob,n=1000,yr=[1.e-6,1],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,int1,/prob,n=1000,yr=[1.e-4,10],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
print,intm
;     0.223028
intatstm=aver(intatst)
print,intatstm
;     0.223028
intsstm=aver(intsst)
print,intsstm
;     0.223028
inthinm=aver(inthin)
print,inthinm
;     0.223028
tvscl,congrid(inthin,800,800)
ploghist,int1,/prob,n=1000,yr=[1.e-4,10],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,intsst/intm,/prob,n=1000,oplt=2
ploghist,intatst/intm,/prob,n=1000,oplt=3
ploghist,int1,/prob,n=1000,yr=[1.e-4,10],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,intatst/intm,/prob,n=1000,oplt=2
ploghist,intnst/intm,/prob,n=1000,oplt=3
ploghist,intsst/intm,/prob,n=1000,oplt=4
ploghist,inthin/intm,/prob,n=1000,oplt=1
label,x=.3,y=.5,lin=0,'!5Simulation'
label,x=.4,y=.5,lin=0,'!5Simulation'
label,lin=2,'!5ATST'
label,lin=3,'!5NST'
label,lin=4,'!5SST'
label,lin=1,'!5Hinode'
ps
ploghist,int1,/prob,n=1000,yr=[1.e-4,10],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,intatst/intm,/prob,n=1000,oplt=2
ploghist,intnst/intm,/prob,n=1000,oplt=3
ploghist,intsst/intm,/prob,n=1000,oplt=4
ploghist,inthin/intm,/prob,n=1000,oplt=1
label,x=.4,y=.5,lin=0,'!5Simulation'
label,lin=2,'!5ATST'
label,lin=3,'!5NST'
label,lin=4,'!5SST'
label,lin=1,'!5Hinode'
psend
x
ploghist,int1,/prob,n=500,yr=[1.e-4,10],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,intatst/intm,/prob,n=500,oplt=2
ploghist,int1,/prob,n=200,yr=[1.e-4,10],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,intatst/intm,/prob,n=200,oplt=2
ploghist,intnst/intm,/prob,n=200,oplt=3
ploghist,intsst/intm,/prob,n=200,oplt=4
ploghist,inthin/intm,/prob,n=200,oplt=1
ps
ploghist,int1,/prob,n=200,yr=[1.e-4,10],xr=[.5,2.1],chars=2,xtit='!5I/<I>',ytit='!5Probability'
ploghist,intatst/intm,/prob,n=200,oplt=2
ploghist,intnst/intm,/prob,n=200,oplt=3
ploghist,intsst/intm,/prob,n=200,oplt=4
ploghist,inthin/intm,/prob,n=200,oplt=1
label,x=.45,y=.6,lin=0,'!5Simulation'
label,lin=2,'!5ATST'
label,lin=3,'!5NST'
label,lin=4,'!5SST'
label,lin=1,'!5Hinode'
psend
x
$mv idl.ps Int_hist_12v_6km.ps