;$Id: pstress_time_more.pro,v 1.9 2022/02/15 10:37:30 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=17,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=1.9
!x.margin=[8.8,0.8]
!y.margin=[3.2,0.4]
tilde='!9!s!AA!n!r!6'
dot='!s!9!A .!n!6!r'
siz=.9
si2=.9
si3=.8
xlab1=.53
xlab2=.005
;
dir1='2048_1e2_1e4_4e8_1em6aD_k10_bes'
dir2='2048_1em1_2e2_16e4_1em4a_tet'
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'
  power,'_GWs','_GWs',k=k1,spec1=grav1,spec2=grav1,i=n,tt=tstr1,/noplot,/lks,datatopdir='../'+dir1+'/data'
  power,'_GWs','_GWs',k=k2,spec1=grav2,spec2=grav2,i=n,tt=tstr2,/noplot,/lks,datatopdir='../'+dir2+'/data'
  power,'_GWh','_GWh',k=k1,spec1=grah1,spec2=grah1,i=n,tt=tstr1,/noplot,/lks,datatopdir='../'+dir1+'/data'
  power,'_GWh','_GWh',k=k2,spec1=grah2,spec2=grah2,i=n,tt=tstr2,/noplot,/lks,datatopdir='../'+dir2+'/data'
  lread=0
endif
;
!p.multi=[0,2,3]
li=[0,1,2,0,0,0,0]
col=[col0,col0,col0,55,122,155,122]
lct=[0,0,0,0,1,2,0]
circ_sym,.9,1
ps=8
;
kk1=[1,5,20,100,500]
kk2=[1,4,20,100,400]
;
kk1=[1,3,10,40,100,200,500]
kk2=[1,3,10,40,100,200,400]
;kk2=[1,2,5,10,20]
nkk1=n_elements(kk1)
nkk2=n_elements(kk2)
;=================================================================================
!x.title='!6'
!y.title='!8T!6!drms!n(!8k!6,!8t!6)'
;
xx=.8 & dx=.14 & yy=10.^(-8.5-.50*findgen(nkk1))
xr1=[.5,3.2]
yr_mag=[5d-13,4d-5]
plot_oo,xr1,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.,(tstr1/6.)*sqrt(stre1[kk1[ikk],*])/k1[kk1[ikk]],col=col[ikk],li=li[ikk]
  legend,xx,dx,yy[ikk],li[ikk],'!8k!6='+str(fix(k1[kk1[ikk]])),col=col[ikk],siz=si3
  print,k1[kk1[ikk]]
endfor
yy=1.15e-13
xyouts,0.47,yy,'!60.5',siz=si2
xyouts,1.95,yy,'!62',siz=si2
xyouts,xlab1,4e-6,'!6(a)',siz=siz
;xx=grange(2.20,3.8,200) & oplot,xx-1.,2e-6*xx^(1.),col=122
xx=grange(1.20,2.8,200) & oplot,xx,4e-6*xx^(1.)
xyouts,1.3,1e-5,siz=siz,'!9A!8t!6-1'
;
;---------------------------------------------------------------------------------
yr_mag=[1d-12,2d-5]
xx=.02 & dx=.014 & yy=10.^(-8.2-.48*findgen(nkk1))
xr2=[.004,6.]
xr2b=xr2
plot_oo,xr2,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.,(tstr2/6.)*sqrt(stre2[kk2[ikk],*])/k2[kk2[ikk]],col=col[ikk],li=li[ikk]
  legend,xx,dx,yy[ikk],li[ikk],'!8k!6='+str(nint(k2[kk2[ikk]])),col=col[ikk],siz=si3
  print,k2[kk2[ikk]]
endfor
xyouts,xlab2,2e-6,'!6(b)',siz=siz
;xx=grange(1.08,2.5,200) & oplot,xx-1.,3e-6*xx^(1./3.),col=122
xx=grange(0.08,1.5,200) & oplot,xx,3e-6*xx^(1./3.)
xyouts,.1,3e-6,siz=siz,'!9A!6(!8t!6-1)!u1/3!n'
;
;=================================================================================
!y.title=dot+'!8h!6!drms!n(!8k!6,!8t!6)'
;
xx=.9 & dx=.5 & yy=10.^(-9.2-.48*findgen(nkk1))
xr=[.5,3.2]
yr_mag=[5d-18,2d-6]
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(grav1[kk1[ikk],*])/k1[kk1[ikk]],col=col[ikk],li=li[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=4.5e-19
xyouts,0.47,yy,'!60.5',siz=si2
xyouts,1.95,yy,'!62',siz=si2
xyouts,xlab1,1e-7,'!6(c)',siz=siz
;xyouts,.310,2.8e-13,'!6.',siz=siz
;
;---------------------------------------------------------------------------------
yr_mag=[5d-11,4d-6]
xx=.06 & dx=.04 & yy=10.^(-8.2-.44*findgen(nkk1))
xr=[.02,5.]
plot_oo,xr2b,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(grav2[kk2[ikk],*])/k2[kk2[ikk]],col=col[ikk],li=li[ikk]
  ;legend,xx,dx,yy[ikk],0,'!8k!6='+str(nint(k2[kk2[ikk]])),col=col[ikk],siz=siz
  print,k2[kk2[ikk]]
endfor
xyouts,xlab2,1.3e-6,'!6(d)',siz=siz
arrow,.1,1e-10,.064,2e-10,/data
xyouts,.112,8e-11,'!610x lower cadence',siz=si3
;xyouts,6e-4,5e-9,'!6.',siz=siz
;
;=================================================================================
!x.title='!8t!6-1'
!y.title='!8h!6!drms!n(!8k!6,!8t!6)'
;
xx=.9 & dx=.5 & yy=10.^(-9.2-.48*findgen(nkk1))
xr=[.5,3.2]
yr_mag=[5d-20,2d-7]
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(grah1[kk1[ikk],*])/k1[kk1[ikk]],col=col[ikk],li=li[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=4.6e-21
xyouts,0.47,yy,'!60.5',siz=si2
xyouts,1.95,yy,'!62',siz=si2
xyouts,xlab1,1e-8,'!6(e)',siz=siz
;
;---------------------------------------------------------------------------------
yr_mag=[5d-14,1d-6]
yr_mag=[5d-15,1d-6]
xx=.06 & dx=.03 & yy=10.^(-8.2-.44*findgen(nkk1))
xr=[.02,5.]
plot_oo,xr2,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(grah2[kk2[ikk],*])/k2[kk2[ikk]],col=col[ikk],li=li[ikk]
  ;legend,xx,dx,yy[ikk],0,'!8k!6='+str(nint(k2[kk2[ikk]])),col=col[ikk],siz=siz
  print,k2[kk2[ikk]]
endfor
xyouts,xlab2,2e-7,'!6(f)',siz=siz
;
!p.multi=0
print,'$mv idl.ps ~/GitHub/Axel/GWspectra_from_CME/fig/pstress_time_more.eps'
END