;$Id: pSM_kpq.pro,v 1.1 2014/07/19 06:28:01 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=17,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/pSM_kpq_short2304pm1_kf60b_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pSM_kpq_short288pm1_kf60a_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pSM_kpq_1024_4x32x32_Pm1.ps
;
!p.charsize=1.7
!x.margin=[7.0,0.5]
!y.margin=[3.2,2.5]
xx=2.3 & yy=19.
col=255
siz=1.8
;
;  read & rescale
;
file='SM_kpq.sav'
restore,file & SM_kpq=SM_kpq/max(abs(SM_kpq))
;
m=24
ki=1+findgen(m)
print,ki
;
!p.multi=[0,3,3]
!x.range=[1,24]
!y.range=[-36,10]
plot,ki,total(total(SM_kpq,3),2),xtit='!8k!6',ytit='!7R!8!dp!n!7R!8!dq!n!8S!d8kpq!n!6'
oplot,ki,ki*0 & xyouts,siz=siz,xx,-30,'!6(a)'
;
plot,ki,total(total(SM_kpq,3),1),xtit='!8p!6',ytit='!7R!8!dk!n!7R!8!dq!n!8S!d8kpq!n!6'
oplot,ki,ki*0 & xyouts,siz=siz,xx,-30,'!6(d)'
;
plot,ki,total(total(SM_kpq,2),1),xtit='!8q!6',ytit='!7R!8!dk!n!7R!8!dp!n!8S!d8kpq!n!6'
oplot,ki,ki*0 & xyouts,siz=siz,18,-30,'!6(g)'
;
lev=.4*grange(-1,1,21)
mm=minmax(lev)
!x.range=[1,24]
!y.range=[1,24]
;
SM_3pq=reform(SM_kpq(4,*,*))
print,minmax(SM_3pq)
contour,clip(SM_3pq,mm),ki,ki,lev=lev,/fil,xtit='!8p!6',ytit='!8q!6',tit='!8k!6=4'
xyouts,siz=siz,xx,yy,'!6(b)',col=col
;
SM_3pq=reform(SM_kpq(*,4,*))
print,minmax(SM_3pq)
contour,clip(SM_3pq,mm),ki,ki,lev=lev,/fil,xtit='!8k!6',ytit='!8q!6',tit='!8p!6=4'
xyouts,siz=siz,xx,yy,'!6(e)',col=col
;
SM_3pq=reform(SM_kpq(*,*,4))
print,minmax(SM_3pq)
contour,clip(SM_3pq,mm),ki,ki,lev=lev,/fil,xtit='!8k!6',ytit='!8p!6',tit='!8q!6=4'
xyouts,siz=siz,xx,yy,'!6(h)',col=col
;
;  second
;
SM_3pq=reform(SM_kpq(8,*,*))
print,minmax(SM_3pq)
contour,clip(SM_3pq,mm),ki,ki,lev=lev,/fil,xtit='!8p!6',ytit='!8q!6',tit='!8k!6=8'
xyouts,siz=siz,xx,yy,'!6(c)',col=col
;
SM_3pq=reform(SM_kpq(*,8,*))
print,minmax(SM_3pq)
contour,clip(SM_3pq,mm),ki,ki,lev=lev,/fil,xtit='!8k!6',ytit='!8q!6',tit='!8p!6=8'
xyouts,siz=siz,xx,yy,'!6(f)',col=col
;
SM_3pq=reform(SM_kpq(*,*,8))
print,minmax(SM_3pq)
contour,clip(SM_3pq,mm),ki,ki,lev=lev,/fil,xtit='!8k!6',ytit='!8p!6',tit='!8q!6=8'
xyouts,siz=siz,xx,yy,'!6(i)',col=col
;
bcharsize=1.8 & divbar=2
xx0=.3 & dxx=.02
yy1=.41 & dyy=.19
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev),/left,/vert, $
  form='(f4.1)',charsize=bcharsize,div=divbar,ytit='!6',col=255
;
END