shock = fltarr(mx,my,mz)*one
 uu = fltarr(mx,my,mz,3)*one
 lnrho = fltarr(mx,my,mz)*one
 ss = fltarr(mx,my,mz)*one
 cooling = fltarr(mx,my,mz)*one