Total pressure not set generally

[BigTailWolfLSL]

In lecture 7 "HLLC_Flux.f90", "PtotR" and "PtotL" variables may be set problem based (only for that example problem).

I guass the Ptot stands for "total pressure" which Ptot=P(from ideal plasma EoS)+0.5*\mathbf{B}^2. In our case, take "PtotR = gamma1*(ER - 0.5d0rhoR(u1R2 + u2R2 + u3R**2))" for example. PtotR=gamma1*(\rho*\epsilon+0.5*\mathbf{B}^2)=P+gamma1 (1/2) \mathbf{B}^2≠Ptot.

However, in shock tube test (our example), gamma=2.0 and gamma1=1.0 just make it works in this case PtotR=Ptot. But it is not general.