;$Id: pstress_strain4.pro,v 1.2 2021/02/27 07:30:43 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=4,yoffset=3
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
end
;
;  mv idl.ps ~/tex/roper/GW/fig/pspecm_scl2.ps
;  mv idl.ps ~/tex/prop/computing/pdc/18/fig/pspecm_scl3n.ps
;
siz=1.0
!p.charsize=1.0
;!p.charsize=1.8
!x.margin=[7.9,0.5]
!y.margin=[3.2,0.4]
!p.multi=[0,2,2]
!x.title='!8k!6'
;
dir0='../'
dire='/data/'
dir1='K512irrok6c'
dir2='K512sig0_k6_ramp1a'
dir3='M512sig0_k6_ramp1b'
dir4='M512sig1_k6_ramp1b'
;
xout=14e3
yout=5e-3
yr=[1e-10,3e-2]
power,'_Str','_GWh',k=k,spec1=spec1 ,spec2=spec2 ,i=n,tt=t1,/noplot,/lkscale,datatopdir=dir0+dir1+dire
i2=n_elements(t1)-1
!y.title='!6Sp!6(!7s!6)  and  Sp(!8k!6!u2!n!8h!6)'
plot_oo,k[1:*],spec1[1:*,99],yr=yr,xtit='!6'
for i=99,(i2<189) do oplot,k[1:*],spec2[1:*,i]*k[1:*]^4,li=1,col=122
xyouts,xout,yout,siz=siz,'!6(a)'
xyouts,6e3,8.0e-5,siz=siz,'!6Sp(!8k!6!u2!n!8h!6)',col=122
xyouts,8e3,1.2e-8,siz=siz,'!6Sp(!7s!6)'
;
power,'_Str','_GWh',k=k,spec1=spec1 ,spec2=spec2 ,i=n,tt=t1,/noplot,/lkscale,datatopdir=dir0+dir2+dire
i2=n_elements(t1)-1
plot_oo,k[1:*],spec1[1:*,99],yr=yr,xtit='!6'
for i=99,(i2<189) do oplot,k[1:*],spec2[1:*,i]*k[1:*]^4,li=1,col=122
xyouts,xout,yout,siz=siz,'!6(b)'
;
power,'_Str','_GWh',k=k,spec1=spec1 ,spec2=spec2 ,i=n,tt=t1,/noplot,/lkscale,datatopdir=dir0+dir3+dire
i2=n_elements(t1)-1
!y.title='!6Sp!6(!7s!6)  and  Sp(!8k!6!u2!n!8h!6)'
plot_oo,k[1:*],spec1[1:*,99],yr=yr,xtit='!8k!6'
for i=99,(i2<189) do oplot,k[1:*],spec2[1:*,i]*k[1:*]^4,li=1,col=122
xyouts,xout,yout,siz=siz,'!6(c)'
;
power,'_Str','_GWh',k=k,spec1=spec1 ,spec2=spec2 ,i=n,tt=t1,/noplot,/lkscale,datatopdir=dir0+dir4+dire
i2=n_elements(t1)-1
plot_oo,k[1:*],spec1[1:*,99],yr=yr,xtit='!8k!6'
for i=99,(i2<189) do oplot,k[1:*],spec2[1:*,i]*k[1:*]^4,li=1,col=122
xyouts,xout,yout,siz=siz,'!6(d)'
;
print,'$mv idl.ps ~/tex/sayan/SVT/fig/pstress_strain4.ps'
END