file='VAR27'
file='VAR26'
file='VAR28'
file='VAR30'
;ixs=1230 & ixe=ixs+255
;iys= 315 & iye=iys+255
ixs=0 & ixe=2047
iys=0 & iye=2047
izs=1024 & ize=1024
izs=500 & ize=500
izs=1000 & ize=1000
izs=1023 & ize=1023
pc_read_subvol_raw,varfile=file,obj=var,xs=ixs,xe=ixe, ys=iys,ye=iye,zs=izs,ze=ize, allprocs=0,/single,/addghost
;divu=(div(var[*,*,*,0:2]))[3:2050,3:2050,3]
;save,file=file+'_divu.sav',divu,x,y
;oz=(curl(var[*,*,*,0:2]))[3:2050,3:2050,3,2]
;save,file=file+'_oz.sav',oz,x,y
lnrho=(var[*,*,*,3])[3:2050,3:2050,3]
save,file=file+'_lnrho.sav',lnrho,x,y
;
length=4*!pi
nx=2048
dx=length/nx
x=dx*(.5+findgen(nx)-nx/2)
y=x
nlev=17 & color,/sgi,steps=nlev
;contour,divu,lev=grange(-1,1,30)*5,/iso,/fil,x,y 
;levo=grange(-1,1,30)*.5
;contour,clip(oz,minmax(levo)),lev=levo,/iso,/fil,x,y 
;levo=grange(-1,1,30)*.5
;contour,clip(oz,minmax(levo)),lev=levo,/iso,/fil,x,y 
levr=grange(-1,1,30)*6.
contour,clip(lnrho,minmax(levr)),lev=levr,/iso,/fil,x,y 
END