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/pumzm1_256_512a_Om02.ps.gz
;  scp $dcsc/pencil-code/axel/conv-slab/chit/256_512a_Om02/idl.ps.gz fig/256_512a_Om02.ps.gz
;
nlev=20
bar='!20!s!u$!n!r!6'
!p.charsize=2
!x.margin=[10.6,1.2]
;
loadct,0
!p.multi=[0,1,2]
!x.title='!8z!6'
!y.title=bar+'u!8!dy!n!6'
!p.title='!6'
error_range,xya.TTmz,2,fm,fml,fmu,/plot,x=zzz
TTm=fm & TTml=fml & TTmu=fmu
;
error_range,xya.rhomz,2,fm,fml,fmu,/plot,x=zzz
rhom=fm & rhoml=fml & rhomu=fmu
;
save,file='TTmz.sav',zzz,TTm,TTml,TTmu,rhom,rhoml,rhomu
END