















































  






































































 


































































































































































 

















































































































































































































































































































































































































































 






































































































 










 








































































































































































































































 










                                
                                

                                
                                

                                
                                

                                
















































 
subroutine num_x_fluxM(dens,sx,sy,sz,tau,Bx,By,Bz,&
     densf,sxf,syf,szf,tauf,Bxf,Byf,Bzf,vxt,vyt,vzt,pressstar,&
     bsubx,bsuby,bsubz,ab0,w,det,alp,beta)

  implicit none

  REAL*8 :: dens,sx,sy,sz,tau,Bx,By,Bz
  REAL*8 :: densf,sxf,syf,szf,tauf,Bxf,Byf,Bzf
  REAL*8 :: vxt,vyt,vzt,bsubx,bsuby,bsubz,ab0,w
  REAL*8 :: det,alp,beta,pressstar
  REAL*8 :: velm
  REAL*8 :: sdet,psipstar

  sdet=sqrt(det)
  psipstar=pressstar*sdet


  velm=vxt+beta/alp



  densf = dens * vxt

  sxf = sx*vxt+psipstar-bsubx*Bx/w

  syf = sy*vxt-bsuby*Bx/w

  szf = sz*vxt-bsubz*Bx/w


  tauf = tau*vxt + psipstar*velm - ab0*Bx/w


  bxf = 0.0
  byf = By * vxt - Bx*vyt
  bzf = Bz * vxt - Bx*vzt

end subroutine num_x_fluxM

