;$Id: pcomp_seriesB.pro,v 1.6 2021/02/08 17:23:42 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=25,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=2.8
!x.margin=[8.8,0.5]
!y.margin=[3.2,0.3]
!x.title='!6'
!p.multi=[0,1,3]
si2=1.0
siz=1.3
;
;  sample to compare time series form different directories
;
dir1='512_1e2_1e4_4e8_1em6a'
dir2='512_1e2_1e4_4e8_2em6a'
dir3='512_1e2_1e4_4e8_5em6a'
dir4='512_1e2_1e4_4e8_1em5a'
dir5='512_1e2_1e4_4e8_2em5a'
dir6='512_1e2_1e4_4e8_5em5a'
dir7='512_1e2_1e4_4e8_1em4a'
dir8='512_1e2_1e4_4e8_2em4a'
dir9='512_1e2_1e4_4e8_5em4a'
dir0='512_1e2_1e4_4e8_1em3a'
;
mu=1d4
mu2=mu^2
default,iread,0
if iread eq 0 then begin
  pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
  pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
  pc_read_ts,o=ts3,datadir='../'+dir3+'/data'
  pc_read_ts,o=ts4,datadir='../'+dir4+'/data'
  pc_read_ts,o=ts5,datadir='../'+dir5+'/data'
  pc_read_ts,o=ts6,datadir='../'+dir6+'/data'
  pc_read_ts,o=ts7,datadir='../'+dir7+'/data'
  pc_read_ts,o=ts8,datadir='../'+dir8+'/data'
  pc_read_ts,o=ts9,datadir='../'+dir9+'/data'
  pc_read_ts,o=ts0,datadir='../'+dir0+'/data'
  iread=1
