;
; This script will compare the alphas for different runs and plot
; them in a common plot.
;
size=2
if !d.name eq 'PS' then begin
   device,xsize=20,ysize=30,yoffset=3
   !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
   !p.charsize=size
endif else begin
   !p.charsize=size
   device,decompose=0
   loadct,5
end
varfile='data/ppar_real.dat'
;
; Read data from the different run-directories with forcing at k=5
;
maindir='./forcing_at_k5'
;
dir='/128a3/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_240_128=kf
St_240_128=Stokes
alpha_240_128=alpha
alpha2_240_128=alpha2
St_240_128[0]=Stokes[1]
St_240_128[1]=Stokes[0]
alpha_240_128[0]=alpha[1]
alpha_240_128[1]=alpha[0]
alpha2_240_128[0]=alpha2[1]
alpha2_240_128[1]=alpha2[0]
;
dir='/128a1/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_240_128=[kf_240_128,kf]
St_240_128=[St_240_128,Stokes]
alpha_240_128=[alpha_240_128,alpha]
alpha2_240_128=[alpha2_240_128,alpha2]
;
dir='/128a4/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_800_128=kf
St_800_128=Stokes
alpha_800_128=alpha
alpha2_800_128=alpha2
;
dir='/128a2/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_800_128    =[kf_800_128,kf]
St_800_128    =[St_800_128,Stokes]
alpha_800_128 =[alpha_800_128,alpha]
alpha2_800_128=[alpha2_800_128,alpha2]
;
dir='/128a6/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_80_128    =kf
St_80_128    =Stokes
alpha_80_128 =alpha
alpha2_80_128=alpha2
;
dir='/128a5/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_80_128    =[kf_80_128,kf]
St_80_128    =[St_80_128,Stokes]
alpha_80_128 =[alpha_80_128,alpha]
alpha2_80_128=[alpha2_80_128,alpha2]
;
dir='/128a9/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_240_128_mu    =kf
St_240_128_mu    =Stokes
alpha_240_128_mu =alpha
alpha2_240_128_mu=alpha2
;
dir='/128a8/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_240_128_mu    =[kf_240_128_mu,kf]
St_240_128_mu    =[St_240_128_mu,Stokes]
alpha_240_128_mu =[alpha_240_128_mu,alpha]
alpha2_240_128_mu=[alpha2_240_128_mu,alpha2]
;
dir='/128a10/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_240_128_musqrt    =kf
St_240_128_musqrt    =Stokes
alpha_240_128_musqrt =alpha
alpha2_240_128_musqrt=alpha2
;
dir='/128a14/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_500_128_mu_lowmach    =kf
St_500_128_mu_lowmach    =Stokes
alpha_500_128_mu_lowmach =alpha
alpha2_500_128_mu_lowmach=alpha2
;
dir='/128a13/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_500_128_mu_lowmach    =[kf_500_128_mu_lowmach,kf]
St_500_128_mu_lowmach    =[St_500_128_mu_lowmach,Stokes]
alpha_500_128_mu_lowmach =[alpha_500_128_mu_lowmach,alpha]
alpha2_500_128_mu_lowmach=[alpha2_500_128_mu_lowmach,alpha2]
;
; Forcing at k=15
;
maindir='./forcing_at_k15'
;
 dir='/256a3/'
 print,'Reading: ',dir
 file=maindir+dir+varfile
 restore,file=file
 kf_80_256_15=kf
 St_80_256_15=Stokes
 alpha_80_256_15=alpha
 alpha2_80_256_15=alpha2
;
 dir='/256a2/'
 print,'Reading: ',dir
 file=maindir+dir+varfile
 restore,file=file
 kf_80_256_15=[kf_80_256_15,kf]
 St_80_256_15=[St_80_256_15,Stokes]
 alpha_80_256_15=[alpha_80_256_15,alpha]
 alpha2_80_256_15=[alpha2_80_256_15,alpha2]
;
; Plot results
;
!p.multi=[0,1,2]
plot_oi,St_80_128/sqrt(80),alpha_80_128,yr=[-0.3,1.2],xr=[0.003,20.],$
  xtit='St',ytit='!4a!6',li=0,xstyle=1,ystyle=1;,ps=-0
oplot,St_240_128/sqrt(240),alpha_240_128,li=5;,ps=-1
oplot,St_800_128/sqrt(800),alpha_800_128,li=2;,ps=-2
oplot,St_80_256_15/sqrt(80),alpha_80_256_15,li=3;,ps=-3
oplot,St_240_128_mu/sqrt(240),alpha_240_128_mu,li=4;,ps=-4
oplot,St_240_128_mu/sqrt(240),alpha_240_128_mu*0.0,li=1,thick=2
;oplot,St_240_128_musqrt,alpha_240_128_musqrt,li=5;,ps=-5
;
; Legends
;
xx=0.007
dx=0.0065
yy=0.6
dy=-0.12
legend,xx,dx,yy+dy*0.,0,'Re= 80, k!df!n=5 ';, ps=-0
legend,xx,dx,yy+dy*1.,5,'Re=240, k!df!n=5 ';, ps=-1
legend,xx,dx,yy+dy*2.,2,'Re=800, k!df!n=5 ';, ps=-2
legend,xx,dx,yy+dy*3.,3,'Re= 80, k!df!n=15';, ps=-3
legend,xx,dx,yy+dy*4.,4,'Re=240, k!df!n=5, !4l!6=const';, ps=-4
;legend,xx,dx,yy+dy*5.,5,'Re=240, k!df!n=5, sqrt(!4q!6)!4m!6=const';, ps=-5
;
plot_oi,St_80_128/sqrt(80),alpha2_80_128,yr=[-0.5,2.],xr=[0.003,20.],$
  xtit='St',ytit='!4a!6!d2!n',li=0,xstyle=1,ystyle=1;,ps=-0
oplot,St_240_128/sqrt(240),alpha2_240_128,li=5;,ps=-1
oplot,St_800_128/sqrt(800),alpha2_800_128,li=2;,ps=-2
oplot,St_80_256_15/sqrt(80),alpha2_80_256_15,li=3;,ps=-3
oplot,St_240_128_mu/sqrt(240),alpha2_240_128_mu,li=4;,ps=-4
;oplot,St_240_128_musqrt,alpha2_240_128_musqrt,li=5;,ps=-5
oplot,St_240_128_mu/sqrt(240),alpha2_240_128_mu*0.0,li=1,thick=1;,ps=-5
;
;
print,'mv idl.ps ~/tex/axel/nathan_igor/nils/fig/fig1.ps'
;
!p.multi=[0,1,1]
;
END