;$Id: pcomp_stress.pro,v 1.5 2022/10/05 01:39:18 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=19,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
  col0=255
endif else begin
  col0=0
endelse
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
dir0='P1024_k1_kf10d_stress1b'
dir1='P1024_k1_kf10d_stress07'
dir2='P1024_k1_kf10d_stress05b'
dir2b='P1024_k1_kf10d_stress05_nophase'
dir3='P1024_k1_kf10d_stress03'
dir4='P1024_k1_kf10d_stress0'
;
pc_read_ts,o=ts0,datadir='../'+dir0+'/data'
pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
pc_read_ts,o=ts2b,datadir='../'+dir2b+'/data'
pc_read_ts,o=ts3,datadir='../'+dir3+'/data'
pc_read_ts,o=ts4,datadir='../'+dir4+'/data'
;
!p.multi=[0,1,2]
!x.range=[1,120]
!y.range=[3e-6,.02]
!y.title='rms of stress, Trms'
plot_oo,ts1.t,ts1.nlin1
oplot,ts2.t,ts2.nlin1,col=122
oplot,ts2.t,ts2b.nlin1,col=122,li=2,th=5 & loadct,6
oplot,ts4.t,ts4.nlin1,col=122 & loadct,5
oplot,ts3.t,ts3.nlin1,col=55
oplot,ts0.t,ts0.nlin1,col=155
xx=[16.,50.] & oplot,xx,1e-2/xx^1.2,col=122
;
;  GW spectra
;
!y.range=[3e-8,4e-4]
!y.title='!13E!6!dGW!n'
plot_oo,ts1.t,ts1.EEGW
;plot,ts1.t,ts1.EEGW
oplot,ts2.t,ts2.EEGW,col=122
oplot,ts2b.t,ts2b.EEGW,col=122,li=2 & loadct,6
oplot,ts4.t,ts4.EEGW,col=122 & loadct,5
oplot,ts3.t,ts3.EEGW,col=55
oplot,ts0.t,ts0.EEGW,col=155
;
xx=3. & dx=xx & yy=1e-6 & dy=3.
legend,/log,xx,dx,yy,0,col=155,'n=1'
legend,/log,xx,dx,yy/dy,0,'n=0.7'
legend,/log,xx,dx,yy/dy^2,0,col=122,'n=0.5'
;
xx=13. & dx=xx & yy=1e-6 & dy=3.
legend,/log,xx,dx,yy,0,col=55,'n=0.3' & loadct,6
legend,/log,xx,dx,yy/dy,0,col=122,'n=0' & loadct,5
;legend,/log,xx,dx,yy/dy^2,0,col=55,'n=0.3'
xx=[16.,50.] & oplot,xx,1e-5*xx^.4,col=122
;
print,'$mv idl.ps ~/GitHub/Axel/VortexGen/fig/pcomp_stress.eps'
END