endif
s1=1e-6*mu2
s2=2e-6*mu2
s3=5e-6*mu2
s4=1e-5*mu2
s5=2e-5*mu2
s6=5e-5*mu2
s7=1e-4*mu2
s8=2e-4*mu2
s9=5e-4*mu2
s0=1e-3*mu2
t=ts7.t*1d0 & nt=n_elements(t) & t7=(grange(t[0],t[nt-1],nt)-1.)*s7
t=ts8.t*1d0 & nt=n_elements(t) & t8=(grange(t[0],t[nt-1],nt)-1.)*s8
t=ts9.t*1d0 & nt=n_elements(t) & t9=(grange(t[0],t[nt-1],nt)-1.)*s9
t=ts0.t*1d0 & nt=n_elements(t) & t0=(grange(t[0],t[nt-1],nt)-1.)*s0
xx2=462.
xx1=80.
x1=480.
;
!y.title='!13E!6!dM!n(!8t!6)'
;xr=[0.,90.] & yr=[2d-7,2d-2] & yr_EEGW=[2d-42,2d-8]
xr=[70.,490.]
yr=[2d-3,1d-1]
yr_EEGW=[2d-12,8d-9]
yr_mu=[1e2,1.2e4]
y1=.023 & y2=.036 & y3=1.1*yr[0] & y4=.013
xbox=[x1,xr[1],xr[1],x1,x1]
plot_io,xr,yr,/nodata & loadct,0
polyfill,xbox,[y1,y1,y2,y2,y1],col=199
polyfill,xbox,[y3,y3,y4,y4,y3],col=233 & loadct,5
oplot,s1*(ts1.t-1.),ts1.EEM,li=1
oplot,s2*(ts2.t-1.),ts2.EEM,li=2
oplot,s3*(ts3.t-1.),ts3.EEM
oplot,s4*(ts4.t-1.),ts4.EEM,col=55 & loadct,6
oplot,s5*(ts5.t-1.),ts5.EEM,col=122 & loadct,5
oplot,s6*(ts6.t-1.),ts6.EEM,col=155
oplot,t7,ts7.EEM,col=122
oplot,t8,ts8.EEM,col=122,li=3
oplot,t9,ts9.EEM,col=122,li=2
oplot,t0,ts0.EEM,li=1,thick=8
;xx=[10.,80.] & oplot,xx,1e-20*exp(.5*xx)
xyouts,170,.0042,siz=si2,'!6B1'
xyouts,220,.0063,siz=si2,'!6B2'
xyouts,275,.0064,siz=si2,'!6B3'
xyouts,300,.0112,siz=si2,'!6B4',col=55 & loadct,6
xyouts,320,.016,siz=si2,'!6B5',col=122 & loadct,5
xyouts,340,.023,siz=si2,'!6B6',col=155
xyouts,360,.026,siz=si2,'!6B7',col=122
xyouts,380,.038,siz=si2,'!6B10'
xyouts,430,4.0e-2,siz=si2,'!6regime II'
xyouts,430,4.2e-3,siz=si2,'!6regime I'
xyouts,xx2,7.0e-2,siz=siz,'!6(a)'
;
f=1.
!y.title='!13E!6!dGW!n(!8t!6)'
y1=7e-10 & y2=6e-9 & y3=1.2*yr_EEGW[0] & y4=4e-10
plot_io,xr,yr_EEGW,/nodata & loadct,0
polyfill,xbox,[y1,y1,y2,y2,y1],col=199
polyfill,xbox,[y3,y3,y4,y4,y3],col=233 & loadct,5
oplot,s1*(ts1.t-1.),f*ts1.EEGW,li=1
oplot,s2*(ts2.t-1.),f*ts2.EEGW,li=2
oplot,s3*(ts3.t-1.),f*ts3.EEGW
oplot,s4*(ts4.t-1.),f*ts4.EEGW,col=55 & loadct,6
oplot,s5*(ts5.t-1.),f*ts5.EEGW,col=122 & loadct,5
oplot,s6*(ts6.t-1.),f*ts6.EEGW,col=155
oplot,t7,f*ts7.EEGW,col=122
oplot,t8,f*ts8.EEGW,col=122,li=3
oplot,t9,f*ts9.EEGW,col=122,li=2
oplot,t0,f*ts0.EEGW,li=1,thick=8
;xx=[10.,80.] & oplot,xx,1d-45*exp(xx)
xyouts,150,7.0e-12,siz=si2,'!6B1'
xyouts,215,1.6e-11,siz=si2,'!6B2'
xyouts,280,9.0e-11,siz=si2,'!6B3'
xyouts,260,4e-10,siz=si2,'!6B7',col=122
xyouts,240,8e-10,siz=si2,'!6B8',col=122
xyouts,220,1.7e-9,siz=si2,'!6B9',col=122
xyouts,160,2.0e-9,siz=si2,'!6B10'
xyouts,430,1.2e-9,siz=si2,'!6regime II'
xyouts,430,1.1e-11,siz=si2,'!6regime I'
xyouts,xx1,3.6e-9,siz=siz,'!6(b)'
;
f=1.
!y.title='!6<!7l!6>(!8t!6)'
!x.title='!64!7c!6!d0!n!8t!6'
y1=1e3 & y2=2e3 & y3=1.1*yr_mu[0] & y4=6e2
plot_io,xr,yr_mu,/nodata & loadct,0
polyfill,xbox,[y1,y1,y2,y2,y1],col=199
polyfill,xbox,[y3,y3,y4,y4,y3],col=233 & loadct,5
oplot,s1*(ts1.t-1.),f*ts1.mu5m,li=1
oplot,s2*(ts2.t-1.),f*ts2.mu5m,li=2
oplot,s3*(ts3.t-1.),f*ts3.mu5m
oplot,s4*(ts4.t-1.),f*ts4.mu5m,col=55 & loadct,6
oplot,s5*(ts5.t-1.),f*ts5.mu5m,col=122 & loadct,5
oplot,s6*(ts6.t-1.),f*ts6.mu5m,col=155
oplot,t7,f*ts7.mu5m,col=122
oplot,t8,f*ts8.mu5m,col=122,li=3
oplot,t9,f*ts9.mu5m,col=122,li=2
oplot,t0,f*ts0.mu5m,li=1,thick=8
xyouts,210,  170,siz=si2,'!6B1'
xyouts,240,  270,siz=si2,'!6B2'
xyouts,280,  350,siz=si2,'!6B3'
xyouts,300,  480,siz=si2,'!6B4',col=55 & loadct,6
xyouts,320,  700,siz=si2,'!6B5',col=122 & loadct,5
xyouts,340,  810,siz=si2,'!6B6',col=155
xyouts,360, 1150,siz=si2,'!6B7',col=122
xyouts,380, 1950,siz=si2,'!6B10'
xyouts,430,2.2e3,siz=si2,'!6regime II'
xyouts,430,1.6e2,siz=si2,'!6regime I'
xyouts,xx2,7.6e3,siz=siz,'!6(c)'
;
cwd,run
!x.range=0
!p.multi=0
print,'$mv idl.ps ~/tex/rei/GW/fig/pcomp_seriesB.ps'
END