pro show_slab
  !P.MULTI=[0,3,3]

  gamma=5.0/3.0
  gam1 = (gamma-1)

  rho_ic=0.125
  P_ic = 1.0
  B_ic = 1.0
  csnd = sqrt(gamma * P_ic/rho_ic)
  valp = sqrt(B_ic^2/rho_ic)
  
  for i=0,30 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"
     readnew,name,b,"BFLD"

     bb=SQRT(b[0,*]^2+b[1,*]^2+b[2,*]^2)
     vv=SQRT(v[0,*]^2+v[1,*]^2+v[2,*]^2)
     
     p = (gamma-1)*u*rho
     e=u*rho+0.5*bb^2+0.5*vv^2*rho
     u=gam1*u*rho
     
     x0=50
     dx=20
     
     plot,x[0,*],rho[*],psym=3,xtit='x',ytit='rho',xth=2,yth=2,charsize=2,charth=2,yr=[-0.1,1.2],yst=1,xr=x0+dx*[-0.6,1.4],tit=snr

     plot,x[0,*],p[*],psym=3,xtit='x',ytit='P',xth=2,yth=2,charsize=2,charth=2,yr=[-0.1,1.2],yst=1,xr=x0+dx*[-0.6,1.4],tit=string(h.time)
     
     plot,x[0,*],e[*],psym=3,xtit='x',ytit='E',xth=2,yth=2,charsize=2,charth=2,yr=[0.6,2.6],yst=1,xr=x0+dx*[-0.6,1.4]

     plot,x[0,*],v[0,*],psym=3,xtit='x',ytit='vx',xth=2,yth=2,charsize=2,charth=2,yr=[-0.4,0.8],yst=1,xr=x0+dx*[-0.6,1.4]

     plot,x[0,*],v[1,*],psym=3,xtit='x',ytit='vy',xth=2,yth=2,charsize=2,charth=2,yr=[-2.0,0.4],yst=1,xr=x0+dx*[-0.6,1.4]
     
     plot,x[0,*],v[2,*],psym=3,xtit='x',ytit='vz',xth=2,yth=2,charsize=2,charth=2,yr=[-1.2,1.2],yst=1,xr=x0+dx*[-0.6,1.4]
     
     plot,x[0,*],b[1,*],psym=3,xtit='x',ytit='By',xth=2,yth=2,charsize=2,charth=2,yr=[-1.2,1.2],yst=1,xr=x0+dx*[-0.6,1.4]
     
     plot,x[0,*],b[2,*],psym=3,xtit='x',ytit='Bz',xth=2,yth=2,charsize=2,charth=2,yr=[-1.2,1.2],yst=1,xr=x0+dx*[-0.6,1.4]
     
     plot,x[0,*],b[0,*],psym=3,xtit='x',ytit='Bx',xth=2,yth=2,charsize=2,charth=2,yr=[-1.2,1.2],yst=1,xr=x0+dx*[-0.6,1.4]
     
     save_screen,"frame_"+snr
  end
  stop
end