if !d.name eq 'PS' then begin
  device,xsize=22,ysize=8,yoffset=3
  !p.charthick=4 & !p.thick=4 & !x.thick=4 & !y.thick=4
end
;
;  plots only the mean uy velocity
;
;  scp $dcsc/pencil-code/axel/conv-slab/chit/256_512a_Om02/idl.ps.gz fig/puyuzmz_256_512a_Om02.ps.gz
;
;  need to read structure xya first, using pumt.pro
;
nlev=20
bar='!20!s!u$!n!r!6'
!p.charsize=1.6
!x.margin=[10.6,1.2]
;
loadct,0
!p.title='!6'
!x.title='!8z!6'
!y.title=bar+'u!8!dy!n!6'
error_range,xya.uyuzmz,2,fm,fml,fmu,/plot,x=zzz
;
save,file='uyuzmz.sav',fm,fml,fmu,zzz
;
END