pro show_galaxy
  !P.MULTI=[0,2,1]

  loadct,5
  
  for i=0,100 do begin
     snr=string(i,form='(i03)')
     name="snap_"+snr
     
     readnew,name,h,"HEAD"
     
     readnew,name,x_bulge,"POS",part=3
     readnew,name,x_disk,"POS",part=2
     readnew,name,x_gas,"POS",part=0
     IF(h.npart[4] GT 0) THEN readnew,name,x_stars,"POS",part=4

     dx=30

     plot,[1],[1],xtit='x [kpx]',ytit='y [kpc]',xth=2,yth=2,charsize=2,charth=2,yr=dx*[-1,1],yst=1,xr=dx*[-1,1],xst=1,tit=snr,/nodata

     oplot,x_disk[0,*],x_disk[1,*],psym=3,col=160
     oplot,x_bulge[0,*],x_bulge[1,*],psym=3,col=140
     oplot,x_gas[0,*],x_gas[1,*],psym=3,col=50
     IF(h.npart[4] GT 0) THEN oplot,x_stars[0,*],x_stars[1,*],psym=3,col=100
     
     plot,[1],[1],xtit='x [kpx]',ytit='z [kpc]',xth=2,yth=2,charsize=2,charth=2,yr=dx*[-1,1],yst=1,xr=dx*[-1,1],xst=1,tit=snr,/nodata
     oplot,x_disk[0,*],x_disk[2,*],psym=3,col=160
     oplot,x_bulge[0,*],x_bulge[2,*],psym=3,col=140
     oplot,x_gas[0,*],x_gas[2,*],psym=3,col=50
     IF(h.npart[4] GT 0) THEN oplot,x_stars[0,*],x_stars[2,*],psym=3,col=100

     save_screen,"frame_"+snr
  end
  stop
end