rho=exp(pc_noghost(lnrho))
ntot=1d0*nx*ny*nz
print,n_elements(where(rho gt 5.2))/ntot
print,n_elements(where(rho lt 5.0 and rho gt 1.))/ntot
print,n_elements(where(rho lt 0.9))/ntot
;
print
Rgas=8.314d7
mu=.62
unit_rho=1d-24
unit_pp=1d-14
unit_TT=1.
Rgas_code=Rgas*unit_rho*unit_TT/(unit_pp*mu)
cp=2.5*Rgas_code
cs2=exp(2./3.*lnrho+1.66667*ss/cp)
TT=cs2/(1.6667*Rgas_code/mu)
TTT=TT(l1:l2,3,3)
print,n_elements(where(TTT lt 420.))/ntot
END