@parameters
;
default,ivar,-1
default,lread,1
if lread then begin
  pc_read_param,/param2,obj=param
  if ivar ge 0 then pc_read_var,/bb,/jj,/oo,/trimall,obj=var,ivar=ivar else pc_read_var,/bb,/jj,/oo,/trimall,obj=var
  x=var.x
  y=var.y
  z=var.z
  t=var.t
  ;B0=param.B_ext(2)
  lread=0
endif
nz=n_elements(z)
;
ab=dot(var.aa,var.bb)
azbz=var.aa[*,*,2]*var.bb[*,*,2]
;
;!p.multi=[0,2,2]
!p.multi=0
;tvscl,var.ss(*,*,nz-1),0
;tvscl,var.ss(*,*,nz*3/4),1
;plot,zzz,haver(var.ss)
;tvscl,var.bb(*,*,nz-1,2)
END