left state
rhLL = uL[i,1]
uuLL = uL[i,2]/rhLL
eeLL = uL[i,3]/rhLL
ppLL = gm*(eeLL*rhLL - 0.5*rhLL*(uuLL*uuLL))
```

```
aaLL = sqrt(abs(gamma*ppLL/rhLL))