;$Id: pTM_kpq_comp.pro,v 1.5 2014/07/27 17:35:58 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=6,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/pTM_kpq_comp.ps
;
k0=60.
dir1='short2304pm1_kf60b_noforce'
dir2='1024_4x32x32_Pm1'
;
!p.multi=[0,3,1]
!p.charsize=1.7
!x.margin=[7.0,0.5]
!y.margin=[3.2,1.9]
xx=2.3 & yy=21
col=255
siz=1.3
siz2=.9
;
;  read & rescale
;
file='TM_kpq.sav'
restore,'../'+dir1+'/'+file & TM_kpq1=TM_kpq/max(abs(TM_kpq))
restore,'../'+dir2+'/'+file & TM_kpq2=TM_kpq/max(abs(TM_kpq))
;
m=24
ki=1+findgen(m)
print,ki
;
;  panel (a)
;
!x.range=[1,24]/k0
!y.range=[-25,20]
plot,ki/k0,total(total(TM_kpq1,3),2),xtit='!8k!6/!8k!6!d0!n',ytit='!7R!8!dp!n!7R!8!dq!n!8T!d8kpq!n!6'
oplot,ki/k0,total(total(TM_kpq2,3),2),li=2
oplot,ki/k0,ki*0 & xyouts,siz=siz,19.5/k0,13,'!6(a)'
xyouts,.11,15.,siz=siz2,'with mag.!chelicity'
xyouts,.23,-8.,siz=siz2,'no mag.!chelicity'
;
;  panel (b)
;
lev=.4*grange(-1,1,21)
mm=minmax(lev)
!x.range=[1,24]/k0
!y.range=[1,24]/k0
;
TMx=reform(TM_kpq1(4,*,*))
contour,clip(TMx,mm),ki/k0,ki/k0,lev=lev,/fil,xtit='!8p!6/!8k!6!d0!n',ytit='!8q!6/!8k!6!d0!n',tit='!8k!6/!8k!6!d0!n=0.07'
xyouts,siz=siz,xx/k0,yy/k0,'!6(b)',col=col
;
;  panel (c)
;
!x.range=[1,23.9]/k0
TMx=reform(TM_kpq1(*,4,*))
contour,clip(TMx,mm),ki/k0,ki/k0,lev=lev,/fil,xtit='!8k!6/!8k!6!d0!n',ytit='!8q!6/!8k!6!d0!n',tit='!8p!6/!8k!6!d0!n=0.07'
xyouts,siz=siz,xx/k0,yy/k0,'!6(c)',col=col
;
;  color bar
;
bcharsize=1.5 & divbar=2
xx0=.975 & dxx=.01
yy1=.20 & dyy=.18
colorbar,pos=[xx0,yy1,xx0+dxx,yy1+dyy], range=minmax(lev),/left,/vert, $
  form='(f4.1)',charsize=bcharsize,div=divbar,ytit='!6',col=255
;
;  inset to panel (b)
;
!x.range=[1,14]/k0
!y.range=[1,14]/k0
yy1=.60 & yy2=.88 & xx1=.55 & dxx=.105
!p.position=[xx1,yy1,xx1+dxx,yy2]
TMx=reform(TM_kpq2(4,*,*))
contour,clip(TMx,mm),ki/k0,ki/k0,lev=lev,/fil,xtit='!6',ytit='!6',/noera,col=255,xticks=1,yticks=1,xtickv=[.1,.2],ytickv=[.1,.2]
;
;  inset to panel (c)
;
xx1=.884
!p.position=[xx1,yy1,xx1+dxx,yy2]
TMx=reform(TM_kpq2(*,4,*))
contour,clip(TMx,mm),ki/k0,ki/k0,lev=lev,/fil,xtit='!6',ytit='!6',/noera,col=255,xticks=1,yticks=1,xtickv=[.1,.2],ytickv=[.1,.2]
;
!p.position=0
END