pro show_slab
  !P.MULTI=[0,1,4]

  gamma=5.0/3.0
  
  for i=11,12 do begin
     snr=string(i,form='(i03)')
     name="snap_"+snr
     
     readnew,name,h,"HEAD"
     readnew,name,x,"POS"
     readnew,name,v,"VEL"
     readnew,name,rho,"RHO"
     readnew,name,u,"U"

     p = (gamma-1)*u*rho
     
     plot,x[0,*],p[*],psym=4,xtit='x',ytit='P',xth=2,yth=2,charsize=2,charth=2,yr=[0.4,1.1],yst=1,xr=[10,30]
     plot,x[0,*],u[*],psym=4,xtit='x',ytit='e',xth=2,yth=2,charsize=2,charth=2,yr=[0.4,1.6],yst=1,xr=[10,30]
     plot,x[0,*],rho[*],psym=4,xtit='x',ytit='rho',xth=2,yth=2,charsize=2,charth=2,yr=[0.5,2.5],yst=1,xr=[10,30]
     plot,x[0,*],v[0,*],psym=4,xtit='x',ytit='vx',xth=2,yth=2,charsize=2,charth=2,yr=[0.0,0.5],yst=1,xr=[10,30]
     
     save_screen,"frame_"+snr
;     stop
  end
  stop
end