;$Id: pstress_time.pro,v 1.4 2022/01/23 12:33:32 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=7,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
  col0=0
endif else begin
  col0=255
end
;
!p.charsize=.8
!x.margin=[8.0,0.3]
!y.margin=[3.2,0.3]
tilde='!9!s!aA!n!r!6'
siz=.9
;
dir1='2048_1e2_1e4_4e8_1em6aD_k10_bes'
dir2='2048_1e1_1e4_4e8_1em4a'
;
default,lread,1
if lread then begin
  power,'_Str','_Str',k=k1,spec1=stre1,spec2=stre1,i=n,tt=tstr1,/noplot,/lks,datatopdir='../'+dir1+'/data'
  power,'_Str','_Str',k=k2,spec1=stre2,spec2=stre2,i=n,tt=tstr2,/noplot,/lks,datatopdir='../'+dir2+'/data'
  lread=0
endif
;
!p.multi=[0,2,1]
col=[col0,55,122,155,122]
lct=[0,0,1,2,0]
dot='!s!9!A .!n!6!r'
circ_sym,.9,1
ps=8
;
kk1=[1,5,20,100,500]
kk2=[1,4,20,100,400]
;
kk1=[1,50,100,200,500]
kk2=[1,50,100,200,400]
nkk1=n_elements(kk1)
nkk2=n_elements(kk2)
;---------------------------------------------------------------------------------
!x.title='!8t!6-1'
!y.title=tilde+'!8T!6(!8k!6,!8t!6)'
;
xx=.9 & dx=.5 & yy=10.^(-9.2-.48*findgen(nkk1))
xr=[.5,3.2]
yr_mag=[5d-13,2d-5]
plot_oo,xr,yr_mag,/nodata
for ikk=0,nkk1-1 do begin
  if lct[ikk] eq 1 then loadct,6
  if lct[ikk] eq 2 then loadct,5
  oplot,tstr1-1.,sqrt(stre1[kk1[ikk],*])/k1[kk1[ikk]],col=col[ikk]
  legend,xx,dx,yy[ikk],0,'!8k!6='+str(fix(k1[kk1[ikk]])),col=col[ikk],siz=siz
  print,k1[kk1[ikk]]
endfor
yy=1.23e-13
xyouts,0.47,yy,'!60.5'
xyouts,1.95,yy,'!62'
;oplot,tstr1-1.,1./tstr1^1.5
oplot,tstr2-1.,1e-7/(tstr2-1.)^3,li=3
;
;---------------------------------------------------------------------------------
yr_mag=[5d-12,1d-5]
xx=.06 & dx=.03 & yy=10.^(-8.2-.44*findgen(nkk1))
xr=[.02,5.]
plot_oo,xr,yr_mag,/nodata
for ikk=0,nkk2-1 do begin
  if lct[ikk] eq 1 then loadct,6
  if lct[ikk] eq 2 then loadct,5
  oplot,tstr2-1.,sqrt(stre2[kk2[ikk],*])/k2[kk2[ikk]],col=col[ikk]
  legend,xx,dx,yy[ikk],0,'!8k!6='+str(nint(k2[kk2[ikk]])),col=col[ikk],siz=siz
  print,k2[kk2[ikk]]
endfor
oplot,tstr2-1.,1e-9/(tstr2-1.)^3,li=3
;
!p.multi=0
print,'$mv idl.ps ~/GitHub/Axel/GWspectra_from_CME/fig/pstress_time.eps'
END