;$Id: pstress_strain.pro,v 1.1 2021/02/26 20:48:37 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=8,yoffset=3
  device,xsize=18,ysize=11,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
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.7
@parameters
!p.charsize=1.7
;!p.charsize=1.8
!x.margin=[6.6,0.5]
!y.margin=[3.2,0.4]
!x.title='!8k!6'
;
power,'_Str','_GWh',k=k,spec1=spec1 ,spec2=spec2 ,i=n,tt=t1,/noplot,/lkscale
i2=n_elements(t1)-1
plot_oo,k[1:*],spec1[1:*,99],yr=[1e-12,2e-2],xtit='!8k!6'
for i=99,(i2<189) do oplot,k[1:*],spec2[1:*,i]*k[1:*]^4,li=1
;
cwd,run
print,'$mv idl.ps ~/tex/sayan/SVT/fig/pstress_strain_'+run+'.ps'
END