@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,/oo,/jj,/bb,/trimall,obj=var,ivar=ivar else pc_read_var,/oo,/jj,/bb,/trimall,obj=var
  if ivar ge 0 then pc_read_var,/oo,/trimall,obj=var,ivar=ivar else pc_read_var,/oo,/trimall,obj=var
  ;if ivar ge 0 then pc_read_var,/jj,/bb,/oo,/trimall,obj=var,ivar=ivar else pc_read_var,/jj,/bb,/oo,/trimall,obj=var
  t=var.t
  x=var.x
  y=var.y
  z=var.z
  lread=0
endif
nz=n_elements(zzz)
;
;!p.multi=[0,2,2]
!p.multi=0
;stop
;tvscl,var.ss(*,*,nz-1),0
;tvscl,var.ss(*,*,nz*3/4),1
;plot,zzz,haver(var.ss)
;tvscl,var.bb(*,*,nz-1,2)
;
;
;  compute 2-d pdf of lnrho and lnB:
;
!p.multi=[0,2,1]
nx=500 & ny=400
lr=var.lnrho
xr=minmax(lr)
;
;lb=.5*alog(dot2(var.bb))
;yr=minmax(lb)
;pdf2d,lr,lb,xx,yy,ww,nx,ny,xr=xr,yr=yr,/unnorm
;contour,alog(ww>1.),xx,yy,nlev=30,/fil
;save,file='Larger_pdf2d_lr_lb_'+str(ivar)+'.sav',xx,yy,ww,t
;
lo=.5*alog(dot2(var.oo))
yr=minmax(lo)
pdf2d,lr,lo,xx,yy,ww,nx,ny,xr=xr,yr=yr,/unnorm
contour,alog(ww>1.),xx,yy,nlev=30,/fil
save,file='Larger_pdf2d_lr_lo_'+str(ivar)+'.sav',xx,yy,ww,t
;
!p.multi=0
END