; $Id$
;
;  expansion shock, to be compared with Fig 1 of S.A.E.G. Falle (2002)
;  ApJL 577, L123-L126 "Rarefaction shocks, shock errors, and low
;  order accuracy in ZEUS"
;
exact=[[0.52,  0.4993], $
[0.5828,  0.4993], $
[0.5828,  0.5414], $
[0.7126,  0.5414], $
[0.7216,  0.5299], $
[0.7843,  0.5299], $
[0.7843,  0.4318], $
[0.8298,  0.4318], $
[0.8298,  0.1961], $
[0.9649,  0.1961], $
[0.9649,  0.1019], $
[1.020,  0.1019]]

tex=texsyms()

!p.multi=0
!p.charsize=1.5
!p.title='t=30'
!x.title='x'
!y.title='!7q!6'
circ_sym,.4,0
;
;  in Falle's paper shock was initially at x=?
;
x0=0

;
;  exclude ghost zones
;
pc_read_var,obj=data,variables=['rho','tt','uu','pp'],/magic,/trimall
;
plot,data.x,data.rho,yst=3,ps=8,back=255,col=1,ytitle=tex.rho
oplot,exact[0,*]*1E3-800,exact[1,*],col=122
;
print,'import riemann_bfield.png'
END