;
;  compute 2-d pdf of lnrho and lnomega:
;
nx=500 & ny=400
lr=var.lnrho
lo=.5*alog(dot2(var.oo))
xr=minmax(lr)
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
;
END