;$Id: pstress_etc_B7post.pro,v 1.1 2022/05/09 00:36:11 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=10,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
;
;  this is to be run in: 2048_1e1_1e4_4e8_1em4a
;
!p.charsize=0.9
!x.margin=[10.5,1.8]
!y.margin=[3.2,0.3]
siz=.9
;
@parameters
default,lread,1
if lread then begin
  ;power,'_mag','_Str',k=k,spec1=spec1,spec2=stre1,i=n,tt=t,/noplot,/lks
  ;power,'_Str','_Str',k=k,spec1=stre1,spec2=stre1,i=n,tt=tstr,/noplot,/lks
  ;power,'_mag','_mag',k=k,spec1=spec1,spec2=spec1,i=n,tt=tmag,/noplot,/lks
  power,'_GWs','_GWh',k=k,spec1=grav1,spec2=grah1,i=n,tt=t,/noplot,/lks
  pc_read_ts,obj=ts
  lread=0
endif
;
;  select times from t array
;
ntt=n_elements(tt)
it_ts=intarr(ntt)
it_str=intarr(ntt)
it_mag=intarr(ntt)
it_gws=intarr(ntt)
for itt=0,ntt-1 do begin
  it_gws[itt]=findex(tt[itt],t)
  it_ts[itt]=findex(tt[itt],ts.t)
endfor
print,it_str
print,it_mag
print,it_gws
;
xlab1=6.2
xlab1b=2500.
xlab2a=9e-3
xlab2b=2.5
xrt=[0e0,80.]
xr=[k[1],max(k)]
!p.multi=[0,2,2]
col=[col0,col0,col0,col0,55,122,155]
dot='!s!9!A .!n!6!r'
circ_sym,.9,1
ps=8
;
;---------------------------------------------------------------------------------
!y.title='!6Sp('+dot+'!8h!6)'
;!y.title='!8E!6!dM!n(!8k!6)'
yr_GWs=[1d-14,7d-11]
plot_oo,xr,yr_GWs,/nodata
for itt=0,ntt-1 do oplot,k[1:*],grav1[1:*,it_GWs[itt]],col=col[itt]
oplot,[1,1]*5000.,yr_GWs,li=3
;xyouts,siz=siz,56.,1e-20,'!8k!6!d0!n'
;xyouts,siz=siz,3000.,5e-12,'!8k!7!dl!n!6'
xyouts,xlab1b,1.7e-11,siz=siz,'!6(e)'
xx=[14.,80.] & oplot,xx,1.5e-10/xx^0.5,li=3
xyouts,60.,2.3e-11,siz=siz,'!9A!8k!6!u-0.5!n'
;
yr_GWs_ts=[5.4d-10,5.8d-10]
!y.title='!13E!6!dGW!n'
tot_grav=total(grav1,1)/6.
;plot_oo,t-1.,tot_grav,xr=xrt,yr=[5.4e-10,5.6e-10]
plot,t-1.,tot_grav,xr=xrt,yr=[5.58e-10,5.7e-10]
for itt=0,ntt-1 do oplot,[1,1]*t[it_GWs[itt]]-1.,[1,1]*tot_grav[it_GWs[itt]],col=col[itt],ps=ps
xyouts,xlab2a,5.5e-10,siz=siz,'!6(f)'
;
;---------------------------------------------------------------------------------
!x.title='!8k!6'
!y.title='!6Sp(!8h!6)'
yr_GWh=[1d-20,3d-12]
plot_oo,xr,yr_GWh,/nodata
for itt=0,ntt-1 do oplot,k[1:*],grah1[1:*,it_GWs[itt]],col=col[itt]
oplot,[1,1]*5000.,yr_GWh,li=3
;xyouts,siz=siz,56.,1e-24,'!8k!6!d0!n'
;xyouts,siz=siz,3000.,1e-14,'!8k!7!dl!n!6'
xyouts,xlab1b,15e-14,siz=siz,'!6(g)'
xx=[12.,150.] & oplot,xx,4e-10/xx^2.3,li=3
xyouts,80.,3.3e-14,siz=siz,'!9A!8k!6!u-2.3!n'
;
!x.title='!8t!6-1'
!y.title='!8h!6!s!u2!n!r!drms!n'
plot,ts.t-1.,ts.hrms^2,xr=xrt,yr=[1e-12,2e-12]
for itt=0,ntt-1 do oplot,[1,1]*ts.t[it_ts[itt]]-1.,[1,1]*ts.hrms[it_ts[itt]]^2,col=col[itt],ps=ps
xyouts,xlab2a,2e-12,siz=siz,'!6(h)'
;
print,'$mv idl.ps ~/GitHub/Axel/GWspectra_from_CME/fig/pstress_etc_B7post.eps'
END