;$Id: pSCL.pro,v 1.2 2019/12/21 12:16:53 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=8,yoffset=3
  device,xsize=18,ysize=16,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.3
@parameters
!p.charsize=1.7
;!p.charsize=1.8
!x.margin=[7.6,0.5]
!y.margin=[3.2,0.2]
;!p.multi=[0,1,2]
;
;power,'_GWh','_GWh',k=k,spec1=spec1 ,spec2=spec1 ,i=n,tt=t,/noplot
;power,'_SCL','_SCL',k=k,spec1=spec1 ,spec2=spec1 ,i=n,tt=t,/noplot
;power,'_SCL','_SCL',k=k,spec1=spec2 ,spec2=spec2 ,i=n,tt=t,/noplot
;power,'_VCT','_VCT',k=k,spec1=spec2 ,spec2=spec2 ,i=n,tt=t,/noplot
;power,'_VCT','_VCT',k=k,spec1=spec2 ,spec2=spec2 ,i=n,tt=t,/noplot
power,'_GWh','_GWh',k=k,spec1=spec0 ,spec2=spec0 ,i=n,tt=t,/noplot,/lkscale
power,'_Str','_Str',k=k,spec1=spec1 ,spec2=spec1 ,i=n,tt=t,/noplot,/lkscale
power,'_Tpq','_Tpq',k=k,spec1=spec2 ,spec2=spec2 ,i=n,tt=t,/noplot,/lkscale
power,'_SCL','_SCL',k=k,spec1=spec3 ,spec2=spec3 ,i=n,tt=t,/noplot,/lkscale
power,'_VCT','_VCT',k=k,spec1=spec4 ,spec2=spec4 ,i=n,tt=t,/noplot,/lkscale
;spec2=reform(spec2,576,1)
s0=size(spec0) & nt0=s0[2]
s1=size(spec1) & nt1=s1[2]
s2=size(spec2) & nt2=s2[2]
s3=size(spec3) & nt3=s3[2]
s4=size(spec4) & nt4=s4[2]
print,'nt0,nt1,nt2=',nt0,nt1,nt2,nt3,nt4
!x.title='!8k!6'
!y.title='!6'
;
plot_oo,yr=[1e-22,1e-3],k[1:*],spec0[1:*,nt0-1]*k[1:*]^2,li=1
;plot_oo,yr=[1e-22,1e10],k[1:*],spec0[1:*,nt0-1]*k[1:*]^2,li=1
;plot_oo,yr=[1e-7,1e-4],k[1:*],spec0[1:*,nt0-1]*k[1:*]^2,li=1
oplot,k[1:*],spec1[1:*,nt1-1],col=122
oplot,k[1:*],spec2[1:*,nt2-1],col=55,li=1,thick=5
oplot,k[1:*],spec3[1:*,nt3-1]*k[1:*]^4,col=155
loadct,6
oplot,k[1:*],spec4[1:*,nt4-1]*k[1:*]^2,col=122,li=2
loadct,5
;
xx=140. & dx=100.
legend,xx,dx,1e-14,0,col=155,'scalar'
legend,xx,dx,1e-16,1,col= 55,'full  '
loadct,6
legend,xx,dx,1e-18,1,col=122,'vector'
loadct,5
;
xx=1000. & dx=400.
legend,xx,dx,1e-16,0,col=122,'projected'
legend,xx,dx,1e-18,1,        '!6d!8h!6/d!8t!6'
;
cwd,run
print,"$convert idl.ps /D/Print/GW_prog2019/pSCL_"+run+".png"

END