;
; This script will compare the alphas for different runs and plot
; them in a common plot.
;
default,png,0
size=2
if !d.name eq 'PS' then begin
   device,xsize=20,ysize=16,yoffset=3
   !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
   !p.charsize=size
   thick2=8
endif else begin
   !p.charsize=size
   device,decompose=0
   loadct,5
   thick2=2
   if (png) then begin
      !p.charthick=1 & !p.thick=4 & !x.thick=4 & !y.thick=4
   endif
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='/256a3/'
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]
;
dir='/128a16/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_240_128_mu_temp    =kf
St_240_128_mu_temp    =Stokes
alpha_240_128_mu_temp =alpha
alpha2_240_128_mu_temp=alpha2
;
dir='/128a15/'
print,'Reading: ',dir
file=maindir+dir+varfile
restore,file=file
kf_240_128_mu_temp    =[kf_240_128_mu_temp,kf]
St_240_128_mu_temp    =[St_240_128_mu_temp,Stokes]
alpha_240_128_mu_temp =[alpha_240_128_mu_temp,alpha]
alpha2_240_128_mu_temp=[alpha2_240_128_mu_temp,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]


col=intarr(7)
if png then begin
   xmax=200.
   backgroundcol=255
   for i=0,6 do begin
      col[i]=10+i*30
   end
endif else begin
   xmax=50.
   if !d.name eq 'PS' then begin
      col[*]=1
      backgroundcol=255
   endif else begin
      col[*]=255
      backgroundcol=0
   end
end

;
; Plot results
;
!p.multi=[0,1,1]
plot_oi,St_80_128,alpha_80_128,yr=[-0.3,1.2],xr=[0.05,xmax],$
  xtit='St',ytit='!4a!6',li=0,xstyle=1,ystyle=1, background=backgroundcol,$
   col=col[0]     
oplot,St_240_128,alpha_240_128,li=5,col=col[1] ;,ps=-1
oplot,St_800_128,alpha_800_128,li=2,col=col[2];,ps=-2
oplot,St_80_256_15,alpha_80_256_15,li=3,col=col[3];,ps=-3
oplot,St_240_128_mu,alpha_240_128_mu,li=4,col=col[4];,ps=-4
oplot,St_240_128_mu_temp,alpha_240_128_mu_temp,li=1,col=col[0];,thick=thick2
;oplot,St_240_128_mu,alpha_240_128_mu*0.0,li=1,thick=2,col=col[0]
;oplot,St_240_128_musqrt,alpha_240_128_musqrt,li=5;,ps=-5
;
; Legends
;
xx=0.07
dx=0.065
yy=0.65
dy=-0.12
legend,xx,dx,yy+dy*0.,0,'Re= 80, k!df!n=5 ',col=col[0];, ps=-0
legend,xx,dx,yy+dy*1.,5,'Re=240, k!df!n=5 ',col=col[1];, ps=-1
legend,xx,dx,yy+dy*2.,2,'Re=800, k!df!n=5 ',col=col[2];, ps=-2
legend,xx,dx,yy+dy*3.,3,'Re= 80, k!df!n=15',col=col[3];, ps=-3
legend,xx,dx,yy+dy*4.,4,'Re=240, k!df!n=5, !4l!6=const',col=col[4];, ps=-4
legend,xx,dx,yy+dy*5.,1,'Re=240, k!df!n=5, temp',col=col[4];,thick=2;, 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,alpha2_80_128,yr=[-0.5,2.],xr=[0.05,200.],$
;;   xtit='St',ytit='!4a!6!d2!n',li=0,xstyle=1,ystyle=1;,ps=-0
;; oplot,St_240_128,alpha2_240_128,li=5;,ps=-1
;; oplot,St_800_128,alpha2_800_128,li=2;,ps=-2
;; oplot,St_80_256_15,alpha2_80_256_15,li=3;,ps=-3
;; oplot,St_240_128_mu,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,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]

if (png) then begin
   write_png,'fig1.png',tvrd(/TRUE)
endif

;
END