;$Id: pTM_kpq.pro,v 1.7 2014/03/03 07:26:52 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/pTM_kpq_short2304pm1_kf60b_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pTM_kpq_short288pm1_kf60a_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pTM_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='TM_kpq.sav'
restore,file & TM_kpq=TM_kpq/max(abs(TM_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(TM_kpq,3),2),xtit='!8k!6',ytit='!7R!8!dp!n!7R!8!dq!n!8T!d8kpq!n!6'
oplot,ki,ki*0 & xyouts,siz=siz,xx,-30,'!6(a)'
;
plot,ki,total(total(TM_kpq,3),1),xtit='!8p!6',ytit='!7R!8!dk!n!7R!8!dq!n!8T!d8kpq!n!6'
oplot,ki,ki*0 & xyouts,siz=siz,xx,-30,'!6(d)'
;
plot,ki,total(total(TM_kpq,2),1),xtit='!8q!6',ytit='!7R!8!dk!n!7R!8!dp!n!8T!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]
;
TM_3pq=reform(TM_kpq(4,*,*))
print,minmax(TM_3pq)
contour,clip(TM_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
;
TM_3pq=reform(TM_kpq(*,4,*))
print,minmax(TM_3pq)
contour,clip(TM_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
;
TM_3pq=reform(TM_kpq(*,*,4))
print,minmax(TM_3pq)
contour,clip(TM_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
;
TM_3pq=reform(TM_kpq(8,*,*))
print,minmax(TM_3pq)
contour,clip(TM_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
;
TM_3pq=reform(TM_kpq(*,8,*))
print,minmax(TM_3pq)
contour,clip(TM_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
;
TM_3pq=reform(TM_kpq(*,*,8))
print,minmax(TM_3pq)
contour,clip(TM_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