if !d.name eq 'PS' then begin
  device,xsize=18,ysize=8,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  Note: before running this routine, you must run pumt.pro before.
;  Unlike pumzt.pro, this routine does not plot the 3 stress components.
;
;  scp $dcsc/pencil-code/axel/conv-slab/chit/256_512a_Om02/idl.ps.gz fig/pumzt1_256_512a_Om02.ps.gz
;  scp $dcsc/pencil-code/axel/conv-slab/chit/256_128a_Om02_theta90/idl.ps.gz fig/256_128a_Om02_theta90_umzt.ps.gz
;
nlev=20
bar='!20!s!u$!n!r!6'
!p.multi=[0,3,1]
!p.charsize=1.7
!x.title='!6'
!x.margin=[6.6,.5]
;
!x.title='!8t!6'
!y.title='!8z!6'
!p.title=bar+'u!8!dx!n!6'
contour,transpose(xya.uxmz),/fill,nlev=nlev,tt,zzz
;
!y.title='!6'
!p.title=bar+'u!8!dy!n!6'
contour,transpose(xya.uymz),/fill,nlev=nlev,tt,zzz
;
!p.title=bar+'U!8!dz!n!6'
contour,transpose(xya.uzmz),/fill,nlev=nlev,tt,zzz
;
!p.multi=0
END